Semileptonic decay of $B^{-}_c$ into $X(3930)$, $X(3940)$, $X(4160)$

We study the semileptonic decay of $B^{-}_c$ meson into $\bar{\nu} l^-$ and the isospin zero $X(3930)~(2^{++})$, $X(3940)~(0^{++})$, $X(4160)~(2^{++})$ resonances. We look at the reaction from the perspective that these resonaces appear as dynamically generated from the vector-vector interaction in the charm sector, and couple strongly to $D^{*}\bar{D^{*}}$ and $D^{*}_{s}\bar{D_{s}^{*}}$. We also look into the $B^{-}_c \rightarrow \bar{\nu}_{l} l^- D^{*}\bar{D^{*}}$ and $B^{-}_c \rightarrow \bar{\nu}_{l} l^- D^{*}_{s}\bar{D_{s}^{*}}$ reactions close to threshold and relate the $D^{*}\bar{D^{*}}$ and $D^{*}_{s}\bar{D_{s}^{*}}$ mass distribution to the rate of production of the $X$ resonances.


I. INTRODUCTION
The X, Y, Z states, that challenge the constituent quark model picture of meson [1,2] have been one of the most spectacular findings in hadron spectroscopy recently [3][4][5][6]. Their advent has stimulated much theoretical work aimed at unravelling their structure. Tetraquark pictures have been proposed [7,8] as well as molecular pictures stemming from the interaction of more elementary mesons [6,9,10]. One of these pictures deals with the interaction of vector mesons with charm, leading to hidden charm quasibound meson states [10]. In that work, the channels D * D * , D * sD * s , K * K * , ρρ, ωω, φφ, J/ψJ/ψ, ωJ/ψ, φJ/ψ, ωφ were considered and the interaction between them was obtained using an extention of the local hidden gauge approach [11][12][13], exchanging vector mesons, and through contact terms provided by the theory. Some quasibound states were found which could be associated to known resonances. These states were: one state around 3943 MeV with I G [J P C ] = 0 + [0 ++ ], which was associated to the X(3940) [14,15]; another state around 3922 MeV with 0 + [2 ++ ], which was associated to the X(3930) [16] (now classified in the PDG [17] as the the χ c2 (2P ) ), which could also correspond to the X(3915) [18,19], and a third one at 4169 MeV with 0 + [2 ++ ] that was associated to the X(4160) [20].
The X(3940) was found to couple mostly to D * D * in [10], the X(3930) also had its strongest coupling to D * D * and the X(4160) had its strongest coupling to D * sD * s . The light vector-vector channels couple weakly to those states, but given the large space available, they are the biggest source of the width, which in the theoretical work is also found in reasonable agreement with experiment. It is interesting to mention that there is a large list of works suggesting a bound D * sD * s state [21][22][23][24][25][26]. QCD * Electronic address: ikeno@rs.tottori-u.ac.jp † Electronic address: melahat.bayar@kocaeli.edu.tr ‡ Electronic address: oset@ific.uv.es sum rules, although with its usual large uncertainties, have also speculated on this possibility [27,28]. The curious thing is that all these work aimed at reproducing the X(4140) resonance not the X(4160). One can think that the fact that light vector channels were not included as coupled channel in these studies had as a consequence a small width for the resonance which made it more appealing to have it associated to the X(4140). Yet, the quantum numbers 0 + [1 ++ ] for this resonance determined lately make the association of the state found in these works to the X(4140) inappropriate, while the association to the X(4160) is more natural.
The discussion on these states becomes more actual when one recalls the experimental work [29,30] in the B + → J/ψφK + reaction, where the data analysis brought the surprising result that the X(4140) has a width of around Γ ≃ 83 MeV, while former experiments give a width around 19 MeV [31][32][33][34][35][36][37][38]. This puzzle found a reasonable explanation in a recent work [39], where the J/ψφ invariant mass distribution at low invariant masses was analyzed in terms of the X(4140) and X(4160) and the mass distribution was better reproduced.
The result of this new analysis was that the X(4140) has a width compatible with 19 MeV, and it is the X(4160) resonance the one that fills the strength in that region. The striking thing is that, since the X(4160) in that work is supposed to be a D * sD * s bound state, but with a relatively large coupling to J/ψφ in the coupled channels study of [10], when the J/ψφ mass distribution is studied, a large cusp structure develops in this distributions at the D * sD * s threshold and such cusp is present in the experiment.
In view of this puzzling situation, any other reaction that brings light into these issues should be most welcome. This is the purpose of the present work, where we propose to measure the semileptonic decay of B − c →ν l l − X i , with X i any of the three resonances, X(3930), X(3940), X(4160). Actually, this reaction has been studied recently [41] from the perspective that the X(3940) and X(4160) resonances are radial high excitations of the charmonium states, corresponding to the η c (3S) and η c (4S) respectively. Our picture, where these resonances are generated from the interaction of vector meson with charm is quite different and the study of the B c → ν e e + X i decays from this perspective is worth pursuing.
The use of semileptonic weak decays aiming at determining the structure of resonances has been exploited before in different cases. In [42] the B s and B semileptonic decays,B 0 s → D * s0 (2317)ν l l − ,B 0 → D * 0 (2400) +ν l l − were studied and compared to related reactions likē B 0 s → (DK) +ν l l − . In [43] the production of light scalar mesons and light vector mesons was also investigated in semileptonic decays of D and D s mesons. In [44] the Λ c → Λ(1405)e + ν e was studied, looking at the decays of Λ(1405) into π + Σ − , π − Σ + , π 0 Σ 0 andKN production. In this case the weak decay filters I = 0 in the final meson baryon system, which makes this reaction special to investigate the properties of the Λ(1405). Related to these works, but with a different aim, one has the work of [45] where the Λ b → Λ c (2595)ν l l and Λ b → Λ c (2625)ν l l decays are studied in order to test the pseudoscalarbaryon and vector-baryon components of the Λ c (2595) and Λ c (2625) resonances [46].
In the present work we take advantage of these previous studies and calculate the B c → ν l l + X i decay rates and compare them with the B c → ν l l + D * D * , ν l l + D * sD * s decays. We establish a link between these processes which is tied to the molecular nature of these resonances, making predictions to be tested in future experiments from where much valuable information concerning the nature of these states is to be expected.

