format: text

text:
# QCD factorization for 𝐵→𝑃𝑃 *footnote **footnote *Supported in part by National Natural Science Foundation of China and State Commission of Science and Technology of China

## I introduction

The study of $`B`$ decays plays an important role in understanding the origin of $`CP`$ violation and physics of heavy flavor. We expect that the parameters of the Cabibbo-Kobayashi-Maskawa (CKM) matrix in the standard model, for instance, the three angles $`\alpha `$, $`\beta `$ and $`\gamma `$ in the unitary triangle, can be well-determined from $`B`$ decays, especially from the charmless non-leptonic two-body B decays. Experimentally, many $`B`$ experiment projects have been running (CLEO, BaBar, Belle etc.), or will run in forthcoming years (BTeV, CERN LHCb, DESY HeraB etc.). With the accumulation of the data, the theorists will be urged to gain a deeper sight into $`B`$ decays, and to reduce the theoretical errors in determining the CKM parameters from the experimental data.

In the theoretical frame, the standard approach to deal with such decays is based on the low-energy effective Hamiltonian which is obtained by the Wilson operator product expansion method (OPE). In this effective Hamiltonian, the short-distance contributions from the scale above $`\mu m_b`$ have been absorbed into the Wilson coefficients with the perturbative theory and renormalization group method. The Wilson coefficients have been evaluated to next-to-leading order. Then the main task in studying non-leptonic two-body $`B`$ decays is to calculate the hadronic matrix elements of the effective operators. However, we do not have a reliable approach to evaluate them from the first principles of QCD dynamics up to now.

Generally, we must resort to the factorization assumption to calculate the hadronic matrix elements for non-leptonic $`B`$ decays, in which the hadronic matrix element of the effective operator (in general, which is in the form of current-current four-quark operator) can be approximated as a product of two single current hadronic matrix elements; then it is parameterized into meson decay constant and meson-meson transition form factor. The most popular factorization model is the Bauer-Stech-Wirbel (BSW) model. In many cases, BSW model achieves great success, which can predict the branching ratios of many modes of non-leptonic $`B`$ decays in correct order of magnitude. This factorization assumption does hold in the limit that the soft interactions in the initial and final states can be ignored. It seems that the argument of color-transparency can give reasonable support to the above limit. Because $`b`$ quark is heavy, the quarks from $`b`$ quark decay move so fast that a pair of quarks in a small color-singlet object decouple from the soft interactions. But the shortcomings of this simple model are obvious. First, the renormalization scheme and scale dependence in the hadronic matrix elements of the effective operators are apparently missed. Then the full decay amplitude predicted by BSW model remains dependent on the renormalization scheme and scale, which are mainly from Wilson coefficients. In past years, many researchers improved the simple factorization scheme and made many remarkable progresses, such as scheme and scale independent effective Wilson coefficients, effective color number which is introduced to compensate the ‘non-factorizable’ contributions, etc. Furthermore, some progresses in nonperturbative methods, such as lattice QCD, QCD sum rule etc., allow us to compute many non-perturbative parameters in $`B`$ decays, such as the meson decay constants and meson-meson transition form factors. Every improvement allows us to have a closer look at the $`B`$ nonleptonic decays.

Except for the factorization approximation, another important approach has been applied to study many $`B`$ exclusive hadronic decay channels, such as $`BD\pi `$, $`\pi \pi `$, $`\pi K`$ etc. This is PQCD method. In this method, people assumes that $`B`$ exclusive hadronic decay is dominated by hard gluon exchange. It is analogous to the framework of perturbative factorization for exclusive processes in QCD at large momentum transfer, such as the calculation of the electromagnetic form factor of the pion . The decay amplitude for $`B`$ decay can be written as a convolution of a hard-scattering kernel with light-cone wave functions of the participating mesons. Furthermore, in Ref. the Sudakov suppression has been taken into account.

Two year ago, Beneke, Buchalla, Neubert, and Sachrajda (BBNS) gave a QCD factorization formula in the heavy quark limit for the decays $`B\pi \pi `$. They pointed out that the radiative corrections from hard gluon exchange can be calculated by use of the perturbative QCD method if one neglects the power contributions of $`\mathrm{\Lambda }_{QCD}/m_b`$. This factorization formula can be justified in case that the ejected meson from the $`b`$ quark decay is a light meson or an onium, no matter whether the other recoiling meson which absorbs the spectator quark in $`B`$ meson is light or heavy. But for the case that the ejected meson is in an extremely asymmetric configuration, such as D meson, this factorization formula does not hold. The contributions from the hard scattering with the spectator quark in B meson are also involved in their formula. This kind of contribution cannot be contained in the naive factorization. But it appears in the order of $`\alpha _s`$. So they said that the naive factorization can be recovered if one neglects the radiative corrections and power $`\mathrm{\Lambda }_{QCD}/m_b`$ suppressed contributions in the QCD factorization, and the ‘non-factorizable’ contributions in the naive factorization can be calculated perturbatively, then we do not need a phenomenological parameter $`N_c^{eff}`$ to compensate the ‘non-factorizable’ effects any more.

This QCD factorization (BBNS approach) has been applied to study many $`B`$ meson decay modes, such as $`BD^{()}\pi ^{}`$, $`\pi \pi `$, $`\pi K`$ and other interesting channels. Some theoretical generalizations of BBNS approach have also been made, such as the chirally enhanced power corrections from the twist-3 light-cone distribution amplitudes of the light pseudoscalar mesons. In this work, we will take a closer look at this issue. This work is organized as follows: Sect. II is devoted to a sketch of the low energy effective Hamiltonian; in Sect.III, we will give a detailed overview of QCD factorization, in which some elaborate calculations are shown, especially for the chirally enhanced power corrections; Sect. IV is for some detailed discussions and comparison of BBNS approach to the generalized factorization and PQCD method; we conclude in Sect.V with a summary.

## II Effective Hamiltonian — First Step Factorization

$`B`$ decays involve three characteristic scales which are strongly ordered: $`m_Wm_b\mathrm{\Lambda }_{QCD}`$. How to separate or factorize these three scales is the most essential question in B hadronic decays.

With the operator product expansion method (OPE), the relevant $`|\mathrm{\Delta }B|=1`$ effective Hamiltonian is given by :

$`_{eff}`$ $`=`$ $`{\displaystyle \frac{G_F}{\sqrt{2}}}[{\displaystyle \underset{q=u,c}{}}v_q(C_1(\mu )Q_1^q(\mu )+C_2(\mu )Q_2^q(\mu )+{\displaystyle \underset{k=3}{\overset{10}{}}}C_k(\mu )Q_k(\mu ))`$ (2)
$`v_t(C_{7\gamma }(\mu )Q_{7\gamma }(\mu )+C_{8G}(\mu )Q_{8G}(\mu ))]+h.c.,`$
where $`v_q=V_{qb}V_{qd}^{}`$(for $`bd`$ transition) or $`v_q=V_{qb}V_{qs}^{}`$(for $`bs`$ transition) and $`C_i(\mu )`$ are the Wilson coefficients which have been evaluated to next-to-leading order approximation with the perturbative theory and renormalization group method.

In the Eq.(1), the four-quark operators $`Q_i`$ are given by

$$\begin{array}{c}\begin{array}{cc}Q_1^u=(\overline{u}_\alpha b_\alpha )_{VA}(\overline{q}_\beta u_\beta )_{VA}\hfill & Q_1^c=(\overline{c}_\alpha b_\alpha )_{VA}(\overline{q}_\beta c_\beta )_{VA}\hfill \\ Q_2^u=(\overline{u}_\alpha b_\beta )_{VA}(\overline{q}_\beta u_\alpha )_{VA}\hfill & Q_2^c=(\overline{c}_\alpha b_\beta )_{VA}(\overline{q}_\beta c_\alpha )_{VA}\hfill \\ Q_3=(\overline{q}_\alpha b_\alpha )_{VA}\underset{q^{}}{}(\overline{q}_\beta ^{^{}}q_\beta ^{^{}})_{VA}\hfill & Q_4=(\overline{q}_\beta b_\alpha )_{VA}\underset{q^{}}{}(\overline{q}_\alpha ^{^{}}q_\beta ^{^{}})_{VA}\hfill \\ Q_5=(\overline{q}_\alpha b_\alpha )_{VA}\underset{q^{}}{}(\overline{q}_\beta ^{^{}}q_\beta ^{^{}})_{V+A}\hfill & Q_6=(\overline{q}_\beta b_\alpha )_{VA}\underset{q^{}}{}(\overline{q}_\alpha ^{^{}}q_\beta ^{^{}})_{V+A}\hfill \\ Q_7=\frac{3}{2}(\overline{q}_\alpha b_\alpha )_{VA}\underset{q^{}}{}e_q^{}(\overline{q}_\beta ^{^{}}q_\beta ^{^{}})_{V+A}\hfill & Q_8=\frac{3}{2}(\overline{q}_\beta b_\alpha )_{VA}\underset{q^{}}{}e_q^{}(\overline{q}_\alpha ^{^{}}q_\beta ^{^{}})_{V+A}\hfill \\ Q_9=\frac{3}{2}(\overline{q}_\alpha b_\alpha )_{VA}\underset{q^{}}{}e_q^{}(\overline{q}_\beta ^{^{}}q_\beta ^{^{}})_{VA}\hfill & Q_{10}=\frac{3}{2}(\overline{q}_\beta b_\alpha )_{VA}\underset{q^{}}{}e_q^{}(\overline{q}_\alpha ^{^{}}q_\beta ^{^{}})_{VA}\hfill \end{array}\hfill \end{array}$$
(3)
and

$$Q_{7\gamma }=\frac{e}{8\pi ^2}m_b\overline{q}_\alpha \sigma ^{\mu \nu }(1+\gamma _5)b_\alpha F_{\mu \nu },Q_{8G}=\frac{g}{8\pi ^2}m_b\overline{q}_\alpha \sigma ^{\mu \nu }t_{\alpha \beta }^a(1+\gamma _5)b_\beta G_{\mu \nu }^a,(q=d\mathrm{or}s).$$
(4)
with $`Q_1^q`$ and $`Q_2^q`$ being the tree operators, $`Q_3Q_6`$ the QCD penguin operators, $`Q_7Q_{10}`$ the electroweak penguin operators, and $`Q_{7\gamma }`$, $`Q_{8G}`$ the magnetic-penguin operators.

In this effective Hamiltonian for $`B`$ decays, the contributions from large virtual momenta of the loop corrections from scale $`\mu =𝒪(m_b)`$ to $`m_W`$ are attributed to the Wilson coefficients, and the low energy contributions are fully incorporated into the matrix elements of the operators. So the derivation of the effective Hamiltonian can be called “the first step factorization”.

To evaluate the Wilson coefficients, we must extract them at a large renormalization scale (for example $`\mu =𝒪(m_W)`$ in the standard model) by matching the amplitude of the effective Hamiltonian ($`A_{eff}`$) to that of the full theory ($`A_{full}`$), then evolve them by the renormalization group equations from the scale $`\mu =𝒪(m_W)`$ to the scale $`\mu =𝒪(m_b)`$. It should be noted that the extraction of the Wilson coefficients $`C_i`$ by matching does not depend on the choice of the external states, if we regularize the infrared (and mass) singularities properly. All dependence on the choice of external states only appears in the matrix elements $`Q_i`$, and is not contained in $`C_i`$. So $`C_i`$ only contains the short-distance contributions from the region where the perturbative theory can be applied. But for the matrix elements $`Q_i`$, the long-distance contributions appear, and are process-dependent.

Several years ago, the perturbative corrections to the Wilson coefficients in SM have been evaluated to next-to-leading order with renormalization group method. As we know, the Wilson coefficients are generally renormalization scheme and scale dependent. So, in order to cancel such dependence, we must calculate the hadronic matrix elements of the effective operators to the corresponding perturbative order with the same renormalization scheme and at the same scale, then we can obtain a complete decay amplitude which is free from those unphysical dependences.

## III QCD Factorization For $`BPP`$

After “the first step factorization”, the decay amplitude for $`Bh_1h_2`$ can be written as

$$𝒜(Bh_1h_2)=\underset{i}{}v_iC_i(\mu )h_1h_2|Q_i(\mu )|\text{B},$$
(5)
in which, as mentioned in the previous section, the contributions from the large scale $`m_W`$ down to $`m_b`$ has been separated into the Wilson coefficients $`C_i(\mu )`$. The remaining task is to calculate the hadronic matrix elements of the effective operators. But for the complexity of QCD dynamics, it is difficult to calculate these matrix elements reliably from first principles. The most popular approximation is factorization hypothesis, in which the matrix element of the current-current operator is approximated to a product of two matrix elements of single current operator:

$$h_1h_2|Q_i|\text{B}h_2|J_2|0h_1|J_1|\text{B}.$$
(6)
Obviously, under this approximation, the original hadronic matrix element $`Q_i(\mu )`$ misses the dependence of the renormalization scheme and scale which should be used to cancel the corresponding dependence in the Wilson coefficients $`C_i(\mu )`$. A plausible solution to recover this scale and scheme dependence of $`Q_i`$ is to calculate the radiative corrections. In one-loop level, they can be written as :

$$𝐐=[\widehat{\mathrm{𝟏}}+\frac{\alpha _s}{4\pi }\widehat{𝐦}_s+\frac{\alpha _{em}}{4\pi }\widehat{𝐦}_e]𝐐_{\mathrm{Tree}}.$$
(7)
Here $`\widehat{𝐦}_s`$ and $`\widehat{𝐦}_e`$ represent the one loop corrections of QCD and QED respectively. Then one takes

$$h_1h_2|Q_i|\text{B}_{\mathrm{Tree}}h_2|J_2|0h_1|J_1|\text{B}.$$
(8)
Therefore, the scheme and scale dependence of $`Q_i`$ which are expressed in the form of $`\widehat{𝐦}_s`$ and $`\widehat{𝐦}_e`$ is recovered. But in quark level, $`\widehat{𝐦}_s`$ and $`\widehat{𝐦}_e`$ usually contain infrared divergences if we take the external quarks on-shell. To remove or regularize the infrared divergence, the conventional treatment is to assume that external quarks are off-shell by $`p^2`$. But this introduction of the infrared cutoff $`p^2`$ results in a gauge dependence of one-loop corrections. So how to factorize the infrared part of the matrix elements is a very subtle question. But maybe this question would get an important simplification in the case that the final states of B meson decay are two light mesons.

