Disclosing D∗D¯∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^*\bar{D}^*$$\end{document} molecular states in the Bc-→π-J/ψω\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_c^- \rightarrow \pi ^- J/\psi \omega $$\end{document} decay

We study the Bc-→π-J/ψω\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_c^- \rightarrow \pi ^- J/\psi \omega $$\end{document} and Bc-→π-D∗D¯∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_c^- \rightarrow \pi ^- D^* \bar{D}^*$$\end{document} 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∗D¯∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^* \bar{D}^*$$\end{document}, but also to J/ψω\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J/\psi \omega $$\end{document}. Because of that, in the J/ψω\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J/\psi \omega $$\end{document} mass distribution we find a cusp with large strength at the D∗D¯∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^* \bar{D}^*$$\end{document} threshold and predict the ratio of strengths between the peak of the cusp and the maximum of the D∗D¯∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^* \bar{D}^*$$\end{document} distribution close to D∗D¯∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^* \bar{D}^*$$\end{document} threshold, which are distinct features of the molecular nature of these two resonances.


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, multimesons states, not just meson-meson molecules, have also been advocated, like multi-rho states in Ref. [4], K * -multirho 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. Even then, according to [4][5][6], states up to 6 a e-mail: dailr@lnnu.edu.cn b e-mail: jdias@if.usp.br c e-mail: oset@ific.uv.es vector mesons can be detected and the f 6 (2510) qualifies as a six-rho meson state [4].
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 *  [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 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 light into these issues. 1 The state at 3943 MeV can be associated with the X (3940) of [24,25] and the X (3922) with the X (3930) [27] now classified in the PDG as the χ c2 (2P). The width of the X (3943) MeV obtained in [18] of 15-17 MeV looks small compared to = 87 ± 22 ± 26 MeV of [24]. These results are improved in [25] by measuring the DD and DD * decay channels with = 15.4 ± 10.1 MeV, in good agreement with the results of [18]. A further experiment for the X (3940) in [26] from the measurement of the D * D * channel close to threshold, provides = 37 +26 −15 ± 8 MeV, which is also compatible with the results of [18] within the large errors. The width of about 50 MeV of [18] for the X (3922) is also in fair agreement with = 29 ± 10 ± 2 MeV of [27].cc

Formalism
We look into the B − c decay mechanism at quark level depicted in Fig. 1a.
The mechanism qualifies as external emission [28] 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 spectator c quark from the B − c hadronize, with the incorporation ofqq with the quantum numbers of the vacuum (see Fig. 1b) to give two mesons. The resulting two mesons are easily obtained by writing where M i j 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 Tree level contribution corresponding to the hadronization depicted in Fig. 1b and we get The intrinsic phase convention for isospin multiplets in 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. (3) accounts only for the flavor composition. We need to take into account the spinangular momentum structure of the vertices. If we produce a 0 + [0 ++ ] D * D * state we have 0 − → 0 − 0 + transition and we adopt the common choice of taking the lowest possible angular momentum in the vertex, L = 0. The s-wave and the J P = 1 − of the D * leads us to a vertex of the type 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 Figs. 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. 4a) and rescattering (Fig.  4b).
Analytically, we have where and while for D * D * production we have with and where I , I I stand for the 0 ++ and 2 ++ states, respectively. Since the · and · k · k − 1 3 |k| 2 · structures filter spin 0 and 2 respectively, the structure is kept in the iterations implicit in Eqs.  [18] using dimensional regularization with the subtraction constant a = −2.07 and μ = 1000 MeV. Here, we follow the prescription of Refs. [17,29] and we use the cutoff method with q max fixed to reproduce the results of Ref. [18]. In the former equations A and B are functions (we take them as constants in the limited range of invariant mass studied) which have to do with the weight of the weak process and hadronization before the final state interaction is taken into account. We shall vary A and B within a reasonable range to see the results.
With the amplitudes of Eqs. (7) and (10) the mass distributions, summing |t| 2 over the final vector polarizations, given by 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, with λ(x, y, z) = We get the amplitudes t I and t I I 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 1.
The amplitudes are given by with i = I, I I , and j = J/ψω or D * D * . We also have where the width is taken as andp ω given by Eq. (17) withp D * asp ω in Eq.
Note that in Eq. (22), 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.
In [18] only the vector-vector channels were used in the coupled channel calculation. The D * D * channel can also couple to DD via pion exchange. This contribution was taken into account in [18], including box diagrams with DD in the intermediate states. As shown there, the inclusion of these channels barely changes the mass of the states and modifies the width in only a few MeV. We thus ignore it here. A possible η c η decay channel is even more suppressed because if involves the D exchange, rather than the pion, and the heavy mass reduces the strength of the propagator.

Results
We will present the invariant mass distribution in arbitrary units, but ization. For this purpose, we take A = 1 and look at the results for different values of B. Since A and B have been normalized to have the same dimensions, providing similar strength for the two terms for A = B, we will take values of B close to 0.5, 1, 1.5 and 2. The values chosen for this ratio are not arbitrary. Using the 3 P 0 model to implement the qq hadronization [30] and Racah algebra one can relate the states for scalar and vector or tensor production. This exercise is done in [31] to relate pseudoscalar and vector pro- 2625) decays and the difference results in a factor of √ 3. Furthermore, we have empirical information more closely related to the problem discussed here in the D + → π + f 0 (1370) and  (1270). Indeed, the f 0 (1370) and f 2 (1270) are obtained from the interaction of ρρ in [10][11][12]. They are the equivalent in the light sector to the X (3940) and X (3930) studied here in the heavy quark sector. The branching ratios are (8 ± 4) × 10 −5 and (4.8 ± 0.8) × 10 −4 , respectively, with the resonances decaying to π + π − [16]. Correcting by 3 2 to account for π 0 π 0 decay of the I = 0 state, and dividing by the branching fractions to ππ from the PDG [16], we obtain the rates (4.6 ± 2.3) × 10 −4 and (8.5 ± 1.4) × 10 −4 for the production of the two resonances irrespectively of their decay mode. As we can see these rates differ by about a factor of two. This also indicates a preference for value of B/A bigger than 1 in our case. We show in Fig. 5  ues. The absolute normalization is arbitrary and the shape changes a bit since one gives more strength to one or another resonance changing B.
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 We show these numbers in Table 2 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. It is interesting to further investigate the role of the resonance in the D * D * cusp. Certainly, the fact that a cusp in D * D * appears observing the J/ψω invariant mass indicates that these channels are coupled, otherwise there would be no cusp. Even more, for the particular case, since J/ψω does not appear in the contribution of primary formed states of Eq. (2), the presence of a cusp at the D * D * threshold implies a transition from D * D * or D * sD s * to J/ψω. This information is already helpful by itself, indicating the presence of coupled channels. On the other hand, we can simulate a situation where we have coupled channels but we have no resonances. Without changing the formalism we can artificially place the positions of the resonances, at 3800 MeV, for instance, and see what happens. We show the results in Fig. 7 and we can see that the resonances do not show up in the range of energies of Fig. 5, which we also consider in Fig. 7. We still see a cusp at the D * D * threshold, which is a necessary consequence of the coupled channels, with or without resonances. However, for the same strength of B/A, the strength of the cusp is much weaker than in Fig. 5, when we have the presence of the resonances nearby. Also the shape of the cusp is sharper in the presence of the resonances, as seen in Fig. 5.

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 relative 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.