II. FORMALISM
The B − c →νe − X process proceeds at the quark level through a first step shown in Fig. 1. The process involves the bc weak transition, which is the same one as in the B decays studied in Ref. [42].
There is, however, a novelty in the present process. Indeed, if we want to see two mesons, the cc quarks of Fig. 1 must hadronize into two mesons components. This is easily done for mesons since one introduces an extraqq pair with vacuum quantum number,ūu+dd+ss+cc, and then the two quarks after the weak process participate in the formation of the two mesons. With two quarks after the weak vertex, as in Fig. 2, the newqq pair can be placed in between these cc quarks.
The procedure followed here is inspired in the approach of Ref. [47] where the basic mechanisms at the quark level are investigated, then pairs of hadrons are produced after implementing hadronization, and finally these hadrons are allowed to undergo final state interaction.
The hadronization of cc, introducing the qq pair is done as follows [42,47]. We take the qq matrix M , uū ud us uc dū dd ds dc sū sd ss sc cū cd cs cc We now write M in terms of vector mesons and we have the vector matrix V , Then M 2 becomes V 2 and We neglect the J/ψJ/ψ channel since it has too high energy relative to other channels. Only an I = 0 state is produced from the cc component since the hadronization is a strong interaction and does not change isospin. We can write the D * D * combination in terms of the isospin doublets (D * + , −D * 0 ) and (D * 0 ,D * − ) and then the production vertex is written as [47] (V · V ) 44  The final state interaction of D * D * is depicted in Fig. 3, and we use the interaction of Ref. [10] where, using an extension of the local hidden gauge approach, the interaction of D * D * generates several resonances, and some XY Z states were dynamically generated. As suggested in Refs. [10,47], the resonances most strongly coupled to the D * D * channel correspond to the experimental states Y (3940), Z(3930) and the D * sD * s channel corresponds to the X(4160).

