Disclosing $D^*\bar{D}^*$ molecular states in the $B_c^- \to \pi^- J/\psi \omega$ decay

We study the $B_c^- \to \pi^- J/\psi \omega$ and $B_c^- \to \pi^- D^* \bar{D}^*$ reactions and show that they are related by the presence of two resonances, the $X(3940)$ and $X(3930)$, that are of molecular nature and couple most strongly to $D^* \bar{D}^*$, but also to $J/\psi\omega$. Because of that, in the $J/\psi\omega$ mass distribution we find a cusp with large strength at the $D^* \bar{D}^*$ threshold and predict the ratio of strengths between the peak of the cusp and the maximum of the $D^* \bar{D}^*$ distribution close to $D^* \bar{D}^*$ threshold, which are distinct features of the molecular nature of these two resonances.


I. INTRODUCTION
Molecular states of mesons have long been the subject of study in hadron physics. Detailed recent reviews can be seen in Refs. [1,2]. As commented in Ref. [3] the support for hadron molecules is quite obvious once we realize that baryon molecules exist in the form of nuclei. In fact, multi-mesons states, not just meson-meson molecules, have also been advocated, like multi-rho states in Ref. [4], K * -multi-rho states in Ref. [5], D * -multi-rho states in Ref. [6], two mesons and a baryon states [7,8] and many others (see a recent review in Ref. [9]). Actually, the interaction between mesons, particularly vector mesons in spin two, is very strong [10][11][12], even stronger than between nucleons, and the only limit to the formation of multi-meson states is that we do not have the meson number conservation, unlike baryon number conservation for the nucleons forming nuclei. This allows the multi-meson states to decay in states of fewer, or lighter mesons, the width increases with the number of mesons of the cluster, and at some point they are no longer identifiable experimentally.
The identification of states as being of molecular nature is not an easy task, and in general standard quark structures, or multiquark states are competing in the interpretation [1,3]. Yet, there are several experimental features that reveal the molecular structure [2] and ultimately it is the systematic and correct description of experimental features and the accuracy of the predictions what builds up in favor of this structure for many states.
The weak decay of heavy mesons and baryons has turned out into one important tool to identify states of molecular type [13]. Curiously, an interaction that does not respect parity and isospin, has shown itself as a great tool to identify molecular states because certain decays filter good quantum numbers due to selection rules, like Cabibbo and color enhancement in some topologies of decay modes.
One of the features attached to the molecular states that couple to several hadron-hadron channels, is that by looking at one of the channels with relatively small strength one finds a strong and unexpected cusp in the threshold of the channels corresponding to the main component of the molecule. One recent example of this was found in the B + → J/ψφK + reaction measured at LHCb [14,15]. The reaction was analyzed in [14,15] and at low invariant masses only the X(4140) state was included, concluding that its width had to be considerably larger than the average of the PDG [16] from other experiments. A different interpretation, with a better fit to the data, was given in [17], where, in addition to the X(4140), the X(4160) was included in the fit, assuming that this state is the D * sD * s state predicted in [18] as a 0 + [2 ++ ] state. It is worth noting that other works have also suggested a bound state of D * sD * s [19][20][21][22], although it was originally associated to the X(4140). This bound D * sD * s state also couples to other light vector states and to J/ψφ, hence, it can be observed in this latter channel. However, the fact that the resonance couples most strongly to D * sD * s has the consequence that the J/ψφ mass spectrum develops a strong cusp at the D * sD * s threshold, something visible in the experimental spectra with an increased strength in that region. It is also worth mentioning that a similar enhancement is seen, although with poor statistics, in the recent BESIII work on the e + e − → γJ/ψφ reaction [23].
In the present work we want to continue along this line of research and present results for a reaction that should reveal the D * D * nature of two states found in [18] as 0 + [0 ++ ] and 0 + [2 ++ ] at 3943 MeV and 3922 MeV, respectively, which can be identified with some experimental states in that region [16,18] 1 . In this case we note that the states found, mostly D * D * bound states, also couple to J/ψω in the second place, and J/ψφ with smaller strength. So we choose the J/ψω observation channel looking into the necessary cusp that should develop at the D * D * threshold. For this purpose we look into the B − c → π − J/ψω decay and then into the J/ψω invariant mass distribution. The choice of this reaction is that in a first stage of the reaction the D * D * state is formed with a dominant weak decay mechanism, but the J/ψω state is not formed at this level. Then the J/ψω is finally produced via rescattering of the D * D * component with the other components that make up the two molecular states. This stresses the role of the resonance since there is no tree level J/ψω contribution. Thus, we obtain two peaks in the J/ψω mass distribution corresponding to the molecular states and a strong cusp at the D * D * threshold. In addition we also look at the D * D * mass distribution in the B − c → π − D * D * reaction and evaluate its strength above the D * D * threshold, which is closely connected to the strength of the J/ψω mass distribution.
The D * D * cusp feature, together with the relative strength of the D * D * compared to the one of J/ψω, are two magnitudes which are tied to the molecular structure of these two resonances, and we encourage the performance of the experiment that should bring valuable 1 The state at 3943 MeV can be associated with the X(3940) of [24,25] and the X(3922) with the X(3930) [26] now classified in the PDG as the χ c2 (2P ). light into these issues.

