Shedding light on the X ( 3930 ) and X ( 3960 ) states with the B − → K − J /ψω reaction

We have studied the contribution of the state X ( 3930 ) , coming from the interaction of the DD and D + s D − s channels, to the B − → K − J /ψω decay. The purpose of this work is to offer a complementary tool to see if the X ( 3930 ) state observed in the D + D − channel is the same or not as the X ( 3960 ) resonance claimed by the LHCb Collaboration from a peak in the D + s D − s mass distribution around threshold. We present results for what we expect in the J /ψω mass distribution in the B − → K − J /ψω decay and conclude that a clear signal should be seen around 3930 MeV. At the same time, ﬁnding no extra resonance signal at 3960 MeV would be a clear indication that there is not a new state at 3960 MeV, supporting the hypothesis that the near-threshold


Introduction
In Refs.[1,2] the LHCb collaboration analyzed the  + →  +   −   + decay and reported a peak in the  +   −  mass distribution close to the threshold, which was associated to a new resonance,  (3960).Soon after disclosing the experimental results in Ref. [3], it was suggested that this peak would be a necessary consequence of a weakly bound  +   −  state and there was no need to invoke a new resonance if such a state existed [4,5].Actually such a state with scalar quantum numbers is predicted in the lattice QCD (LQCD) calculation of Ref. [6].Indeed, a state of molecular nature coupling to  +   −  and  is found close to  +   −  threshold, with a strong coupling to  +   −  and a weaker one to .At the same time another state coupling to  below its threshold is also found in Ref. [6].Such a  bound state had been previously predicted in Ref. [7][8][9].Although no  +   −  bound state was obtained in Ref. [7], its existence with a mass in the 3915 − 3935 MeV range was predicted in the single-channel calculation carried out in Ref. [9].The coupled-channels scheme of Ref. [7] was reviewed in Ref. [4], where it was shown that a small decrease of the  →  +   −  transition potential reverted into the appearance of a  +   −  bound state close to the  +   −  threshold.It was explained there that this decrease actually was naturally obtained by considering the full propagator ( 2 −  2  * ) −1 in the  * exchange instead of the  2 → 0 approximation used in Ref. [7].The results of Ref. [4] agree then with those of Ref. [6], and a bound state in the  +   −  - coupled-channel amplitude is obtained, coupling more strongly to  +   −  than to , and hence it approximately qualifies as a  +   −  bound state.However, given the fact that this state also couples to , it would be expected to be visible in the  +  − spectrum of some decay.Actually, the LHCb collaboration also finds a scalar state in the  +  − spectrum in the  + →  +  −  + decay [10,11] branded as  0 (3930) with   = 0 ++ and with mass and width given by:  = 3924 ± 2 MeV, Γ = 17 ± 5 MeV.
The  +   −  threshold is at 3936.7 MeV, and hence this state is barely 10 MeV below the  +   −  threshold, so it could as well correspond to the  +   −  bound state found in Refs.[4][5][6]9].This is the suggestion made in Refs.[4,5], where, in addition, it was shown that the existence of this state produces an enhancement above the  +   −  threshold in the mass distribution in the  + →  +   −   + decay, compatible with the experimental results.In Ref. [4] it was also reproduced the experimental ratio of the  →  +  − to  →  +   −  decay widths, of the order of 0.3, assuming  to be the same state responsible for the two peaks.In turn, the authors of Ref. [5] employed an effective field theory based on Heavy Quark Spin Symmetry (HQSS) to conclude that the bump in the  +   −  mass distribution could be described by a bound or virtual state below threshold.Using the pole (either virtual or bound) position and the existing information on the  (3872) and   (3900) resonances, a complete spectrum of the -wave hadronic molecules formed by a pair of ground state charmed and anticharmed mesons was established in that work.In a following paper [12] the available data in  and  +   −  channels from both  decays and  fusion reaction are analyzed, firmly concluding that the  (3930) and  (3960) should correspond to the same state.
The molecular picture to explain the  +   −  bump is also supported using QCD sum rules in Refs.[13,14].The boson exchange model is used in Ref. [15] to obtain also a molecular state with the  +   −  , ,  * +   * −  channels, but a resonance rather than a bound state is found.The molecular picture is also assumed in Ref. [16] where production rates are evaluated.The molecular picture is not the only suggested possibility.In fact, in Ref. [17] and using also QCD sum rules, a scalar diquark-antidiquark state is studied to explain the  +   −  peak.On the other hand, a state of tetraquark nature is also proposed in Ref. [18] using a chromomagnetic interaction model.In Ref. [19], the state is supposed to be the ordinary  0 (2) state of the constituent quark model.In Ref. [20], the  (3915) and  (3960) are considered as four-quark states.
With different views on the subject, it is important to come with new ideas that can help us gain further insight into the origin and meaning of the observed peak and the nature of whatever state is responsible for it.Knowing that: (i) the  (3915) was observed by the Babar collaboration in the  → / reaction [21] with quantum numbers   = 0 ++ [25]; (ii) a supposedly equivalent state was reported by the Belle collaboration at 3943 MeV in the  → / reaction [26]; (iii) the  0 (3930) was found in the  +  − mass distribution of  + →  +  −  + [10,11]; and (iv) the  0 (3960) was found in the  +   −  mass distribution in we propose to look into the / mass distribution of the  + → / + reaction.Actually, data for the  + → / + decay are already available [27].Although promising, it could be improved in future LHCb runs.A peak around 3930 MeV in the / mass distribution in  + → / + is also seen in a BaBar experiment [28], but as in the  +  − mass distribution of the LHCb, the peak could be due to a superposition of a  = 0 and a  = 2 states.
According to the pictures in Refs.[4,12] a peak should be seen around 3930 MeV and no extra peak at 3960 MeV, since we assume that the  +  − and  +   −  peaks are due to the same state  0 (3930).We can make predictions on the position and shape of the peak and the strength compared to those seen in the  +  − and  +   −  mass distributions.Alternatively, and due to a destructive interference with a background, one could also have a dip.Cases like the latter one are seen in hadron physics, for instance in the  cross section in the region of the  0 (980) [29] (see also the general discussion in [30]).
On the theoretical side, Ref. [31] deals with the  + → / + reaction from a different perspective to the one in this manuscript.Indeed, that work looks for signals of  *  * bound states in that reaction.In the 3920 − 3940 MeV region there are predictions of  *  * bound states in Ref. [32], where a state with     = 0 + (0 ++ ) is found around 3940 MeV and another one with 0 + (2 ++ ) around 3920 MeV, both of them coupling to / in -wave.The results of Ref. [31] indicated that a broad peak around 3920 MeV in Ref. [27] could be attributed mostly to the  = 2  *  * state and a cusp like sharp peak of smaller strength appeared at the  *  * threshold.The present work It has been argued [22] that the helicity-2 dominance hypothesis, which is reasonable for the coupling of a 2 ++  c state to two photons [23] and is supported by experimental measurements [24], adopted in the BaBar analysis [25], was not reliable since the  (3915) may not be a purely  c state.If such an assumption was removed, the data appear more consistent with the assignment of 2 ++ to the  (3915) [22].A  *  * bound state with 0 + (0 ++ ) quantum numbers well below the corresponding threshold is possible according also to Ref. [5].A 0 + (2 ++ )  *  * state is also predicted and studied as a HQSS partner of  (3872) in Refs.[8,9,33], but only very close to the  *  * thresold.is also complementary to the one of Ref. [31], looking for the contribution of the  and     states.Given the small mixing of pseudoscalar-pseudoscalar and vector-vector components [34], the contributions that we find here are additional to those found in Ref. [31].

