Cabibbo-Kobayashi-Maskawa-favored B decays to a scalar meson and a D meson.

In this work, we attempt to study the Cabibbo-Kobayashi-Maskawa-favored [Formula: see text] ("S" denoting the scalar meson) decays within the perturbative QCD approach at the leading order and the leading power. Although the light scalar mesons are widely perceived as primarily the four-quark bound states, in practice it is hard for us to make quantitative predictions based on the four-quark picture for light scalars. Hence, we calculate the decays with light scalars in the two-quark model. For the decays with scalar mesons above 1 GeV, we have explored two possible scenarios, depending on whether the light scalars are treated as the lowest lying [Formula: see text] states or four-quark particles. In total, we calculated the branching fractions of 72 decay modes, and most of them are in the range [Formula: see text]-[Formula: see text], which are measurable in the on-going LHCb experiment and the forthcoming Belle-II experiment. Moreover, since in the standard model these decays occur only through tree operators and have no CP asymmetries, any deviation will be a signal of new physics beyond the standard model. Despite large uncertainties induced by nonperturbative parameters and corrections of high order and high power, our results and discussions will be useful for the on-going LHCb and the forthcoming Belle-II experiments.

Introduction

Even though the quark–antiquark model works well for the pseudoscalar mesons and vector mesons, the study of the inner substructure of the scalar mesons is quite non-trivial, because the conventional quark–antiquark model cannot explain the properties, such as the decay rates and the mass spectrum, especially for ones below 1 GeV. Therefore, the understanding of the internal structure of the scalar mesons is one of the most interesting topics in hadron physics.

The scalar mesons reported by experiments include the isosinglet (), , , and , the isodoublet () and , and the isovector and [1]. Studies of the mass spectrum of scalar mesons and their strong as well as electromagnetic decays suggest that the scalar mesons with the mass below 1 GeV constitute one nonet, while those near 1.5 GeV form another one [2-6]. Irrespective of the existence of the and the mesons, in the literature, the scalar mesons have been identified as ordinary states, four-quark states or meson–meson bound states or even those supplemented with a scalar glueball. Unfortunately, we have not reached a definite conclusion yet till now, due to the unknown nonperturbative properties of QCD, though many efforts have been made to interpret the quark contents of the scalar mesons [7-11]. Now, a consistent picture [2] provided by the data suggests that the scalar meson states above 1 GeV can be identified as a conventional nonet with some possible glue content, which has widely been accepted. However, the quark structure of the light scalar mesons below or near 1 GeV has been quite controversial, though they are widely perceived as primarily four-quark bound states. For example, has been treated as a traditional state [12-14], as a four-quark state [15, 16], and even as a bound state of hadrons [17, 18]. In fact, even in the two-quark picture, the quark component is still unclear; for example, the observation of decay introduced the probability of the component of , while indicated the existence of the non-strange components [19, 20]. Therefore, if the light scalar mesons are dominated by two-quark states, the isoscalars and perhaps should be mixing states like the – system [21, 22]. In the literature, according to the category that the light mesons belong to, there are two typical scenarios for describing the scalar mesons [5]. The scenario 1 (S1) is the naive two-quark model: the nonet mesons below 1 GeV, such as , , , and , are treated as the lowest lying states, and those near 1.5 GeV, such as , , , are the first orbitally excited states. In scenario 2 (S2), the nonet mesons near 1.5 GeV are viewed as the lowest lying states, while the mesons below 1 GeV may be the exotic states beyond the quark model such as four-quark bound states.

In hadron physics, most studies of the light scalar mesons are concentrated on the decay properties of the scalar mesons and the production of the scalar mesons in (or np) collisions or the radiative decays [23]. After the first B decay into a scalar meson, , was observed by Belle [24] and confirmed by BaBar [25], the studies of the scalar mesons in hadronic B decays have attracted more attention because of the large phase space of the B decays. Theoretically, from the studies of QCD sum rules, one finds that the nonperturbative parameters, such as decay constants and the distribution amplitudes, are related to the certain scenarios seriously, which affect the experimental observables noticeably, such as branching fractions and CP asymmetries. So, by comparing the experimental data and the theoretical predictions, one could deduce which scenario is favorable. Motivated by this viewpoint, many charmless B decays with scalars have been studied extensively [5, 6, 26–47].

Recently, the LHCb collaboration reported their first measurements of the decays and [48, 49],Note that in these three decay modes, not only the scalars but also the charmed mesons are involved, and these decays are induced by the transition. Very recently, in Ref. [50], we have attempted to study the decays induced by the transition within the perturbative QCD (PQCD) approach, which are suppressed by the Cabibbo–Kobayashi–Maskawa (CKM) matrix element , but which evade the suppression by the vector decay constants of the scalar mesons. Now, since some experimental data are available, it is worthwhile for us to extend our study to the CKM favored decays induced by the transition.

