Two-body nonleptonic decays of the heavy mesons in the factorization approach

Abstract

In the framework of the factorization approach we calculate the branching fractions of 100 two-body nonleptonic decay channels in total, including 44 channels of the charm meson decays and 56 channels of the bottom meson decays. For charm meson decays, we test and confirm the previous observation that taking the limit for the number of colors N\to \mathrm{\infty } significantly improves theoretical predictions. For bottom meson decays, the penguin contributions are included in addition. As an essential input, we employ the weak decay form factors obtained in the framework of the relativistic quark model based on the quasi-potential approach. These form factors have well been tested by calculating observables in the semileptonic D and B meson decays and confronting obtained results with experimental data. In general, the predictions for the nonleptonic decay branching fractions are acceptable. However, for a quantitative calculation it is necessary to account for a more subtle effects of the final-state interaction.

Keywords

form factor / quark model / nonleptonic decay / factorization method

1 Introduction

Nonleptonic decays of the heavy mesons offer an environment to understand the nature of quantum chromodynamics (QCD). Experimentally many such decays have been measured. Theoretically nonleptonic decays involve more complex mechanism than the leptonic and semileptonic ones due to the local four-quark operators. A usual treatment is the factorization approach, where the decay amplitude is factorized into the product of the meson decay constant and weak transition form factors. An intuitive justification for the factorization approximation comes from the so-called colour transparency proposed by Bjorken [1], where for the energetic B decay, the final light meson flies very fast in the opposite direction to the other meson, thus almost escaping the color field of the parent particle. As a result, the factorization holds. However, the strong and complicated final state interaction does challenge the factorization approximation. And this part is very hard to be quantified, therefore, in this paper we test how factorization works for the nonleptonic decays.

The form factors embodying the dynamics of the meson weak transitions are an essential input. As a motivation and also a new point of this paper, we adopt the form factors which are derived from the relativistic quark model based on the quasi-potential approach. The numerical values of form factor parameters can be found in Refs. [2, 3], containing the results for the weak D and B decays to the pseudoscalar and/or vector mesons in the final states. In this relativistic quark model, the meson wave functions are explicitly obtained as numerical solutions of the relativistic Schrödinger-like bound-state equation and not assumed to be an empirical Gaussian function. Moreover, no free parameters are involved since they have been fixed by the previous studies of hadron spectroscopy. All relativistic effects including the transformation of meson wave functions from the rest reference frame to the moving one, and the contribution of intermediate negative energy states are included. It is important to note that the form factors are predicted in the whole kinematically allowed region. The values of these form factors have been well tested confronting with experiment by a series of the calculated semileptonic decay observables, e.g., the branching fractions, forward-backward asymmetries and polarizations. Similar work concerning the application of those form factors to the nonleptonic decays can be found in Refs. [4-6]. A very recent study of the charmless two-body B meson decays is performed in the perturbative QCD factorization approach [7] as a more advanced tool. See also Refs. [8-10] for some earlier works.

In this paper we calculate the branching fractions of the charm and bottom meson nonleptonic decays in the framework of the factorization approach based on the effective weak Hamiltonian. Experimental results from PDG and other theoretical predictions are also compiled for a direct comparison. For charm meson decays, the penguin contributions are highly suppressed, thus we neglect them. We have considered the cases of the color numbers N=3 as in reality and N\to \mathrm{\infty } . In the latter case, the discrepancy between the theory and experimental results are expected to be significantly reduced due to the experience in 1980s [14-18]. For the decay process of the B mesons, we consider both tree-level and penguin loop-level contributions. The latter can be as large as the former, or even dominant.

This paper is organized as follows. In Section 2 and Section 3 we briefly describe the effective Hamiltonian governing the D and B weak decays. In Section 4 we collect the input values. In Section 5 we show our numerical results and discuss them. Some care should also be taken for the conventions for the definitions of the decay constants and form factors. Conclusions are given in Section 6.

2 Factorization in the charm meson two-body decays

In the standard model, the effective Hamiltonian of the weak charm meson decay process reads

(1)\begin{array}{rl}{H}_{\text{eff}}=& \frac{G_F}{\sqrt{2}}{V}_{c{q}_1}^{\ast }{V}_{u{q}_2}[c_1(\mu )({\overline{q}}_{1\alpha }{c}_{\alpha }{)}_{V-A}({\overline{u}}_{\beta }{q}_{2\beta }{)}_{V-A}\\ & +c_2(\mu )({\overline{q}}_{1\alpha }{c}_{\beta }{)}_{V-A}({\overline{u}}_{\beta }{q}_{2\alpha }{)}_{V-A}],\end{array}

where both q_{1,2} can be either s or d quarks, ( {\overline{q}}_1{q}_2{)}_{V-A}= {\overline{q}}_1{\gamma }_{\mu }(1-{\gamma }_5)q_2 and G_F=1.166\times {10}^{-5} {\text{GeV}}^{-2} is the Fermi coupling constant; \alpha and \beta are color indexes; c_1(\mu ) and c_2(\mu ) are Wilson coefficients for which we use the values c_1=1.26 ,c_2= -0.51 [19] at the scale \mu = m_c. The penguin contributions are tiny and thus can be ignored.

Using the Fierz identity (i,j,k ,l are color indices)

(2)\begin{array}{r}{\delta }_{lj}{\delta }_{ik}=\frac{1}{N}{\delta }_{ij}{\delta }_{lk}+2T_{lk}^a{T}_{ij}^a,\end{array}

with T^a=\frac{{\lambda }^a}{2}, and {\lambda }^a (a=1,2 ,\cdots ,8) being the Gell-Mann matrices, and N being the number of colors, one has

(3)\begin{array}{rl}(\overline{s}s{)}_{V-A}(\overline{u}c{)}_{V-A}=& \frac{1}{N}(\overline{s}c{)}_{V-A}(\overline{u}s{)}_{V-A}\\ & +2({\overline{s}}_{\alpha }{T}_{\beta \rho }{s}_{\rho }{)}_{V-A}({\overline{u}}_{\beta }{T}_{\alpha \sigma }{c}_{\sigma }{)}_{V-A}.\end{array}

In the factorization approximation the second term of Eq. (3), which contributes to the nonfactorizable part, is neglected. Then we can write the effective Hamiltonian corresponding to the color-favored decay process

(4)\begin{array}{r}{H}_{\text{cf}}=\frac{G_F}{\sqrt{2}}{V}_{c{q}_1}^{\ast }{V}_{u{q}_2}{a}_1({\overline{q}}_{1}c{)}_{V-A}(\overline{u}{q}_2{)}_{V-A},\end{array}

and for the color-suppressed case

(5)\begin{array}{r}{H}_{\text{cs}}=\frac{G_F}{\sqrt{2}}{V}_{c{q}_1}^{\ast }{V}_{u{q}_2}{a}_2({\overline{q}}_1{q}_2{)}_{V-A}(\overline{u}c{)}_{V-A},\end{array}

with a_1=c_1+\frac{c_2}{N},a_2=c_2+\frac{c_1}{N}. Empirically, taking the choice of N\to \mathrm{\infty } will generally improve theoretical predictions for the charm meson decays, as it was already mentioned in Introduction. In this sense, part of the nonfactorizable effects have been compensated by the choice of N\to \mathrm{\infty }.

In the factorization approach, the hadronic matrix element can be expressed by the product of decay constant and the invariant form factors. The decay constant is defined as the matrix element of the weak current between the vacuum and a pseudoscalar (P) or a vector (V) meson:

(6)\begin{array}{rl}& ⟨P(p_{\mu })|({\overline{q}}_1{q}_2{)}_{V-A}|0⟩=-\mathrm{i}{f}_P{p}_{\mu },\\ & ⟨V|({\overline{q}}_1{q}_2{)}_{V-A}|0⟩=\mathrm{i}{f}_V{m}_V{ϵ}_{\mu }^{\ast },\end{array}

where f_P and f_V are the decay constants of the pseudoscalar and vector meson, respectively; m_V and {ϵ}_{\mu } are the mass and polarization vector of the vector meson.

For the D\to P transition (with the momenta p_D,p_P and masses m_D,m_P of the initial and final mesons, respectively), the matrix element of the weak current is parameterized as

(7)\begin{array}{rl}& ⟨P|(\overline{q}{\gamma }^{\mu }c)|D⟩=f_{+}(q^2)\left(p_D^{\mu }+p_P^{\mu }-\frac{m_D^2-m_P^2}{q^2}{q}^{\mu }\right)\\ & +f_0(q^2)\frac{m_D^2-m_P^2}{q^2}{q}^{\mu },\\ & ⟨P|(\overline{q}{\gamma }^{\mu }{\gamma }^{5}c)|D⟩=0.\end{array}

For the D\to V transition,

(8)\begin{array}{rl}& ⟨V|(\overline{q}{\gamma }^{\mu }c)|D⟩=\frac{2\mathrm{i}V(q^2)}{m_D+m_V}{\epsilon }^{\mu \nu \rho \sigma }{ϵ}_{\nu }^{\ast }{p}_{D\rho }{p}_{V\sigma },\\ & ⟨V|(\overline{q}{\gamma }^{\mu }{\gamma }_{5}c)|D⟩=2m_V{A}_0(q^2)\frac{{ϵ}^{\ast }\cdot q}{q^2}{q}^{\mu }\\ & +(m_D+m_V)A_1(q^2)\left({ϵ}^{\ast \mu }-\frac{{ϵ}^{\ast }\cdot q}{q^2}{q}^{\mu }\right)\\ & -A_2(q^2)\frac{{ϵ}^{\ast }\cdot q}{m_D+m_V}\\ & \times \left(p_D^{\mu }+p_V^{\mu }-\frac{m_D^2-m_V^2}{q^2}{q}^{\mu }\right).\end{array}

In these equations, q=p_D-p_{P,V} is the four-momentum transfer between the initial D and final P or V mesons. For the case of the nonleptonic two-body decay considered below, q is just the on-shell momentum of the meson created from vacuum.

The matrix element X_{D\to M_1,M_2}=⟨M_1|(\overline{q}c{)}_{V-A}|D⟩⟨ M_2|(\overline{u}q{)}_{V-A}|0⟩ is then simplified as

(9)\begin{array}{rl}& M_2=P:X_{D\to M_1,P}=-\mathrm{i}{f}_P{m}_P{ϵ}_{\lambda }^{†\mu }⟨M_1|(\overline{q}c{)}_{V-A}|D⟩,\\ & M_2=V:X_{D\to M_1,V}=\mathrm{i}{f}_V{m}_V{ϵ}_{\lambda }^{†\mu }⟨M_1|(\overline{q}c{)}_{V-A}|D⟩,\end{array}

where we adopt the convention where the final meson M_2 (P or V) after the comma in the subscript of X is generated from vacuum.

In the rest frame of the initial D meson, one has the explicit representations of the momentum and polarization vectors:

(10)\begin{array}{rl}& p_D^{\mu }=(m_D,0,0,0),p_{M_1}^{\mu }=(E_1,0,0,|\mathit{p}|),\\ & q^{\mu }=(E_2,0,0,-|\mathit{p}|),{ϵ}_t^{†\mu }=\frac{1}{\sqrt{q^2}}(E_2,0,0,-|\mathit{p}|),\\ & {ϵ}_{\pm }^{†\mu }=\frac{1}{\sqrt{2}}(0,\pm 1,\mathrm{i},0),{ϵ}_0^{†\mu }=\frac{1}{\sqrt{q^2}}(|\mathit{p}|,0,0,-E_2),\end{array}

where E_1 is the energy of M_1 and E_1+E_2=m_D, |\mathit{p}|={\lambda }^{1/2}( m_D^2,m_1^2,m_2^2)/(2m_D ) is the momentum of the daughter meson with \lambda (x,y,z)=x^2 +y^2+ z^2-2(xy +yz+x z). For convenience, we define the helicity amplitudes

(11)H_{\lambda }\equiv {ϵ}_{\lambda }^{†\mu }⟨M_1|(\overline{q}c{)}_{V-A}|D⟩,

with \lambda =t for M_2=P and \lambda = \pm ,0 for M_2=V.

For the process D\to P_1,P_2, one has

(12)\begin{array}{rl}{X}_{D\to P_1,P_2}=& -\mathrm{i}{f}_{P_2}{m}_{P_2}{ϵ}_t^{†\mu }⟨P_1|(\overline{q}c{)}_{V-A}|D⟩\\ =& -\mathrm{i}{f}_{P_2}(m_D^2-m_{P_1}^2)f_0(m_{P_2}^2)\\ =& -\mathrm{i}{f}_{P_2}{m}_{P_2}{H}_t;\end{array}

for D\to P,V,

(13)\begin{array}{rl}{X}_{D\to P,V}=& \mathrm{i}{f}_V{m}_V{ϵ}_{\lambda }^{†\mu }⟨P|(\overline{q}c{)}_{V-A}|D⟩\\ =& \mathrm{i}2f_V{f}_{+}(m_V^2)m_D|\mathit{p}|\\ =& \mathrm{i}{f}_V{m}_V{H}_0.\end{array}