Formalism
In order to work with  quarks rather than b we study the charge conjugate reaction,  − → /   − .We can see in Fig. 1 the mechanism at the quark level for the direct production of the final state without the need of rescattering (tree level).The hadronization is achieved including a q pair with the quantum numbers of the vacuum.Concerning flavor space, we denote q ≡  q   ( = {, , }).Then, since we want to get a pseudoscalar meson and a vector, we obtain where , are the  q matrices in SU(4) flavor space written in terms of pseudoscalar or vector mesons: In the matrix  we have considered the standard  −  mixing of Ref. [35].We could as well have the ( ) 31 combination, but since we want a final state with a  − , the  product is the appropriate one.Thus, We obtain then the /   − component for the tree level background in the  − decay.
Next, we look for the resonant contribution.The ,     components are obtained by hadronizing the  c pair in Fig. 1(a) as shown in Fig. 2. We have then: Hadronization of the  c quark pair.We remark that the     component does not play any role due to its large mass.The non-resonant (tree-level) and resonant parts of the reaction  − → /   − are produced according to the diagrams shown in Fig. 3. Our isospin phase convention of the  and  doublets is the following:  ≡ ( + , − 0 ) and  ≡ ( 0 ,  − ).Then, it follows that the  = 0 state is given by We associate the weights  and  to the mechanisms depicted in Figs.1(b) and 2, respectively.Since both of them come from internal emission and hadronization, we will assume that  and  are of the same order of magnitude, but the relative sign is not fixed.In addition, we consider the process to proceed via -wave, which is allowed for  − → /   − and  − →    − /  +   −   − decays.Therefore, the production amplitude is given by where the ì   ( = {/, }) are the polarization vector of the vector mesons,   ( inv ) and   +   −  ( inv ) are respectively the loop functions of the  and  +   −  intermediate states, depicted in Fig. 3(b), and  inv is the invariant mass of the  or  +   −  system.In addition,  ,  represent the elements of the unitarized transition matrix between ,  +   −  states, obtained in Ref. [4] from J/ψ(p 1 ) where  here represents the interaction potential matrix.We use the same  and  matrices as in Ref. [4].It is worth mentioning that even if throughout the manuscript we are apparently using a SU(4) formalism, one is actually only making use of SU(3) symmetry and the  q character of the meson states [36].
On the other hand, the functions   (  =  +  − ,  0  0 ,  +   −  ) correspond to the triangle loops in Fig. 3(b), with  +  − ,  0  0 , and  +   −  being the intermediate states.We shall see that the   also contains the ì   / • ì   factor.Before the evaluation of their magnitudes, it is easy to see that the  +   −  channel does not contribute.Indeed, it would necessarily involve a  +   −   vertex which is zero in the scheme of Ref. [4], since the  meson has only light-flavour quarks and has no overlap with the charm and strange content of the  +   −  pair.As a consequence, t/  in Eq. ( 6) gets simplified and can be rewritten as In order to calculate the triangle loop function, we show in Fig. 4 the diagram with the explicit momenta.The analytical expression for diagram (a) reads: where  1 , 2 correspond to / +  − or  +  − vertices, which are evaluated through the following Lagrangian: with , denoting the matrices of Eqs.[37,38].Here, we take  max = 750 MeV as in Ref. [4].The factor (  * /  * ) 2 in Eq. ( 9) is a normalization stemming from the different weight factors 1/