To handle the hadronic effects for decays, we use the factorization formalism, called perturbative QCD approach [51-54], which is based on the factorization and the transition matrix element is described by the convolution of hadron wave functions and the hard kernel. In the limit of heavy quarks, in order to guarantee that color transparency mechanism is satisfied, i.e., no soft gluon exchange occurring between the final states, we argue that the following hierarchy must be postulated [55, 56]:with . The relation justifies the perturbative analysis of the form factors at large recoil and the definition of light-cone meson wave functions. The relation justifies the power expansion in the parameter . The small ratio is regarded as being of higher power. Because of the inclusion of parton transverse degrees of freedom, large double logarithmic corrections appear and should be summed to all orders. It turns out that the resultant Sudakov factor for an energetic meson is similar to that for a B meson. Including the Sudakov effects from resummation and from threshold resummation for hard amplitudes, the end-point singularities do not exist, and soft contributions can be suppressed effectively. Although the applicability PQCD is still in controversy [57-60], it has been employed for studying the two-body charmed B decays, such as decays [56, 61–66, 97], where P, V, A, T denote the pseudoscalar, vector, axial-vector, and tensor mesons, respectively. Most of the predictions were in good agreement with the present experimental data. Compared with decays, for some decays, the factorizable amplitude will vanish or will be heavily suppressed due to the vanishing or tiny vector decay constants of the scalar mesons. However, the hard-scattering emission diagrams and annihilation type diagrams perhaps provide sizable contributions; such cases are similar to the decays [56, 61–66, 97]. We thus expect that the branching fractions of some decays are large enough to be measured in the current LHC experiment and/or the forthcoming Belle-II in the future. It is worth pointing out that the annihilation type diagrams can be perturbatively calculated in the PQCD approach without end-point singularity, and the predictions of some pure annihilation decays are well in agreement with data, such as the and [67-69].

It should be stressed that although many experimental data indicate that the light scalar mesons, such as and , are predominately four-quark states, in practice it is very difficult for us to make quantitative predictions on based on the four-quark picture for S, because both the decay constants and the distribution amplitudes of S are beyond the conventional quark model. Hence, in practice we shall assume the two-quark scenario for light scalar mesons in the current work.

This paper is organized as follows: we will give a brief review of the formalism of the PQCD approach and specify the wave functions of the initial and final states in Sect. 2. The perturbative calculations and the analytic formulas are given in Sect. 3. The numerical results and phenomenological discussions will be presented in Sect. 4. The final section is reserved for the summary.

Formalism and wave function

As aforementioned, based on the factorization, PQCD approach can effectively avoid the end-point singularity by keeping the intrinsic transverse momenta of inner quarks. The kept transverse momenta will introduce the additional energy scale and lead to the double logarithms appearing in the QCD radiative corrections, which can be resummed into the Sudakov factor. As a result, the Sudakov factor will suppress the soft region contribution and make the calculation in the PQCD approach reliable and consistent.

The effective Hamiltonian related to the decays is given as [70]with the four-quark tree-diagram operatorswhere and are the color indices, and . and are the CKM matrix elements. are the so-called Wilson coefficients at the renormalization scale .

To deal with the hadronic B decays with multiple scales, the factorization hypothesis is usually adopted. The physics regime higher than the scale of the W boson mass is electroweak and can be calculated perturbatively. Using renormalization group techniques, we can evaluate the dynamical effects and obtain the Wilson coefficients from the scale to the b quark mass () scale. The physics between scale and the factorization scale (t) can be calculated perturbatively, which is the so-called hard kernel in the PQCD approach. The dynamics below the factorizable scale is soft and nonperturbative but universal, which can be described by the hadronic wave function. So, in the PQCD approach, the decay amplitude can be written as the convolution of the Wilson coefficients C(t), the hard kernel , and the hadronic wave functions [71, 72],where denotes the trace over Dirac and color indices, the and are the longitudinal momentum fractions and conjugate variables of of the valence quarks in each meson, respectively. The threshold resummation of the double logarithms lead to the jet function [73]. The aforementioned Sudakov factor , coming from the resummation of the double logarithms, can suppress the soft dynamics effectively, i.e. the long distance contributions in the small region [74, 75].

In PQCD, the most important inputs are the wave functions. For the scalar meson, in the two-quark picture, the wave function can be defined aswith the lightlike vectors and . and are the leading-twist and twist-3 light-cone distribution amplitudes respectively, where x is the momentum fraction of the “quark”. The leading-twist light-cone distribution amplitude of the scalar meson has the general form [5, 6]with the Gegenbauer moments and the Gegenbauer polynomials . For the twist-3 distribution amplitudes, we adopt the asymptotic forms for simplicity,The and are the vector decay constant and scalar decay constant of the scalar meson, respectively. For the neutral scalar mesons (, , and ), the vector decay constant vanishes due to the conservation of the vector current. However, the scalar decay constant , related by the equationremains finite. Note that the above parameters, , , and , depend on certain scenarios, for the numerical results of which one is referred to Refs. [5, 6].