The expressions for H_0 and H_t coincide with the ones given in Ref. [20] for the D\to P transition.

For D\to V,P, one has

(14)\begin{array}{rl}{X}_{D\to V,P}=& -\mathrm{i}{f}_P{m}_P{ϵ}_t^{†\mu }⟨V|(\overline{q}c{)}_{V-A}|D⟩\\ =& 2\mathrm{i}{f}_P{A}_0(m_P^2)m_D|\mathit{p}|\\ =& -\mathrm{i}{f}_P{m}_P{H}_t;\end{array}

for D\to V_1,V_2,

(15)\begin{array}{rl}& |X_{D\to V_1,V_2}{|}^2=f_{V_2}^2{m}_{V_2}^2(|H_{+}{|}^2+|H_{-}{|}^2+|H_0{|}^2),\\ & H_{\pm }=-(m_D+m_{V_1})A_1(m_{V_2}^2)\pm \frac{2m_D|\mathit{p}|}{m_D+m_{V_1}}V(m_{V_2}^2),\\ & H_0=-(m_D+m_{V_1})A_1(m_{V_2}^2)\frac{m_D^2-m_{V_1}^2-m_{V_2}^2}{2m_{V_1}{m}_{V_2}}\\ & +\frac{2m_D^2|\mathit{p}{|}^2}{(m_D+m_{V_1})m_{V_1}{m}_{V_2}}{A}_2(m_{V_2}^2),\end{array}

where the definitions of H_t,H_{\pm },H_0 coincide with the ones given in Ref. [20] for the D\to V transition.

3 Factorization of the bottom meson two-body decay amplitudes

We classify the bottom meson decay channels into two classes according to the effective Hamiltonian. For \mathrm{\Delta } B=\pm 1,\mathrm{\Delta }C=\pm 1 transitions, e.g., the processes \overline{b}\to \overline{c}u \overline{d}, \overline{b}\to \overline{c}u\overline{s}, \overline{b}\to \overline{u}c\overline{d}, and \overline{b}\to \overline{u}c\overline{s}, the effective Hamiltonian reads

(16)H_{\text{eff}}= {\displaystyle \frac{G_F}{\sqrt{2}}}[V_{qb}^{\ast }{V}_{q_1{q}_2}(c_1 (\mu )O_1+ c_2(\mu )O_2)].

Then such category is similar to charm decays described above.

For \mathrm{\Delta }B=\pm 1,\mathrm{\Delta }C=0 transitions, e.g., \overline{b}\to \overline{c}c \overline{d}, \overline{b}\to \overline{c}c\overline{s}, \overline{b}\to \overline{u}u\overline{d} and \overline{b}\to \overline{u}u\overline{s}, the effective Hamiltonian reads

(17)\begin{array}{rl}{H}_{\text{eff}}=& \frac{G_F}{\sqrt{2}}\{V_{q^{\prime }b}^{\ast }{V}_{q^{\prime }q}[c_1(\mu )O_1+c_2(\mu )O_2]\\ & -V_{tb}^{\ast }{V}_{tq}\sum _{i=3}^{10}{c}_i(\mu )O_i\}.\end{array}

where q(q^{\prime })=s ,d,c;

(18)\begin{array}{rl}& O_1=(\overline{b}{q}^{\prime }{)}_{V-A}({\overline{q}}^{\prime }q{)}_{V-A},O_2=({\overline{b}}_{\alpha }{q}_{\beta }^{\prime }{)}_{V-A}({\overline{q}}_{\beta }^{\prime }{q}_{\alpha }{)}_{V-A},\\ & O_3=(\overline{b}q{)}_{V-A}\sum _{q^{\prime }}({\overline{q}}^{\prime }{q}^{\prime }{)}_{V-A},O_4=({\overline{b}}_{\alpha }{q}_{\beta }{)}_{V-A}\sum _{q^{\prime }}({\overline{q}}_{\beta }^{\prime }{q}_{\alpha }^{\prime }{)}_{V-A},\\ & O_5=(\overline{b}q{)}_{V-A}\sum _{q^{\prime }}({\overline{q}}^{\prime }{q}^{\prime }{)}_{V+A},O_6=({\overline{b}}_{\alpha }{q}_{\beta }{)}_{V-A}\sum _{q^{\prime }}({\overline{q}}_{\beta }^{\prime }{q}_{\alpha }^{\prime }{)}_{V+A},\\ & O_7=\frac{3}{2}(\overline{b}q{)}_{V-A}\sum _{q^{\prime }}{e}_{q^{\prime }}({\overline{q}}^{\prime }{q}^{\prime }{)}_{V+A},O_8=\frac{3}{2}({\overline{b}}_{\alpha }{q}_{\beta }{)}_{V-A}\sum _{q^{\prime }}{e}_{q^{\prime }}({\overline{q}}_{\beta }^{\prime }{q}_{\alpha }^{\prime }{)}_{V+A},\\ & O_9=\frac{3}{2}(\overline{b}q{)}_{V-A}\sum _{q^{\prime }}{e}_{q^{\prime }}({\overline{q}}^{\prime }{q}^{\prime }{)}_{V-A},O_{10}=\frac{3}{2}({\overline{b}}_{\alpha }{q}_{\beta }{)}_{V-A}\sum _{q^{\prime }}{e}_{q^{\prime }}({\overline{q}}_{\beta }^{\prime }{q}_{\alpha }^{\prime }{)}_{V-A},\end{array}

and e_{q^{\prime }} is the charge of the q^{\prime } quark.

The even operators O_{2-10} can be rearranged to a color singlet form by the Fierz transformation

(19)\begin{array}{r}({\overline{\psi }}_1{O}^i{\psi }_2)({\overline{\psi }}_3{O}_i{\psi }_4)=\sum _j{C}_{ij}({\overline{\psi }}_1{O}^j{\psi }_4)({\overline{\psi }}_3{O}_j{\psi }_2),\end{array}

where C_{ij} are the Fierz coefficients that are presented in Tab.1. In this way, we have

(20)\begin{array}{rl}{O}_2=& (\overline{b}q{)}_{V-A}({\overline{q}}^{\prime }{q}^{\prime }{)}_{V-A},\\ O_4=& \sum _{q^{\prime }}(\overline{b}{q}^{\prime }{)}_{V-A}({\overline{q}}^{\prime }q{)}_{V-A},\\ O_6=& -2\sum _{q^{\prime }}(\overline{b}{q}^{\prime }{)}_{S+P}({\overline{q}}^{\prime }q{)}_{S-P},\\ O_8=& -2\sum _{q^{\prime }}\frac{3}{2}(\overline{b}{q}^{\prime }{)}_{S+P}({\overline{q}}^{\prime }q{)}_{S-P},\\ O_{10}=& \sum _{q^{\prime }}\frac{3}{2}(\overline{b}{q}^{\prime }{)}_{V-A}{e}_{q^{\prime }}({\overline{q}}^{\prime }q{)}_{V-A}.\end{array}

In the above equations, ({\overline{q}}_1{q}_2{)}_{V+A}\equiv {\overline{q}}_1{\gamma }_{\mu }( 1+{\gamma }_5)q_2 and ({\overline{q}}_1{q}_2{)}_{S\pm P}\equiv {\overline{q}}_1(1\pm {\gamma }_5)q_2.

Here we take B\to P_1(q_s{q}^{\prime })P_2(q^{\prime }q) (q_s is the spectator quark), as an example, to demonstrate the calculation of the penguin contribution, with P_2 representing a charged pseudoscalar meson. The contribution of O_4 therein is proportional to a_4 and the contribution of O_{10} is proportional to \frac{3}{2}{e}_{q^{\prime }}{a}_{10}. The coefficient a_i is related to the Wilson one:

(21)\begin{array}{rl}& \text{for}{a}_{\text{odd}},a_i=c_i+\frac{c_{i+1}}{N},\\ & \text{for}{a}_{\text{even}},a_i=c_i+\frac{c_{i-1}}{N}.\end{array}

The operator O_6 can be further written as

(22)\begin{array}{rl}{O}_6=& -2\sum _{q^{\prime }}(\overline{b}{q}^{\prime }{)}_{S+P}({\overline{q}}^{\prime }q{)}_{S-P}\\ =& -2[(\overline{b}{q}^{\prime })({\overline{q}}^{\prime }q)+(\overline{b}{\gamma }^5{q}^{\prime })({\overline{q}}^{\prime }q)\\ & -(\overline{b}{q}^{\prime })({\overline{q}}^{\prime }{\gamma }^{5}q)-(\overline{b}{\gamma }^5{q}^{\prime })({\overline{q}}^{\prime }{\gamma }^{5}q)].\end{array}

Parity conservation leads to

(23)\begin{array}{l}⟨P|({\overline{q}}_1{\gamma }^{\mu }{q}_2)|0⟩=0,⟨P|(\overline{q}{\gamma }^{\mu }{\gamma }^{5}b)|B⟩=0,\end{array}

while equations of motion read (with m_{1,2} being the masses of quarks q_{1,2})

(24)\begin{array}{rl}& ({\overline{q}}_1{\gamma }^5{q}_2)=\frac{-\mathrm{i}}{m_1+m_2}{\mathrm{\partial }}_{\mu }({\overline{q}}_1{\gamma }^{\mu }{\gamma }^5{q}_2),\\ & ({\overline{q}}_1{q}_2)=\frac{-\mathrm{i}}{m_1-m_2}{\mathrm{\partial }}_{\mu }({\overline{q}}_1{\gamma }^{\mu }{q}_2),\end{array}

and thus only the third term of Eq. (22) survives. Then according to Eq. (24),

(25)\begin{array}{rl}& ⟨P_1|(\overline{b}{q}^{\prime })|B⟩=\frac{-\mathrm{i}}{m_b-m_{q^{\prime }}}(-\mathrm{i}{q}_{\mu })⟨P_1|(\overline{b}{\gamma }^{\mu }{q}^{\prime })|B⟩,\\ & ⟨P_2|({\overline{q}}^{\prime }{\gamma }^{5}q)|0⟩=\frac{-\mathrm{i}}{m_q+m_{q^{\prime }}}(\mathrm{i}{q}^{\mu })⟨P_2|(\overline{b}{\gamma }^{\mu }{\gamma }^5{q}^{\prime })|0⟩,\end{array}

and the product is given by

(26)\begin{array}{rl}& ⟨P_1|(\overline{b}{q}^{\prime })|B⟩⟨P_2|({\overline{q}}^{\prime }{\gamma }^{5}q)|0⟩\\ & =-\frac{m_{P_2}^2}{(m_q+m_{q^{\prime }})(m_b-m_{q^{\prime }})}\\ & \times ⟨P_1|(\overline{b}{\gamma }^{\mu }{q}^{\prime })|B⟩⟨P_2|(\overline{b}{\gamma }^{\mu }{\gamma }^5{q}^{\prime })|0⟩\\ & =\frac{m_{P_2}^2}{(m_q+m_{q^{\prime }})(m_b-m_{q^{\prime }})}{X}_{B\to P_1,P_2}.\end{array}

Therefore, the contribution of O_6 is proportional to \frac{2m_{P_2}^2}{(m_q+m_{q^{\prime }})(m_b-m_{q^{\prime }})}{a}_6. And similarly, the contribution of O_8 is proportional to \frac{3}{2}{e}_{q^{\prime }}\frac{2m_{P_2}^2}{(m_q+m_{q^{\prime }})(m_b-m_{q^{\prime }})}. In Tab.2 we summarize the total penguin contributions in various processes. In our convention, the second meson (M_2) corresponds to the one generated from vacuum; and in the lower half of this table, M_2 is flavor neutral. The values of Wilson coefficients at the scale \mu =m_b used in our calculation are c_1=1.105, c_2=-0.228, c_3=0.013, c_4=- 0.029, c_5=0.009, c_6=- 0.033, c_7/\alpha =0.005, c_8/\alpha =0.060, c_9/\alpha =-1.283, c_{10}/\alpha =0.266 [21], where \alpha is the fine structure constant.

From Tab.2 one can easily read out the amplitude for a given decay process. However, some decay amplitudes contain more than one class. We take B^{+}\to {\pi }^{+}\eta as an example. In this decay either \eta or \pi can be produced from vacuum. For the case when \eta is produced from vacuum, the decay amplitude can be written as follows:

(27)\begin{array}{rl}& q^{\prime }=u,q=d,e_{q^{\prime }}=\frac{2}{3}:A_1=\frac{G_F}{\sqrt{2}}\{V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}[a_3-a_5+a_9-a_7]\}{X}_{B^{+}\to {\pi }^{+},{\eta }_u},\\ & q^{\prime }=d,q=d,e_{q^{\prime }}=-\frac{1}{3}:A_2=\frac{G_F}{\sqrt{2}}\{-V_{tb}^{\ast }{V}_{td}[a_3-a_5-\frac{1}{2}{a}_9+\frac{1}{2}{a}_7+a_4-\frac{1}{2}{a}_{10}\\ & +\frac{m_{\eta }^2}{m_s(m_b-m_d)}(a_6-\frac{1}{2}{a}_8)(\frac{f_{\eta }^s}{f_{\eta }^u}-1)r_{\eta }]\}{X}_{B^{+}\to {\pi }^{+},{\eta }_u},\\ & q^{\prime }=s,q=d,e_{q^{\prime }}=-\frac{1}{3}:A_3=\frac{G_F}{\sqrt{2}}\{-V_{tb}^{\ast }{V}_{td}[a_3-a_5-\frac{1}{2}(a_9-a_7)]\}{X}_{B^{+}\to {\pi }^{+},{\eta }_s},\\ & X_{B^{+}\to {\pi }^{+},{\eta }_u}=⟨{\pi }^{+}|(\overline{b}d{)}_{V-A}|B^{+}⟩⟨\eta |(\overline{u}u{)}_{V-A}|0⟩=⟨{\pi }^{+}|(\overline{b}d{)}_{V-A}|B^{+}⟩f_{\eta }^u,\\ & X_{B^{+}\to {\pi }^{+},{\eta }_s}=⟨{\pi }^{+}|(\overline{b}d{)}_{V-A}|B^{+}⟩⟨\eta |(\overline{s}s{)}_{V-A}|0⟩=⟨{\pi }^{+}|(\overline{b}d{)}_{V-A}|B^{+}⟩f_{\eta }^s.\end{array}

The definition of r_{\eta } is given later in Eq. (33). For the case when \pi is produced from vacuum, the decay amplitude reads

(28)\begin{array}{rl}& q^{\prime }=u,q=d,e_{q^{\prime }}=\frac{2}{3}:\\ & A_4=\frac{G_F}{\sqrt{2}}\{V_{ub}^{\ast }{V}_{ud}{a}_1-V_{tb}^{\ast }{V}_{td}[a_4+a_{10}+\frac{2m_{\pi }^2}{(m_u+m_d)(m_b-m_u)}(a_6+a_8)]\}{X}_{B^{+}\to \eta ,{\pi }^{+}},\\ & X_{B^{+}\to \eta ,{\pi }^{+}}=⟨\pi |(\overline{u}d{)}_{V-A}|0⟩⟨\eta |(\overline{b}u{)}_{V-A}|B⟩.\end{array}

The total amplitude for B^{+}\to {\pi }^{+}\eta is then given by the sum of these amplitudes

(29)\begin{array}{r}A(B\to {\pi }^{+}\eta )=A_1+A_2+A_3+A_4.\end{array}

4 The input

In our consideration we use the following quark compositions of the light mesons

(30)\begin{array}{rl}& K^{+}=u\overline{s},K^0=d\overline{s},K^{-}=s\overline{u},\\ & {\pi }^{+}({\rho }^{+})=u\overline{d},{\pi }^0({\rho }^0)=\frac{u\overline{u}-d\overline{d}}{\sqrt{2}},{\pi }^{-}({\rho }^{-})=d\overline{u},\\ & {\eta }_0=\frac{d\overline{d}+u\overline{u}+s\overline{s}}{\sqrt{3}},{\eta }_8=\frac{d\overline{d}+u\overline{u}-2s\overline{s}}{\sqrt{6}},\\ & \eta ={\eta }_8\cos \theta -{\eta }_0\sin \theta ,{\eta }^{\prime }={\eta }_8\sin \theta +{\eta }_0\cos \theta ,\end{array}

with \theta =-{15.4}^{\circ }, which corresponds to the mixing angle \varphi = {39.3}^{\circ } [22]. Note that such value of \varphi was previously used in Refs. [23,24]. This value of \varphi agrees with the CLEO measurement {42}^{\circ }\pm {2.8}^{\circ } [25] and also with the recent BESIII measurement {40.1}^{\circ }\pm {2.1}^{\circ }\pm {0.7}^{\circ } [26]. Reference [22] presents a nice analysis of the \eta - {\eta }^{\prime } mixing both from the theoretical and phenomenological standpoint, where in the former only the masses of pseudoscalar mesons are involved as inputs and in the latter the experimental measurements of branching fractions are used. The relation between the decay constants f_{{\eta }^{{(}^{\prime })}}^u and f_{{\eta }^{{(}^{\prime })}}^s in the singlet-octet mixing scheme and the ones f_q and f_s in the quark flavor basis q\overline{q}= \frac{1}{\sqrt{2}}(u\overline{u}+d\overline{d}) and s\overline{s} is given by

(31)\begin{array}{rl}{f}_{\eta }^u& =\frac{1}{\sqrt{2}}{f}_q\cos \varphi ,f_{\eta }^s=-f_s\sin \varphi ,\\ f_{{\eta }^{\prime }}^u& =\frac{1}{\sqrt{2}}{f}_q\sin \varphi ,f_{{\eta }^{\prime }}^s=f_s\cos \varphi ,\end{array}

with f_q/f_{\pi }=1.07,f_s/f_{\pi }= 1.34. The values of the decay constants used in our calculations [27-33] are as follows (in MeV)

(32)\begin{array}{rl}& f_{\pi }=130.2,f_K=155.6,f_{K^{\ast }}=217,\\ & f_{\eta }^u=78,f_{\eta }^s=-112,f_{{\eta }^{\prime }}^u=63,f_{{\eta }^{\prime }}^s=137,\\ & f_{D^{+}}=212.7,f_{D^0}=211.6,f_{D_s}=249.9,\\ & f_{\rho }=205,f_{\omega }=187,f_{\varphi }=215.\end{array}

Calculating the matrix elements of the scalar and pseudoscalar currents, one needs to use the equations of motion, Eq. (24). When {\eta }^{{(}^{\prime })} is generated from vacuum, the hadron matrix element is treated differently due to the SU(3) breaking [34,35]:

(33)\begin{array}{rl}& ⟨{\eta }^{{(}^{\prime })}|\overline{s}{\gamma }^{5}s|0⟩=-\mathrm{i}\frac{m_{{\eta }^{{(}^{\prime })}}}{2m_s}(f_{{\eta }^{{(}^{\prime })}}^s-f_{{\eta }^{{(}^{\prime })}}^u),\\ & ⟨{\eta }^{{(}^{\prime })}|\overline{u}{\gamma }^{5}u|0⟩=⟨{\eta }^{{(}^{\prime })}|\overline{d}{\gamma }^{5}d|0⟩=r_{{\eta }^{{(}^{\prime })}}⟨{\eta }^{{(}^{\prime })}|\overline{s}{\gamma }^{5}s|0⟩,\\ & r_{{\eta }^{{}^{\prime }}}=\frac{\sqrt{2f_0^2-f_8^2}}{\sqrt{2f_8^2-f_0^2}}\left(\frac{\cos \theta +{\displaystyle \frac{1}{\sqrt{2}}}\sin \theta }{\cos \theta -\sqrt{2}\sin \theta }\right),\\ & r_{\eta }=-\frac{1}{2}\frac{\sqrt{2f_0^2-f_8^2}}{\sqrt{2f_8^2-f_0^2}}\left(\frac{\cos \theta -\sqrt{2}\sin \theta }{\cos \theta +{\displaystyle \frac{1}{\sqrt{2}}}\sin \theta }\right),\end{array}

with f_0/f_{\pi }=1.17 and f_8/f_{\pi }=1.26. The axial-vector anomaly effect has been incorporated into this equation in order to ensure the correct behavior in the chiral limit. By using Eq. (2.12) and Eq. (2.18) from Ref. [22], we have

(34)\begin{array}{rl}⟨0|\overline{u}{\gamma }^{5}u|\eta ⟩& =-\frac{\mathrm{i}}{\sqrt{2}}\frac{m_{\pi }^2}{2m_u}{f}_q\cos \varphi ,\\ ⟨0|\overline{s}{\gamma }^{5}s|\eta ⟩& =-\mathrm{i}\frac{2m_K^2-m_{\pi }^2}{2m_s}{f}_s\sin \varphi .\end{array}

Considering the fact that in the chiral limit

(35){\displaystyle \frac{m_{\pi }^2}{2m_u}}={\displaystyle \frac{2m_K^2-m_{\pi }^2}{2m_s}},

we arrive at Eq. (33). Note also that in the limit of f_0=f_8, Eq. (33) reproduces Eq. (19) of Ref. [36].

The running quark masses at the scale \mu =m_b have the following values [37]

(36)m_u=1.86 ,m_d=4.22, m_c=901,m_s =80,m_b =4200,

in units of MeV. For the CKM matrix we use the Wolfenstein parameterization

(37)\left(\begin{array}{ccc}1-{\displaystyle \frac{{\lambda }^2}{2}}& \lambda & A{\lambda }^3(\rho -\mathrm{i}\eta )\\ -\lambda & 1-{\displaystyle \frac{{\lambda }^2}{2}}& A{\lambda }^2\\ A{\lambda }^3(1-\rho -\mathrm{i}\eta )& -A{\lambda }^2& 1\end{array}\right),

with central values \lambda =0.2265, A=0.790, \overline{\rho } =0.141 and \overline{\eta }=0.357 taken from PDG [38].

We employ the form factor values from Refs. [2, 3] calculated in RQM, which have been well tested in the semileptonic decays. These form factors are in agreement with lattice determination, and the resulting observables (not only the branching fractions but also the forward-backward asymmetries, polarizations of the leptons or the vector mesons), agree with lattice and experimental results. As the function of momentum transfer squared, the relevant form factors are expressed by

• f_{+}(q^2), V(q^2), A_0(q^2):

(38)F(q^2)={\displaystyle \frac{F(0)}{\left(1-{\displaystyle \frac{q^2}{M^2}}\right)\left[1-{\sigma }_1{\displaystyle \frac{q^2}{M_1^2}}+{\sigma }_2{\left({\displaystyle \frac{q^2}{M_1^2}}\right)}^2\right]}},

• f_0(q^2), A_1(q^2 ),A_2(q^2 ):

(39)F(q^2)= {\displaystyle \frac{F(0)}{1-{\sigma }_1{\displaystyle \frac{q^2}{M_1^2}}+{\sigma }_2{\left({\displaystyle \frac{q^2}{M_1^2}}\right)}^2}}.

For the c\to s transition, M=M_{D_s^{\ast }}=2.112 GeV for the form factors f_{+}(q^2), V(q^2), and M=M_{D_s}=1.968 GeV for the form factor A_0(q^2). For the c\to d transition, M=M_{D^{\ast }}=2.010 GeV for the form factors f_{+}(q^2), V(q^2), and M=M_D=1.870 GeV for the form factor A_0(q^2). For the b\to c transition, M=M_{B_c^{\ast }}=6.332 GeV for the form factors f_{+}(q^2), V(q^2), and M=M_{B_c}=6.227 GeV for the form factor A_0(q^2). For the b\to u transition, M=M_{B^{\ast }}=5.325 GeV for the form factors f_{+}(q^2),V(q^2), and M=M_B=5.280 GeV for the form factor A_0(q^2). The mass M_1 is always taken as the pole mass between the active quarks: M_1=M_{D_s^{\ast }} for c\to s transition, M_1=M_{D^{\ast }} for c\to d transition, M_1=M_{B_c^{\ast }} for b\to c transition, M_1=M_{B^{\ast }} for b\to u transition. For convenience, we compile these mass parameters for the charm meson decays in Tab.3 while for the bottom meson case one refers to Tab.1 in Ref. [3]. The values of F(0), {\sigma }_1, {\sigma }_2 are easily found in Refs. [2,3].

5 Results and discussion

5.1 Branching fractions

The decay branching fractions can be calculated by the equation

(40)$$\begin{gathered} \mathcal{B}=\tau {\displaystyle \frac{|\mathit{p}|}{8\pi m^2}}|A{|}^2, \\ \end{gathered}$$

where for the two-body nonleptonic decays, \tau and m are the lifetime and mass of the parent particle, respectively, |\mathit{p}| is the magnitude of the three-momentum of the final mesons in the rest frame of the decaying heavy meson, and expressions for the amplitudes A are given in Appendices A and B. The involved expressions for X are given in Eqs. (12)−(15) for PP,PV ,VP, and VV modes, respectively. For VV modes, the helicities H_{+},H_{-} and H_0 are involved. The results for the branching fractions are shown in Tab.4 and Tab.5 for charm and bottom meson decays, respectively.