√
2 of the light and heavy mesons [39].On the other hand, it can also be seen that the contributions of diagrams in Fig. 4 [(a) and (b)] are the same, and that   +  − =   0  0 .Hence, we have with One can take advantage of the fact that the -meson exchanged between the / and  in Fig. 4 is far off-shell and we can factor it out of the integral.Then, we perform the  0 integration analytically in Eq. ( 12) via Cauchy's theorem, taking the pole of the  line, and thus  is placed on-shell.Therefore, where we have considered that / and the intermediate  on-shell come with relative small momentum.As a result, we obtain where  0 =  inv and   () ≡ In Appendix A we discuss the integral appearing in Eq. ( 14).So, by using we can now write The mass distribution is obtained from the expression: where the factor 3 comes from summing over the / and  polarizations, and 3 Partial decay width of the  0 (3930) to /  In Ref. [4] the  0 (3930) state came from the coupled channels  and  +   −  , mostly  +   −  .The width originating from its decay to the open  channel was found to be 13 MeV.Here we benefit from this previous evaluation to obtain the decay width of  0 (3930) to / .

We can take the diagrams of
where  has been defined in Eqs. ( 14) and ( 15), and  , is the coupling of the resonance to  found in Ref. [4],  , = −1.00− 2.71GeV .The  0 (3930) width is then given by where the factor 3 comes again from the sum of /,  polarizations.For the sake of coherence, we use the mass   = 3932.72MeV from Ref. [4].