For the neutral light scalars and , in the two-quark model, much experimental evidence indicates that there is a mixing between and , which is similar to the – system,with and . For the mixing angle , various experimental measurements have provided different values [76-78]. Recently, the LHCb has proposed the upper limit by the process [79]. Analyzing the present experimental implications, we prefer to adopt the two possible ranges of and [80-85]. It is noted that, for the and mesons, there are other interpretations, for example, the generalized distribution amplitudes [86, 87].

For the – system, according to Ref. [88] and neglecting the tiny contribution from the scalar glueball [89-91], the mixing form can be simplified asFor the initial B meson, neglecting the numerically suppressed Lorentz structure, the remaining leading order wave function can be decomposed as [92, 93]withThe contribution from starts from the next-to-leading-power , which should be included together with other next-to-leading-power contributions in order to form a complete analysis. So, we here neglected it in the current work. The light-cone distribution amplitude can be written as [92-94]with the normalization constant , which can be determined through the following normalization condition:For the shape parameter and the decay constant , we will take GeV and GeV for the B meson, respectively, and take GeV and GeV for the meson, due to the SU(3) breaking effects [51-54].

In terms of the heavy quark limit, the two-parton light-cone distribution amplitudes of the meson will be taken as [95-98] with the distribution amplitudes [96-98]with the shape parameter GeV. Note that the high-twist distribution amplitudes are not included either, because they are suppressed by . We choose and MeV for the D meson, and and MeV for the meson [99]. The parameters are fitted from the and decays [96-98]. For , the decay constants can be obtained through the relation based on the heavy quark effective theory, which can be found in Refs. [61–64, 97].

Perturbative calculation

In this section, within the PQCD approach, we specifically calculate the decay amplitudes without the Wilson coefficients in Eq. (5) for each Feynman diagram, and we express the calculated amplitudes as a convolution of the hard kernel and the mesons’ wave functions. It is noted that there are two kinds of diagrams contributing to the considered decays at the leading order. The diagrams with a meson emitted are presented in Fig. 1, and those with a scalar meson emitted are listed in Fig. 2.

Fig. 1

Leading order Feynman diagrams contributing to the decays in PQCD approach

Fig. 2

Leading order Feynman diagrams contributing to the decays with a scalar meson emitted in PQCD

For the two factorizable emission diagrams (a) and (b) in Fig. 1, the amplitudes can be written aswhere , , and the color factor for B decays. The expressions for the scale t, the Sudakov factor E, and the hard functions h from the denominator of the propagators can be found in Appendix A of Ref. [63].

The two diagrams (c) and (d) in Fig. 1 are the so-called hard-scattering emission diagrams. Compared to the previous two, each decay amplitude involves three meson wave functions. After integrating out with the function , the amplitudes for these two diagrams can be expressed byThe four diagrams in the second row are the annihilation type diagrams, which can be perturbatively calculated in the PQCD approach. Diagrams (e) and (f) are the factorizable diagrams with the B meson factorized out, and the amplitudes can be written asFor the nonfactorizable diagrams (g) and (h), the corresponding amplitudes are given as follows:

For these diagrams with a scalar meson emitted in Fig. 2, the decay amplitudes are expressed as From Eq. (23), one finds that the amplitudes of factorizable emission diagrams with a scalar meson emitted are proportional to the vector decay constant of the scalar, so it will be highly suppressed by the tiny vector decay constant of the scalar meson, or even vanishes for the decays emitting a neutral scalar meson, because the neutral scalar meson cannot be produced through current.

For the decays, there are only the longitudinal polarization contributions in the light of the conservation of angular momentum. The expressions of the factorizable emission contributions can be obtained by the following substitutions in Eq. (19):For the hard-scattering emission diagrams, the decay amplitudes can be expressed asSimilarly, the annihilation type contributions can be written as For those diagrams with a scalar meson emitted in decays, the factorizable emission contributions can be obtained from Eq. (23) directly by adopting the same substitutions as Eq. (25), and the hard-scattering emission contributions can be expressed asAs for the Wilson coefficients of each decay modes, they are the same as those of decays, and can be found in Ref. [63].

Numerical results and discussions

In this section, we will present the numerical results and some phenomenological analyses of the considered decays. First of all, we should list the input parameters in our numerical calculations [1]:For the decay constants of the scalar mesons calculated within the QCD sum rules, we adopt the values in Ref. [6].