A. Coalescence
From the present perspective, we have studied the semileptonic decay process in Ref. [42]. First, we study the coalescence process which produces the resonances after rescattering, as shown in Fig. 4.
This process has a three-body final state with a lepton, its neutrino and the resonance R. The resonance R stands for the X, Y, Z resonances. The hadronization factor V had can be obtained as for the resonance R in J = 0 which requires L = 0. For the constant C, we use the value C = 7.22 of the semileptonic B decays as established in Ref. [42]. G D * D * and G D * sD * s are the two meson loop functions, and g R,D * D * and g R,D * sD * s are the couplings of the resonance to these channels. We use the values reported in Ref. [10]. The two meson loop function G i for each channel i is where P is the total four-momentum of the two mesons, and m 1 and m 2 are the masses of the two mesons in channel i. We use cut off regularization as done in [39] to avoid potential problems of dimensional regularization pointed out in [40]. The G function has the form where q max stands for the cutoff in the three momentum, the square of center of mass energy (P 0 ) 2 = s and ω i = q 2 i + m 2 i . As in [39] we use q max = 690 MeV. The whole amplitude T Bc for the semileptonic decay of the B c meson is written as, where By following the steps of Refs. [42,43], we find for the sum and average over the polarization of the fermions Further steps are done in Ref. [42] to perform the angular integrations of the resonance in the B c rest frame and the lepton in the νe rest frame, and finally one obtains a formula of the decay widths Γ coal for the coalescence of the X, Y, Z resonances by Here, P cm R is the momentum of the resonance R in the B c rest frame, andp ν is the momentum of the neutrino in the νe rest frame, with the Källen function defined as, The energiesẼ Bc andẼ R are calculated in the νe rest frame,Ẽ andp Bc is the momentum of the B c in the νe rest frame, The integral range of M (νe) inv is [m e + m ν , m Bc − m R ]. We take the Fermi coupling constant G F = 1.166 × 10 −5 GeV −2 and the Cabibbo-Kobayashi-Maskawa matrix element V bc = 0.0411.
For the other resonance of 2 ++ , we can not use the formula in Eq. (12) because we need L = 2 state and hence the matrix element would be different. For the L = 2 case we replace the P cm R with (P cm R ) 5 in Eq. 12, but we do not obtained an absolute value for the width. Yet, we can obtain the ratio of rates for the two 2 ++ resonances.
where G is the two meson loop function for each channel in Eq. (8). The factor C is not the same for J = 0 and J = 2, and we will come back to that. On the other hand, in the other case of resonance coupled to the D * sD * s state for I = 0, J = 2, the hadronization amplitude V ′ had is written as, The scattering amplitudes t D  The differential decay widths dΓ dM inv are given by [42][43][44] where i corresponds to the D * D * or D * sD * s states. In the case of the L = 2 state, we replace P cm by (P cm ) 5 as done before in the integrand of Eq. (12) for the coalescence case. Here, P cm is the momentum of the νe system in the B c rest frame,p ν is the momentum of the neutrino in the νe rest frame as defined in Eq. (14), andp i is the relative momentum of the final mesons in their rest frame, with m i , m ′ i the two meson masses of the final state. The energiesẼ Bc andẼ i are calculated in the νe rest frame, andp Bc is given by Eq. (18).

III. RESULTS
First, we show the results of the coalescence process. We consider X(3934) as the resonance R, namely the B − c → X(3934)νl − process. We evaluate the decay widths Γ coal in Eq. (12) for the resonance L = 0 and J = 0 state using the mass m R = 3943 MeV, and find Γ coal = 2.6 × 10 −14 MeV.
In Fig. 5 we show the integrand of Eq. (12). The mean life of the B c is 0.507 × 10 −12 s, and then the branching ratio is evaluated as For the other resonance of 2 ++ , we can not use the formula of the decay widths in Eq. (12) because we need an L = 2 state and hence the matrix element would be different. Thus, we calculated the ratio of the two L = 2. We consider two options. In the first one the ratio is evaluated as, where (2) By replacing the integral in Eq. (12) by Eq. (32), the other ratio is evaluated as This latter result is more realistic and we take it. We also show the integrand of Eq. (32) in Fig. 6 for the two tensor resonances. Next, we show the results for the rescattering process. In order to study the rescattering process the scattering amplitudes t D * D * ,D * D * , t D *     As important as the shape, showing the presence of a resonance below threshold, the values in the scale, corresponding to ratios, are absolute values of our predictions, tied to the molecular nature of these resonances and their strong coupling to D * D * and D * sD * s . We should note that we have compared the D * D * or D * sD * s production with the production of each particular resonance. This implies that in the experiment the S-wave is separated from the D-wave in each case, something that is at reach in present partial wave analysis of data [30,48] IV. CONCLUSIONS We have studied the semileptonic decay of B c in the reaction B − c →ν e e − X, with X any of the resonances X(3930) (2 ++ ), X(3940) (0 ++ ) and X(4160) (2 ++ ). The main point of the approach is that we treat these resonances as dinamically generated from the vector-vector interaction in the charm sector.
The X(3940) and X(3930) states are basically D * D * molecules in that approach, although they also couple to other channels with a smaller intensity. To produce these X states one proceeds in three steps. The first one looks into the elementary process B − c →ν e e − cc. In the second step the cc pair hadronizes producing an extrā qq with the vacuum quantum numbers, which leads to D * D * and D * sD * s pairs. In the last step these mesons are allowed to undergo final state interaction from where the three resonances appear. By analogy with other reactions producing scalar mesons in the final state, we make an estimation of the rate of B − c →ν e e − X(3940). For the production of the two tensor states we can not obtain the absolute rate of production, but we can obtain the ratio for the X(3930) and X(4160) states. We also look at the production of D * D * and D * sD * s close to threshold and we can make predictions of the ratio of this differential mass distribution to the rate of resonance production, which are tied to the nature of these resonances as dynamically generated from the vector-vector interaction in the charm sector. As more decay modes of B c become available, it would be interesting to look into these modes which will provide good information on the nature of these resonances.