Further considerations
As mentioned above, in Ref. [31] the same reaction  − → /   − was studied but looking for states of  *  * nature.The same quark topologies as in Figs. 1 and 2 were considered but, in addition, a diagram corresponding to external emission, shown here in Fig. 5, was also taken into account.The hadronization of the  c pair gives rise to: Among all these terms, there is one that can lead to our resonance, namely  −  0 .Together with the  0 meson of the  ū pair in Fig. 5 we obtain the term Considering Eqs. ( 5) and (6) thus will lead to an extra term,

𝑉
Hence, the term mentioned above can be incorporated by performing the following replacement in Eq. ( 6), Since the external emission topology is color favored, we can expect that  3 .As a consequence, the new mechanism can give about the same contribution as the former intermediate  state contribution.Yet, since the resonance that we consider couples mostly to  +   −  and is close to the  +   −  threshold, both the   +   −

𝑠 and 𝑇 𝐷 +
−  , have larger strength than   and  , , and the new  contribution, as well as the former one, will be small compared to the contribution of the intermediate  +   −  state.In Ref. [31] the largest contribution was obtained from a  *  * (  = 2 ++ ) state.On experimental grounds, in order to isolate the  ,  +   −  state, the 2 ++ contribution can be identified either by a partial wave analysis or by using the moments method [40][41][42][43], projecting the mass distribution as 3.88 3.90 3.92 3.94 3.96 3.98 4.00 Real and imaginary parts of triangle loop contribution  given in Eqs. ( 14) and ( 15).
where  is the angle between the  − and  in the /  rest frame, which is the same as the angle between  in the /  rest frame and  − in the  − rest frame; Ω is the solid angle of /  in their rest frame.