In Tables 1, 2, 3 and 4, we tabulate our numerical results of the branching fractions with uncertainties. Considering the four-quark picture beyond the conventional quark model, we cannot calculate the decay constants and the distribution amplitudes within the QCD sum rules and the lattice QCD approach. So, for the light scalars below 1 GeV, we only list the results of S1, i.e., the two-quark scenario. For the heavy scalars, the numerical results in the two scenarios are both presented. Indeed, previous studies [5, 6] indicated that S2 is favorable. Of course, there are many uncertainties in our calculations, and here we mainly evaluate three kinds. The first errors are caused by nonperturbative parameters, such as the decay constants , , , the shape parameters in distribution amplitudes of B / D mesons, the Gegenbauer moments in the distribution amplitudes of the scalar mesons. The second ones are from the unknown next-to-leading order (NLO) corrections, characterized by the choice of the and the variations of the factorization scales t (). NLO contributions in PQCD are not available. In Refs. [100, 101], the authors have estimated the NLO effects in decays by including the NLO form factors, NLO Wilson coefficients, vertex corrections, the quark loops, and the chromomagnetic penguin. The results showed that the NLO effects are modest and the uncertainties due to the scale variation in the LO are reasonable. The last errors come form the uncertainties of the CKM matrix elements listed in Eq. (30). In our calculations, we admit that the next-leading power () corrections are not included, which have been estimated to be no more than 10% [55, 56]. Moreover, the next-leading power corrections of are even larger. Yet now, all next-leading power corrections are not available, so we will not include them here and leave them as our future work. From the tables, it is apparent that the most significant theoretical uncertainties are from the nonperturbative parameters. In these tables, in order to indicate the dominant contributions, we also mark each channel by the symbols “T” (color-allowed tree contributions), “C” (color-suppressed tree contributions), and “E” (W exchange type contributions). Because all of these decays only occur through tree operators, the CP asymmetries of these decays are null in SM.

Table 1

Branching fractions of decays calculated in the PQCD approach in S1

Decay modes	Class	BRs(10^-6)
B^+→D¯^0a0+	C	483-215-68-12+244+64+26
B^0→D^-a0+	T	17.6-7.6-7.7-0.4+9.8+11.4+0.9
B^0→D¯^0a00	C	160-75-31-4+88+35+9
B^0→D¯^0σ(f_n)	C	134-55-39-3+65+32+7
B^0→D¯^0f_0(f_n)	C	78.4-36.2-36.8-1.9+42.8+33.6+4.3
B^0→Ds-κ^+	E	72.6-22.1-8.8-1.8+24.3+8.0+4.0
B_s→D¯^0κ¯	C	262-131-66-7+154+57+14
B_s→Ds-a0+	T	64.3-28.7-10.2-1.6+42.3+11.0+3.5
B^+→D¯^0κ^+	C	10.8-6.4-1.9-0.3+7.9+1.9+0.7
B^0→D^-κ^+	T	4.83-1.64-0.63-0.14+1.92+0.73+0.29
B^0→D¯^0κ	C	6.89-4.28-2.51-0.21+5.26+2.36+0.41
B_s→D^-a0+	E	3.42-1.17-0.44-0.10+1.31+0.41+0.21
B_s→D¯^0a_0	E	1.70-0.58-0.21-0.05+0.65+0.20+0.10
B_s→D¯^0σ(f_n)	E	1.13-0.40-0.14-0.04+0.42+0.17+0.07
B_s→D¯^0σ(f_s)	C	14.2-6.7-2.7-0.4+7.9+2.5+0.9
B_s→D¯^0f_0(f_n)	E	1.36-0.45-0.15-0.04+0.51+0.15+0.08
B_s→D¯^0f_0(f_s)	C	10.6-5.3-2.6-0.3+6.1-2.5-0.6
B_s→Ds-κ^+	T	1.02-0.47-0.29-0.03+0.59+0.29+0.06

Table 2

Branching fractions of , and calculated in the PQCD approach in S1 and S2, respectively