Two years ago, Beneke, Buchalla, Neubert and Sachrajda proposed a promising QCD factorization method for $`B\pi \pi `$. They pointed out that in the heavy quark limit $`m_b\mathrm{\Lambda }_{QCD}`$, the hadronic matrix elements for $`B\pi \pi `$ can be written in the form

$$\pi \pi |Q|\text{B}=\pi |J_2|0\pi |J_1|B[1+r_n\alpha _s^n+𝒪(\mathrm{\Lambda }_{QCD}/m_b)].$$
(9)
Obviously, the above formula reduces to the naive factorization if we neglect the power corrections in $`\mathrm{\Lambda }_{QCD}/m_b`$ and the radiative corrections in $`\alpha _s`$. They find that the radiative corrections, which are dominated by hard gluon exchange, can be calculated systematically with the perturbative theory in the limit $`m_b\mathrm{}`$, in terms of the convolution of the hard scattering kernel and the light-cone distribution amplitudes of the mesons. This is also similar to the framework of perturbative factorization for exclusive processes in QCD at large momentum transfer, such as the calculation of the electromagnetic form factor of the pion . Then a factorization formula for $`B\pi \pi `$ can be written as :

$$\pi (p^{})\pi (q)|Q_i|B(p)=F^{B\pi }(q^2)_0^1𝑑xT_i^I(x)\mathrm{\Phi }_\pi (x)+_0^1𝑑\xi 𝑑x𝑑yT^{II}(\xi ,x,y)\mathrm{\Phi }_B(\xi )\mathrm{\Phi }_\pi (x)\mathrm{\Phi }_\pi (y).$$
(10)
We call this factorization formalism as QCD factorization or the BBNS approach. In the above formula, $`\mathrm{\Phi }_B(\xi )`$ and $`\mathrm{\Phi }_\pi (x)`$ are the leading-twist wave functions of $`B`$ and pion mesons respectively, and the $`T_i^{I,II}`$ denote hard-scattering kernels which are calculable in perturbative theory. At the order of $`\alpha _s`$, the hard kernels $`T^{I,II}`$ can be depicted by Fig.1. Figures 1(a)-1(d) represent vertex corrections, Figs 1(e) and 1(f) penguin corrections, and Figs 1(g) and 1(h) hard spectator scattering.

In the heavy quark limit, both pions are energetic. The pion ejected from $`b`$ quark decay moves so fast that it can be described by its leading-twist light-cone distribution amplitude. The $`q\overline{q}`$ pair in the ejected pion is produced as a small-size color dipole. Consequently, the ejected pion decouples from the soft gluons at leading order of $`\mathrm{\Lambda }_{QCD}/m_b`$. Of course, only the cancellation of soft gluons is not enough to make the factorization hold, it is necessary that the $`q\overline{q}`$ pair also decouples from the collinear gluons. Both the cancellations of soft gluons and collinear gluons guarantee that the hard kernel $`T_i^I`$ is of infrared finiteness. Contrast to the pion ejected from b quark weak decay, the recoiling pion which picks up the spectator in B meson can not be described by its leading-twist light-cone distribution amplitude(LCDA), because the spectator is transferred to the recoiling pion as a soft quark. Here Beneke et al. take the point of view that the form factor $`F^{B\pi }`$ cannot be calculated perturbatively. If we attempt to calculate the form factor within the perturbative framework, by the naive power counting, we find that the leading twist LCDA of pion does not fall fast enough to suppress the singularity at the endpoint where the quark from $`b`$ decay carries almost all momentum of the pion. It indicates that the contributions to form factor are dominated by the soft gluon exchange. This point of view can be justified also from the calculation of the form factor $`F^{B\pi }`$ by using light-cone sum rule (LCSR), in which the dominated contribution to $`F^{B\pi }`$ comes from the region where the the spectator quark is transferred as a soft quark to the pion. So the transition form factor survives in the factorization formula as a nonperturbative parameter. However, when the spectator quark in $`B`$ meson interacts with a hard gluon from the ejected pion, the recoiling pion can be also described by its light-cone distribution amplitude. This hard spectator scattering is missed in the naive factorization, but calculable in the perturbative QCD at the leading power in $`\mathrm{\Lambda }_{QCD}/m_b`$. So with this factorization formula, the remaining hard part of the hadronic matrix element $`Q_i`$ from the scale about $`m_b`$ has been factorized into the hard scattering kernel, and the long distance contributions are absorbed into the transition form factors and the light-cone wave functions of the participating mesons. Thus this is the “final factorization” for the two-body nonleptonic charmless $`B`$ decays.

An explicitly technical demonstration of the above argument has been presented in one-loop level in Refs.. For $`BD\pi `$, this QCD factorization has been proved to two-loop order. In the literature, the ejected pion is represented by its leading twist light-cone distribution amplitude(LCDA). However, since the mass of $`b`$ quark is not asymptotically large, in particular, some power corrections might be enhanced by certain factors, such as the scale of chiral symmetry breaking $`\mu _\pi =m_\pi ^2/(m_u+m_d)1.5`$ GeV, and have significant effects in studying $`B`$ two-body nonleptonic charmless decays. So, in this manner, the chirally enhanced power corrections must be taken into account. Accordingly, describing the ejected pion by its leading twist LCDA is not enough, the two-particle twist-3 LCDAs must be taken into account. Below, we will show the elaborate results of QCD factorization in these two cases. For illustration, we take $`\overline{B}_d^0\pi ^+\pi ^{}`$ as an example, but the result is easily generalized to the cases that the final states are the other light pseudoscalars.

### A Leading-twist Distribution Amplitude Insertion

When inserting leading-twist LCDA of the light pseudoscalar, in the heavy quark limit, the quark constituents of the ejected pion can be treated as a pair of collinear massless quark and antiquark with the momentum $`uq`$ and $`\overline{u}q`$ respectively ($`q`$ is the momentum of the ejected pion and we take $`q`$ as a hard light-cone momentum in calculation, $`\overline{u}=1u`$), because that the contributions from the transverse momenta of the quarks in ejected pion are power suppressed .

#### 1 Vertex Corrections

Now we move on to the explicit one-loop calculation of the diagram Figs. 1(a)-1(d) for $`B\pi \pi `$. For illustration, we write down the one-gluon exchange contribution to the $`\overline{B}_d^0\pi ^+\pi ^{}`$ matrix element of the operator $`Q_2^u=(\overline{u}_\alpha b_\beta )_{VA}(\overline{d}_\beta u_\alpha )_{VA}=(\overline{d}_\alpha b_\alpha )_{VA}(\overline{u}_\beta u_\beta )_{VA}`$.

$`Q_2^u_{(a)}`$ $`=`$ $`g_s^2{\displaystyle \frac{f_\pi }{4}}{\displaystyle \frac{C_F}{N}}{\displaystyle _0^1}𝑑u\varphi (u){\displaystyle \frac{d^4k}{(2\pi )^4}\frac{1}{k^2(uqk)^2[(pk)^2m_b^2]}}`$ (12)
$`\times \pi ^+|\overline{u}_i\gamma ^\mu (1\gamma _5)/q\gamma _5\gamma ^\alpha (u/q/k)\gamma _\mu (1\gamma _5)(/p/k+m_b)\gamma _\alpha b_i|\overline{B}_d^0,`$
$`Q_2^u_{(b)}`$ $`=`$ $`g_s^2{\displaystyle \frac{f_\pi }{4}}{\displaystyle \frac{C_F}{N}}{\displaystyle _0^1}𝑑u\varphi (u){\displaystyle \frac{d^4k}{(2\pi )^4}\frac{1}{k^2(\overline{u}qk)^2[(pk)^2m_b^2]}}`$ (14)
$`\times \pi ^+|\overline{u}_i\gamma ^\mu (1\gamma _5)(\overline{u}/q/k)\gamma ^\alpha /q\gamma _5\gamma _\mu (1\gamma _5)(/p/k+m_b)\gamma _\alpha b_i|\overline{B}_d^0,`$
$`Q_2^u_{(c)}`$ $`=`$ $`g_s^2{\displaystyle \frac{f_\pi }{4}}{\displaystyle \frac{C_F}{N}}{\displaystyle _0^1}𝑑u\varphi (u){\displaystyle \frac{d^4k}{(2\pi )^4}\frac{1}{k^2(uq+k)^2(pqk)^2}}`$ (16)
$`\times \pi ^+|\overline{u}_i\gamma _\alpha (/p/q/k)\gamma ^\mu (1\gamma _5)/q\gamma _5\gamma _\alpha (u/q+/k)\gamma _\mu (1\gamma _5)b_i|\overline{B}_d^0,`$
$`Q_2^u_{(d)}`$ $`=`$ $`g_s^2{\displaystyle \frac{f_\pi }{4}}{\displaystyle \frac{C_F}{N}}{\displaystyle _0^1}𝑑u\varphi (u){\displaystyle \frac{d^4k}{(2\pi )^4}\frac{1}{k^2(\overline{u}q+k)^2(pqk)^2}}`$ (18)
$`\times \pi ^+|\overline{u}_i\gamma _\alpha (/p/q/k)\gamma ^\mu (1\gamma _5)(\overline{u}/q+/k)\gamma _\alpha /q\gamma _5\gamma _\mu (1\gamma _5)b_i|\overline{B}_d^0,`$
When we calculate the vertex corrections in the leading power of $`\mathrm{\Lambda }_{QCD}/m_b`$, not only ultraviolet divergence emerges but infrared divergence does also. Infrared divergence arises from two regions where the virtuality of the loop $`k`$ is soft or collinear to the momentum of the pions. In Ref., the authors gave an explicit cancellation of soft and collinear divergence in vertex corrections for $`BD\pi `$ in eikonal approximation. Figures 1(a),1(b) and 1(c),1(d) cancel the soft divergence; 1(a),1(c) and 1(b),1(d) cancel the collinear divergence. For $`B\pi \pi `$, the cancellation is similar except that the collinear divergence also arises from the region where $`k`$ is collinear to the momentum of the recoiling pion. So Figs. 1(c),1(d) cancel not only part of soft divergence but also part of collinear divergence. Below, we give an explicit calculation of the Feynman diagrams Figs. 1(a)-1(d) to show the cancellation of the infrared divergences. In order to regularize the infrared divergence, there are two choices for us. One is the dimensional regularization (DR) scheme, in which the infrared divergence can be regularized into the pole terms $`1/(d4)`$. In contrast to the dimensional regularization of ultraviolet divergence, the infrared divergence arises when $`d4`$, instead of $`d4`$ in the case of the ultraviolet divergence. So the dimension $`d`$ in regularization for infrared divergence must be set to be greater than 4. This is a subtle point, but it will not cause any ambiguity in our calculation because the infrared part and ultraviolet part can be safely separated. The other method to regularize the infrared divergence is the well-known massive gluon (MG) scheme, in which the infrared divergence is handled by replacing $`1/k^2`$ by $`1/(k^2m_g^2)`$ in the gluon propagator. Similar scheme has been applied in earlier computation of the radiative corrections for $`\mu ^{}e^{}\overline{\nu }_e\nu _\mu `$, in which the massless photon is replaced by a massive photon. In addition, in our latter calculation, there are also several schemes in treating $`\gamma _5`$, the most popular two are the naive dimensional regularization (NDR) scheme and the ’t Hooft-Veltman renormalization (HV) scheme. Both have been applied to calculate the Wilson coefficients . In this work, if there is no specification, the NDR scheme is always applied in our calculations for its simplicity <sup>§</sup><sup>§</sup>§Such choice does cause a scheme dependence in the matrix elements. However, when we choose the Wilson coefficients in the same scheme as for the matrix elements, the final full decay amplitude is free of scheme-dependence..

After a straightforward calculation in DR scheme and using the corresponding Feynman parameter integrals listed in Appendix C, we obtain