Results
In Ref. [4] the dimensional regularization was used, and accordingly two subtraction constants were introduced: one in each of the  and  +   −  channels.Notwithstanding this, one can find the equivalent cutoff values by demanding that the loop  functions are the same at their respective thresholds.In this way, we obtain  max = 550 MeV for  and 1070 MeV for  +   −  .Since we are concerned about the state coupling mostly to  +   −  , the use of a different  max for the  channel does not change much the results for that state.We need max for the triangle loop function and we shall take an intermediate value, i.e.  max = 750 MeV.We shall quantify the uncertainties concerning small changes of that parameter.
Let us start the discussion of the results by evaluating the triangle contribution  .As we can see in Eqs. ( 14) and ( 15),  looks like the ordinary meson-meson  function, except for the extra ì  2 factor in the numerator (see also Appendix A).Since the threshold of  is at 3730 MeV, in the region of 3870 MeV to 4000 MeV it should behave like a  function above threshold where the real part increases from a negative number and the imaginary part is negative and grows in magnitude when the energy increases.This is what we observe in Fig. 6.The magnitude of  around 3930 MeV is of the order of (2 -3) × 10 −1 .
For completeness, in Fig. 7 we plot the squared modulus of the elements of the unitarized transition matrix  that are relevant in the estimation of the t/  in Eq. ( 8).As we discussed previously, in the amplitude of Eq. ( 8) we needed the  , and   +   −  , transitions, and in Fig. 7 it can be clearly observed that We plot in Fig. 8 the mass distribution Γ/ inv for / production in the reaction  − →  − / .We consider different sets of weights {,  ,  }, present in t/ , based on our estimates of the strength of internal and external emissions.We take  and  of the order of 1 an  of the order of 3, but play with possible relative signs of the magnitudes.Further experimental works will give us information on which set should be better to reproduce experimental data.However, these estimations suggest interesting findings.In the region of 3920 − 3950 MeV an interference appears among the tree-level and resonant terms present in the amplitude t/  , which might be   constructive or destructive, depending on the set of parameters utilized.One should recall that interferences leading to dips in amplitudes or mass distributions are common when dealing with coupled channels [30].At any rate, a clear signal of the resonance is found.Obviously, our calculations do not show any resonant signal around 3960 MeV, since a state at that energy has not been included in our amplitudes.The line shapes shown here could therefore serve as a baseline of future theoretical or experimental analysis to detect or discard the presence of an eventual additional state around 3960 MeV.In Fig. 9 we compare the calculated mass distribution Γ/ inv (/) taking the normalization ( = 1,  = −1,  = 3), with those studied in Ref. [4], considering  inv ( +  − ),  inv ( +   −  ) for the reactions  + →  +  +  − and  + →  +  +   −  .The scattering amplitudes for these two reactions correspond to the mechanism of Fig. 3 removing the triangle loop and is then given by (see also [4]) From Fig. 9 we conclude that the expected signal is of the same order of that observed for  − →  −  +  − and is reachable with the statistics of present runs and should be even clearer in future runs of LHCb.
Finally, in Fig. 10 we plot the width Γ  0 (3930) from the /  decay channel given by Eq. ( 22), as a function of the /  invariant mass implicit in the  function, defined in the coupling   0 (3930)  /  .Taking the mass of the resonance considered here, one can conclude that the contribution to the width Γ  0 (3930) from the / channel is around 1.7 MeV.This value is about one order of magnitude smaller than the width coming from the decay to .Even if the coupling of the resonance to  is small compared to that to  +   −  , the width for the decay into the  channel is still larger than the one into / .Then it might be surprising to see that the signals in Fig. 8 are of comparable order.This is due to the fact that the signal for /  comes from the interference with a large tree level, and hence is linear, not quadratic in the resonance amplitude.
We have checked uncertainties coming from changes in  max , by evaluating the results with the use of  max = 700 MeV and  max = 800 MeV.The changes in the results of Fig. 8 are small, at the level of 7% and do not change qualitatively what is found there.The modifications in the  0 (3930) width for / are larger, because it is quadratic in the resonant amplitude.On the other hand, the results for Γ/ inv in Fig. 8 are linearly dependent on the resonance signal.Numerically we find   Hence, the changes are around 15%.Yet, we still find that this width is reasonably smaller than the one for  0 (3930) decay to  .
The data in Refs.[27,28] are indicative and promising.A clear peak is seen in the / mass distribution in the  + →  + / decay.It is unclear, however, how much of the peak can be due to a / in -wave with  = 0 or in -wave with  = 2.This is actually the case in the  +  − mass distribution in the LHCb experiment [10,11], where the peak in this region is found to come from the  0 (3930) and  2 (3930).Disentangling the content of the peak along the lines used in the LHCb analysis [10,11] is necessary to allow a comparison with the predictions in the present work.In any case, it is already illustrative to see that in this experiment one cannot find a trace of a possible peak around 3960 MeV.This feature is even clearer in the  *  * → / spectrum of Ref. [25].These experiments, so far, do not support the hypothesis that there is an extra state around 3960 MeV.