Decay modes	Class	BRs(10^-5)S_1	BRs(10^-5)S_2
B^+→D¯^0a0+(1450)	C	72.1-31.6-12.5-1.8+39.2+13.9+3.9	123-64-14-4+73+12+5
B^0→D^-a0+(1450)	T,E	4.09-2.09-0.86-0.10+2.63+1.59+0.23	0.92-0.43-0.24-0.02+0.57+0.41+0.05
B^0→D¯^0a_0(1450)	C	31.3-14.3-5.6-0.8+17.4+4.8+1.7	66.2-32.8-7.1-1.7+37.6+6.4+3.6
B^0→D¯^0f_0(1370)(f_n)	C	28.6-12.0-3.4-0.8+13.7+2.5+1.5	16.3-7.6-4.0-0.4+9.1+4.9+0.8
B^0→D¯^0f_0(1500)(f_n)	C	27.2-11.1-2.8-0.7+13.2+2.6+1.5	13.2-5.9-2.6-0.3+7.4+3.9+0.7
B^0→Ds-K0∗+(1430)	E	1.41-0.46-0.46-0.03+0.53+0.46+0.08	9.08-3.76-0.58-0.22+4.21+0.63+0.50
B_s→D¯^0K¯0∗(1430)	C	53.9-23.0-6.2-1.3+26.8+3.8+2.9	68.8-35.4-8.5-1.7+38.3+7.6+3.8
B_s→Ds-a0+(1450)	T	10.3-4.9-1.8-0.3+5.8+1.9+0.5	4.35-2.43-0.97-0.11+3.28+1.21+0.23
B^+→D¯^0K0∗+(1430)	C	4.72-1.83-0.94-0.15+2.22+1.03+0.28	4.96-3.15-0.39-0.15+4.02+0.98+0.29
B^0→D^-K0∗+(1430)	T	0.97-0.36-0.12-0.03+0.45+0.14+0.06	0.79-0.38-0.11-0.03+0.46+0.14+0.05
B^0→D¯^0K0∗0(1430)	C	3.39-1.27-0.46-0.10+1.47+0.43+0.21	3.19-2.34-0.40-0.09+3.21+0.44+0.20
B_s→D^-a0+(1450)	E	0.14-0.07-0.01-0.01+0.06+0.01+0.01	0.43-0.18-0.03-0.01+0.21+0.02+0.02
B_s→D¯^0a00(1450)	E	0.07-0.03-0.01-0.01+0.03+0.01+0.01	0.21-0.09-0.01-0.01+0.11+0.02+0.02
B_s→D¯^0f_0(1370)(f_n)	E	0.05-0.02-0.01-0.01+0.04+0.01+0.01	0.17-0.07-0.01-0.01+0.09+0.01+0.01
B_s→D¯^0f_0(1370)(f_s)	C	2.30-1.12-0.25-0.07+1.46+0.15+0.14	2.97-1.53-0.32-0.09+2.99+0.30+0.18
B_s→D¯^0f_0(1500)(f_n)	E	0.05-0.02-0.01-0.01+0.04+0.01+0.01	0.17-0.07-0.01-0.01+0.10+0.01+0.01
B_s→D¯^0f_0(1500)(f_s)	C	2.24-1.08-0.23-0.07+1.39+0.18+0.13	2.71-2.04-0.22-0.08+2.81+0.26+0.17
B_s→Ds-K0∗+(1430)	T,E	0.36-0.18-0.08-0.01+0.24+0.13+0.02	0.39-0.16-0.04-0.01+0.21+0.04+0.02

Table 3

Branching fractions of decays calculated in the PQCD approach in S1

Decay modes	Class	BRs(10^-6)
B^+→D¯^∗0a0+	C	520-188-160-13+215+127+29
B^0→D^∗-a0+	T	250-84-41-7+91+63+13
B^0→D¯^∗0a_0	C	128-67-40-3+76+34+7
B^0→D¯^∗0σ(f_n)	C	171-70-54-4+78+45+9
B^0→D¯^∗0f_0(f_n)	C	119-51-50-4+57+48+6
B^0→Ds∗-κ^+	E	13.0-4.3-2.5-0.4+4.5+1.8+0.7
B_s→D¯^∗0κ¯	C	320-153-80-8+178+69+17
B_s→Ds∗-a0+	T	169-63-52-4+73+56+10
B^+→D¯^∗0κ^+	C	8.80-4.41-3.67-0.27+5.78+3.09+0.53
B^0→D^∗-κ^+	T	3.42-1.32-0.98-0.10+1.58+1.44+0.21
B^0→D¯^∗0κ	C	9.25-5.03-3.55-0.27+6.54+3.03+0.56
B_s→D^∗-a0+	E	0.76-0.28-0.15-0.03+0.30+0.14+0.04
B_s→D¯^∗0a_0	E	0.38-0.15-0.08-0.02+0.14+0.06+0.02
B_s→D¯^∗0σ(f_n)	E	0.19-0.08-0.04-0.01+0.07+0.03+0.01
B_s→D¯^∗0σ(f_s)	C	16.4-7.6-3.4-0.5+8.8+2.9+1.0
B_s→D¯^∗0f_0(f_n)	E	0.27-0.10-0.06-0.01+0.10+0.04+0.01
B_s→D¯^∗0f_0(f_s)	C	12.7-6.0-3.4-0.3+7.1+3.1+0.8
B_s→Ds∗-κ^+	T	6.83-2.92-1.80-0.20+3.43+2.08+0.41