$`Q_2^u_{(a)}`$ $`=`$ $`{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}\pi ^{}|\overline{d}_\alpha \gamma ^\mu (1\gamma _5)u_\alpha |0\pi ^+|\overline{u}_\beta \gamma _\mu (1\gamma _5)b_\beta |\overline{B}_d^0`$ (19)
$`\times `$ $`{\displaystyle _0^1}du\varphi (u)\{[{\displaystyle \frac{1}{ϵ}}\gamma _E+\mathrm{ln}4\pi +2\mathrm{ln}{\displaystyle \frac{\mu }{m_b}}+1+{\displaystyle \frac{u}{1u}}\mathrm{ln}u]`$ (20)
$``$ $`{\displaystyle \frac{\mathrm{\Gamma }(1a)}{(4\pi )^a}}\left({\displaystyle \frac{m_b}{\mu }}\right)^{2a}[{\displaystyle \frac{1}{a^2}}+{\displaystyle \frac{2(\mathrm{ln}u1)}{a}}+\mathrm{ln}^2u2\mathrm{L}\mathrm{i}_2(1{\displaystyle \frac{1}{u}})4\mathrm{ln}u+5+{\displaystyle \frac{2\mathrm{ln}u}{1u}}]\},`$ (21)
$`Q_2^u_{(b)}`$ $`=`$ $`{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}\pi ^{}|\overline{d}_\alpha \gamma ^\mu (1\gamma _5)u_\alpha |0\pi ^+|\overline{u}_\beta \gamma _\mu (1\gamma _5)b_\beta |\overline{B}_d^0`$ (22)
$`\times `$ $`{\displaystyle _0^1}du\varphi (u)\{4[{\displaystyle \frac{1}{ϵ}}\gamma _E+\mathrm{ln}4\pi +2\mathrm{ln}{\displaystyle \frac{\mu }{m_b}}+{\displaystyle \frac{11}{4}}+{\displaystyle \frac{\overline{u}}{1\overline{u}}}\mathrm{ln}\overline{u}]`$ (23)
$`+`$ $`{\displaystyle \frac{\mathrm{\Gamma }(1a)}{(4\pi )^a}}\left({\displaystyle \frac{m_b}{\mu }}\right)^{2a}[{\displaystyle \frac{1}{a^2}}+{\displaystyle \frac{2(\mathrm{ln}\overline{u}1)}{a}}+\mathrm{ln}^2\overline{u}2\mathrm{L}\mathrm{i}_2(1{\displaystyle \frac{1}{\overline{u}}})4\mathrm{ln}\overline{u}+6+{\displaystyle \frac{2\mathrm{ln}\overline{u}}{1\overline{u}}}]\},`$ (24)
$`Q_2^u_{(c)}`$ $`=`$ $`{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}\pi ^{}|\overline{d}_\alpha \gamma ^\mu (1\gamma _5)u_\alpha |0\pi ^+|\overline{u}_\beta \gamma _\mu (1\gamma _5)b_\beta |\overline{B}_d^0`$ (25)
$`\times `$ $`{\displaystyle _0^1}du\varphi (u)\{4[{\displaystyle \frac{1}{ϵ}}\gamma _E+\mathrm{ln}4\pi +2\mathrm{ln}{\displaystyle \frac{\mu }{m_b}}+{\displaystyle \frac{11}{4}}\mathrm{ln}(u)]`$ (26)
$`+`$ $`{\displaystyle \frac{\mathrm{\Gamma }(1a)}{(4\pi )^a}}\left({\displaystyle \frac{m_b}{\mu }}\right)^{2a}[{\displaystyle \frac{2}{a^2}}+{\displaystyle \frac{2(\mathrm{ln}(u)2)}{a}}+10{\displaystyle \frac{\pi ^2}{3}}4\mathrm{ln}(u)+\mathrm{ln}^2(u)]\},`$ (27)
$`Q_2^u_{(d)}`$ $`=`$ $`{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}\pi ^{}|\overline{d}_\alpha \gamma ^\mu (1\gamma _5)u_\alpha |0\pi ^+|\overline{u}_\beta \gamma _\mu (1\gamma _5)b_\beta |\overline{B}_d^0`$ (28)
$`\times `$ $`{\displaystyle _0^1}du\varphi (u)\{[{\displaystyle \frac{1}{ϵ}}\gamma _E+\mathrm{ln}4\pi +2\mathrm{ln}{\displaystyle \frac{\mu }{m_b}}+1\mathrm{ln}(\overline{u})]`$ (29)
$``$ $`{\displaystyle \frac{\mathrm{\Gamma }(1a)}{(4\pi )^a}}\left({\displaystyle \frac{m_b}{\mu }}\right)^{2a}[{\displaystyle \frac{2}{a^2}}+{\displaystyle \frac{2(\mathrm{ln}(\overline{u})2)}{a}}+10{\displaystyle \frac{\pi ^2}{3}}4\mathrm{ln}(\overline{u})+\mathrm{ln}^2(\overline{u})]\}.`$ (30)
In above, we have set $`d=4+2a`$ ($`a>0`$) in regularizing the infrared divergence. Then, after summing over all four diagrams, we find that all pole terms in $`1/a`$ are really cancelled before we integrate over the momentum fraction variable $`u`$. So after modified minimal subtraction ($`\overline{MS}`$), we get

$`\overline{Q_2^u}_{(a)+(b)+(c)+(d)}`$ (31)
$`=`$ $`{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}\pi ^{}|\overline{d}_\alpha \gamma ^\mu (1\gamma _5)u_\alpha |0\pi ^+|\overline{u}_\beta \gamma _\mu (1\gamma _5)b_\beta |\overline{B}_d^0`$ (32)
$`\times `$ $`{\displaystyle _0^1}du\varphi (u)[1812\mathrm{ln}{\displaystyle \frac{\mu }{m_b}}+{\displaystyle \frac{u}{1u}}\mathrm{ln}u{\displaystyle \frac{4\overline{u}}{1\overline{u}}}\mathrm{ln}\overline{u}+4\mathrm{ln}(u)\mathrm{ln}(\overline{u})`$ (33)
$``$ $`\left(\mathrm{ln}^2u\mathrm{ln}^2\overline{u}\right)+2\left(\mathrm{Li}_2(1{\displaystyle \frac{1}{u}})\mathrm{Li}_2(1{\displaystyle \frac{1}{\overline{u}}})\right)`$ (34)
$``$ $`({\displaystyle \frac{2\mathrm{ln}u}{1u}}{\displaystyle \frac{2\mathrm{ln}\overline{u}}{1\overline{u}}})+(\mathrm{ln}^2(u)\mathrm{ln}^2(\overline{u}))].`$ (35)
Assuming that the light-cone distribution amplitude $`\varphi (u)`$ is symmetric, then the above equation can be simplified as follows:

$`\overline{Q_2^u}_{(a)+(b)+(c)+(d)}`$ (36)
$`=`$ $`{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}\pi ^{}|\overline{d}_\alpha \gamma ^\mu (1\gamma _5)u_\alpha |0\pi ^+|\overline{u}_\beta \gamma _\mu (1\gamma _5)b_\beta |\overline{B}_d^0`$ (37)
$`\times `$ $`{\displaystyle _0^1}𝑑u\varphi (u)\left[1812\mathrm{ln}{\displaystyle \frac{\mu }{m_b}}+3{\displaystyle \frac{12u}{1u}}\mathrm{ln}u3i\pi \right].`$ (38)
It is easy to check that the above equations are consistent with the results in previous works. Actually, with the MG scheme, we get the same results as that by using the DR scheme.

With Eqs.(35,38), we can compute the vertex corrections no matter the LCDA $`\varphi (u)`$ is symmetric or asymmetric. This is very important in principle. For instance, when kaon is ejected from b quark decay, we must take the contributions from the asymmetric part of LCDA of kaon into account, although the contributions from the asymmetric part are very small numerically.

#### 2 Penguin Corrections

There are two kinds of penguin corrections. One is the four quark operators insertion \[Fig. 1(e)\]; the other is the magnetic penguin insertion \[Fig. 1(f)\]. The first kind is generally called BSS mechanism. In generalized factorization, BSS mechanism plays a very important role in $`CP`$ violation because it is the unique source of strong phases. But in generalized factorization, the virtuality of the gluon or photon is ambiguous; usually one varies $`k^2`$ around $`m_B^2/2`$. This variation does not have an important effect on the branching ratios, but it does for $`CP`$ asymmetries. In QCD factorization, this ambiguity is rendered by taking the virtuality of the gluon as $`k^2=(puq)^2=\overline{u}m_b^2/2`$. When treating penguin contractions, one should be careful that Fig. 1(e) contains two kinds of topology, which is depicted in Fig. 2. They are equivalent in 4 dimensions according to Fierz rearrangement. However, since penguin corrections contain ultraviolet divergence, we must do calculations in d dimensions where these two kinds of topology are not equivalent . Below we give an explicit calculation of $`Q_4`$ or $`Q_6`$ penguin insertions for $`\overline{B}_d^0\pi ^+\pi ^{}`$ which belong to the second topology, Fig. 2(b):

$`Q_{4,6}_{(e)}^{\mathrm{twist}2}`$ $`=`$ $`{\displaystyle \frac{f_\pi }{4}}g_s^2\mu ^{2ϵ}{\displaystyle \frac{C_F}{N}}{\displaystyle _0^1}𝑑u\varphi (u)\pi ^+|\overline{u}_i\gamma _\alpha /q\gamma _5\gamma _\mu (1\gamma _5)b_i|\overline{B}_d^0`$ (39)
$`\times `$ $`{\displaystyle \underset{q}{}}{\displaystyle \frac{d^dk}{(2\pi )^d}\frac{\mathrm{Tr}[(/l/km_q)\gamma ^\alpha (/k+m_q)\gamma ^\mu (1\gamma _5)]}{[(lk)^2m_q^2][k^2m_q^2]l^2}}|_{l=puq}`$ (40)
$`=`$ $`2if_\pi {\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}{\displaystyle _0^1}𝑑u\varphi (u)\pi ^+|\overline{u}_i\gamma _\alpha /q\gamma _5\gamma _\mu (1\gamma _5)b_i|\overline{B}_d^0\left[{\displaystyle \frac{l_\alpha l_\mu }{l^2}}g_{\alpha \mu }\right]`$ (41)
$`\times `$ $`{\displaystyle \underset{q}{}}\left[{\displaystyle \frac{1}{6}}({\displaystyle \frac{1}{ϵ}}\gamma _E+\mathrm{ln}4\pi )+{\displaystyle _0^1}𝑑tt(1t)\mathrm{ln}{\displaystyle \frac{\mu ^2}{m_q^2t(1t)l^2iϵ}}\right]|_{l=puq}.`$ (42)
After $`\overline{MS}`$ subtraction and using the equations of motions, we get the finite result

$`\overline{Q_{4,6}}_{(e)}^{\mathrm{twist}2}`$ $`=`$ $`{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}\pi ^{}|\overline{d}_i\gamma ^\mu (1\gamma _5)u_i|0\pi ^+|\overline{u}_j\gamma _\mu (1\gamma _5)b_j|\overline{B}_d^0`$ (43)
$`\times `$ $`{\displaystyle \underset{q}{}}\left[{\displaystyle \frac{4}{3}}\mathrm{ln}{\displaystyle \frac{\mu }{m_b}}4{\displaystyle _0^1}𝑑u\varphi (u){\displaystyle _0^1}𝑑tt(1t)\mathrm{ln}(s_qt(1t)\overline{u}iϵ)\right],`$ (44)
where $`s_q=m_q^2/m_b^2`$. The first topology, Fig. 2(a), for example, $`Q_1^c`$ penguin insertion for $`\overline{B}_d^0\pi ^+\pi ^{}`$, is similar to the results of the second topology, Fig.2(b), except that there is an extra factor $`2/3`$:

$`\overline{Q_1^c}_{(e)}^{\mathrm{twist}2}`$ $`=`$ $`{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}\pi ^{}|\overline{d}_i\gamma ^\mu (1\gamma _5)u_i|0\pi ^+|\overline{u}_j\gamma _\mu (1\gamma _5)b_j|\overline{B}_d^0`$ (45)
$`\times `$ $`\left[{\displaystyle \frac{2}{3}}+{\displaystyle \frac{4}{3}}\mathrm{ln}{\displaystyle \frac{\mu }{m_b}}4{\displaystyle _0^1}𝑑u\varphi (u){\displaystyle _0^1}𝑑tt(1t)\mathrm{ln}(s_ct(1t)\overline{u}iϵ)\right].`$ (46)
For the magnetic penguin insertion, it is the easiest calculation of the radiative corrections. The result of $`Q_{8G}`$ insertion for $`\overline{B}_d^0\pi ^+\pi ^{}`$ is

$`Q_{8G}_{(f)}^{\mathrm{twist}2}`$ $`=`$ $`{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}f_\pi m_b{\displaystyle _0^1}𝑑u\varphi (u){\displaystyle \frac{1}{k^2}}\pi ^+|\overline{u}_i\gamma ^\alpha /q\gamma _5\sigma _{\beta \alpha }k^\beta (1+\gamma _5)b_i|\overline{B}_d^0|_{k=puq}`$ (47)
$`=`$ $`{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}{\displaystyle _0^1}𝑑u{\displaystyle \frac{2\varphi (u)}{\overline{u}}}\pi ^{}|\overline{d}_i\gamma ^\mu (1\gamma _5)|0\pi ^+|\overline{u}_i\gamma _\mu (1\gamma _5)|\overline{B}_d^0`$ (48)
#### 3 Hard scattering with the spectator

Hard spectator scattering \[Fig.1(g) and (h)\] is completely missing in the naive factorization. But in QCD factorization, it can be calculated in the perturbative QCD, and expressed by a convolution of the hard kernel $`T^{II}`$ and the LCDAs of mesons. At the leading power of $`\mathrm{\Lambda }_{QCD}/m_b`$, both of the light pseudoscalars from the $`B`$ meson decay can be represented by their leading twist LCDAs. So after a straightforward calculation, we obtain this contribution for $`\overline{B}_d^0\pi ^+\pi ^{}`$ from the operator $`Q_2^u`$ insertion,

$`Q_2^u_{(g)+(h)}^{\mathrm{twist}2}`$ $`=`$ $`{\displaystyle \frac{if_Bf_\pi ^2}{64}}{\displaystyle \frac{C_F}{N^2}}g_s^2{\displaystyle _0^1}𝑑\xi 𝑑u𝑑v\varphi _B(\xi )\varphi (u)\varphi (v)`$ (49)
$`\times `$ $`\{{\displaystyle \frac{\mathrm{Tr}[(/p+m_B)\gamma _5\gamma ^\alpha /q_1\gamma _5\gamma _\rho (1\gamma _5)/q_2\gamma _5\gamma _\alpha /l_d\gamma _\rho (1\gamma _5)]}{k^2l_d^2}}|_{l_d=uq_2k}^{k=\xi p\overline{v}q_1}`$ (51)
$`+{\displaystyle \frac{\mathrm{Tr}[(/p+m_B)\gamma _5\gamma ^\alpha /q_1\gamma _5\gamma _\rho (1\gamma _5)/l_u\gamma _\alpha /q_2\gamma _5\gamma _\rho (1\gamma _5)]}{k^2l_u^2}}|_{l_u=k\overline{u}q_2}^{k=\xi p\overline{v}q_1}\}`$
$`=`$ $`i\pi \alpha _sf_Bf_\pi ^2{\displaystyle \frac{C_F}{N^2}}{\displaystyle _0^1}𝑑\xi {\displaystyle \frac{\varphi _B(\xi )}{\xi }}{\displaystyle _0^1}𝑑u{\displaystyle \frac{\varphi (u)}{u}}{\displaystyle _0^1}𝑑v{\displaystyle \frac{\varphi (v)}{\overline{v}}}.`$ (52)
### B Chirally enhanced corrections — Twist-3 LCDAs insertion