Conclusions
We have studied the  − →  − / decay, paying attention to the formation of the state  0 (3930), which in our approach comes from the interaction of the coupled channels  and  +   −  .Prior to the present work we had looked at the  +  − and  +   −  mass distributions in the  − →  −  +  − and  − →  −  +   −  reactions, reaching the conclusion that the  (3960) state, claimed from the peak around the  +   −  threshold in the  − →  −  +   −  reaction was actually the same state  0 (3930) observed in the  +  − mass distribution in the  − →  −  +  − decay.Indeed, the existence of a resonance below the threshold of some channel induces an enhanced mass distribution around the threshold of that channel if the state couples to it.Such a state below the  +   −  threshold is found in Refs.[4,5,12] and also in the lattice results of Ref. [6].The purpose of the present work is to offer additional information to corroborate, or refute, this reasonable hypothesis.For this purpose we have used the same formalism and input as in Ref. [4] and implemented the transition from the  and  +   −  components to / via a triangle loop diagram with  exchange.In addition, we have also considered the unavoidable tree level contribution to  − →  − /, which we have added to the resonant mechanisms.
We have some unknown magnitudes in the approach, but, up to a normalization not relevant in the study, we control the order of magnitude of these magnitudes and we study different possibilities within the freedom that we have.We find important consequences of our study.The first one is that the tree level and the resonant terms add coherently and we find clear interference of the amplitudes, constructive some times and destructive other times, but the signal of the resonance is clearly seen either way.This should serve as a warning for experimental analyses, where commonly the non resonant part, background, is added incoherently to the resonant contributions.The other relevant consequence is that, since the threshold of / is at 3880 MeV, far below the  +   −  threshold, then the decay to this channel is open for 3930 MeV and 3960 MeV.With our hypothesis that there is only one state at 3930 MeV, we obviously observe only the resonant signal around this energy, but the shape of this signal, whether constructive or destructive, can be approximately controlled and, furthermore, is relatively narrow.Should there be an extra state at around 3960 MeV or, more generally, above the  +   −  threshold, we expect to find an additional peak at that energy.Hence, the performance of the experiment will tell us whether there is one or two states.
So far there are data on this reaction from the Thesis work of Guido Andreassi [27] in 2014.However, the data are not yet published.The statistics of the data do not allow us to make conclusions at present, but, with increased number of data in present and future runs of LHCb, the spectra of this reaction should be precise enough to provide an answer to the questions raised in the present work.
As for the BaBar analysis of Ref. [28], where a peak is clearly seen in the / mass distribution of the  + →  + / decay, as well as for the work of Ref. [27], it is also important to perform a separation of the contribution of  = 0 and  = 2, once it has been shown that the peak at this energy seen in the  +  − mass distribution in the LHCb experiment in  + →  +  +  − decay comes from the contribution of the  0 (3930) and the  2 (3930).The techniques to do this separation are, thus, available and their use in further analysis of the  + →  + / reaction is most welcome.Moreover, from the data of these works, and also from  *  * → / [25], one cannot find any signature of the existence of a state at 3960 MeV.
We subtract and add the on-shell momentum in the numerator, and we find: being Ω Λ just a number.Taking a sharp cutoff regulator, as pointed out above, we have indeed: In addition, an explicit, algebraic expression for  Λ () can be found in (the erratum to) Ref. [48].Physically, this reduction of a three-to a two-point loop function means that, because the exchanged -meson is far off-shell and its propagator is factored out, we have reduced the three-point loop integral with a -meson exchange to a two-point loop function and a contact interaction for the amplitude  →  .

Fig. 1
Fig. 1 Mechanism at the quark level for the reaction  − →  /   − .(a) Primary quark production with internal emission.(b) Hadronization of the  ū component into two mesons.

Fig. 4
Fig. 4 Triangle diagram with  +  − intermediate state.The panels (a) and (a) are essentially the same but with  / ad  exchanged.
(2)-(3) respectively, • • • the trace over the flavor space, and  the coupling  =   /(2   ) (  = 800 MeV,   = 93 MeV).The factor  ( max − | ì |) appears due to the sharp cutoff regularization used in the    matrix, which takes the form Fig 4 and together with the  0  0 intermediate diagram we find an effective vertex of the resonance and /  as

Fig. 7
Fig. 7 Squared modulus of the transition matrix elements   D, D and   +   −  , D appearing in t/  in Eq. (8).

Fig. 8
Fig. 8 Differential decay for the reaction  + →  +  / , as a function of the  /  mass distribution, denoted as √ , with different values of the weights ,  ,  appearing in t/ .

Fig. 10
Fig. 10 Contribution to the width Γ  0 (3930) from the  /  channel, as a function of the  /  invariant mass.The vertical lines delimit the region of the resonance.