Table 4

Branching fractions of , and calculated in the PQCD approach in S1 and S2, respectively

Decay modes	Class	BRs(10^-5)(S_1)	BRs(10^-5)(S_2)
B^+→D¯^∗0a0+(1450)	C	207-81-23-5+90+24+11	98.2-56.1-29.2-2.5+70.1+28.9+5.3
B^0→D^∗-a0+(1450)	T	26.8-11.1-7.3-0.6+12.7+7.5+1.5	11.3-5.9-2.7-0.2+7.4+2.1+0.7
B^0→D¯^∗0a_0(1450)	C	40.1-17.3-5.3-1.0+21.2+5.2+2.1	58.9-31.2-8.4-1.4+36.7+7.7+3.3
B^0→D¯^∗0f_0(1370)(f_n)	C	44.8-18.4-4.0-1.1+21.3+2.1+2.4	27.8-21.4-6.0-0.7+32.0+6.7+1.6
B^0→D¯^∗0f_0(1500)(f_n)	C	44.5-18.1-3.4-1.0+21.2+2.6+2.5	25.1-20.1-4.6-0.6+30.0+5.7+1.4
B^0→Ds∗-K0∗+(1430)	E	0.53-0.17-0.17-0.01+0.22+0.19+0.03	0.99-0.49-0.21-0.03+0.56+0.28+0.05
B_s→D¯^∗0K¯0∗(1430)	C	73.0-30.7-5.9-1.8+36.5+4.0+4.0	79.0-55.6-11.6-2.0+72.0+9.6+4.3
B_s→Ds∗-a0+(1450)	T	23.0-10.1-7.8-0.6+11.1+8.9+1.2	6.20-3.62-1.64-0.16+4.58+1.46+0.33
B^+→D¯^∗0K0∗+(1430)	C	10.1-3.3-1.1-0.3+3.9+1.3+0.6	3.08-2.25-0.71-0.08+3.62+0.82+0.18
B^0→D^∗-K0∗+(1430)	T	0.75-0.24-0.27-0.02+0.26+0.32+0.05	0.09-0.07-0.02-0.01+0.14+0.02+0.01
B^0→D¯^∗0K0∗0(1430)	C	4.90-1.81-0.51-0.14+2.08+0.44+0.30	3.80-2.79-0.60-0.11+3.80+0.61+0.23
B_s→D^∗-a0+(1450)	E	0.06-0.03-0.02-0.01+0.03+0.01+0.01	0.05-0.02-0.01-0.01+0.03+0.01+0.01
B_s→D¯^∗0a00(1450)	E	0.03-0.01-0.01-0.01+0.01+0.01+0.01	0.03-0.01-0.01-0.01+0.01+0.01+0.01
B_s→D¯^∗0f_0(1370)(f_n)	E	0.02-0.01-0.01-0.01+0.01+0.01+0.01	0.02-0.01-0.01-0.01+0.01+0.01+0.01
B_s→D¯^∗0f_0(1370)(f_s)	C	3.29-1.53-0.21-0.10+1.94+0.18+0.20	3.43-2.56-0.45-0.11+3.50+0.39+0.20
B_s→D¯^∗0f_0(1500)(f_n)	E	0.02-0.01-0.01-0.01+0.01+0.01+0.01	0.02-0.01-0.01-0.01+0.01+0.01+0.01
B_s→D¯^∗0f_0(1500)(f_s)	C	3.25-1.49-0.19-0.09+1.90+0.23+0.20	3.21-2.43-0.36-0.10+3.32+0.34+0.19
B_s→Ds∗-K0∗+(1430)	T	1.40-0.50-0.41-0.04+0.55-0.44-0.09	0.24-0.20-0.08-0.01+0.27+0.06+0.01

From the tables, one can find that, compared with the processes, the processes are all suppressed by the CKM matrix elements . For these T-type decays with a scalar meson emitted, the contributions from factorizable emission diagrams are either suppressed by the tiny vector decay constant of the scalar meson or even vanish for the neutral scalar mesons, though they have large Wilson coefficients. For the two hard-scattering diagrams ((c) and (d) in Fig. 2), because the light-cone distribution amplitude of the scalar meson is antisymmetric, the contributions of these two diagrams no longer cancel but are enhanced, which can be seen from Eqs. (24) and (29). So, although the hard-scattering diagrams are suppressed by the Wilson coefficient , they also provide sizable contributions and even dominate the decay amplitudes in some decay modes. In addition, we note that, for the T-type decays with emission, the contributions from factorizable emission diagrams are still sizable, because the vector decay constants of are not too small due to the mass difference of the up and strange quarks.