For charm meson decays, both N=3 (the number of colors in reality) and the limit N\to \mathrm{\infty } are considered in Tab.4. As mentioned above, the case of N\to \mathrm{\infty } compensates the nonfactorizable effects to some extent and is expected to improve the theoretical predictions empirically. We confirm this point, e.g., for the channels D^{+}\to {\pi }^0{\pi }^{+}, D^0\to K^{+}{K}^{-}, D^0\to K^{-} {\rho }^{+} the results for N\to \mathrm{\infty } are altered by a factor of about 2 compared to the ones for N=3, improving the agreement with the experimental values. For the channels D^{+}\to {\pi }^{+}\varphi, D^0\to \eta \eta, D_s\to K^{+}{\overline{K}}^0 and D_s\to K^{+}{\pi }^0, the effect is even more pronounced, the results are changed by one or two orders of magnitude compared to the ones for N=3 bringing them closer to the measured values. In Ref. [39], a more elaborate phenomenological analysis is performed, where the annihilation and exchange contributions as well as the resonant final-state interaction (FSI) are considered. As a result, the branching fractions for D\to KK and D\to \pi \pi as a long-standing puzzle get correctly treated in Ref. [39] compared to the experimental values. We find in our simple treatment that only the value for the D^{+}\to {\pi }^0{\pi }^{+} channel agrees with the experimental value within 2 standard deviation while for the D^0\to {\pi }^{-} {\pi }^{+},{\pi }^0{\pi }^0 ones the branching fractions differ from experiments by a factor of 2−3. This is in line with the observation of Ref. [39] showing the importance of the nonfactorizable effects.

Here we discuss the rule of N\to \mathrm{\infty } in more detail. The phenomenon that this rule greatly improves predictions for branching fractions of the nonleptonic two-body D decays was realized by the community in 1980s, as shown in Refs. [14-16]. In Ref. [16], Buras et al. made a more complete analysis of charm decays, where the effectiveness of N\to \mathrm{\infty } is clearly demonstrated compared to the case of N=3, and also the result of the 1/ N expansion is phrased much better in terms of simple diagrammatical rules. But we stress that this rule is purely empirical. As mentioned in Ref. [14], it is not clear whether this rule is just a coincidence or has a deeper meaning. Note that the generalization of the N\to \mathrm{\infty } to the B decays will lead to predictions in contradiction with experiment. Also, the 1/N suppression varies in different channels and is rather of a dynamical origin.

In cases where the rule of N\to \mathrm{\infty } works well, we can understand what happens for the factorizable and nonfactorizable contributions in the spirit of the large N QCD [16]. In the usual procedure, the 1/N term in Eq. (21), being part of the factorizable term, is kept while the nonfactorizable term is not considered since there is no reliable way to calculate it. In such a situation the leading and nonleading 1 /N contributions mix up. The nonfactorizable one, e.g., the final state interaction effect, is nonleading in the 1/N expansion. By dropping the 1/ N term in Eq. (21), one will work in a self-consistent expansion of 1/ N. Or we can say that the 1/N term in factorizable part is almost compensated by the (unknown) nonfactorizable one. There is an explicit calculation to demonstrate this point [17], where the author shows that the soft gluon exchange mechanism (a type of nonfactorizable contribution) tends to cancel the 1/N term in the factorized amplitude by using the light cone sum rule. In a more physical picture, we can say that the quarks belonging to different color singlet currents do not easily form a meson and thus the 1 /N term is highly suppressed.

The results for the B decays are shown in Tab.5. The theoretical predictions should be better consistent with the experimental data than in the D meson case. Indeed, the factorization assumption works better for the heavier B meson since the final mesons carry larger momenta. And for some decay channels, such as B^{+}\to {\rho }^{+}{\eta }^{{(}^{\prime })} and B_s\to D_s^{\ast -}{\rho }^{+}, the results for N=3 perfectly match the experimental values within 1.5\sigma uncertainty. We have calculated branching fractions for the three sets of color number N. Those results constitute a range of branching fractions varying with the choices of N, which may be understood as an error estimate in some sense. However, there is an exception, for the penguin governed decay B^{+}\to {\pi }^{+}\varphi the result for N=2 deviates from the one for N=3 by two orders of magnitude. We should compare our results to the ones given in Ref. [40] since we work in the same framework, considering the tree-level as well as penguin contributions. However, in Ref. [40] a different set of Wilson coefficients (known as the generalized factorization) is used. Besides, we employ the form factor values predicted by our relativistic quark model, as a more advanced tool from today's perspective compared to their BSW ones. Our results for N=3 are of similar magnitude with Ref. [40] under the same condition N_c^{\text{eff}}(LL )=N_c^{\text{eff}}(LR) =3. For most of channels, our ranges of branching fractions formed by N=2,3,\mathrm{\infty } are close to the corresponding experimental values within 2\sigma uncertainty. But there are a few channels where results differ from experimental ones by larger or around factor of 5. Then we also compare with the predictions of Refs. [40,47], and find such deviations also happens in their results. In general, for the color-suppressed decay channels (involving a_2 terms), such as B^0\to {\pi }^0{\pi }^0, B^0\to {\pi }^0{\rho }^0, B^0\to \omega \omega and B^0\to {\rho }^0{\rho }^0, the predicted branching fractions are lower than the experiment values. One of the reasons is due to the smallness of a_2, but more importantly, the strong FSI effects should play an essential role, as has been explicitly demonstrated in Ref. [40] in detail. In fact, as we know, the interaction between pseudoscalar octet, e.g., the \pi \pi -K\overline{K} system, is very strong, for which some of our recent investigations can be found in Refs. [48-52]. That is, the B\to \pi \pi decay will receive large contributions from the intermediate states K\overline{K} and \eta \eta etc. On another hand, it has been found in Ref. [53] that the spacelike penguin contributions may be sizable in B\to PP decays, where the authors showed that such corrections to the branching fraction for B\to \pi \pi may be more than 100%. However, in Ref. [54] the authors assume that such contributions in B\to PV,V V decays are not as severe as in B\to PP. Reference [40] provides a careful examination but those effects of FSI and spacelike penguins can not be reliably determined yet. So conservatively speaking, the branching fraction can be trusted by its order of magnitude.

In this paragraph we give a few comments on comparison of our results with the ones in Ref. [40] by Cheng et al. In this reference many sets of numbers for the branching fraction values are calculated, and these numbers constitute an interval. Such an interval may contain the experimental value, which is very encouraging. But in some cases their ranges span two orders of magnitude. For example, for the B\to \omega \omega decay it reads 7\times {10}^{-8}-2\times {10}^{-6}, and for the B\to {\rho }^0{\rho }^0 5\times {10}^{-8}-2.57 \times {10}^{-6}. Their preferred values correspond to N_c^{\text{eff}}(LL) =2 and N_c^{\text{eff}}(LR) =5. For the same value for the color number, the differences between results of Ref. [40] and ours mainly come from the different inputs. Especially, in Ref. [40] complex-valued numbers for the set of the Wilson coefficients are used while we use the real-valued ones from Buchalla et al. in Ref. [21]. Our main goal is not to reproduce the experimental values exactly or to match them well. We want to test the factorization hypothesis by using our most recent form factor values calculated from an advanced relativistic quark model. To get a more quantitative calculation, the nonfactorizable contribution should be included anyway.

As is known and also mentioned earlier, the nonfactorizable effects may dominate in a specific decay, and there is currently no method to calculate them beforehand. However, in literature there are important works dedicated to the analysis of such nonfactorizable effects by confronting with experimental data. One typical example is the factorization-assisted topological (FAT) approach [41-44] which combines the naive factorization hypothesis and the topological diagram approach [45,46]. In these papers the authors identify the possible sources of nonfactorizable contributions and then parameterize them, in order to fit to the existing experimental data. It is found that with the inclusion of the factorization, FAT generally works better than the topological diagram approach, with less parameters and better {\chi }^2 per degree of freedom. Specifically, in Ref. [41] for the analysis of D\to PP decays, the authors assign a nonfactorizable term (magnitude and phase) to each of the color-suppressed, W-exchange and W-annihilation amplitudes, and the Glauber phase is additionally associated to a pion (which is important to resolve the {\pi }^{+}{\pi }^{-} and K^{+}{K}^{-} branching fraction puzzles). Then 12 parameters are used to fit 28 D\to PP branching fractions and good results are achieved. For the D\to PV decays [42], two more parameters are involved compared to the D\to PP ones, with 33 experimental numbers of branching fractions in total. Once the parameters are determined, the authors predict the CP asymmetries for D decays. The results in Ref. [43] are very impressive. The 4 universal parameters are associated to the color-suppressed amplitude and the W-exchange amplitude, i.e., parameterizing their sizes and phases, which are used to describe the 31 decay branching fractions induced by the b\to c transition in the B\to D^{\ast } M decays with M denoting a light pseudoscalar or a vector meson. If available in experiment, the predicted values are consistent with them. Then other 120 decay branching fractions are predicted. The similar analysis of the charmless two-body non-leptonic B decays B\to PP,P V is done in Ref. [44]. In brief, the nonfactorizable contributions require a fine analysis which is essential for a quantitative prediction of the branching fractions, and it is worth working in this direction in the future.

Here we stress again that we use the most recent form factor values. This is one of our important motivations and improvements. It is known that the form factors, which encode the underlying dynamics, play a significant role in calculations of the nonleptonic decays, as also noted in Ref. [55]. In the earlier works [40,47,55], the authors use the form factor values from the sum rule calculations, which are more appropriate for the small values of the momentum-transfer to leptons, or use the older predictions from the BSW model [40]. In our case, the RQM includes all sources of relativistic effects, and the form factors are obtained in the whole kinematically allowed region without any extrapolations. Transitions like B\to D and B\to D^{\ast } are also considered without using the heavy quark limit.

Moreover, we have calculated as many channels as possible. Previous papers studied only some of them (although the more advanced tools in a formal perspective were used). For example, in Ref. [56] only two decay modes are discussed. In Ref. [55] the B\to PP, PV channels are calculated, but not the B\to VV case. We have performed a complete calculation of the B_{(s)}, D_{(s)} decay to PP, PV and VV. In this way, we could show how the form factors influence the results from a holistic point of view based on such framework. So our calculations should be, at least, a useful complement and an important update for the previous ones.

5.2 A note on the conventions for the definitions of form factors and decay constants

In some references [57,58], the following quark compositions for the octet mesons are used

(41)\begin{array}{rl}& K^{+}=u\overline{s},K^0=d\overline{s},K^{-}=-s\overline{u},\\ & {\pi }^{+}({\rho }^{+})=u\overline{d},{\pi }^0({\rho }^0)=\frac{d\overline{d}-u\overline{u}}{\sqrt{2}},{\pi }^{-}({\rho }^{-})=-d\overline{u},\end{array}

which are different from Eq. (30) for the K^{-},{\pi }^0 ,{\rho }^0,{\pi }^{-},{\rho }^{-} cases. Then the definitions of the decay constants as well as the corresponding transition form factors will change by an overall sign. Any physical result is not affected.

We have also checked different conventions on definitions of the form factors and decay constants, which differ by factors of (-1) and/or i. Note that this detail may influence the calculation if using an inappropriate/incosistent convention. In the factorization scheme we are treating the product of ⟨M_1|J^{\mu }|B⟩⟨ M_2|J_{\mu }|0⟩, and surely an overall sign does not matter. However, for e.g., the channel B\to \pi \rho has the subprocess ⟨\pi |J^{\mu }|B⟩⟨ \rho |J_{\mu } | 0⟩ and ⟨\rho |J^{\mu }|B⟩⟨ \pi |J_{\mu } | 0⟩ and thus their interference occurs. Note also the convention difference {ϵ}^{0123}\equiv -1 and {ϵ}^{0123}\equiv +1, where the former is used in Refs. [56,59] and the latter is used in Refs. [4,60]. We have checked that the final results in Refs. [35,40,54-56] agree with each other just up to an overall factor of (-\mathrm{i}),\mathrm{i} or -1, which have no influence for branching fraction of a two-body decay. As a result, the vector decay constant should be defined by ⟨V(ϵ ,q)|V_{\mu }|0⟩=-\mathrm{i}{f}_V{m}_V {ϵ}_{\mu }^{\ast } in Ref. [4]. Then the amplitude for {\overline{B}}^0\to {\pi }^0{\rho }^0 follows:

(42)\begin{array}{rl}A({\overline{B}}^0\to {\pi }^0{\rho }^0)=& \frac{G_F}{\sqrt{2}}[-\mathrm{i}{f}_{\rho }{m}_{\rho }{ϵ}^{\ast }\cdot p_{\pi }{F}_1(m_{\rho }^2)\\ & -\mathrm{i}{f}_{\pi }{m}_{\rho }{ϵ}^{\ast }\cdot p_{\pi }{A}_0(m_{\pi }^2)].\end{array}

It is also important to mention that all authors use real and positive form-factor and decay-constant values.

6 Conclusions

Based on the form factors computed in the relativistic quark model, we calculate the branching fraction of 100 nonleptonic decay channels of charm and bottom mesons. We provide the detailed derivation for the decay amplitudes and branching fractions. The numerical results are shown in the Section 5.

For the nonleptonic decay process of the D mesons, we consider only the tree-level contributions and use the number of colors N=3 and N\to \mathrm{\infty } for demonstration. Taking the value of N different from 3 is a way to parametrize the nonfactorizable effects. And indeed we find that the limit N\to \mathrm{\infty } works much better than other numbers of N generally. Some typical examples are D^{+}\to {\pi }^0{\pi }^{+}, D^0\to K^{-} K^{+}, D^0\to \eta \eta, D^0\to K^{-} {\rho }^{+}, and D_s\to K^{+}\omega.