II. FORMALISM
We look into the B − c decay mechanism at quark level depicted in Fig. 1(a). The mechanism qualifies as external emission [27] and is both Cabibbo favored in W ud vertex, and color favored (the W bc vertex is also the least Cabibbo suppressed of the b decays). Then the final c quark from b decay and the spectatorc quark from the B − c hadronize, with the incorporation ofqq with the quantum numbers of the vacuum (see Fig. 1(b)) to give two mesons. The resulting two mesons are easily obtained by writing where M ij is the qq matrix with the u, d, s, c quarks. However, it is convenient to write the qq matrix in terms of physical mesons, in this case vector mesons as and we get Tree level contribution corresponding to the hadronization depicted in Fig. 1(b).
The intrinsic phase convention for isospin multiplets in (D * + , −D * 0 ), (D * 0 , D * − ) indicates that the isospin combination of H is I = 0, as it should be since it comes from cc. Thus, we can write where we have neglected the J/ψJ/ψ component which is far beyond in energy from our range of concern. In addition, the coupling of the resonances found in [18] to J/ψJ/ψ is negligibly small.
The combination of |H in Eq.
with , the polarization vertices of D * ,D * . Note that we shall work in the rest frame of the resonances produced, where D * ,D * momenta are small with respect to their masses and then we neglect the 0 component. On the other hand, if we produce a 2 ++ state, the 0 − → 0 − 2 + requires L = 2 and we shall then take the D-wave structure where k is the momentum of the pion. Hence, the tree level amplitude for B − s → π − D * D * shown in Fig. 2 is given by where we have substituted A of Eq. (4) by A | k av | 2 , with k av , an average value of k, just to make A and B have the same dimension. We take | k av | = 1000 MeV.
After the first step for D * D * and D * sD * s production, these mesons undergo final state interaction, as depicted in Fig. 3 and 4, to produce J/ψω and D * D * in the final state. In the case of J/ψω production shown in Fig. 3, since this state is not primarily produced in |H , it is produced through rescattering via the resonances X(3922) and X(3943). In the case of D * D * production, shown in Fig. 4, it proceeds via tree level (primary production, Fig. 4(a)) and rescattering ( Fig. 4(b)).
Analytically, we have where and while for D * D * production we have with and where I, II stand for the 0 ++ and 2 ++ states, respectively. Since the · and · k · k − where k is the π momentum in the B − c rest frame,p ω the ω momentum in the J/ψω rest frame and k the pion momentum in the J/ψω rest frame for the J/ψω final state, We get the amplitudes t I and t II from Ref. [18]. We take them using the Flatté form of the amplitude in terms of the couplings obtained in Ref. [18] and the width. The couplings are given in Table I.
The amplitudes are given by with i = I, II, and j = J/ψω or D * D * . We also have where the width is taken as with andp ω given by Eq. (17) withp D * asp ω in Eq. (17) Note that in Eq. (22), Γ (i) D * D * only starts above the D * D * threshold, but since the coupling of the resonance to this channel is so large, it grows fast above threshold giving rise to the Flatté effect.

III. RESULTS
We will present the invariant mass distribution in arbitrary units, but dΓ In Fig. 5 we see that due to the proximity of the two resonances, and the fact that both of them can be produced in this reaction, the two peaks actually merge into a broader  one, although a precise measurement could maybe allow a separation of the two peaks, particularly if a partial wave analysis is done that separates the two different spin resonances.
Interesting, however, is the fact that the cusp appears always at the same place, the D * D * threshold. The other relevant feature is that its strength grows with increasing weight of the tensor resonance, indicating that the cusp is basically tied to the 2 ++ X(3930) state.
In Fig. 6 we show the D * D * mass distribution in the B − c → π − D * D * decay. We observe a distribution quite different from ordinary phase space, sticking close to threshold, indicating that it is influenced by a resonance below threshold. Its strength also grows with increasing strength of the tensor resonance, which establishes a link between this state and the D * D * distribution.
Very interesting is the ratio of the strengths of dΓ dM J/ψω inv at the peak of the cusp of the D * D * threshold versus the strength at the peak of dΓ dM D * D * inv . We show these numbers in Table II for different values of B. As we can see, this ratio is relatively stable and tied to the dynamically generated nature of the two resonances discussed.  The fact that the ratio R is essentially independent on the strength B of the tensor resonance indicates that it is this resonance in practice the one that is responsible for both the cusp in the J/ψω and the D * D * mass distributions in the B − c → π − D * D * reaction.

IV. CONCLUSIONS
We have looked into the B − c → J/ψω decay and in particular in the J/ψω mass distribution. We find that this observable is much influenced by the role of the X(3940) and X(3930) resonances, which in Ref. [18] appear dynamically generated from the vector-vector meson interaction in the charm sector. These resonances couple mostly to D * D * in 2 ++ and 0 ++ , respectively. In order to find support for this nature of the resonances we stress two particular features: the first one is to observe that J/ψω is not the main channel for this resonances, but D * D * . As a consequence, one finds a strong cusp at the D * D * threshold in the J/ψω mass distribution. The other feature is that since the resonances are tied to D * D * , they should influence the D * D * mass distribution close to threshold in the B − c → π − D * D * reaction. What we find is that, within uncertainties tied to our ignorance of the weight by which the X(3940) and X(3930) resonances are produced, the ratio of the strength at the cusp peak and the strength at the maximum of the D * D * mass distribution are related and quite independent of the relative weight of these two resonances. This is because the D * D * mass distribution is more influenced by the X(3930) resonance that has a larger width.
In addition we observe also a peak around 3930−3940 MeV in the J/ψω mass distribution, corresponding to the excitation of these two resonances, and show that the cusp at the D * D * threshold has similar strength as the peak. All these features, when observed, should serve to support the molecular nature of these resonances and we can only encourage the performance of the experiments.