We now discuss the C-type decays with a meson emitted. The factorizable emission diagrams are suppressed by the small Wilson coefficients . Since the cancellation between the hard-scattering emission diagrams (c and d in Fig. 1) is suppressed by the mass difference between the quark and the “light” quark in the emitted meson, the contributions of the hard-scattering emission diagrams with the large Wilson coefficient are no longer negligible; even they dominate the decay amplitudes. Therefore, their branching fractions are expected to be large enough to be detected at on-going experiments, especially for these processes.

As is well known, the annihilation type diagrams are power suppressed in the PQCD approach, and the branching fractions of the E-type decays are expected to be much smaller than the others. However, for the decays, because there is a large mass difference between the meson and the scalar meson, which will weaken the cancellation between the two nonfactorizable annihilation type diagrams (g) and (h) in Fig. 1, the contributions of the annihilation diagrams might be sizable. As a result, the branching fractions of these E-type decays are not too small, as usual, especially for these processes. For example, enhanced by the CKM matrix elements, the branching fractions of the even reach , the order of which is measurable in the on-going experiments. When the experimental data is available, it will provide another platform to study the dynamical mechanism of the annihilation diagrams in two-body hadronic B decays.

Specially, and decays have both a T-type contribution with a scalar meson emitted and a C-type contribution with emitted. For decays, the constructive interference between those two contributions makes the branching fractions larger than the pure C-type decays, such as decay. Similarly, the constructive (destructive) interferences also lead the branching fractions of larger (smaller) than the pure C-type decays. In particular, because the vector decay constant of is not tiny, the T-type contributions with a emitted are sizable. From the Table 1, one can also find that and , which can be understood by the interference between T-type and E-type contributions. The relations and in Table 3 can also be explained in the same manner.

From Tables 2 and 4, it is found that, for these C-type decays, such as and decays, the branching fractions in S1 are roughly equal to those in S2. This can be explained by the fact that the two dominant nonfactorizable diagrams (c) and (d) in Fig. 1 will be cancelled out by each other. Thus, the effects caused by the wave functions of scalar mesons are inconspicuous. In fact, such cases also occur in the color-suppressed decays in Ref. [50]. We also note that the branching fraction of in S2 is larger than that in S1. That is because the constructive interference between the C-type contributions with emitted and the T-type contributions with the scalar meson emitted in S2 is much larger than that in S1. However, for the decays and , the destructive interference causes their branching fractions in S2 to be smaller than those in S1. As for , although the interference between the two type contributions is also constructive in S2, its branching fraction in S2 is only slightly larger than that in S1, because the contributions from factorizable diagrams and the nonfactorizable ones are canceled out by each other, especially when the vector decay constant of is no longer as small as the others. We also note that the decays can be used to identify different scenarios. Similarly, for the decays and , the branching fraction differences between the two scenarios are ascribed to the interference between the emission contributions and the annihilation ones.

The T-type decays, which are pure emission processes with emitted, are dominated by the hard-scattering emission diagrams, since the factorizable diagrams are highly suppressed by the small vector decay constant of the . The ratio of branching fractions between S1 and S2 is about 2 and 4, for and , respectively, which indicates that the branching fractions are sensitive to the scenarios. From Eq. (20), it is found that the contributions from the two hard-scattering diagrams are enhanced by each other. When we switch from S1 to S2, the changes induced by the distribution amplitudes in S2 will overlap with each other, which makes the branching fractions different from those in S1. As for decays, the ratio between S1 and S2 is about 4 for and about 2 for , which is contrary to the cases in decays. This is caused by the interferences between the emission diagrams and the annihilation ones. In S2, the interference is destructive for decays, but constructive for . Unlike the cases of the above T-type decays with the , the branching fractions of the and decays in the two scenarios are roughly equal. However, for and decays, the branching fractions in S1 are much larger (about seven–eight times larger) than those in S2. For the above four decays, the color-allowed factorizable emission contributions are sizable in S2, because the vector decay constant of the in S2 is larger than in S1. For and decays, the interference between the above contributions and the ones of the hard-scattering emission diagrams is constructive, so their branching fractions in S2 are roughly equal to those in S1. However, this kind of interference is destructive for and decays, and their branching fractions in S2 are smaller than those in S1.

Now, we turn to a discussion of the pure annihilation decays, which are dominated by the nonfactorizable annihilation diagrams. From Table 2, one can find that the branching fractions of pure E-type decays in S2 are much larger than those in S1. As is well known, the cancellation between two nonfactorizable annihilation diagrams is suppressed by the large mass difference between the b quark and the light quark. So, the changes induced by the distribution amplitudes of the scalars become important, which leads to the fact that the branching fractions are obviously dependent on the scenarios. Taking for illustration, the branching ratios in S2 are about three times larger than that in S1. However, from Table 4, we find that the situation is reversed for the pure annihilation decays, the discrepancies of branching fractions in different scenarios are quit small. Comparing the Eq. (21) with Eq. (27), we notice that the two factorizable annihilation diagrams are cancelled by each other in decays but enhanced in decays. So, in decays, the contributions from two factorizable diagrams are comparable with those from nonfactorizable ones. Moreover, the interference between the factorizale annihilation diagrams and the nonfactorizable ones is destructive (or constructive) in S1 (S2), which causes that the branching fractions in S2 are almost equal to or even larger than those in S1.