For the nonleptonic decay process of the bottom mesons, we consider both the tree-level and penguin contributions. The results for branching fractions are in agreement with the experimental data for most of decay channels. However, for some decays, e.g., B^{+}\to {\pi }^{+}{K}^0 ,{\rho }^{+}{K}^0 and B^0\to {\pi }^0{\eta }^{\prime },{\pi }^0 {\pi }^0,{\rho }^{-}{K}^{+} our results are too small, and as it has been demonstrated in Ref. [40], the final-state interaction effects play an indispensable role to get the quantitatively correct values.

7 Appendix A: Decay amplitudes of the charm mesons

Here we list the amplitudes for the charm meson decays. The definition of X is given in Section 2.

(A1)\begin{array}{rl}A(D^{+}\to {\pi }^0{\pi }^{+})=& \frac{G_F}{\sqrt{2}}{V}_{cd}^{\ast }{V}_{ud}[a_1{X}_{D^{+}\to {\pi }^0,{\pi }^{+}}\\ & +a_2{X}_{D^{+}\to {\pi }^{+},{\pi }^0}],\end{array}

(A2)A(D^{+}\to {\pi }^0{K}^{+})={\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cd}^{\ast }{V}_{us}{a}_1{X}_{D^{+}\to {\pi }^0,K^{+}},

(A3)A(D^{+}\to {\eta }^{{(}^{\prime })}{K}^{+})={\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cd}^{\ast }{V}_{us}{a}_1{X}_{D^{+}\to {\eta }^{{(}^{\prime })},K^{+}},

(A4)\begin{array}{rl}A(D^{+}\to {\eta }^{{(}^{\prime })}{\pi }^{+})=& \frac{G_F}{\sqrt{2}}[V_{cd}^{\ast }{V}_{ud}{a}_1{X}_{D^{+}\to {\eta }^{{(}^{\prime })},{\pi }^{+}}\\ & +V_{cd}^{\ast }{V}_{ud}{a}_2{X}_{D^{+}\to {\pi }^{+},{\eta }_u^{{(}^{\prime })}}\\ & +V_{cs}^{\ast }{V}_{us}{a}_2{X}_{D^{+}\to {\pi }^{+},{\eta }_s^{{(}^{\prime })}}],\end{array}

(A5)\begin{array}{rl}A(D^{+}\to {\pi }^{+}{\rho }^0)=& \frac{G_F}{\sqrt{2}}{V}_{cd}^{\ast }{V}_{ud}(a_1{X}_{D^{+}\to {\rho }^0,{\pi }^{+}}\\ & +a_2{X}_{D^{+}\to {\pi }^{+},{\rho }^0}),\end{array}

(A6)A(D^{+}\to {\pi }^{+}\varphi )={\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cs}^{\ast }{V}_{us}{a}_2{X}_{D^{+}\to {\pi }^{+},\varphi },

(A7)\begin{array}{rl}A(D^{+}\to {\pi }^{+}\omega )=& \frac{G_F}{\sqrt{2}}{V}_{cd}^{\ast }{V}_{ud}(a_1{X}_{D^{+}\to \omega ,{\pi }^{+}}\\ & +a_2{X}_{D^{+}\to {\pi }^{+},\omega }),\end{array}

(A8)A(D^{+}\to {\rho }^0{K}^{+})={\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cd}^{\ast }{V}_{us}{a}_1{X}_{D^{+}\to {\rho }^0,K^{+}},

(A9)A(D^{+}\to {\rho }^{+}\varphi )={\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cs}^{\ast }{V}_{us}{a}_2{X}_{D^{+}\to {\rho }^{+},\varphi },

(A10)A(D^0\to K^{-}{\pi }^{+})= {\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cs}^{\ast }{V}_{ud} a_1{X}_{D^0\to K^{-},{\pi }^{+}},

(A11)A(D^0\to {\pi }^{-}{\pi }^{+})= {\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cd}^{\ast }{V}_{ud} a_1{X}_{D^0\to {\pi }^{-},{\pi }^{+}},

(A12)A(D^0\to {\pi }^0{\pi }^0)=2{\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cd}^{\ast }{V}_{ud} a_2{X}_{D^0\to {\pi }^0,{\pi }^0},

(A13)A(D^0\to K^{-}{K}^{+})= {\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cs}^{\ast }{V}_{us} a_1{X}_{D^0\to K^{-},K^{+}},

(A14)\begin{array}{rl}A(D^0\to \eta \eta )=& 2\frac{G_F}{\sqrt{2}}(V_{cs}^{\ast }{V}_{us}{a}_2{X}_{D^0\to \eta ,{\eta }_s}\\ & +V_{cd}^{\ast }{V}_{ud}{a}_2{X}_{D^0\to \eta ,{\eta }_u}),\end{array}

(A15)A(D^0\to {\pi }^{-}{K}^{+})= {\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cd}^{\ast }{V}_{us} a_1{X}_{D^0\to {\pi }^{-},K^{+}},

(A16)\begin{array}{rl}A(D^0\to {\eta }^{{(}^{\prime })}{\pi }^0)=& \frac{G_F}{\sqrt{2}}(V_{cd}^{\ast }{V}_{ud}{a}_2{X}_{D^0\to {\eta }^{{(}^{\prime })},{\pi }^0}\\ & +V_{cd}^{\ast }{V}_{ud}{a}_2{X}_{D^0\to {\pi }^0,{\eta }_u^{{(}^{\prime })}}\\ & +V_{cs}^{\ast }{V}_{us}{a}_2{X}_{D^0\to {\pi }^0,{\eta }_s^{{(}^{\prime })}}),\end{array}

(A17)\begin{array}{rl}A(D^0\to \eta {\eta }^{\prime })=& \frac{G_F}{\sqrt{2}}(V_{cd}^{\ast }{V}_{ud}{a}_2{X}_{D^0\to \eta ,{\eta }_u^{\prime }}\\ & +V_{cs}^{\ast }{V}_{us}{a}_2{X}_{D^0\to \eta ,{\eta }_s^{\prime }}\\ & +V_{cd}^{\ast }{V}_{ud}{a}_2{X}_{D^0\to {\eta }^{\prime },{\eta }_u}\\ & +V_{cs}^{\ast }{V}_{us}{a}_2{X}_{D^0\to {\eta }^{\prime },{\eta }_s}),\end{array}

(A18)\begin{array}{rl}A(D^0\to {\pi }^0\omega )=& \frac{G_F}{\sqrt{2}}(V_{cd}^{\ast }{V}_{ud}{a}_2{X}_{D^0\to {\pi }^0,\omega }\\ & +V_{cd}^{\ast }{V}_{ud}{a}_2{X}_{D^0\to \omega ,{\pi }^0}),\end{array}

(A19)\begin{array}{rl}A(D^0\to \eta \omega )=& \frac{G_F}{\sqrt{2}}(V_{cd}^{\ast }{V}_{ud}{a}_2{X}_{D^0\to \eta ,\omega }\\ & +V_{cd}^{\ast }{V}_{ud}{a}_2{X}_{D^0\to \omega ,{\eta }_u}\\ & +V_{cs}^{\ast }{V}_{us}{a}_2{X}_{D^0\to \omega ,{\eta }_s}),\end{array}

(A20)\begin{array}{rl}A(D^0\to {\rho }^0{\pi }^0)=& \frac{G_F}{\sqrt{2}}(V_{cd}^{\ast }{V}_{ud}{a}_2{X}_{D^0\to {\pi }^0,{\rho }^0}\\ & +V_{cd}^{\ast }{V}_{ud}{a}_2{X}_{D^0\to {\rho }^0,{\pi }^0}),\end{array}

(A21)A(D^0\to {\pi }^{-}{\rho }^{+})= {\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cd}^{\ast }{V}_{ud} a_1{X}_{D^0\to {\pi }^{-},{\rho }^{+}},

(A22)A(D^0\to {\pi }^0\varphi )={\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cs}^{\ast }{V}_{us} a_2{X}_{D^0\to {\pi }^0,\varphi },

(A23)A(D^0\to {\rho }^{-}{\pi }^{+})= {\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cd}^{\ast }{V}_{ud} a_1{X}_{D^0\to {\rho }^{-},{\pi }^{+}},

(A24)A(D^0\to \eta \varphi )= {\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cs}^{\ast }{V}_{us} a_2{X}_{D^0\to \eta ,\varphi },

(A25)A(D^0\to K^{-}{\rho }^{+})= {\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cs}^{\ast }{V}_{ud} a_1{X}_{D^0\to K^{-},{\rho }^{+}},

(A26)A(D^0\to {\eta }^{{(}^{\prime })}{\overline{K}}^{\ast 0})={\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cs}^{\ast }{V}_{ud} a_2{X}_{D^0\to {\eta }^{{(}^{\prime })},{\overline{K}}^{\ast 0}},

(A27)A(D^0\to {\rho }^0{\rho }^0)= 2{\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cd}^{\ast }{V}_{ud} a_2{X}_{D^0\to {\rho }^0,{\rho }^0},

(A28)A(D^0\to \omega \varphi )= {\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cs}^{\ast }{V}_{us} a_2{X}_{D^0\to \omega ,\varphi },

(A29)A(D_s\to K^{+}{\overline{K}}^0)={\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cs}^{\ast }{V}_{ud} a_2{X}_{D_s\to K^{+},{\overline{K}}^0},

(A30)A(D_s\to {\eta }^{{(}^{\prime })}{\pi }^{+})= {\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cs}^{\ast }{V}_{ud} a_1{X}_{D_s\to {\eta }^{{(}^{\prime })},{\pi }^{+}},

(A31)A(D_s\to K^{+}{\pi }^0)= {\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cd}^{\ast }{V}_{ud} a_2{X}_{D_s\to K^{+},{\pi }^0},

(A32)\begin{array}{rl}A(D_s\to {\eta }^{{(}^{\prime })}{K}^{+})=& \frac{G_F}{\sqrt{2}}(V_{cd}^{\ast }{V}_{ud}{a}_2{X}_{D_s\to K^{+},{\eta }_u^{{(}^{\prime })}}\\ & +V_{cs}^{\ast }{V}_{us}{a}_2{X}_{D_s\to K^{+},{\eta }_s^{{(}^{\prime })}}\\ & +V_{cs}^{\ast }{V}_{us}{a}_1{X}_{D_s\to {\eta }^{{(}^{\prime })},K^{+}}),\end{array}

(A33)A(D_s\to {\eta }^{{(}^{\prime })}{\rho }^{+})= {\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cs}^{\ast }{V}_{ud} a_1{X}_{D_s\to {\eta }^{{(}^{\prime })},{\rho }^{+}},

(A34)A(D_s\to K^{+}\omega )={\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cd}^{\ast }{V}_{ud} a_2{X}_{D_s\to K^{+},\omega },

(A35)A(D_s\to K^{+}{\rho }^0)={\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cd}^{\ast }{V}_{ud}{a}_2{X}_{D_s\to K^{+},{\rho }^0},

(A36)A(D_s\to \varphi {\pi }^{+})={\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cs}^{\ast }{V}_{ud}{a}_1{X}_{D_s\to \varphi ,{\pi }^{+}},

(A37)A(D_s\to \varphi {\rho }^{+})={\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cs}^{\ast }{V}_{ud} a_1{X}_{D_s\to \varphi ,{\rho }^{+}}.

8 Appendix B: Decay amplitudes of the bottom mesons

As in Appendix A, here the amplitude for the bottom meson decays are provided.

(B1)\begin{array}{rl}A(B^{+}\to {\pi }^{+}{\eta }^{{(}^{\prime })})=& \frac{G_F}{\sqrt{2}}\{[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}(2a_3-2a_5-\frac{1}{2}{a}_7+\frac{1}{2}{a}_9+a_4\\ & -\frac{1}{2}{a}_{10}+\frac{m_{{\eta }^{{(}^{\prime })}}^2}{m_s(m_b-m_d)}(a_6-\frac{1}{2}{a}_8)(\frac{f_{{\eta }^{{(}^{\prime })}}^s}{f_{{\eta }^{{(}^{\prime })}}^u}-1)r_{{\eta }^{{(}^{\prime })}})]X_{B^{+}\to {\pi }^{+},{\eta }_u^{{(}^{\prime })}}\\ & -V_{tb}^{\ast }{V}_{td}\left(a_3-a_5+\frac{1}{2}{a}_7-\frac{1}{2}{a}_9\right)X_{B^{+}\to {\pi }^{+},{\eta }_s^{{(}^{\prime })}}+[V_{ub}^{\ast }{V}_{ud}{a}_1\\ & -V_{tb}^{\ast }{V}_{td}\left(a_4+a_{10}+\frac{2m_{{\pi }^{+}}^2}{(m_u+m_d)(m_b-m_u)}(a_6+a_8)\right)]X_{B^{+}\to {\eta }^{{(}^{\prime })},{\pi }^{+}}\},\end{array}