It has been observed that QCD factorization is demonstrated only in the strict heavy quark limit. This means that any generalization of QCD factorization to include or partly include power corrections of $`\mathrm{\Lambda }_{QCD}/m_b`$ should redemonstrate that factorization still holds. There are a variety of sources which may contribute to power corrections in $`1/m_b`$; examples are higher twist distribution amplitudes, transverse momenta of quarks in the light meson, annihilation diagrams, etc. Unfortunately, there is no known systematic way to evaluate these power corrections for exclusive decays. Moreover, factorization might break down when these power corrections, for instance, transverse momenta effects, are considered. This indicates that one might have to give up such an ambitious plan that all power corrections could be, at least in principle, incorporated into QCD factorization order by order. One might argue that power corrections in $`B`$ decays are numerically unimportant because these corrections are expanded in order of a small number $`\mathrm{\Lambda }_{QCD}/m_b1/15`$. But this is not true. For instance, the contributions of operator $`Q_6`$ to decay amplitudes would formally vanish in the strict heavy quark limit. However it is numerically very important in penguin-dominated $`B`$ rare decays, such as interesting channels $`B\pi K`$, etc. This is because $`Q_6`$ is always multiplied by a formally power suppressed but chirally enhanced factor $`r_\chi =2m_P^2/m_b(m_1+m_2)𝒪(1)`$, where $`m_1`$ and $`m_2`$ are current quark masses. So power suppression might probably fail at least in this case. Therefore phenomenological applicability of QCD factorization in $`B`$ rare decays requires at least a consistent inclusion of chirally enhanced corrections.

Chirally enhanced corrections arise from the two particle twist-3 light-cone distribution amplitudes, generally called $`\varphi _p(x)`$ and $`\varphi _\sigma (x)`$. So when chirally enhanced corrections are concerned, the final light mesons should be described by leading twist and twist-3 distribution amplitudes. Then it is crucial to show that factorization really holds when considering twist-3 distribution amplitudes. The most difficult part is to demonstrate the infrared finiteness of the hard scattering kernels $`T_i^I`$. In addition, possible chirally enhanced power corrections can also appear in the hard spectator scattering. So, for consistency, we must involve these corrections.

#### 1 Vertex corrections

When we calculate the chirally enhanced power corrections, contrast to the leading-twist light-cone wave function insertion, the cancellation of the infrared divergences in the vertex corrections to $`(VA)(V+A)`$ operator (here it is $`Q_5`$ or $`Q_7`$) can not be shown simply by the eikonal approximation similar to what has been done at the leading power of $`\mathrm{\Lambda }_{QCD}/m_b`$, because the Dirac structure or spin structure of twist-3 light-cone wave functions of the light pseudoscalar makes the “on-shell” condition for the external quarks invalid. Thus, to justify the cancellation of the infrared divergence in $`(VA)(V+A)`$ vertex corrections, we must give the explicit calculation. As mentioned in the previous subsections, we have two choices to regularize the infrared divergence in one-loop calculation. One is the DR scheme; the other is MG scheme. Generally, these two schemes are equivalent, for instance, similar to what has been done in $`(VA)(VA)`$ vertex corrections. However, in the DR scheme, it is difficult to extrapolate the twist-3 wave functions of the light pseudoscalar to d dimensions properly, although they are well-defined in 4 dimensions. Therefore, we prefer to use the MG scheme in our calculation for chirally enhanced corrections to avoid the above possible problems.

In addition, generally we calculate the Feynman diagrams in the momentum space, so the correct projection of the light-cone wave functions of the light pseudoscalar in the momentum space is necessary. From Appendix B, we find that it is easy to find the proper momentum space projection of the leading twist and $`\varphi _p`$ type twist-3 wave function, but for $`\varphi _\sigma `$, the projection is not very clear. Note that the coordinate $`x^\mu `$ in the definition of $`\varphi _\sigma `$ by the non-local matrix element must be transformed into a partial derivative of a certain momentum in the projection of momentum space. However, it is difficult to find the derivative which makes the projection only depend on the structure of the light pseudoscalar itself. Generally, the momentum which the partial derivative acts on is dependent on the hard kernel. Therefore, we prefer to compute the Feynman diagrams of the twist-3 wave functions insertion, especially $`\varphi _\sigma `$ insertion in the coordinate space. We think that such calculation can avoid the ambiguity about how to project the coordinate $`x^\mu `$ into the momentum space. We recalculate the leading twist insertion by using the same method, and obtain the same results as those in the previous sections. Below, we will show how to perform this trick in calculation of $`\varphi _\sigma `$ insertion. For example, let us consider Fig. 3. In coordinate space, we have

$`\mathrm{Fig}.3`$ $`=`$ $`{\displaystyle \frac{f_P\mu _P}{4}}{\displaystyle \frac{C_F}{N}}g_s^2{\displaystyle 𝑑u\frac{\varphi _\sigma (u)}{6}d^4x_1d^4x_2\frac{d^4k}{(2\pi )^4}\frac{d^4l_{\overline{u}}}{(2\pi )^4}\frac{d^4l_b}{(2\pi )^4}\frac{e^{i(\overline{u}qk+l_{\overline{u}})x_2}e^{i(k+l_bp)x_1}}{(k^2m_g^2)l_{\overline{u}}^2(l_b^2m_b^2)}}`$ (53)
$`\times `$ $`\pi ^+|\overline{u}_i\gamma ^\rho (1+\gamma _5)/l_{\overline{u}}\gamma ^\alpha \gamma _5\sigma _{\mu \nu }q^\mu x_2^\nu \gamma _\rho (1\gamma _5)(/l_b+m_b)\gamma _\alpha b_i|\overline{B}_d^0`$ (54)
$`=`$ $`{\displaystyle \frac{f_P\mu _P}{4}}{\displaystyle \frac{C_F}{N}}g_s^2{\displaystyle 𝑑u\frac{\varphi _\sigma (u)}{6}\frac{d^4k}{(2\pi )^4}\frac{d^4l_{\overline{u}}d^4x_2}{(2\pi )^4}\frac{e^{i(\overline{u}qk+l_{\overline{u}})x_2}q^\mu x_2^\nu }{(k^2m_g^2)l_{\overline{u}}^2(l_b^2m_b^2)}}`$ (55)
$`\times `$ $`\pi ^+|\overline{u}_i\gamma ^\rho (1+\gamma _5)/l_{\overline{u}}\gamma ^\alpha \gamma _5\sigma _{\mu \nu }\gamma _\rho (1\gamma _5)(/l_b+m_b)\gamma _\alpha b_i|\overline{B}_d^0|_{l_b=pk}`$ (56)
$`=`$ $`{\displaystyle \frac{if_P\mu _P}{4}}{\displaystyle \frac{C_F}{N}}g_s^2{\displaystyle 𝑑u\frac{\varphi _\sigma (u)}{6}\frac{d^4k}{(2\pi )^4}\frac{q^\mu }{(k^2m_g^2)(l_b^2m_b^2)}}`$ (57)
$`\times `$ $`{\displaystyle \frac{}{l_{\overline{u}\nu }}}\left\{\pi ^+|\overline{u}_i\gamma ^\rho (1+\gamma _5){\displaystyle \frac{/l_{\overline{u}}}{l_{\overline{u}}^2}}\gamma ^\alpha \gamma _5\sigma _{\mu \nu }\gamma _\rho (1\gamma _5)(/l_b+m_b)\gamma _\alpha b_i|\overline{B}_d^0\right\}|_{l_b=pk}^{l_{\overline{u}}=k\overline{u}q}.`$ (58)
The above trick has been applied in the calculation of the proper correlation function to extract the transition form factor $`F^{B\pi }`$ within the frame of the light-cone sum rule. So within the same frame, we obtain

$`Q_5_{(a)}`$ $`=`$ $`g_s^2{\displaystyle \frac{f_\pi \mu _\pi }{4}}{\displaystyle \frac{C_F}{N}}{\displaystyle _0^1}𝑑u{\displaystyle \frac{d^4k}{(2\pi )^4}\frac{1}{(k^2m_g^2)l_d^2(l_b^2m_b^2)}}`$ (62)
$`\times \{\varphi _p(u)\pi ^+|\overline{u}_i\gamma ^\mu (1+\gamma _5)\gamma _5\gamma ^\alpha /l_d\gamma _\mu (1\gamma _5)(/l_b+m_b)\gamma _\alpha b_i|\overline{B}_d^0`$
$`+i{\displaystyle \frac{\varphi _\sigma (u)}{6}}\pi ^+|\overline{u}_i\gamma ^\rho (1+\gamma _5)\gamma _5\sigma _{\mu \nu }q^\mu \gamma ^\alpha \left(\gamma ^\nu {\displaystyle \frac{2l_d^\nu /l_d}{l_d^2}}\right)`$
$`\gamma _\rho (1\gamma _5)(/l_b+m_b)\gamma _\alpha b_i|\overline{B}_d^0\left\}\right|^{l_d=uqk}_{l_b=pk},`$
$`Q_5_{(b)}`$ $`=`$ $`g_s^2{\displaystyle \frac{f_\pi \mu _\pi }{4}}{\displaystyle \frac{C_F}{N}}{\displaystyle _0^1}𝑑u{\displaystyle \frac{d^4k}{(2\pi )^4}\frac{1}{(k^2m_g^2)l_{\overline{u}}^2(l_b^2m_b^2)}}`$ (66)
$`\times \{\varphi _p(u)\pi ^+|\overline{u}_i\gamma ^\mu (1+\gamma _5)/l_{\overline{u}}\gamma ^\alpha \gamma _5\gamma _\mu (1\gamma _5)(/l_b+m_b)\gamma _\alpha b_i|\overline{B}_d^0`$
$`+i{\displaystyle \frac{\varphi _\sigma (u)}{6}}\pi ^+|\overline{u}_i\gamma ^\rho (1+\gamma _5)\left(\gamma ^\nu {\displaystyle \frac{2l_{\overline{u}}^\nu /l_{\overline{u}}}{l_{\overline{u}}^2}}\right)\gamma ^\alpha \gamma _5\sigma _{\mu \nu }q^\mu `$
$`\gamma _\rho (1\gamma _5)(/l_b+m_b)\gamma _\alpha b_i|\overline{B}_d^0\left\}\right|^{l_{\overline{u}}=k\overline{u}q}_{l_b=pk},`$
$`Q_5_{(c)}`$ $`=`$ $`g_s^2{\displaystyle \frac{f_\pi \mu _\pi }{4}}{\displaystyle \frac{C_F}{N}}{\displaystyle _0^1}𝑑u{\displaystyle \frac{d^4k}{(2\pi )^4}\frac{1}{(k^2m_g^2)l_d^2l_u^2}}`$ (70)
$`\times \{\varphi _p(u)\pi ^+|\overline{u}_i\gamma _\alpha /l_u\gamma ^\mu (1+\gamma _5)\gamma _5\gamma _\alpha /l_d\gamma _\mu (1\gamma _5)b_i|\overline{B}_d^0`$
$`+i{\displaystyle \frac{\varphi _\sigma (u)}{6}}\pi ^+|\overline{u}_i\gamma ^\alpha /l_u\gamma ^\mu (1+\gamma _5)\gamma _5\sigma _{\mu \nu }q^\mu \gamma _\alpha `$
$`(\gamma ^\nu {\displaystyle \frac{2l_d^\nu /l_d}{l_d^2}})\gamma _\rho (1\gamma _5)b_i|\overline{B}_d^0\left\}\right|^{l_d=uqk}_{l_u=pq+k},`$
$`Q_5_{(d)}`$ $`=`$ $`g_s^2{\displaystyle \frac{f_\pi \mu _\pi }{4}}{\displaystyle \frac{C_F}{N}}{\displaystyle _0^1}𝑑u{\displaystyle \frac{d^4k}{(2\pi )^4}\frac{1}{(k^2m_g^2)l_{\overline{u}}^2l_u^2}}`$ (74)
$`\times \{\varphi _p(u)\pi ^+|\overline{u}_i\gamma ^\alpha /l_u\gamma ^\mu (1+\gamma _5)/l_{\overline{u}}\gamma _\alpha \gamma _5\gamma _\mu (1\gamma _5)b_i|\overline{B}_d^0`$
$`+i{\displaystyle \frac{\varphi _\sigma (u)}{6}}\pi ^+|\overline{u}_i\gamma ^\alpha /l_u\gamma ^\rho (1+\gamma _5)\left(\gamma ^\nu {\displaystyle \frac{2l_{\overline{u}}^\nu /l_{\overline{u}}}{l_{\overline{u}}^2}}\right)`$
$`\gamma _\alpha \gamma _5\sigma _{\mu \nu }q^\mu \gamma _\rho (1\gamma _5)b_i|\overline{B}_d^0\left\}\right|^{l_{\overline{u}}=k\overline{u}q}_{l_u=pq+k}.`$
Perform the one-loop integrations:

$`Q_5_{(a)}`$ $`=`$ $`2{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}\pi ^{}|\overline{d}_i(1+\gamma _5)u_i|0\pi ^+|\overline{u}_j(1\gamma _5)b_j|\overline{B}_d^0`$ (75)
$`\times `$ $`{\displaystyle _0^1}du\{\varphi _p(u)[({\displaystyle \frac{1}{ϵ}}\gamma _E+\mathrm{ln}4\pi +2\mathrm{ln}{\displaystyle \frac{\mu }{m_b}})`$ (76)
$`+`$ $`{\displaystyle \frac{1}{4}}\mathrm{ln}^2\lambda +\mathrm{ln}(u)\mathrm{ln}\lambda 2\mathrm{ln}u\mathrm{ln}\lambda +{\displaystyle \frac{1}{2}}\mathrm{ln}^2u\mathrm{Li}_2(1{\displaystyle \frac{1}{u}})+{\displaystyle \frac{1}{2}}+{\displaystyle \frac{5}{4}}\pi ^2]`$ (77)
$`+`$ $`{\displaystyle \frac{\varphi _\sigma (u)}{6u}}[{\displaystyle \frac{1}{2}}\mathrm{ln}^2\lambda +2\mathrm{ln}(u)\mathrm{ln}\lambda 4\mathrm{ln}u\mathrm{ln}\lambda +\mathrm{ln}\lambda `$ (78)
$`+`$ $`4\mathrm{ln}^2u\mathrm{ln}^2(u)2\mathrm{ln}u\mathrm{ln}(1u)+2\mathrm{L}\mathrm{i}_2({\displaystyle \frac{1}{u}})+{\displaystyle \frac{7}{6}}\pi ^2\mathrm{ln}u{\displaystyle \frac{\mathrm{ln}u}{1u}}]\},`$ (79)
$`Q_5_{(b)}`$ $`=`$ $`2{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}\pi ^{}|\overline{d}_i(1+\gamma _5)u_i|0\pi ^+|\overline{u}_j(1\gamma _5)b_j|\overline{B}_d^0`$ (80)
$`\times `$ $`{\displaystyle _0^1}du\{\varphi _p(u)[({\displaystyle \frac{1}{ϵ}}\gamma _E+\mathrm{ln}4\pi +2\mathrm{ln}{\displaystyle \frac{\mu }{m_b}})`$ (81)
$`+`$ $`{\displaystyle \frac{1}{4}}\mathrm{ln}^2\lambda +\mathrm{ln}(\overline{u})\mathrm{ln}\lambda 2\mathrm{ln}\overline{u}\mathrm{ln}\lambda +{\displaystyle \frac{1}{2}}\mathrm{ln}^2\overline{u}\mathrm{Li}_2(1{\displaystyle \frac{1}{\overline{u}}}){\displaystyle \frac{5}{2}}+{\displaystyle \frac{5}{4}}\pi ^2]`$ (82)
$`+`$ $`{\displaystyle \frac{\varphi _\sigma (u)}{6\overline{u}}}[{\displaystyle \frac{1}{2}}\mathrm{ln}^2\lambda +2\mathrm{ln}(\overline{u})\mathrm{ln}\lambda 4\mathrm{ln}\overline{u}\mathrm{ln}\lambda +\mathrm{ln}\lambda `$ (83)
$`+`$ $`4\mathrm{ln}^2\overline{u}\mathrm{ln}^2(\overline{u})2\mathrm{ln}\overline{u}\mathrm{ln}(1\overline{u})+2\mathrm{L}\mathrm{i}_2({\displaystyle \frac{1}{\overline{u}}})+{\displaystyle \frac{7}{6}}\pi ^2\mathrm{ln}\overline{u}{\displaystyle \frac{\mathrm{ln}\overline{u}}{1\overline{u}}}]\},`$ (84)
$`Q_5_{(c)}`$ $`=`$ $`2{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}\pi ^{}|\overline{d}_i(1+\gamma _5)u_i|0\pi ^+|\overline{u}_j(1\gamma _5)b_j|\overline{B}_d^0`$ (85)
$`\times `$ $`{\displaystyle _0^1}du\{\varphi _p(u)[({\displaystyle \frac{1}{ϵ}}\gamma _E+\mathrm{ln}4\pi +2\mathrm{ln}{\displaystyle \frac{\mu }{m_b}})`$ (86)
$`+`$ $`{\displaystyle \frac{1}{2}}\mathrm{ln}^2\lambda [\mathrm{ln}(u)2]\mathrm{ln}\lambda +{\displaystyle \frac{1}{2}}\mathrm{ln}^2(u)\mathrm{ln}(u){\displaystyle \frac{3}{2}}+{\displaystyle \frac{\pi ^2}{3}}]`$ (87)
$`+`$ $`{\displaystyle \frac{\varphi _\sigma (u)}{6u}}[\mathrm{ln}^2\lambda [2\mathrm{ln}(u)3]\mathrm{ln}\lambda `$ (88)
$`+`$ $`2\mathrm{ln}(u)\mathrm{ln}u\mathrm{ln}^2u3\mathrm{ln}(u)+3{\displaystyle \frac{\pi ^2}{3}}]\},`$ (89)
$`Q_5_{(d)}`$ $`=`$ $`2{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}\pi ^{}|\overline{d}_i(1+\gamma _5)u_i|0\pi ^+|\overline{u}_j(1\gamma _5)b_j|\overline{B}_d^0`$ (90)
$`\times `$ $`{\displaystyle _0^1}du\{\varphi _p(u)[({\displaystyle \frac{1}{ϵ}}\gamma _E+\mathrm{ln}4\pi +2\mathrm{ln}{\displaystyle \frac{\mu }{m_b}})`$ (91)
$`+`$ $`{\displaystyle \frac{1}{2}}\mathrm{ln}^2\lambda [\mathrm{ln}(\overline{u})2]\mathrm{ln}\lambda +{\displaystyle \frac{1}{2}}\mathrm{ln}^2(\overline{u})\mathrm{ln}(\overline{u})+{\displaystyle \frac{3}{2}}+{\displaystyle \frac{\pi ^2}{3}}]`$ (92)
$`+`$ $`{\displaystyle \frac{\varphi _\sigma (u)}{6\overline{u}}}[\mathrm{ln}^2\lambda [2\mathrm{ln}(\overline{u})3]\mathrm{ln}\lambda `$ (93)
$`+`$ $`2\mathrm{ln}(\overline{u})\mathrm{ln}\overline{u}\mathrm{ln}^2\overline{u}3\mathrm{ln}(\overline{u})+3{\displaystyle \frac{\pi ^2}{3}}]\}.`$ (94)
Here $`\lambda =m_g^2/m_b^2`$. From the above equations, it is observed that, in the case of $`\varphi _\sigma `$ distribution amplitudes, the terms with infrared divergence in vertex correction diagrams cannot cancel unless $`\varphi _\sigma (u)`$ is a symmetric function: $`\varphi _\sigma (u)=\varphi _\sigma (\overline{u})`$. This is an unexpected result, which means QCD factorization is violated for asymmetric twist-3 light-cone distribution amplitudes. This indicates that chirally enhanced corrections can be included consistently in the framework of QCD factorization only when twist-3 light-cone distribution amplitudes are symmetric. Therefore, in the following, we will implicitly assume a symmetric twist-3 light-cone distribution amplitude for light pseudoscalar mesons. It is then straightforward to show that vertex corrections of $`(VA)(V+A)`$ operator are completely cancelled after summing over four diagrams in the case of $`\varphi _\sigma `$ distribution amplitude. The final result of $`(VA)(V+A)`$ vertex corrections, in the condition that the twist-3 LCDA is symmetric, is

$$Q_5_{(a)+(b)+(c)+(d)}=12\frac{\alpha _s}{4\pi }\frac{C_F}{N}\pi ^{}|\overline{d}_i(1+\gamma _5)u_i|0\pi ^+|\overline{u}_j(1\gamma _5)b_j|\overline{B}_d^0$$
(95)
#### 2 Penguin corrections

In quark level, usually one decomposes the basic QCD vertex $`ig_s\gamma ^\mu T_{ij}^a`$ in penguin insertion into the two chiral current couplings $`i\frac{1}{2}g_s\gamma ^\mu T_{ij}^a(1+\gamma _5)`$ and $`i\frac{1}{2}g_s\gamma ^\mu T_{ij}^a(1\gamma _5)`$; then the penguin insertions contribute the same magnitude to the $`(VA)(VA)`$ and $`(VA)(V+A)`$ vertex. But in hadron level, this point of view must be examined in elaborate calculation.

For illustration, we give the results of $`Q_4`$ or $`Q_6`$ penguin corrections to $`\overline{B}_d^0\pi ^+\pi ^{}`$, which belong to the second penguin topology Fig. 2(b), when $`\varphi _p(u)`$ is inserted:

$`Q_{4,6}_{(e)}^{\varphi _p}`$ $`=`$ $`{\displaystyle \frac{f_\pi \mu _\pi }{4}}g_s^2\mu ^{2ϵ}{\displaystyle \frac{C_F}{N}}{\displaystyle _0^1}𝑑u\varphi _p(u)\pi ^+|\overline{u}_i\gamma _\alpha \gamma _5\gamma _\mu (1\gamma _5)b_i|\overline{B}_d^0`$ (96)
$`\times `$ $`{\displaystyle \underset{q}{}}{\displaystyle \frac{d^dk}{(2\pi )^d}\frac{\mathrm{Tr}[(/pu/q/km_q)\gamma ^\alpha (/k+m_q)\gamma ^\mu (1\gamma _5)]}{[(puqk)^2m_q^2][k^2m_q^2](puq)^2}}`$ (97)
$`=`$ $`2if_\pi \mu _\pi {\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}{\displaystyle _0^1}𝑑u\varphi _p(u)\pi ^+|\overline{u}_i\gamma _\alpha \gamma _5\gamma _\mu (1\gamma _5)b_i|\overline{B}_d^0\left[{\displaystyle \frac{l_\alpha l_\mu }{l^2}}g_{\alpha \mu }\right]`$ (98)
$`\times `$ $`{\displaystyle \underset{q}{}}\left[{\displaystyle \frac{1}{6}}({\displaystyle \frac{1}{ϵ}}\gamma _E+\mathrm{ln}4\pi )+{\displaystyle _0^1}𝑑tt(1t)\mathrm{ln}{\displaystyle \frac{\mu ^2}{m_q^2t(1t)l^2iϵ}}\right]|_{l=puq}`$ (99)
After $`\overline{MS}`$ subtraction, we obtain

$`\overline{Q_{4,6}}_{(e)}^{\varphi _p}`$ $`=`$ $`2{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}\pi ^{}|\overline{d}_i(1+\gamma _5)u_i|0\pi ^+|\overline{u}_j(1\gamma _5)b_j|\overline{B}_d^0`$ (101)
$`\times `$ $`{\displaystyle \underset{q}{}}\left[\mathrm{ln}{\displaystyle \frac{\mu }{m_b}}3{\displaystyle _0^1}𝑑u\varphi _p(u){\displaystyle _0^1}𝑑tt(1t)\mathrm{ln}(s_qt(1t)\overline{u}iϵ)\right].`$ (102)
For the first penguin insertion topology, Fig. 2(a), the result is

$`\overline{Q_1^c}_{(e)}^{\varphi _p}`$ $`=`$ $`2{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}\pi ^{}|\overline{d}_i(1+\gamma _5)u_i|0\pi ^+|\overline{u}_j(1\gamma _5)b_j|\overline{B}_d^0`$ (103)
$`\times `$ $`\left[\mathrm{ln}{\displaystyle \frac{\mu }{m_b}}{\displaystyle \frac{1}{2}}3{\displaystyle _0^1}𝑑u\varphi _p(u){\displaystyle _0^1}𝑑tt(1t)\mathrm{ln}(s_ct(1t)\overline{u}iϵ)\right].`$ (104)
Similarly, when $`\varphi _\sigma (u)`$ is inserted, by using the method in the previous subsection, we have

$`Q_{4,6}_{(e)}^{\varphi _\sigma }`$ $`=`$ $`{\displaystyle \frac{if_\pi \mu _\pi }{4}}g_s^2\mu ^{2ϵ}{\displaystyle \frac{C_F}{N}}{\displaystyle _0^1}𝑑u{\displaystyle \frac{\varphi _\sigma (u)}{6}}\pi ^+|\overline{u}_i(1\gamma _5)\gamma _\alpha \sigma _{\mu \nu }q^\mu \gamma _\rho b_i|\overline{B}_d^0`$ (105)
$`\times `$ $`{\displaystyle \underset{q}{}}\left[{\displaystyle \frac{}{k_\nu }}{\displaystyle \frac{d^dl_q}{(2\pi )^d}\frac{\mathrm{Tr}[(/l_q/k+m_q)\gamma ^\alpha (/l_q+m_q)\gamma ^\rho (1\gamma _5)]}{[(l_qk)^2m_q^2][l_q^2m_q^2]k^2}}\right]|_{k=puq}.`$ (106)
After integration and subtraction,

$`\overline{Q_{4,6}}_{(e)}^{\varphi _\sigma }`$ $`=`$ $`2{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}\pi ^{}|\overline{d}_i(1+\gamma _5)u_i|0\pi ^+|\overline{u}_j(1\gamma _5)b_j|\overline{B}_d^0`$ (107)
$`\times `$ $`{\displaystyle \underset{q}{}}{\displaystyle _0^1}du{\displaystyle \frac{\varphi _\sigma }{6\overline{u}}}[{\displaystyle \frac{2}{3}}\mathrm{ln}{\displaystyle \frac{\mu }{m_b}}`$ (108)
$``$ $`{\displaystyle _0^1}dt(2t(1t)\mathrm{ln}(s_qt(1t)\overline{u}iϵ)+{\displaystyle \frac{t^2(1t)^2\overline{u}}{s_qt(1t)\overline{u}iϵ}})].`$ (109)
The magnetic penguin insertion is easier; we write the result of $`Q_{8G}`$ insertion as follows:

$`Q_{8G}_{(f)}^{\mathrm{twist}3}`$ $`=`$ $`{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}f_\pi \mu _\pi m_b`$ (110)
$`\times `$ $`{\displaystyle _0^1}du\{\varphi _p(u){\displaystyle \frac{1}{k^2}}\pi ^+|\overline{u}_i\gamma ^\alpha \gamma _5\sigma _{\beta \alpha }k^\beta (1+\gamma _5)b_i|\overline{B}_d^0`$ (111)
$`+`$ $`i{\displaystyle \frac{\varphi _\sigma (u)}{6}}{\displaystyle \frac{1}{4k^2}}[\pi ^+|\overline{u}_i\gamma ^\alpha \gamma _5[/q,\gamma _\sigma ][\gamma ^\sigma ,\gamma _\alpha ](1+\gamma _5)b_i|\overline{B}_d^0`$ (112)
$``$ $`{\displaystyle \frac{2}{k^2}}\pi ^+|\overline{u}_i\gamma ^\alpha \gamma _5[/q,/k][/k,\gamma _\alpha ](1+\gamma _5)b_i|\overline{B}_d^0]\left\}\right|_{k=puq}`$ (113)
$`=`$ $`2{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}({\displaystyle \frac{3}{2}}+{\displaystyle _0^1}𝑑u{\displaystyle \frac{\varphi _\sigma (u)}{6\overline{u}}})\pi ^{}|\overline{d}_i(1+\gamma _5)u_i|0\pi ^+|\overline{u}_j(1\gamma _5)b_j|\overline{B}_d^0.`$ (114)
#### 3 Hard scattering with the spectator