In the conventional two-quark picture for the light scalars, although the LHCb experiment had measured the branching fractions of and [48, 49], the mixing angle cannot be constrained stringently due to the large uncertainties. For convenience, we present individually the branching fractions under the pure and components in the tables. Once the S1 is confirmed and the mixing angle is fixed, one can obtain the branching fractions directly from the two predictions with the and components. For instance, if the popular value ranges and are adopted, we can predict the branching fractions as listed in Table 5. In the same manner, by neglecting the tiny glueball contents and adopting the results of Eq. (11), we also list the branching fractions of the decay modes with or in Table 6. Note that we here only list the center values for simplicity.

Table 5

The calculated branching fractions of and with the mixing in the PQCD approach (unit:)

Decay modes	[25^∘,40^∘]	[140^∘,165^∘]
B^0→D¯^0σ	78.6∼110	78.6∼125
B^0→D¯^0f_0(980)	46.0∼64.4	46.0∼73.1
B_s→D¯^0σ	3.99∼7.22	1.66∼5.87
B_s→D¯^0f_0(980)	5.46∼7.92	8.09∼10.6
B^0→D¯^∗0σ	100∼140	100∼159
B^0→D¯^∗0f_0(980)	70.2∼98.2	70.2∼111
B_s→D¯^∗0σ	2.93∼6.70	1.37∼7.09
B_s→D¯^∗0f_0(980)	7.62∼10.5	7.55∼11.9

Table 6

The calculated branching fractions of and with the mixing in the PQCD approach (unit:)

Decay modes	S1	S2
B^0→D¯^0f_0(1370)	173	99.2
B^0→D¯^0f_0(1500)	79.3	38.5
B_s→D¯^0f_0(1370)	7.40	4.85
B_s→D¯^0f_0(1500)	14.9	24.5
B^0→D¯^∗0f_0(1370)	272	168
B^0→D¯^∗0f_0(1500)	129	73.2
B_s→D¯^∗0f_0(1370)	10.7	7.55
B_s→D¯^∗0f_0(1500)	20.8	24.6

In fact, under the two-quark assumption, only the component contributes to the decay modes and . Thus, we can define the ratioUsing the latest experimental data in Eq. (1), we can obtain , which can constrain the range of the mixing angle asCompared with the results of Ref. [80-85], the obtuse angle solutions agree with each other, but the acute angle we obtained is a bit smaller than the previous results. Using the mixing angle value in Eq. (32) and the results in Table 1, we get the branching fractions of :where the errors are only from the mixing angle. Compared to Eq. (1), one can find that our numerical results can accommodate the experimental data well within the limit of errors.

It should be noted that there are large uncertainties in our results and discussions, especially ones induced by the nonperturbative parameters, which might lower our prediction power and blur the distinction between the two scenarios even in favorable cases. Therefore more reliable nonperturbative approaches and experimental data are needed for describing the nature of the scalar mesons in the future.

Summary

Motivated by recent results of the charmed B decays with a scalar from LHCb experiment, we attempted to investigate the decays induced by the transition within the framework of the PQCD approach. Although the light scalar mesons, especially and , are widely perceived as primarily four-quark bound states, in practice it is difficult to make quantitative predictions based on the four-quark picture for light scalars. Hence, predictions are made in the two-quark model for the decays with the light scalar mesons. For the decays with scalar mesons above 1 GeV, we have explored two possible scenarios, the difference depending on whether the light scalars are treated as the lowest lying states or four-quark particles. By comparing with the experimental data, we can deduce whether the heavy scalars are the ground states or the first excited ones, though the current studies prefer the former. Since the considered decays occur only through tree operators, there are no CP asymmetries. The branching fractions of most decay modes are in the range of –, which can be tested in the LHCb experiment and by Belle-II in the near future. Some decays have large branching fractions, such as the and , which are measurable in the current experiment. For the and decays, our numerical results accommodate the experimental data well within the limit of the errors, even in the two-quark picture. We note that although the nonperturbative parameters, the higher-order and high-power correction, and even the effects of the final states interaction lead rather large uncertainties, the orders of magnitude and some phenomenological discussions will be useful for the on-going LHCb and the forthcoming Belle-II experiments.