(B2)\begin{array}{rl}A(B^{+}\to {\pi }^{+}\omega )=& \frac{G_F}{\sqrt{2}}\{\left[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}\left(2a_3+2a_5+\frac{1}{2}{a}_7+\frac{1}{2}{a}_9+a_4-\frac{1}{2}{a}_{10}\right)\right]X_{B^{+}\to {\pi }^{+},\omega }\\ & +[V_{ub}^{\ast }{V}_{ud}{a}_1-V_{tb}^{\ast }{V}_{td}(a_4+a_{10}\\ & -\frac{2m_{{\pi }^{+}}^2}{(m_u+m_d)(m_b+m_u)}(a_6+a_8))]X_{B^{+}\to \omega ,{\pi }^{+}}\},\end{array}

(B3)\begin{array}{rl}A(B^{+}\to {\rho }^{+}{\eta }^{{(}^{\prime })})=& \frac{G_F}{\sqrt{2}}\{[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}(2a_3-2a_5-\frac{1}{2}{a}_7+\frac{1}{2}{a}_9+a_4\\ & -\frac{1}{2}{a}_{10}-\frac{m_{{\eta }^{{(}^{\prime })}}^2}{m_s(m_b+m_d)}(a_6-\frac{1}{2}{a}_8)(\frac{f_{{\eta }^{{(}^{\prime })}}^s}{f_{{\eta }^{{(}^{\prime })}}^u}-1)r_{{\eta }^{{(}^{\prime })}})]X_{B^{+}\to {\rho }^{+},{\eta }_u^{{(}^{\prime })}}\\ & -V_{tb}^{\ast }{V}_{td}\left(a_3-a_5+\frac{1}{2}{a}_7-\frac{1}{2}{a}_9\right)X_{B^{+}\to {\rho }^{+},{\eta }_s^{{(}^{\prime })}}\\ & +[V_{ub}^{\ast }{V}_{ud}{a}_1-V_{tb}^{\ast }{V}_{td}\left(a_4+a_{10}\right)]X_{B^{+}\to {\eta }^{{(}^{\prime })},{\rho }^{+}}\},\end{array}

(B4)A(B^{+}\to {\pi }^{+}{K}^0)=-{\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{tb}^{\ast }{V}_{ts} \left(a_4-{\displaystyle \frac{1}{2}}{a}_{10}+{\displaystyle \frac{2m_{K^0}^2}{(m_s+m_d)(m_b-m_d)}}(a_6-{\displaystyle \frac{1}{2}}{a}_8)\right)X_{B^{+}\to {\pi }^{+},K^0},

(B5)A(B^{+}\to {\rho }^{+}{K}^0)=-{\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{tb}^{\ast }{V}_{ts} \left(a_4-{\displaystyle \frac{1}{2}}{a}_{10}-{\displaystyle \frac{2m_{K^0}^2}{(m_s+m_d)(m_b+m_d)}}(a_6-{\displaystyle \frac{1}{2}}{a}_8)\right)X_{B^{+}\to {\rho }^{+},K^0},

(B6)\begin{array}{rl}A(B^{+}\to {\pi }^{+}{\pi }^0)=& \frac{G_F}{\sqrt{2}}\{[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}(\frac{3}{2}{a}_9-\frac{3}{2}{a}_7\\ & -a_4+\frac{1}{2}{a}_{10}-\frac{m_{{\pi }^0}^2}{m_d(m_b-m_d)}(a_6-\frac{1}{2}{a}_8))]X_{B^{+}\to {\pi }^{+},{\pi }^0}\\ & +[V_{ub}^{\ast }{V}_{ud}{a}_1-V_{tb}^{\ast }{V}_{td}(a_4+a_{10}\\ & +\frac{2m_{{\pi }^{+}}^2}{(m_u+m_d)(m_b-m_u)}(a_6+a_8))]X_{B^{+}\to {\pi }^0,{\pi }^{+}}\},\end{array}

(B7)\begin{array}{rl}A(B^{+}\to {\pi }^{+}{\rho }^0)=& \frac{G_F}{\sqrt{2}}\{\left[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}\left(-a_4+\frac{1}{2}{a}_{10}+\frac{3}{2}{a}_9+\frac{3}{2}{a}_7\right)\right]X_{B^{+}\to {\pi }^{+},{\rho }^0}\\ & +[V_{ub}^{\ast }{V}_{ud}{a}_1-V_{tb}^{\ast }{V}_{td}(a_4+a_{10}\\ & -\frac{2m_{{\pi }^{+}}^2}{(m_u+m_d)(m_b+m_u)}(a_6+a_8))]X_{B^{+}\to {\rho }^0,{\pi }^{+}}\},\end{array}

(B8)\begin{array}{rl}A(B^{+}\to {\rho }^{+}{\pi }^0)=& \frac{G_F}{\sqrt{2}}\{[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}(\frac{3}{2}{a}_9-\frac{3}{2}{a}_7-a_4+\frac{1}{2}{a}_{10}\\ & +\frac{m_{{\pi }^0}^2}{m_d(m_b+m_d)}(a_6-\frac{1}{2}{a}_8))]X_{B^{+}\to {\rho }^{+},{\pi }^0}\\ & +\left[V_{ub}^{\ast }{V}_{ud}{a}_1-V_{tb}^{\ast }{V}_{td}\left(a_4+a_{10}\right)\right]X_{B^{+}\to {\pi }^0,{\rho }^{+}}\},\end{array}

(B9)A(B^{+}\to {\pi }^{+}\varphi )=-{\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{tb}^{\ast } V_{td}(a_3+a_5-{\displaystyle \frac{1}{2}}{a}_7- {\displaystyle \frac{1}{2}}{a}_9)X_{B^{+}\to {\pi }^{+},\varphi } ,

(B10)\begin{array}{rl}A(B^{+}\to {\rho }^{+}{\rho }^0)=& \frac{G_F}{\sqrt{2}}\{\left[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}\left(\frac{3}{2}{a}_9+\frac{3}{2}{a}_7-a_4+\frac{1}{2}{a}_{10}\right)\right]X_{B^{+}\to {\rho }^{+},{\rho }^0}\\ & +[V_{ub}^{\ast }{V}_{ud}{a}_1-V_{tb}^{\ast }{V}_{td}\left(a_4+a_{10}\right)]X_{B^{+}\to {\rho }^0,{\rho }^{+}}\},\end{array}

(B11)\begin{array}{rl}A(B^{+}\to {\rho }^{+}\omega )=& \frac{G_F}{\sqrt{2}}\{\left[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}\left(2a_3+2a_5+\frac{1}{2}{a}_7+\frac{1}{2}{a}_9+a_4-\frac{1}{2}{a}_{10}\right)\right]X_{B^{+}\to {\rho }^{+},\omega }\\ & +[V_{ub}^{\ast }{V}_{ud}{a}_1-V_{tb}^{\ast }{V}_{td}\left(a_4+a_{10}\right)]X_{B^{+}\to \omega ,{\rho }^{+}}\},\end{array}

(B12)A(B^0\to D^{-}{\pi }^{+})= {\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cb}^{\ast }{V}_{ud} a_1{X}_{B^0\to D^{-},{\pi }^{+}},

(B13)A(B^0\to D^{-}{K}^{+})= {\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cb}^{\ast }{V}_{us} a_1{X}_{B^0\to D^{-},K^{+}},

(B14)\begin{array}{rl}A(B^0\to {\pi }^{-}{K}^{+})=& \frac{G_F}{\sqrt{2}}[V_{ub}^{\ast }{V}_{us}{a}_1-V_{tb}^{\ast }{V}_{ts}(a_4+a_{10}\\ & +\frac{2m_{K^{+}}^2}{(m_u+m_s)(m_b-m_u)}(a_6+a_8))]X_{B^0\to {\pi }^{-},K^{+}},\end{array}

(B15)\begin{array}{rl}A(B^0\to {\pi }^{-}{\pi }^{+})=& \frac{G_F}{\sqrt{2}}[V_{ub}^{\ast }{V}_{ud}{a}_1-V_{tb}^{\ast }{V}_{td}(a_4+a_{10}\\ & +\frac{2m_{{\pi }^{+}}^2}{(m_u+m_d)(m_b-m_u)}(a_6+a_8))]X_{B^0\to {\pi }^{-},{\pi }^{+}},\end{array}

(B16)\begin{array}{rl}A(B^0\to {\pi }^0{\pi }^0)=& 2\frac{G_F}{\sqrt{2}}[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}(\frac{3}{2}{a}_9-\frac{3}{2}{a}_7-a_4+\frac{1}{2}{a}_{10}\\ & -\frac{m_{{\pi }^0}^2}{m_d(m_b-m_d)}(a_6-\frac{1}{2}{a}_8))]X_{B^0\to {\pi }^0,{\pi }^0},\end{array}

(B17)\begin{array}{rl}A(B^0\to {\pi }^0{\eta }^{{(}^{\prime })})=& \frac{G_F}{\sqrt{2}}\{[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}(2a_3-2a_5-\frac{1}{2}{a}_7+\frac{1}{2}{a}_9+a_4-\frac{1}{2}{a}_{10}\\ & +\frac{m_{{\eta }^{{(}^{\prime })}}^2}{m_s(m_b-m_d)}(a_6-\frac{1}{2}{a}_8)(\frac{f_{{\eta }^{{(}^{\prime })}}^s}{f_{{\eta }^{{(}^{\prime })}}^u}-1)r_{{\eta }^{{(}^{\prime })}})]X_{B^0\to {\pi }^0,{\eta }_u^{{(}^{\prime })}}\\ & -V_{tb}^{\ast }{V}_{td}\left(a_3-a_5+\frac{1}{2}{a}_7-\frac{1}{2}{a}_9\right)X_{B^0\to {\pi }^0,{\eta }_s^{{(}^{\prime })}}\\ & +[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}(\frac{3}{2}{a}_9-\frac{3}{2}{a}_7-a_4+\frac{1}{2}{a}_{10}\\ & -\frac{m_{{\pi }^0}^2}{m_d(m_b-m_d)}(a_6-\frac{1}{2}{a}_8))]X_{B^0\to {\eta }^{{(}^{\prime })},{\pi }^0}\},\end{array}

(B18)\begin{array}{rl}A(B^0\to \eta \eta )=& 2\frac{G_F}{\sqrt{2}}\{[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}(2a_3-2a_5-\frac{1}{2}{a}_7+\frac{1}{2}{a}_9+a_4-\frac{1}{2}{a}_{10}\\ & +\frac{2m_{\eta }^2}{(m_s+m_s)(m_b-m_d)}(a_6-\frac{1}{2}{a}_8)(\frac{f_{\eta }^s}{f_{\eta }^u}-1)r_{\eta })]X_{B^0\to \eta ,{\eta }_u}\\ & -V_{tb}^{\ast }{V}_{td}\left(a_3-a_5+\frac{1}{2}{a}_7-\frac{1}{2}{a}_9\right)X_{B^0\to \eta ,{\eta }_s}\},\end{array}

(B19)\begin{array}{rl}A(B^0\to {\eta }^{\prime }{\eta }^{\prime })=& 2\frac{G_F}{\sqrt{2}}\{[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}(2a_3-2a_5-\frac{1}{2}{a}_7+\frac{1}{2}{a}_9+a_4-\frac{1}{2}{a}_{10}\\ & +\frac{m_{{\eta }^{\prime }}^2}{m_s(m_b-m_d)}(a_6-\frac{1}{2}{a}_8)(\frac{f_{{\eta }^{\prime }}^s}{f_{{\eta }^{\prime }}^u}-1)r_{{\eta }^{\prime }})]X_{B^0\to {\eta }^{\prime },{\eta }_u^{\prime }}\\ & -V_{tb}^{\ast }{V}_{td}\left(a_3-a_5+\frac{1}{2}{a}_7-\frac{1}{2}{a}_9\right)X_{B^0\to {\eta }^{\prime },{\eta }_s^{\prime }}\},\end{array}

(B20)A(B^0\to D^{\ast -}{\pi }^{+})= {\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cb}^{\ast }{V}_{ud} a_1{X}_{B^0\to D^{\ast -},{\pi }^{+}},

(B21)A(B^0\to D^{-}{\rho }^{+})= {\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cb}^{\ast }{V}_{ud} a_1{X}_{B^0\to D^{-},{\rho }^{+}},

(B22)\begin{array}{rl}A(B^0\to {\rho }^0{\pi }^0)=& \frac{G_F}{\sqrt{2}}\{\left[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}\left(\frac{3}{2}{a}_7+\frac{3}{2}{a}_9-a_4+\frac{1}{2}{a}_{10}\right)\right]X_{B^0\to {\pi }^0,{\rho }^0}\\ & +[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}(\frac{3}{2}{a}_9-\frac{3}{2}{a}_7-a_4+\frac{1}{2}{a}_{10}\\ & +\frac{m_{{\pi }^0}^2}{m_d(m_b+m_d)}(a_6-\frac{1}{2}{a}_8))]X_{B^0\to {\rho }^0,{\pi }^0}\},\end{array}