The chirally enhanced power corrections in hard spectator scattering not only occurs in the case of $`(VA)(V+A)`$ vertex insertion, but also in the case of $`(VA)(VA)`$ insertion. But in the case of $`(VA)(V+A)`$ insertion, after a straightforward calculation, we find that there will be serious linear divergence at the end points of the LCDAs if the twist-3 LCDAs are not symmetric. Because infrared finiteness of the vertex corrections requires that the twist-3 LCDAs, especially $`\varphi _\sigma (u)`$, must be symmetric, we shall implicitly assume this symmetric condition for the LCDAs in latter computation. So, in this symmetric condition, the hard scattering with the spectator vanishes when $`(VA)(V+A)`$ vertex is inserted. However, even in this strict symmetric condition, there is still a logarithmic divergence from the end point of the recoiling pion in hard spectator scattering when $`(VA)(VA)`$ vertex is inserted. For example,

$$Q_2^u_{(g)+(h)}=i\pi \alpha _sf_Bf_\pi ^2\frac{C_F}{N^2}_0^1𝑑\xi \frac{\varphi _B(\xi )}{\xi }_0^1𝑑u\frac{\varphi (u)}{u}_0^1𝑑v\left[\frac{\varphi (v)}{\overline{v}}+\frac{2\mu _\pi }{m_B}\frac{\varphi _\sigma (v)}{6\overline{v}^2}\right].$$
(115)
This means that QCD factorization is broken down. But we can still give a phenomenological treatment for this hard spectator scattering. By using the asymptotic form of $`\varphi _\sigma (u)`$, we find that there is a divergent integral over $`v`$: $`_0^1𝑑v(1/\overline{v})`$. In Refs., the authors prefer to introduce a phenomenological parameterization for this logarithmic divergence. They take $`_0^1𝑑v(1/\overline{v})=\mathrm{ln}(m_B/\mathrm{\Lambda }_B)+re^{i\theta }`$, where $`r`$ is taken from 3 (realistic) to 6 (conservative), and the phase $`\theta `$ from $`\pi `$ to $`\pi `$. We shall take similar phenomenological treatment in the numerical computation below.

We notice that the above approach of evaluating hard spectator contribution is naive. For instance, the scale of hard spectator contribution should be different from the vertex correction contribution. While it seems reasonable to take the scale $`\mu 𝒪(m_b)`$ for the vertex correction diagrams to avoid large logarithm $`\alpha _s\mathrm{log}(\mu /m_b)`$, a natural choice of the scale of hard spectator contribution may be around $`𝒪(1GeV)`$ because the average momentum squared of the exchanged gluon is about $`1GeV^2`$. Another disturbing feature of hard spectator contribution is that, as pointed out in Refs. , when including the contribution of $`\varphi _\sigma `$, there would appear a divergent integral $`_0^1𝑑v(1/\overline{v})`$ even if the symmetric distribution amplitude is applied. This divergent integral implies that the dominant contribution comes from the end-point region, or, in other words, it is dominated by soft gluon exchange. However, the transverse momentum may not be omitted in the end-point region ; if so, the corresponding divergent integral would then be changed to

$$𝑑v\frac{1}{\overline{v}}𝑑vd^2k_T\frac{\mathrm{\Psi }(v,k_T)}{\overline{v}\xi m_b^2+k_T^2}.$$
(116)
As an illustration, we do not consider the $`k_T`$ dependence of wave functions (though it is certainly not a good approximation); then the above integral is proportional to

$$\frac{dvdk_T^2}{\overline{v}\xi m_b^2+k_T^2}\frac{dxdy}{x+y}.$$
(117)
The above integration converges now; furthermore it is not dominated by end-point contribution. This illustrates that the treatment of hard spectator diagrams may need further discussion.

There exists “annihilation topology” contributions which may belong to chirally enhanced corrections. In Ref. , the authors have discussed this topic and find that a divergent integral $`[(dx/x)]^2`$ will appear. We suspect that this divergence may disappear, similar to the hard spectator term, if the effect of transverse momenta can be included. It is also possible that “annihilation topology” contributions are really dominated by soft interactions and thus violate factorization. Due to its complexity, we do not include “annihilation topology” contributions in this work.

### C Final formulas

With these effective operators, $`BP_1P_2`$ decay amplitudes in QCD factorization can be written as

$$A(BP_1P_2)=\frac{G_F}{\sqrt{2}}\underset{p=u,c}{}\underset{i=1,10}{}v_pa_i^pP_1P_2|Q_i|\text{B}_F,$$
(118)
where $`v_p`$ is CKM factor, $`P_1P_2|Q_i|`$B$`_F`$ is the factorized matrix element and is the same as the definition of the BSW Largrangian. Then as an illustration, the explicit expressions of $`a_i^p`$ ($`i=1`$ to $`10`$) for $`B\pi \pi `$ (using symmetric LCDAs of the pion) are obtained. But it is easy to generalize these formulas to the case that the final states are other light pseudoscalars. Furthermore, we take only part of QED corrections into account in our final formula, in particular the QED penguin insertions. Now $`a_i^p`$ for $`B\pi \pi `$ in NDR $`\gamma _5`$ scheme is listed as follows Because of the tedium, we do not calculate the radiative corrections in the HV scheme. However, generally, the results in the NDR scheme and HV scheme can be related by a constant matrix $`\mathrm{\Delta }\widehat{r}_s=\widehat{r}_{s,HV}\widehat{r}_{s,NDR}`$  which is free from the gauge dependence and infrared structure of the theory. Thus, in principle, we can obtain the results in the HV scheme just by using $`\mathrm{\Delta }\widehat{r}_s`$. In , the constant matrix has been applied to obtain the results in the NDR and HV scheme for the coefficients $`a_i`$ which only contain the current-current vertex corrections. But whether we can obtain the expression of $`a_6`$ or $`a_8`$ in HV scheme in a similar way needs further discussion.:

$`a_1^u`$ $`=`$ $`C_1+{\displaystyle \frac{C_2}{N}}+{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}C_2F,`$ (119)
$`a_2^u`$ $`=`$ $`C_2+{\displaystyle \frac{C_1}{N}}+{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}C_1F,`$ (120)
$`a_3`$ $`=`$ $`C_3+{\displaystyle \frac{C_4}{N}}+{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}C_4F,`$ (121)
$`a_4^p`$ $`=`$ $`C_4+{\displaystyle \frac{C_3}{N}}+{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}C_3F`$ (124)
$`{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}\{C_1({\displaystyle \frac{4}{3}}\mathrm{log}{\displaystyle \frac{\mu }{m_b}}+G(s_p){\displaystyle \frac{2}{3}})+(C_3{\displaystyle \frac{C_9}{2}})({\displaystyle \frac{8}{3}}\mathrm{log}{\displaystyle \frac{\mu }{m_b}}+G(0)+G(1){\displaystyle \frac{4}{3}})`$
$`+{\displaystyle \underset{q=u,d,s,c,b}{}}(C_4+C_6+{\displaystyle \frac{3}{2}}\mathrm{e}_\mathrm{q}C_8+{\displaystyle \frac{3}{2}}\mathrm{e}_\mathrm{q}C_{10})({\displaystyle \frac{4}{3}}\mathrm{log}{\displaystyle \frac{\mu }{m_b}}+G(s_q))+G_8C_{8G}\},`$
$`a_5`$ $`=`$ $`C_5+{\displaystyle \frac{C_6}{N}}+{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}C_6(F12),`$ (125)
$`a_6^p`$ $`=`$ $`C_6+{\displaystyle \frac{C_5}{N}}{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}6C_5`$ (130)
$`{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}\{C_1((1+{\displaystyle \frac{2}{3}}A_\sigma )\mathrm{log}{\displaystyle \frac{\mu }{m_b}}{\displaystyle \frac{1}{2}}{\displaystyle \frac{1}{3}}A_\sigma +G^{}(s_p)+G^\sigma (s_p))`$
$`+{\displaystyle \underset{q=d,b}{}}(C_3{\displaystyle \frac{C_9}{2}})((1+{\displaystyle \frac{2}{3}}A_\sigma )\mathrm{log}{\displaystyle \frac{\mu }{m_b}}{\displaystyle \frac{1}{2}}{\displaystyle \frac{1}{3}}A_\sigma +G^{}(s_q)+G^\sigma (s_q))`$
$`+{\displaystyle \underset{q=u,d,s,c,b}{}}(C_4+C_6+{\displaystyle \frac{3}{2}}\mathrm{e}_\mathrm{q}C_8+{\displaystyle \frac{3}{2}}\mathrm{e}_\mathrm{q}C_{10})\left((1+{\displaystyle \frac{2}{3}}A_\sigma )\mathrm{log}{\displaystyle \frac{\mu }{m_b}}+G^{}(s_q)+G^\sigma (s_q)\right)`$
$`+({\displaystyle \frac{3}{2}}+A_\sigma )C_{8G}\},`$
$`a_7`$ $`=`$ $`C_7+{\displaystyle \frac{C_8}{N}}+{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}C_8(F12),`$ (131)
$`a_8^p`$ $`=`$ $`C_8+{\displaystyle \frac{C_7}{N}}{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}6C_7`$ (136)
$`{\displaystyle \frac{\alpha _{em}}{9\pi }}\{(C_2+{\displaystyle \frac{C_1}{N}})((1+{\displaystyle \frac{2}{3}}A_\sigma )\mathrm{log}{\displaystyle \frac{\mu }{m_b}}{\displaystyle \frac{1}{2}}{\displaystyle \frac{1}{3}}A_\sigma +G^{}(s_p)+G^\sigma (s_p))`$
$`+(C_4+{\displaystyle \frac{C_3}{N}}){\displaystyle \underset{q=d,b}{}}{\displaystyle \frac{3}{2}}\mathrm{e}_\mathrm{q}((1+{\displaystyle \frac{2}{3}}A_\sigma )\mathrm{log}{\displaystyle \frac{\mu }{m_b}}{\displaystyle \frac{1}{2}}{\displaystyle \frac{1}{3}}A_\sigma +G^{}(s_q)+G^\sigma (s_q))`$
$`+(C_3+{\displaystyle \frac{C_4}{N}}+C_5+{\displaystyle \frac{C_6}{N}}){\displaystyle \underset{q=u,d,s,c,b}{}}{\displaystyle \frac{3}{2}}\mathrm{e}_\mathrm{q}\left((1+{\displaystyle \frac{2}{3}}A_\sigma )\mathrm{log}{\displaystyle \frac{\mu }{m_b}}+G^{}(s_q)+G^\sigma (s_q)\right)`$
$`+({\displaystyle \frac{3}{4}}+{\displaystyle \frac{1}{2}}A_\sigma )C_{7\gamma }\},`$
$`a_9`$ $`=`$ $`C_9+{\displaystyle \frac{C_{10}}{N}}+{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}C_{10}F,`$ (137)
$`a_{10}^p`$ $`=`$ $`C_{10}+{\displaystyle \frac{C_9}{N}}+{\displaystyle \frac{\alpha _s}{4\pi }}{\displaystyle \frac{C_F}{N}}C_9F{\displaystyle \frac{\alpha _{em}}{9\pi }}\{(C_2+{\displaystyle \frac{C_1}{N}})({\displaystyle \frac{4}{3}}\mathrm{log}{\displaystyle \frac{\mu }{m_b}}+G(s_p){\displaystyle \frac{2}{3}})`$ (140)
$`+(C_4+{\displaystyle \frac{C_3}{N}}){\displaystyle \underset{q=d,b}{}}{\displaystyle \frac{3}{2}}\mathrm{e}_\mathrm{q}({\displaystyle \frac{4}{3}}\mathrm{log}{\displaystyle \frac{\mu }{m_b}}+G(s_q){\displaystyle \frac{2}{3}})`$
$`+(C_3+{\displaystyle \frac{C_4}{N}}+C_5+{\displaystyle \frac{C_6}{N}}){\displaystyle \underset{q=u,d,s,c,b}{}}{\displaystyle \frac{3}{2}}\mathrm{e}_\mathrm{q}({\displaystyle \frac{4}{3}}\mathrm{log}{\displaystyle \frac{\mu }{m_b}}+G(s_q))+{\displaystyle \frac{1}{2}}G_8C_{7\gamma }\}.`$
Here $`N=3`$ is the number of color, $`C_F=(N^21)/2N`$ is the factor of color, $`s_q=m_q^2/m_b^2`$ and we define the other symbols in the above expressions as

$`F=12\mathrm{ln}{\displaystyle \frac{\mu }{m_b}}18+f^I+f^{II},`$ (141)
$`f^I={\displaystyle \underset{0}{\overset{1}{}}}𝑑xg(x)\varphi (x),G_8={\displaystyle \underset{0}{\overset{1}{}}}𝑑xG_8(x)\varphi (x),`$ (142)
$`G(s)={\displaystyle \underset{0}{\overset{1}{}}}𝑑xG(s,x)\varphi (x),`$ (143)
$`G^{}(s)={\displaystyle \underset{0}{\overset{1}{}}}𝑑xG^{}(s,x)\varphi _p(x),`$ (144)
$`G^\sigma (s)={\displaystyle \underset{0}{\overset{1}{}}}𝑑xG^\sigma (s,x){\displaystyle \frac{\varphi _\sigma (x)}{6(1x)}},A_\sigma ={\displaystyle \underset{0}{\overset{1}{}}}𝑑x{\displaystyle \frac{\varphi _\sigma (x)}{6(1x)}},`$ (145)
where $`\varphi (x)`$ \[$`\varphi _p(x)`$, $`\varphi _\sigma (x)`$\] is leading twist (twist-3) LCDA of the ejected pion, and the hard-scattering functions are