(B23)\begin{array}{rl}A(B^0\to \omega {\pi }^0)=& \frac{G_F}{\sqrt{2}}\{\left[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}\left(2a_3+2a_5+\frac{1}{2}{a}_7+\frac{1}{2}{a}_9+a_4-\frac{1}{2}{a}_{10}\right)\right]X_{B^0\to {\pi }^0,\omega }\\ & +[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}(\frac{3}{2}{a}_9-\frac{3}{2}{a}_7-a_4+\frac{1}{2}{a}_{10}\\ & +\frac{m_{{\pi }^0}^2}{m_d(m_b+m_d)}(a_6-\frac{1}{2}{a}_8))]X_{B^0\to \omega ,{\pi }^0}\},\end{array}

(B24)\begin{array}{rl}A(B^0\to {\rho }^{-}{K}^{+})=& \frac{G_F}{\sqrt{2}}[V_{ub}^{\ast }{V}_{us}{a}_1-V_{tb}^{\ast }{V}_{ts}(a_4+a_{10}\\ & -\frac{2m_{K^{+}}^2}{(m_s+m_u)(m_b+m_u)}(a_6+a_8))]X_{B^0\to {\rho }^{-},K^{+}},\end{array}

(B25)A(B^0\to {\pi }^{-}{K}^{\ast +})= {\displaystyle \frac{G_F}{\sqrt{2}}}[V_{ub}^{\ast }{V}_{us} a_1- V_{tb}^{\ast }{V}_{ts}\left(a_4+a_{10}\right)]X_{B^0\to {\pi }^{-},K^{\ast +}},

(B26)\begin{array}{rl}A(B^0\to {\rho }^{-}{\pi }^{+})=& \frac{G_F}{\sqrt{2}}[V_{ub}^{\ast }{V}_{ud}{a}_1-V_{tb}^{\ast }{V}_{td}(a_4+a_{10}\\ & -\frac{2m_{{\pi }^{+}}^2}{(m_u+m_d)(m_b+m_u)}(a_6+a_8))]X_{B^0\to {\rho }^{-},{\pi }^{+}},\end{array}

(B27)A(B^0\to D^{-}{K}^{\ast +})= {\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cb}^{\ast }{V}_{us} a_1{X}_{B^0\to D^{-},K^{\ast +}},

(B28)A(B^0\to D^{\ast -}{K}^{+})= {\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cb}^{\ast }{V}_{us} a_1{X}_{B^0\to D^{\ast -},K^{+}},

(B29)\begin{array}{rl}A(B^0\to {\eta }^{\prime }\eta )=& \frac{G_F}{\sqrt{2}}\{[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}(2a_3-2a_5-\frac{1}{2}{a}_7+\frac{1}{2}{a}_9+a_4-\frac{1}{2}{a}_{10}\\ & +\frac{m_{\eta }^2}{m_s(m_b-m_d)}(a_6-\frac{1}{2}{a}_8)(\frac{f_{\eta }^s}{f_{\eta }^u}-1)r_{\eta })]X_{B^0\to {\eta }^{{}^{\prime }},{\eta }_u}\\ & -V_{tb}^{\ast }{V}_{td}\left(a_3-a_5+\frac{1}{2}{a}_7-\frac{1}{2}{a}_9\right)X_{B^0\to {\eta }^{\prime },{\eta }_s}\\ & +[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}(2a_3-2a_5-\frac{1}{2}{a}_7+\frac{1}{2}{a}_9+a_4-\frac{1}{2}{a}_{10}\\ & +\frac{m_{{\eta }^{{}^{\prime }}}^2}{m_s(m_b-m_d)}(a_6-\frac{1}{2}{a}_8)(\frac{f_{{\eta }^{{}^{\prime }}}^s}{f_{{\eta }^{{}^{\prime }}}^u}-1)r_{{\eta }^{{}^{\prime }}})]X_{B^0\to \eta ,{\eta }_u^{{}^{\prime }}}\\ & -V_{tb}^{\ast }{V}_{td}\left(a_3-a_5+\frac{1}{2}{a}_7-\frac{1}{2}{a}_9\right)X_{B^0\to \eta ,{\eta }_s^{\prime }}\},\end{array}

(B30)\begin{array}{rl}A(B^0\to {\rho }^0{\eta }^{{(}^{\prime })})=& \frac{G_F}{\sqrt{2}}\{[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}(2a_3-2a_5-\frac{1}{2}{a}_7+\frac{1}{2}{a}_9+a_4-\frac{1}{2}{a}_{10}\\ & -\frac{m_{{\eta }^{{(}^{\prime })}}^2}{m_s(m_b+m_d)}(a_6-\frac{1}{2}{a}_8)(\frac{f_{{\eta }^{{(}^{\prime })}}^s}{f_{{\eta }^{{(}^{\prime })}}^u}-1)r_{{\eta }^{{(}^{\prime })}})]X_{B^0\to {\rho }^0,{\eta }_u^{{(}^{\prime })}}\\ & -V_{tb}^{\ast }{V}_{td}\left(a_3-a_5+\frac{1}{2}{a}_7-\frac{1}{2}{a}_9\right)X_{B^0\to {\rho }^0,{\eta }_s^{{(}^{\prime })}}\\ & +\left[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}\left(\frac{3}{2}{a}_7+\frac{3}{2}{a}_9-a_4+\frac{1}{2}{a}_{10}\right)\right]X_{B^0\to {\eta }^{{(}^{\prime })},{\rho }^0}\},\end{array}

(B31)\begin{array}{rl}A(B^0\to \omega {\eta }^{{(}^{\prime })})=& \frac{G_F}{\sqrt{2}}\{[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}(2a_3-2a_5-\frac{1}{2}{a}_7+\frac{1}{2}{a}_9+a_4-\frac{1}{2}{a}_{10}\\ & -\frac{m_{{\eta }^{{(}^{\prime })}}^2}{m_d(m_b+m_d)}(a_6-\frac{1}{2}{a}_8)(\frac{f_{{\eta }^{{(}^{\prime })}}^s}{f_{{\eta }^{{(}^{\prime })}}^u}-1)r_{{\eta }^{{(}^{\prime })}})]X_{B^0\to \omega ,{\eta }_u^{{(}^{\prime })}}\\ & -V_{tb}^{\ast }{V}_{td}\left(a_3-a_5+\frac{1}{2}{a}_7-\frac{1}{2}{a}_9\right)X_{B^0\to \omega ,{\eta }_s^{{(}^{\prime })}}+[V_{ub}^{\ast }{V}_{ud}{a}_2\\ & -V_{tb}^{\ast }{V}_{td}\left(2a_3+2a_5+\frac{1}{2}{a}_7+\frac{1}{2}{a}_9+a_4-\frac{1}{2}{a}_{10}\right)]X_{B^0\to {\eta }^{{(}^{\prime })},\omega }\},\end{array}

(B32)A(B^0\to {\pi }^{-}{\rho }^{+})= {\displaystyle \frac{G_F}{\sqrt{2}}}[V_{ub}^{\ast }{V}_{ud} a_1- V_{tb}^{\ast }{V}_{td}\left(a_4+a_{10}\right)]X_{B^0\to {\pi }^{-},{\rho }^{+}},

(B33)A(B^0\to {\rho }^{-}{\rho }^{+})= {\displaystyle \frac{G_F}{\sqrt{2}}}[V_{ub}^{\ast }{V}_{ud} a_1- V_{tb}^{\ast }{V}_{td}\left(a_4+a_{10}\right)]X_{B^0\to {\rho }^{-},{\rho }^{+}},

(B34)A(B^0\to {\rho }^0{\rho }^0)= 2{\displaystyle \frac{G_F}{\sqrt{2}}}\left[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}\left({\displaystyle \frac{3}{2}}{a}_9+{\displaystyle \frac{3}{2}}{a}_7-a_4+{\displaystyle \frac{1}{2}}{a}_{10}\right)\right]X_{B^0\to {\rho }^0,{\rho }^0},

(B35)A(B^0\to \omega \omega )=2{\displaystyle \frac{G_F}{\sqrt{2}}}\left[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}\left(2a_3+2a_5+{\displaystyle \frac{1}{2}}{a}_7+{\displaystyle \frac{1}{2}}{a}_9+a_4-{\displaystyle \frac{1}{2}}{a}_{10}\right)\right]X_{B^0\to \omega ,\omega },

(B36)\begin{array}{rl}A(B^0\to \omega {\rho }^0)=& \frac{G_F}{\sqrt{2}}\{\left[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}\left(2a_3+2a_5+\frac{1}{2}{a}_7+\frac{1}{2}{a}_9+a_4-\frac{1}{2}{a}_{10}\right)\right]X_{B^0\to {\rho }^0,\omega }\\ & +\left[V_{ub}^{\ast }{V}_{ud}{a}_2-V_{tb}^{\ast }{V}_{td}\left(\frac{3}{2}{a}_7+\frac{3}{2}{a}_9-a_4+\frac{1}{2}{a}_{10}\right)\right]X_{B^0\to \omega ,{\rho }^0}\},\end{array}

(B37)A(B^0\to D^{\ast -}{\rho }^{+})= {\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cb}^{\ast }{V}_{ud} a_1{X}_{B^0\to D^{\ast -},{\rho }^{+}},

(B38)A(B^0\to D^{\ast -}{K}^{\ast +})={\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cb}^{\ast } V_{us}{a}_1{X}_{B^0\to D^{\ast -},K^{\ast +}},

(B39)A(B_s\to D_s^{-}{\pi }^{+}) ={\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cb}^{\ast }{V}_{ud} a_1{X}_{B_s\to D_s^{-},{\pi }^{+}},

(B40)A(B_s\to D_s^{-}{K}^{+}) ={\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cb}^{\ast }{V}_{us} a_1{X}_{B_s\to D_s^{-},K^{+}},

(B41)\begin{array}{rl}A(B_s\to {\pi }^{+}{K}^{-})=& \frac{G_F}{\sqrt{2}}[V_{ub}^{\ast }{V}_{ud}{a}_1-V_{tb}^{\ast }{V}_{td}(a_4+a_{10}\\ & +\frac{2m_{{\pi }^{+}}^2}{(m_u+m_d)(m_b-m_u)}(a_6+a_8))]X_{B_s\to K^{-},{\pi }^{+}},\end{array}

(B42)\begin{array}{rl}A(B_s\to K^{+}{K}^{-})=& \frac{G_F}{\sqrt{2}}[V_{ub}^{\ast }{V}_{us}{a}_1-V_{tb}^{\ast }{V}_{ts}(a_4+a_{10}\\ & +\frac{2m_{K^{+}}^2}{(m_u+m_s)(m_b-m_u)}(a_6+a_8))]X_{B_s\to K^{-},K^{+}},\end{array}

(B43)A(B_s\to D_s^{-}{\rho }^{+}) ={\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cb}^{\ast }{V}_{ud} a_1{X}_{B_s\to D_s^{-},{\rho }^{+}},

(B44)A(B_s\to D_s^{\ast -}{\pi }^{+})={\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cb}^{\ast }{V}_{ud} a_1{X}_{B_s\to D_s^{\ast -},{\pi }^{+}},

(B45)A(B_s\to D_s^{\ast -}{K}^{+})={\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cb}^{\ast }{V}_{us} a_1{X}_{B_s\to D_s^{\ast -},K^{+}},

(B46)A(B_s\to D_s^{-}{K}^{\ast +})={\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cb}^{\ast }{V}_{us} a_1{X}_{B_s\to D_s^{-},K^{\ast +}},

(B47)\begin{array}{rl}A(B_s\to K^{\ast -}{\pi }^{+})=& \frac{G_F}{\sqrt{2}}[V_{ub}^{\ast }{V}_{ud}{a}_1-V_{tb}^{\ast }{V}_{td}(a_4+a_{10}\\ & -\frac{2m_{{\pi }^{+}}^2}{(m_u+m_d)(m_b+m_u)}(a_6+a_8))]X_{B_s\to K^{\ast -},{\pi }^{+}},\end{array}

(B48)\begin{array}{rl}A(B_s\to K^{\ast -}{K}^{+})=& \frac{G_F}{\sqrt{2}}[V_{ub}^{\ast }{V}_{us}{a}_1-V_{tb}^{\ast }{V}_{ts}(a_4+a_{10}\\ & -\frac{2m_{K^{+}}^2}{(m_u+m_s)(m_b+m_u)}(a_6+a_8))]X_{B_s\to K^{\ast -},K^{+}},\end{array}

(B49)A(B_s\to K^{-}{K}^{\ast +})={\displaystyle \frac{G_F}{\sqrt{2}}}[V_{ub}^{\ast }{V}_{us}{a}_1-V_{tb}^{\ast }{V}_{ts}(a_4+a_{10})]X_{B_s\to K^{-},K^{\ast +}},

(B50)A(B_s\to D_s^{\ast -}{K}^{\ast +})= {\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cb}^{\ast }{V}_{us} a_1{X}_{B_s\to D_s^{\ast -},K^{\ast +}},

(B51)A(B_s\to D_s^{\ast -}{\rho }^{+}) ={\displaystyle \frac{G_F}{\sqrt{2}}}{V}_{cb}^{\ast }{V}_{ud} a_1{X}_{B_s\to D_s^{\ast -},{\rho }^{+}}.