$`g(x)=3{\displaystyle \frac{12x}{1x}}\mathrm{ln}x3i\pi ,G_8(x)={\displaystyle \frac{2}{1x}},`$ (146)
$`G(s,x)=4{\displaystyle \underset{0}{\overset{1}{}}}𝑑uu(1u)\mathrm{ln}(su(1u)(1x)iϵ),`$ (147)
$`G^{}(s,x)=3{\displaystyle \underset{0}{\overset{1}{}}}𝑑uu(1u)\mathrm{ln}(su(1u)(1x)iϵ),`$ (148)
$`G^\sigma (s,x)=2{\displaystyle \underset{0}{\overset{1}{}}}𝑑uu(1u)\mathrm{ln}(su(1u)(1x)iϵ)`$ (149)
$`+{\displaystyle \underset{0}{\overset{1}{}}}𝑑u{\displaystyle \frac{u^2(1u)^2(1x)}{su(1u)(1x)iϵ}}.`$ (150)
The contributions from the hard spectator scattering \[Figs. 1(g),(1)(h)\] are reduced to the factor $`f^{II}`$:

$$f^{II}=\frac{4\pi ^2}{N}\frac{f_\pi f_B}{F_+^{B\pi }(0)m_B^2}\underset{0}{\overset{1}{}}𝑑\xi \frac{\mathrm{\Phi }_B(\xi )}{\xi }\underset{0}{\overset{1}{}}𝑑x\frac{\varphi (x)}{x}\underset{0}{\overset{1}{}}𝑑y\left[\frac{\varphi (y)}{1y}+\frac{2\mu _\pi }{M_B}\frac{\varphi _\sigma (y)}{6(1y)^2}\right].$$
(151)
There contains a divergent integral in $`f^{II}`$. Here we simply assume that $`(dy/y)ln(m_b/\mathrm{\Lambda }_{QCD})`$ (similar to what has been done in Refs. , though our assumption here is certainly an oversimplification). We thus illustrate numerically the scale dependence of $`a_i^p(\pi \pi )`$ in Table 1. Here we use the asymptotic form of the LCDAs of the light pseudoscalar meson which are listed in Appendix A, and the other input parameters are taken as follows: $`F^{B\pi }(0)=0.33`$, $`f_B=0.2GeV`$, $`f_\pi =133MeV`$, the pole masses $`m_b=4.8GeV`$, $`m_c=1.4GeV`$, the $`\overline{MS}`$ masses $`\overline{m}_t(\overline{m}_t)=170GeV`$, $`\overline{m}_b(\overline{m}_b)=4.4GeV`$, $`\overline{m}_u(2GeV)=4.2MeV`$, $`\overline{m}_d(2GeV)=7.6MeV`$ and $`\mathrm{\Lambda }_{QCD}^{(5)}=225MeV`$.

## IV Discussions and General Remarks

### A Color transparency and factorization

Color transparency gives a clear physics picture of QCD factorization. In the argument of the color transparency, a fast-moving small color singlet formed by a pair of $`q\overline{q}`$ decouples from the surrounding soft gluons. Of course, as mentioned in the previous section, only the decoupling with the soft gluons is not enough for a factorization formula, the decoupling from the collinear divergence is also necessary. We really find that both of the requirements can be satisfied in the one-loop calculations. So the QCD factorization is guaranteed. Therefore, the calculations in the above sections seem to be a one-loop technical manifestation (or demonstration) of the color transparency. On the other hand, at the leading power of $`\mathrm{\Lambda }_{QCD}/m_b`$, the soft or collinear gluons only “see” the direction of the light meson, but are “blind” to the spins of the quark constituents. So the soft or collinear gluon cannot distinguish whether the ejected meson from $`b`$ quark decay is a light pseudoscalar or a light longitudinally polarized vector meson. As a consequence, the cancellation of the infrared divergence is universal for $`B`$ decays to two light mesons, no matter whether the meson is a pseudoscalar meson or a vector meson. Therefore, the QCD factorization formula for $`BPP`$ at the leading power of $`\mathrm{\Lambda }_{QCD}/m_b`$ is easy to be generalized to $`BPV`$ and $`VV`$.

Similarly, the color transparency argument can not only be applied to the strong interactions, but also be generalized to the electromagnetic interactions. When the ejected meson is electric neutral, the soft photons also decouple from the fast moving small electric dipole. So QED vertex corrections are also of infrared finiteness. But for the case that the ejected meson is charged, QED corrections are infrared divergent, and the infrared divergence must be cancelled by the soft photon emission mechanism, which is common in the calculation of the radiative corrections for $`\mu ^{}e^{}\nu _\mu \overline{\nu }_e`$. About this, it is easy to be covered in the calculation in the previous section, just replacing the QCD vertex by a QED vertex. This can be called a one-loop demonstration for “charge” transparency.

It should be noted that the above arguments must be based on the condition that the ejected meson is in a very compact configuration, then it, as a small color dipole, is disentangled with the soft gluons. Otherwise, if its size is too large, it is difficult to decouple from the soft gluons. For example, the spectator quark in D meson is very soft, and runs around c quark like a soft quark cloud, which has a large overlap with the $`B`$ meson spectator system. As a consequence, the process in which D meson is ejected from $`b`$ decay is dominated by the soft gluon exchange.

### B The scale dependence

From the expressions of the QCD coefficients $`a_i^p`$ obtained in previous sections, it is found that the renormalization scale dependence of the hadronic matrix elements of the effective operators is recovered. Apparently, we expect this recovered dependence can cancel the scale dependence of the Wilson coefficients $`C_i`$.

With the renormalization group equations for the Wilson coefficients $`C_i(\mu )`$ at leading order logarithm approximation,

$$\mu \frac{d}{d\mu }𝐂(\mu )=\frac{\alpha _s}{4\pi }\widehat{\gamma }^{(0)T}𝐂(\mu )$$
(152)
we do find

$$\mu \frac{d}{d\mu }a_i^p=0(\mathrm{for}\mathrm{i}=15\mathrm{and}7,9,10)$$
(153)
when we neglect the contributions from higher order of $`\alpha _s`$. But for $`a_6`$ or $`a_8`$, some scale dependence at the order of $`\alpha _s`$ still remains. Note that other QCD coefficients ($`a_{i=1,2,3,4,5,7,9,10}`$) are multiplied by the product of the matrix elements of the conserved currents which are independent of the renormalization scale; whereas the coefficient $`a_6`$ or $`a_8`$ is multiplied by a product of the two matrix elements of scalar and pseudoscalar current

$`2P_1|\overline{d}(1+\gamma _5)q|0P_2|\overline{q}(1\gamma _5)`$b$`|`$B$`.`$which is still of scale dependence. This scale dependence is generally represented by the factor

$`r_\chi (\mu )={\displaystyle \frac{2m_{P_1}^2}{\overline{m}_b(\mu )(\overline{m}_1(\mu )+\overline{m}_2(\mu ))}}`$after we apply the equations of motion to transform the $`(S+P)(SP)`$ matrix elements into the type of $`(VA)(VA)`$. Here $`m_1`$ and $`m_2`$ are the current masses of the valence quarks in meson $`P_1`$. With the renormalization group equations for the running mass of the current quark

$$\mu \frac{d}{d\mu }\overline{m}(\mu )=6\frac{\alpha _s}{4\pi }C_F\overline{m}(\mu ),$$
(154)
we have

$$\mu \frac{d}{d\mu }r_\chi (\mu )=12\frac{\alpha _s}{4\pi }C_Fr_\chi (\mu ).$$
(155)
Consequently, we find

$$\mu \frac{d}{d\mu }\left(a_{6,8}^pr_\chi \right)=0$$
(156)
with the asymptotic form of $`\varphi _\sigma (u)=6u(1u)`$. Then, as we expect, the decay amplitude for $`B`$ decays to two light pseudoscalars predicted by QCD factorization is really independent of the renormalization scale within the constraint $`A_\sigma =1/2`$. This also can be obviously seen from the numerical results of $`a_i^p(\pi \pi )`$ listed in Table 1. In particular, if we think that the results of QCD factorization are reliable and really independent of the renormalization scale, maybe $`A_\sigma =1/2`$ is a strict constraint for the form of $`\varphi _\sigma (u)`$.

It should be noted that the imaginary part of QCD coefficients $`a_i`$ only arises at the order of $`\alpha _s`$, and depends on the renormalization scale. This dependence would bring some uncertainties in determining the $`CP`$ asymmetries in $`B`$ decays. Maybe this scale dependence of the imaginary part could be cancelled by the results on higher order of $`\alpha _s`$.

### C Comparison to the generalized factorization and PQCD method

Comparing QCD factorization approach with the generalized factorization and PQCD method, some interesting comments are in order.

(i) At the zeroth order of $`\alpha _s`$, both of QCD factorization (BBNS approach) and generalized factorization can reproduce the results of “naive factorization”; at the higher order of $`\alpha _s`$, the renormalization scheme and scale dependence for the hadronic matrix elements can be recovered from the hard-scattering kernels $`T_i^I`$ in BBNS approach and $`\widehat{𝐦}_s`$ in generalized factorization. However, in generalized factorization, $`\widehat{𝐦}_s`$ is from the one-loop calculations of quark-level matrix elements. According to Buras et al., quark-level matrix elements are accompanied with infrared divergences. To avoid these divergences, one may assume that external quarks are off-shell. Unfortunately, it will introduce gauge dependence which is also unphysical. But in the BBNS approach, it is different because the external states are all physical and can be approximated as on-shell quarks in the leading order of $`\mathrm{\Lambda }_{QCD}/m_b`$. As a consequence, the unphysical gauge dependence does not appear. In the PQCD method, Li et al. claim that their method is based on a six-quark system in which the external quarks are on-shell, so the gauge invariance of PQCD predictions is guaranteed. The scale dependence in PQCD prediction is removed by evolving the Wilson coefficients down to the proper hard scale. So no explicit scale dependence is left.

(ii) In generalized factorization and BBNS approach, the hadronic transition form factors are not calculable, and they are dominated by soft gluon exchange and determined only by experiments or some non-perturbative approaches such as sum rule, lattice QCD, etc. In particular, the above assumption can be also justified in the BBNS approach by naive power counting . However, in the PQCD method, this naive power counting rule may be invalid when the transverse momenta of the quark constituents and Sudakov suppression are taken into account. So, Li et al. thought that the hadronic transition form factors can be calculated in the PQCD method because $`b`$ quark is heavy enough and the soft gluon exchange is suppressed by Sudakov form factors. This is the essential difference between the BBNS approach and the PQCD method.

(iii) The generalized factorization considers “nonfactorizable” contributions as intractable. Therefore, one may introduce one or more effective color numbers $`N_c^{eff}`$ to phenomenologically represent “nonfactorizable” contributions . Furthermore, $`N_c^{eff}`$ is assumed to be universal to maintain predictive power. However, “nonfactorizable” contribution is really process-dependent. In the BBNS approach and the PQCD method, such “nonfactorizable” contributions are indeed calculable in perturbative theory. In consequence, $`N_c^{eff}`$ need not be introduced.

(iv) As mentioned in the above sections, the strong phases predicted by generalized factorization are only from the BSS mechanism which is represented by the penguin insertion. However, the virtuality of the gluon or photon $`k^2`$ in the penguin insertion is ambiguous in generalized factorization, and usually it is approximated around $`m_B^2/2`$. This brings significant uncertainties for predicting the $`CP`$ asymmetries for B decays. A particular interesting result of the BBNS approach is that strong phases are not only from the BSS mechanism but also from the hard scattering, and there are no uncertainties in determining $`k^2`$ of penguin insertion. However, compared with the real part of the decay amplitude, the imaginary part is $`𝒪(\alpha _s)`$ or power $`\mathrm{\Lambda }_{QCD}/m_b`$ suppressed and cannot lead to large $`CP`$ asymmetries, since they come solely from hard scattering processes which are only calculable in the heavy quark limit. In the PQCD method, there is no such $`\alpha _s`$ suppression in the imaginary part of the decay amplitude. Thus $`CP`$ asymmetries predicted by the PQCD method are usually greater than the prediction of BBNS approach and generalized factorization. So maybe these differences of prediction for $`CP`$ asymmetries can be an experimental test for the BBNS and PQCD approaches.

(v) Hard spectator contributions \[Fig. 1(g)and 1(h)\], which are leading power effects in QCD factorization, miss out in “naive factorization” and “generalized factorization”. They are, however, $`𝒪(\alpha _s)`$ suppressed compared with the leading factorized contributions (the hadronic transition form factors). But, in the PQCD method, they are of the same order as the form factors.

(vi) In the PQCD method, penguin contributions receive a dynamical enhancement called “Fat Penguin”. But in generalized and QCD factorization, they are missed. This enhancement in the PQCD method arises from the strong scale dependence of the penguin Wilson coefficients $`C_{4,6}`$, etc.

(vii) Final state interactions (FSI) do not appear in the three methods. In QCD factorization, Beneke et al. point out that the cancellation of the infrared divergences implies that the long distance FSI is power suppressed due to the quark-hadron duality. However, this point of view is controversial , but can be examined by the experimental measurements of $`BKK`$.

### D Limitation of QCD factorization

QCD factorization formula only holds in the heavy quark limit $`m_b\mathrm{}`$. In the real world $`m_b`$ is only about $`4.8`$ GeV, the validity of the power suppression may be questionable. In particular, for several cases, the power suppressed corrections can be numerically large, because the perturbative expansion is in order of $`\alpha _s`$ which is not small at the realistic scale $`𝒪(m_b)`$ compared to $`\mathrm{\Lambda }_{QCD}/m_b`$.

(i) The hard “nonfactorizable” contributions computed by QCD factorization are generally small compared to the leading “factorizable” contribution. But when the leading “factorizable” contributions are color suppressed, the “nonfactorizable” contribution may be larger than the leading results. At the same time, the potentially soft contribution, which is formally power suppressed, may be important. For example, in $`\overline{B}_d^0\pi ^0\pi ^0`$, any perturbative and soft power suppressed contributions can have a significant effect on predicting the branching ratio and $`CP`$ asymmetry. Furthermore, this problem also arises when the entire leading power contribution is suppressed by small Wilson coefficients, for example, in $`BKK`$; or when the leading power contribution is suppressed by the small CKM elements.