References

[1] J. D. Bjorken. Topics in B-physics. Nucl. Phys. B Proc. Suppl., 1989, 11: 325

[2] R.N. FaustovV.O. GalkinX.W. Kang, Semileptonic decays of D and D_s mesons in the relativistic quark model, Phys. Rev. D 101(1), 013004 (2020), arXiv: 1911.08209 [hep-ph]

[3] R.N. FaustovV.O. GalkinX.W. Kang, Relativistic description of the semileptonic decays of bottom mesons, Phys. Rev. D 106(1), 013004 (2022), arXiv: 2206.10277 [hep-ph]

[4] D.EbertR. N. FaustovV.O. Galkin, Exclusive nonleptonic decays of B mesons, Phys. Rev. D 56, 312 (1997), arXiv: hep-ph/9701218

[5] R.N. FaustovV.O. Galkin, Weak decays of B_s mesons to D_s mesons in the relativistic quark model, Phys. Rev. D 87(3), 034033 (2013), arXiv: 1212.3167 [hep-ph]

[6] R.N. FaustovV.O. Galkin, Charmless weak B_s decays in the relativistic quark model, Phys. Rev. D 87(9), 094028 (2013), arXiv: 1304.3255 [hep-ph]

[7] J.ChaiS. ChengY.H. JuD.C. YanC.D. LüZ.J. Xiao, Charmless two-body B meson decays in perturbative QCD factorization approach, arXiv: 2207.04190 [hep-ph] (2022)

[8] A.AliG. KramerY.LiC.D. LuY.L. Shen W.WangY. M. Wang, Charmless non-leptonic B_s decays to PP, PV and VV final states in the pQCD approach, Phys. Rev. D 76, 074018 (2007), arXiv: hep-ph/0703162

[9] Y.YangL. LangX.ZhaoJ.HuangJ.Sun, Reinvestigating the B → PP decays by including the contributions from ϕ_B2, Phys. Rev. D 103(5), 056006 (2021), arXiv: 2012.10581 [hep-ph]

[10] Z.J. XiaoW. F. WangY.Y. Fan, Revisiting the pure annihilation decays B_s → π⁺π⁻ and B₀ → K⁺K⁻: The data and the perturbative QCD predictions, Phys. Rev. D 85(9), 094003 (2012), arXiv: 1111.6264 [hep-ph]

[11] A.V. Manohar, Large N QCD, arXiv: hep-ph/9802419 (1998)

[12] Z.H. GuoJ. A. OllerJ.Ruiz. de Elvira, Chiral dynamics in form factors, spectral-function sum rules, meson−meson scattering and semilocal duality, Phys. Rev. D 86(5), 054006 (2012), arXiv: 1206.4163 [hep-ph]

[13] L.Y. DaiU. G. Meißner, A new method to study the number of colors in the final-state interactions of hadrons, Phys. Lett. B 783, 294 (2018), arXiv: 1706.10123 [hep-ph]

[14] D. Tadić, J. Trampetić. Weak meson decays and the expansion. Phys. Lett. B, 1982, 114(2−3): 179

[15] M. Bauer, B. Stech. Exclusive D-decays. Phys. Lett. B, 1985, 152(5−6): 380

[16] A.J. BurasJ. M. GerardR.Ruckl, expansion for exclusive and inclusive charm decays, Nucl. Phys. B 268(1), 16 (1986)

[17] I.E. Halperin, Soft gluon suppression of 1/N_c contributions in color suppressed heavy meson decays, Phys. Lett. B 349, 548 (1995), arXiv: hep-ph/9411422

[18] H.Y. ChengX. W. KangF.Xu, Singly Cabibbo-suppressed hadronic decays of Λc+, Phys. Rev. D 97(7), 074028 (2018), arXiv: 1801.08625 [hep-ph]

[19] A.N. KamalA. B. SantraT.UppalR.C. Verma, Nonfactorization in hadronic two-body Cabibbo-favored decays of D⁰ and D⁺, Phys. Rev. D 53, 2506 (1996), arXiv: hep-ph/9503390

[20] L.ZhangX. W. KangX.H. GuoL.Y. DaiT.Luo C.Wang, A comprehensive study on the semileptonic decay of heavy flavor mesons, J. High Energy Phys. 02, 179 (2021), arXiv: 2012.04417 [hep-ph]

[21] G.BuchallaA. J. BurasM.E. Lautenbacher, Weak decays beyond leading logarithms, Rev. Mod. Phys. 68, 1125 (1996), arXiv: hep-ph/9512380

[22] T.FeldmannP. KrollB.Stech, Mixing and decay constants of pseudoscalar mesons, Phys. Rev. D 58, 114006 (1998), arXiv: hep-ph/9802409

[23] H.Y. ChengX. W. Kang, Branching fractions of semileptonic D and D_s decays from the covariant light-front quark model, Eur. Phys. J. C 77 (9), 587 (2017) [Erratum: Eur. Phys. J. C 77(12), 863 (2017)], arXiv: 1707.02851 [hep-ph]

[24] H.Y. ChengX. W. Kang, The semileptonic baryonic decay D_s⁺ → p p ¯e⁺ν_e, Phys. Lett. B 780, 100 (2018), arXiv: 1712.00566 [hep-ph]

[25] J.HietalaD. Cronin-HennessyT.Pedlar I.Shipsey, Exclusive D_s semileptonic branching fraction measurements, Phys. Rev. D 92(1), 012009 (2015), arXiv: 1505.04205 [hep-ex]

[26] M.AblikimM. N. AchasovS.AhmedM.AlbrechtM.Alekseev, ., Measurement of the dynamics of the decays D_s⁺ → η(′)e⁺ν_e, Phys. Rev. Lett. 122(12), 121801 (2019), arXiv: 1901.02133 [hep-ex]

[27] J.L. RosnerS. StoneR.S. Van de Water, Leptonic decays of charged pseudoscalar mesons – 2015, arXiv: 1509.02220 [hep-ph] (2015)

[28] D.EbertR. N. FaustovV.O. Galkin, Relativistic treatment of the decay constants of light and heavy mesons, Phys. Lett. B 635, 93 (2006), arXiv: hep-ph/0602110

[29] A.BazavovC. BernardN.BrownC.DetarA.X. El-KhadraE.GámizS.GottliebU.M. HellerJ.KomijaniA.S. KronfeldJ.Laiho P.B. MackenzieE.T. NeilJ.N. SimoneR.L. SugarD.ToussaintR.S. Van de Water, B- and D-meson leptonic decay constants from four-flavor lattice QCD, Phys. Rev. D 98(7), 074512 (2018), arXiv: 1712.09262 [hep-lat]

[30] P.MarisP. C. Tandy, Bethe−Salpeter study of vector meson masses and decay constants, Phys. Rev. C 60, 055214 (1999), arXiv: nucl-th/9905056

[31] P.BallR. Zwicky, B_d,s→ρ, ω, K^∗, ϕ decay form factors from light-cone sum rules revisited, Phys. Rev. D 71, 014029 (2005), arXiv: hep-ph/0412079

[32] P.BallG. W. JonesR.Zwicky, B→V gamma beyond QCD factorisation, Phys. Rev. D 75, 054004 (2007), arXiv: hep-ph/0612081

[33] J.C. R. BlochY.L. KalinovskyC.D. RobertsS.M. Schmidt, Describing a₁ and b₁ decays, Phys. Rev. D 60, 111502 (1999), arXiv: nucl-th/9906038

[34] H.Y. ChengB. Tseng, Nonfactorizable effects in spectator and penguin amplitudes of hadronic charmless B decays, Phys. Rev. D 58, 094005 (1998), arXiv: hep-ph/9803457

[35] A.AliC. Greub, An analysis of two-body non-leptonic B decays involving light mesons in the standard model, Phys. Rev. D 57, 2996 (1998), arXiv: hep-ph/9707251

[36] R.AkhouryJ. M. Frere, η, η′ mixing and anomalies, Phys. Lett. B 220(1–2), 258 (1989)

[37] Z.XingH. ZhangS.Zhou, Updated values of running quark and lepton masses, Phys. Rev. D 77(11), 113016 (2008), arXiv: 0712.1419 [hep-ph]

[38] P. A. Zyla, . [Particle Data Group]. . Review of particle physics. Prog. Theor. Exp. Phys., 2020, 2020(8): 083C01

[39] A.BiswasN. SinhaG.Abbas, Nonleptonic decays of charmed mesons into two pseudoscalars, Phys. Rev. D 92(1), 014032 (2015), arXiv: 1503.08176 [hep-ph]

[40] Y.H. ChenH. Y. ChengB.TsengK.C. Yang, Charmless hadronic two-body decays of B_u and B_d mesons, Phys. Rev. D 60, 094014 (1999), arXiv: hep-ph/9903453

[41] H.LiC.D. LüF.S. Yu, Branching ratios and direct CP asymmetries in D → PP decays, Phys. Rev. D 86(3), 036012 (2012), arXiv: 1203.3120 [hep-ph]

[42] Q.QinH. N. LiC.D. LüF.S. Yu, Branching ratios and direct CP asymmetries in D→PV decays, Phys. Rev. D 89(5), 054006 (2014), arXiv: 1305.7021 [hep-ph]

[43] S.H. ZhouY. B. WeiQ.QinY.LiF.S. Yu C.D. Lü, Analysis of two-body charmed b meson decays in factorization-assisted topological-amplitude approach, Phys. Rev. D 92(9), 094016 (2015), arXiv: 1509.04060 [hep-ph]

[44] S.H. ZhouQ. A. ZhangW.R. LyuC.D. Lü, Analysis of charmless two-body B decays in factorization assisted topological amplitude approach, Eur. Phys. J. C 77(2), 125 (2017), arXiv: 1608.02819 [hep-ph]

[45] H.Y. ChengC. W. Chiang, Two-body hadronic charmed meson decays, Phys. Rev. D 81(7), 074021 (2010), arXiv: 1001.0987 [hep-ph]

[46] H.Y. ChengC. W. ChiangA.L. Kuo, Updating B → PP, VP decays in the framework of flavor symmetry, Phys. Rev. D 91(1), 014011 (2015), arXiv: 1409.5026 [hep-ph]

[47] Y.H. ChenH. Y. ChengB.Tseng, Charmless hadronic two-body decays of the B_s mesons, Phys. Rev. D 59, 074003 (1999), arXiv: hep-ph/9809364

[48] Z.Q. WangX. W. KangJ.A. OllerL.Zhang, Analysis on the composite nature of the light scalar mesons f₀ (980) and a₀ (980), Phys. Rev. D 105(7), 074016 (2022), arXiv: 2201.00492 [hep-ph]

[49] L.Y. DaiX. W. KangU.G. Meißner, Test of semilocal duality in a large NC framework, Phys. Rev. D 98(7), 074033 (2018), arXiv: 1808.05057 [hep-ph]

[50] X.W. KangB. KubisC.HanhartU.G. Meißner, ^l4B decays and the extraction of |V_ub|, Phys. Rev. D 89(5), 053015 (2014), arXiv: 1312.1193 [hep-ph]

[51] L. Y. DaiX. W. Kang, U. G. Meißner, X. Y. Song , D. L. Yao., Amplitude analysis of the anomalous decay η′ → π⁺π⁻γ, Phys. Rev. D 97(3), 036012 (2018), arXiv: 1712.02119 [hep-ph]

[52] L.Y. DaiM. R. Pennington, Comprehensive amplitude analysis of γγ → π⁺π⁻, π₀π₀ and K ¯K below 1.5 GeV, Phys. Rev. D 90(3), 036004 (2014), arXiv: 1404.7524 [hep-ph]

[53] D.S. DuM. Z. YangD.Z. Zhang, Space-like Penguin effects in B→PP decays, Phys. Rev. D 53, 249 (1996), arXiv: hep-ph/9507280

[54] A.AliG. KramerC.D. Lü, Experimental tests of factorization in charmless non-leptonic two-body B decays, Phys. Rev. D 58, 094009 (1998), arXiv: hep-ph/9804363

[55] M.BenekeM. Neubert, QCD factorization for B→PP and B→PV decays, Nucl. Phys. B 675, 333 (2003), arXiv: hep-ph/0308039

[56] M.BenekeG. BuchallaM.NeubertC.T. Sachrajda, QCD factorization for exclusive, non-leptonic B meson decays: General arguments and the case of heavy-light final states, Nucl. Phys. B 591, 313 (2000), arXiv: hep-ph/0006124

[57] C.Amsler, The Quark Structure of Hadrons, Lecture Notes in Physics 949, Springer, 2018

[58] Section “Quark Model” in Ref. [38]

[59] H.Y. ChengC. K. ChuaC.W. Hwang, Covariant light-front approach for s-wave and p-wave mesons: Its application to decay constants and form factors, Phys. Rev. D 69(7), 074025 (2004), arXiv: hep-ph/0310359

[60] W. Jaus. Covariant analysis of the light-front quark model. Phys. Rev. D, 1999, 60(5): 054026

RIGHTS & PERMISSIONS

Higher Education Press