(ii) An important power suppressed contribution is from the higher twist light-cone wave functions of the light mesons. The chirally enhanced power correction from the two-particle twist-3 wave functions is the most important, and has been partly involved in this work except for the annihilation topologies. Other contribution from multiparticle non-valence fock state has been proved to be also power suppressed. However, there is no systematic way to evaluate it. The author of Ref. proposed a way to evaluate the soft gluon exchange contribution from higher twist $`q\overline{q}g`$ wave functions within the frame of the light-cone sum rule(LCSR). But the accuracy of LCSR is limited due to the quark-hadron duality approximation. On the other hand, power correction from transverse momenta needs a subtle treatment in the future. In Ref., the authors point out that the contribution from the transverse momenta might be considered when we evaluate the hadronic matrix elements to two-loop order. In this case, it is possible that Sudakov suppression might be taken into account as well.

In summary, up to now, we do not have a systematic way to evaluate many kinds of power suppressed corrections for exclusive processes. How to evaluate such corrections in a consistent way within the frame of QCD factorization is a potentially interesting work.

## V summary

In this work, we give a detailed discussion for QCD factorization involving the complete chirally enhanced power corrections in the heavy quark limit for $`B`$ decays to two light pseudoscalar mesons, and present some elaborate calculations of radiative corrections at the order of $`\alpha _s`$. We point out that the infrared finiteness of the vertex corrections in the chirally enhanced power corrections requires twist-3 light-cone distribution amplitudes (LCDAs) of the light pseudoscalar symmetric. However, even in the symmetric condition, there is also infrared divergence from the end point of the LCDAs in the hard spectator scattering and annihilation topology. So the transverse momenta and Sudakov suppression should be taken into account. We also point out that the decay amplitude of $`BPP`$ predicted by QCD factorization is really independent of the renormalization scale, at least at the order of $`\alpha _s`$. At last, we briefly compare the QCD factorization to the generalized factorization and PQCD method which are generally used in studying $`B`$ exclusive hadronic decays.

Acknowledgements

We thank Professor Hai-Yang Cheng and Kwei-Chou Yang for pointing out errors in the coefficients of $`C_{8G}`$ and $`C_{7\gamma }`$ in Eqs.(50) and (52). And we also thank Professor Maozhi Yang for helpful discussions about QCD factorization and the PQCD method. This work is supported in part by National Science Foundation of China and State Commission of Science and Technology of China.

Appendix A. Twist-2 and -3 LCDAs of Light Pseudoscalar Meson

Two particle twist-2 and twist-3 light-cone distribution amplitudes of light pseudoscalar mesons are defined by the following nonlocal matrix elements:

$`P(p^{})|\overline{q}(y)\gamma _\mu \gamma _5q(x)|0`$ $`=`$ $`if_Pp_\mu ^{}{\displaystyle _0^1}𝑑ue^{iup^{}y+i\overline{u}p^{}x}\varphi (u),`$ (157)
$`P(p^{})|\overline{q}(y)\gamma _5q(x)|0`$ $`=`$ $`if_P\mu _P{\displaystyle _0^1}𝑑ue^{iup^{}y+i\overline{u}p^{}x}\varphi _p(u),`$ (158)
$`P(p^{})|\overline{q}(y)\sigma _{\mu \nu }\gamma _5q(x)|0`$ $`=`$ $`if_P\mu _P(p_\mu ^{}z_\nu p_\nu ^{}z_\mu ){\displaystyle _0^1}𝑑ue^{iup^{}y+i\overline{u}p^{}x}{\displaystyle \frac{\varphi _\sigma (u)}{6}},`$ (159)
with $`f_P`$ being the decay constant of the light pseudoscalar, $`\mu _P=M_P^2/(m_1+m_2)`$ ($`m_1`$ and $`m_2`$ are the masses of the constituent quarks in the pseudoscalar), and $`z=yx`$. Here $`\varphi (u)`$ is the twist-2 light-cone distribution amplitude; $`\varphi _p(u)`$ and $`\varphi _\sigma (u)`$ are two-particle twist-3 distribution amplitudes. The above definitions can be combined into the below nonlocal matrix element:

$`P(p^{})|\overline{q}_\alpha (y)q_\beta (`$ $`x`$ $`)|0={\displaystyle \frac{if_P}{4}}{\displaystyle }_0^1e^{iup^{}y+i\overline{u}p^{}x}`$ (161)
$`\times \{/p^{}\gamma _5\varphi (u)\mu _P\gamma _5(\varphi _p(u)\sigma _{\mu \nu }p^\mu z^\nu {\displaystyle \frac{\varphi _\sigma (u)}{6}})\}_{\beta \alpha }.`$
Neglecting the three-particle twist-3 light-cone wave function, the asymptotic forms of the above distribution amplitudes are given as

$`\varphi (u)`$ $`=`$ $`6u(1u),`$ (162)
$`\varphi _p(u)`$ $`=`$ $`1,`$ (163)
$`\varphi _\sigma (u)`$ $`=`$ $`6u(1u).`$ (164)
Appendix B. The Evolution of $`C_i(\mu )`$

The renormalization group equation for the Wilson coefficients $`C_i(\mu )`$ is written as follows:

$$\mu \frac{d}{d\mu }𝐂(\mu )=\widehat{\gamma }^T𝐂(\mu )$$
(165)
Here $`\gamma `$ is the anomalous dimension matrix, which can be calculated by the perturbative theory and expanded in order of the coupling constants $`\alpha _s`$ and $`\alpha _{em}`$:

$$\widehat{\gamma }=\frac{\alpha _s}{4\pi }\widehat{\gamma }_s^{(0)}+(\frac{\alpha _s}{4\pi })^2\widehat{\gamma }_s^{(1)}+\frac{\alpha _{em}}{4\pi }\widehat{\gamma }_e^{(0)}+\frac{\alpha _s\alpha _{em}}{(4\pi )^2}\widehat{\gamma }_{se}^{(1)}+\mathrm{}$$
(166)
The LO anomalous dimension matrix $`\gamma _s^{(0)}`$ of the above equations has the explicit form

$$\widehat{\gamma }_s^{(0)}=\left(\begin{array}{cccccccccc}\frac{6}{N}& 6& \frac{2}{3N}& \frac{2}{3}& \frac{2}{3N}& \frac{2}{3}& 0& 0& 0& 0\\ 6& \frac{6}{N}& 0& 0& 0& 0& 0& 0& 0& 0\\ 0& 0& \frac{22}{3N}& \frac{22}{3}& \frac{4}{3N}& \frac{4}{3}& 0& 0& 0& 0\\ 0& 0& 6\frac{2f}{3N}& \frac{6}{N}+\frac{2f}{3}& \frac{2f}{3N}& \frac{2f}{3}& 0& 0& 0& 0\\ 0& 0& 0& 0& \frac{6}{N}& 6& 0& 0& 0& 0\\ 0& 0& \frac{2f}{3N}& \frac{2f}{3}& \frac{2f}{3N}& \frac{6(N^21)}{N}+\frac{2f}{3}& 0& 0& 0& 0\\ 0& 0& 0& 0& 0& 0& \frac{6}{N}& 6& 0& 0\\ 0& 0& \frac{2(ud/2)}{3N}& \frac{2(ud/2)}{3}& \frac{2(ud/2)}{3N}& \frac{2(ud/2)}{3}& 0& \frac{6(N^21)}{N}& 0& 0\\ 0& 0& \frac{2}{3N}& \frac{2}{3}& \frac{2}{3N}& \frac{2}{3}& 0& 0& \frac{6}{N}& 6\\ 0& 0& \frac{2(ud/2)}{3N}& \frac{2(ud/2)}{3}& \frac{2(ud/2)}{3N}& \frac{2(ud/2)}{3}& 0& 0& 6& \frac{6}{N}\end{array}\right)$$
(167)
where $`N`$ is the color number, $`f`$ is the active flavor number, and $`u`$ and $`d`$ denote the number of the active up- and down-type flavors respectively.

Appendix C. Some Useful Feynman Parameter Integrals

In calculation of the perturbative diagrams shown in Fig. 1, one might encounter some Feynman parameter integrals which involve nontrivial infrared divergence. To deal with the infrared divergence, as mentioned in preceding sections, the dimensional regularization (DR) and massive gluon (MG) scheme are applied. Below, we give the explicit calculation of some useful Feynman parameter integrals in the above two regularization schemes.

First, we deal with the integrals in DR scheme. In the DR scheme (here we take $`d=4+2a`$ and $`a>0`$), the integrals involving the infrared divergence are written as follows:

$`{\displaystyle _0^1}𝑑t_1{\displaystyle _0^{1t_1}}𝑑t_2{\displaystyle \frac{1}{(t_1(t_1+t_2u))^{1a}}}`$ $`=`$ $`{\displaystyle \frac{1}{u}}\left[{\displaystyle \frac{1}{2a^2}}+{\displaystyle \frac{\mathrm{ln}u}{a}}+{\displaystyle \frac{1}{2}}\mathrm{ln}^2u\mathrm{Li}_2(1{\displaystyle \frac{1}{u}})\right],`$ (168)
$`{\displaystyle _0^1}𝑑t_1{\displaystyle _0^{1t_1}}𝑑t_2{\displaystyle \frac{t_2}{(t_1(t_1+t_2u))^{1a}}}`$ $`=`$ $`{\displaystyle \frac{1}{u}}\left[{\displaystyle \frac{1}{a}}2+\mathrm{ln}u+{\displaystyle \frac{\mathrm{ln}u}{1u}}\right],`$ (169)
$`{\displaystyle _0^1}𝑑t_1{\displaystyle _0^{1t_1}}𝑑t_2{\displaystyle \frac{(1t_1)(1t_2)}{(t_1t_2u)^{1a}}}`$ $`=`$ $`{\displaystyle \frac{1}{u}}[{\displaystyle \frac{1}{a^2}}+{\displaystyle \frac{\mathrm{ln}(u)2}{a}}+{\displaystyle \frac{27\pi ^2}{6}}`$ (170)
$``$ $`2\mathrm{ln}(u)+{\displaystyle \frac{1}{2}}\mathrm{ln}^2(u)].`$ (171)
Here $`\mathrm{Li}_2(x)`$ is the dilogarithm function. It is defined by

$$\mathrm{Li}_2(x)=_0^x\frac{\mathrm{ln}(1t)}{t}𝑑t.$$
(172)
The Feynman parameter integrals in the MG scheme are listed as follows:

$`{\displaystyle _0^1}𝑑t_1{\displaystyle _0^{1t_1}}𝑑t_2{\displaystyle \frac{1}{(t_1(t_1+t_2u)+(1t_1t_2)\lambda )}}`$ (173)
$`=`$ $`{\displaystyle \frac{1}{u}}\left[{\displaystyle \frac{1}{4}}\mathrm{ln}^2\lambda +\mathrm{ln}(u)\mathrm{ln}\lambda 2\mathrm{ln}u\mathrm{ln}\lambda +{\displaystyle \frac{1}{2}}\mathrm{ln}^2u\mathrm{Li}_2(1{\displaystyle \frac{1}{u}})+{\displaystyle \frac{5}{4}}\pi ^2\right],`$ (175)
$`{\displaystyle _0^1}𝑑t_1{\displaystyle _0^{1t_1}}𝑑t_2{\displaystyle \frac{t_2}{(t_1(t_1+t_2\omega )+(1t_1t_2)\lambda )}}`$
$`=`$ $`{\displaystyle \frac{1}{u}}\left[\mathrm{ln}\lambda 1+\mathrm{ln}u+{\displaystyle \frac{\mathrm{ln}u}{1u}}\right],`$ (177)
$`{\displaystyle _0^1}𝑑t_1{\displaystyle _0^{1t_1}}𝑑t_2{\displaystyle \frac{(1t_1)(1t_2)}{t_1t_2u+(1t_1t_2)\lambda }}`$
$`=`$ $`{\displaystyle \frac{1}{u}}\left[{\displaystyle \frac{1}{2}}\mathrm{ln}^2\lambda (\mathrm{ln}(u)2)\mathrm{ln}\lambda 2\mathrm{ln}(u)+{\displaystyle \frac{1}{2}}\mathrm{ln}^2(u)+{\displaystyle \frac{5}{2}}+{\displaystyle \frac{\pi ^2}{3}}\right].`$ (178)
When we calculate the above integrals in the MG scheme, the following equations about dilogarithm function may be useful:

$`\mathrm{Li}_2(x)+\mathrm{Li}_2({\displaystyle \frac{1}{x}})`$ $`=`$ $`{\displaystyle \frac{\pi ^2}{6}}{\displaystyle \frac{1}{2}}\mathrm{ln}^2x(x>0)`$ (179)
$`\mathrm{Li}_2(x)+\mathrm{Li}_2({\displaystyle \frac{1}{x}})`$ $`=`$ $`{\displaystyle \frac{\pi ^2}{3}}{\displaystyle \frac{1}{2}}\mathrm{ln}^2xi\pi \mathrm{ln}x(x>1)`$ (180)
$`\mathrm{Li}_2(ix)+\mathrm{Li}_2({\displaystyle \frac{i}{x}})`$ $`=`$ $`{\displaystyle \frac{\pi ^2}{24}}{\displaystyle \frac{1}{2}}\mathrm{ln}^2x+{\displaystyle \frac{i}{2}}\pi \mathrm{ln}x(x>0)`$ (181)
$`\mathrm{Li}_2(ix)+\mathrm{Li}_2({\displaystyle \frac{i}{x}})`$ $`=`$ $`{\displaystyle \frac{\pi ^2}{24}}{\displaystyle \frac{1}{2}}\mathrm{ln}^2x{\displaystyle \frac{i}{2}}\pi \mathrm{ln}x(x>0)`$ (182)
$`\mathrm{Li}_2(x)+\mathrm{Li}_2(1x)`$ $`=`$ $`{\displaystyle \frac{\pi ^2}{6}}\mathrm{ln}x\mathrm{ln}(1x)`$ (183)