Molecular Ξbc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Xi _{bc}$$\end{document} states from meson–baryon interaction

We have studied the meson–baryon interaction in coupled channels with the same quantum numbers of Ξbc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Xi _{bc}$$\end{document}. The interaction is attractive in some channels and of sufficient intensity to lead to bound states or resonances. We use a model describing the meson–baryon interaction based on an extrapolation of the local hidden gauge approach to the heavy sector, which has been successfully used in predicting Ωc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Omega _c$$\end{document} and hidden charm states. We obtain many states, some of them narrow or with zero width, as a consequence of the interaction, which qualify as molecular states in those channels. The success in related sectors of the picture used should encourage the experimental search for such states.


Introduction
The spectroscopy of baryons with heavy quarks has raised a wave of intensive theoretical work, with models competing to explain experimental facts and make new predictions. The reporting of several pentaquark states in Ref. [1], updated recently [2], was a main trigger of this wave of works, but other works, as the discovery of several narrow c states [3], and the more recent finding of a ++ cc state [4], have also contributed to keep the flame alive.
Related works include the study of systems of a light pseudoscalar with doubly heavy baryons [56][57][58], studies of triple charm molecular pentaquarks of cc D ( * ) systems with pion exchange [59], the use of Lorentz-invariant baryon chiral perturbation theory to study the ground state of cc [60], the study of electromagnetic form factors of cc , cc [61], and the study of cc , bc and bc masses using a scalar confining potential and one gluon exchange with the Bethe-Salpeter equation [62]. Reviews on these topics can be found in Refs. [63][64][65]. Related to the works of Refs. [56][57][58] is the work of Ref. [66], but with more coupled channels. For instance, in addition to cc π , cc η, cc K , that account for a light pseudoscalar and a heavy baryon, the c D, c D, c D s , c D s channels are considered to produce cc states in J P = 1/2 − , which by themselves give rise to molecular states. Similar states using vector-baryon interaction and mesons with 3/2 + baryons are also considered in Ref. [66]. The approach predicts several states of negative parity between 3837 MeV and 4374 MeV. It is clear that when including coupled channels with charmed mesons one can no longer invoke chiral dynamics, as is the case when dealing with light pseudoscalar in Refs. [56][57][58]. Instead, a method was found in Refs. [67,68] to produce a reliable source of interaction in this case: (i) First one realizes that the chiral Lagrangians in SU (3) can be obtained from the local hidden gauge approach [69][70][71][72] by exchanging vector mesons. This was shown in the case of the pseudoscalar-pseudoscalar interaction in Ref. [73]. (ii) Take a typical channel c D and the direct transition c D → c D. The D + flavor wave function is just cd, and for the c and other baryons we single out the heavy quark and impose flavor-spin symmetry in the remaining light quarks. Thus, explicitly one is not making use of SU(4) symmetry. The direct c D → c D transition is mediated by the exchange of light vectors and c quarks in D and c act as spectators. Then the interaction follows the SU(3) symmetry of the light quarks. A welcome side effect is that, since the heavy quarks are spectators, the interaction does not depend upon them and heavy quark symmetry is automatically implemented. (iii) Some non-diagonal transitions, like c D → cc π , require the exchange of a D * and here the heavy quarks are no longer spectators. Yet, no SU(4) is used in the approach with the wave functions used, and the vertices V P P, V B B (V for vector, P for pseudoscalar and B for baryon) essentially count the number of quarks involved in the exchange. Yet, these terms no longer comply with heavy quark symmetry, as one finds explicitly, but neither should they, since the interaction goes as O( 1 (m Q for the mass of the heavy quark) from the D * propagator, and these terms are subleading in the O( 1 m Q ) counting, and small in practice. (iv) One needs a piece of experimental information to fine tune the regulator of the loops (usually q max for the modulus of the three momentum), which is adjusted to some mass and should be of natural size in the range of 600-800 MeV. Then the masses and widths of many states are predicted by the approach.
This said it is not surprising that the approaches which use explicitly SU(4) for the evaluation of this interaction obtain the same results for the transition matrix elements led by the exchange of light vectors, since automatically they are effectively using the SU(3) subgroup of this group. This is the case of the work of Ref. [74] in the study of c molecular states. There are differences with respect to Ref. [67] in the transitions including the exchange of heavy vectors, but since this interaction is small, it is not surprising that in the end the results of Ref. [74] using SU(4) and those of Ref. [67] where SU(4) is not used, are very similar. The approach described above has been very successful, and in Ref. [67] three of the c states of Ref. [3] were correctly reproduced in mass and width. In Ref. [75] heavy quark spin symmetry (HQSS) was used to find the relationship between the transition matrix elements of the coupled channels,D ( * ) ( * ) c , J/ψ N and others, to describe the recent hidden charm pentaquark states [2]. The strength of the interaction was obtained from the evaluation of the hidden gauge approach described above. Once again, the states found in the experiment were fairly well reproduced and a few more states were predicted. These results are similar to those of Ref. [76], where also HQSS is used to evaluate masses, but in addition the widths are evaluated in Ref. [67].
The approach of Ref. [67] to study the c states is also used in Ref. [66] in the study of the cc states, in Ref. [77] in the study of c and b states and in Ref. [78] in the study of b states.
In the present work we study in detail the bc states that can emerge from the interaction of pseudoscalar-baryon (1/2 + ) interaction, vector-baryon (1/2 + ) interaction and pseudoscalar-baryon (3/2 + ) interaction. Using the same regulator obtained in cases where we could contrast with experiment, and the same source of information, we obtain several states in each sector. We evaluate binding energies and widths, as well as couplings of the resonant states to the different channels. In some cases we can see a striking dominance of one of the channels, which allows us to deem the state as a molecular state of this channel. Since we work in meson-baryon interaction in s-wave, we also evaluate the wave functions at the origin for the different channels, which provide extra information concerning the relevance of the channels in different reactions.

Baryon states
In order to see the coupled channels that we need, we classify the meson-baryon states as (2) Meson-baryon states with b in the baryon and c in the meson. The coupled channels are Next we take the baryon wave functions isolating the heavy quarks and imposing the spin-flavor symmetry on the light quarks. In our approach it is important to specify the spin of the states because, as we shall see below, the interaction is spin independent, which immediately imposes selection rules in the transitions. The wave functions are summarized in Tables 1 and 2. The corresponding states with different charge are trivial using the u, d quarks. In Tables 1 and 2 χ M S , χ M A , χ S are the spin wave functions of the three quarks [79].
In the interaction we exchange vector mesons as shown in Fig. 1. The lower vertex of V B B is of the type γ μ μ , but, with heavy baryons and close to threshold, only the γ 0 1 term is relevant, which means that this vertex is spin independent, and so is the upper vertex that will go as ( p D + p D ) 0 . For V B → V B transitions, the upper vertex has the same structure, but with the additional · factor for the vector polarizations, which is diagonal in the spin of the vectors. With the spin-independent interaction, we can classify the (A) P B channels with χ M S for the baryon: π bc , η bc , We do not study the interaction of vectors with J P = 3/2 + baryons. From previous works the states obtained belong to a region where signals are difficult to see experimentally.

Isospin states
We take the isospin multiplets: We find the following states classified as B P, BV . For consistency with Ref. [66] we write the states as baryonmeson. Then and the rest are trivial.

Evaluation of matrix elements
The evaluation of the upper vertex of Fig. 1, V P P, is done in Ref. [68] using the vector character of the vertex and the quark content of the mesons. Yet, it was found that for practical purpose one can get it from the Lagrangian where · · · stands for the matrix trace, when dealing with charmed mesons, and when dealing with bottomed mesons. The lower vertex is trivial to evaluate. The Lagrangian for the V B B vertex is given by the operator Hence, for instance, ρ 0 0 bc 0 bc involves the vertex Then we get finally transition matrix elements where k 0 , k 0 are the energies of the mesons and D i j the coefficients which are shown in Tables in the next section.
With the potential of Eq. (27) we solve the Bethe-Salpeter equation in coupled channels where G is the meson-baryon loop function and we use q max = 650 MeV as in Refs. [66,67,78]. ω l , E l are the energies of the meson and baryon respectively, ω l = m 2 l + q 2 , E l = M 2 l + q 2 , and m l , M l the meson and baryon masses. In Eq. (29) p 0 is the energy of the incoming baryon.
In order to get poles in the second Riemann sheet we replace G by G I I , and it is given by where √ s th,l is the threshold mass of the l channel, and For the evaluation of the couplings g l of the state to different coupled channels, we find that the T i j matrix can be expressed in the following form close to the poles, z R , which defines the couplings g i up to a global sign of one of them.
The approach followed relies on the exchange of vector mesons, which we justified from the extension of the chiral Lagrangians via the local hidden gauge approach. In mesonbaryon interactions the exchange of pseudoscalar mesons is also sometimes used [80], together with σ -exchange. Unlike in baryon-baryon interaction where the π exchange plays a very important role, in meson-meson or meson-baryon interaction, pion exchange plays a more moderate role because the direct P B → P B through π exchange does not go, since a three-pseudoscalar vertex is forbidden by parity-angularmomentum conservation. It is through intermediate states transitions P B → V B → P B that the π -exchange can have a contribution. Detailed calculations of its effects have also been considered in the study of baryon states with open charm [81] and hidden charm [82] and the main source of interaction remains vector exchange. Also, to some extend, the two step process P B → V B → P B can be incorporated into an effective transition potential δV for P B → P B, adding to the vector exchange potential, and this extra potential can again be effectively accounted for by changes in the cut off that regularizes the loop function G, since [V −1 − G] will be the same with an extra δV and δG = δ(V −1 ). Yet, not all the contributions can be reabsorbed in this way and one may expect some remnant contributions to break the degeneracy between spin-parity 1/2 − and 3/2 − . That these effects are finally small can be seen in the breakup of the degeneracy of the P c (4440) and P c (4457) states recently observed in Ref. [2]. There are several works including pion exchange that break the degeneracy of the 1/2 − , 3/2 − states around this energy [83][84][85][86], but the small difference of mass between these states gives us an idea of the role played by π exchange in these cases. Uncertainties of this type in our predictions must certainly be admitted.
As to σ -exchange, it is empirically accounted for in some of the former works, but the strength is unknown. However, there is one way to make this exchange quantitative by going to the root of the σ meson as dynamically generated from the ππ interaction [87,88]. Indeed, in Ref. [89] the σ -exchange between nucleons was studied from the point of view of the exchange of two interacting mesons in s-wave. The same picture was used in the study of the meson-meson interaction describing the Z c (3980) [90] and Z c (4000) [91] and its effect was found minor. In particular, for meson-baryon interaction, the σ -exchange from this perspective was considered in Ref. [92] in the study of theK N interaction and it was found very much suppressed.

Pseudoscalar-baryon (1/2 + ) states, mixed symmetric sector
In Table 3, we consider the channels in this sector together with their threshold masses. 1 In Table 4 we show the coefficients D i j of Eq. (27). We find that there are some terms which go with λ. These terms involve the exchange of D * and they are reduced.   The value of λ is the reduction factor versus the exchange of a light vector, which in Ref. [67] is found to be around λ = 0.25 as in Ref. [93]. Apart from that, when exchanging B c meson we set the D i j result as ∼0, since it is highly suppressed.
In Table 5 the poles for three states that appear in this sector are shown for different values of the cutoff q max . As in former works we take q max = 650 MeV. In Tables 6 and 7 we take the first two states of Table 5 and show the couplings g i , and the wave functions at the origin g i G i (see Ref. [94]). We Table 8 The coupling constants to various channels and g i G I I i in MeV with q max = 650 MeV      Table 5 for q max = 650 MeV, appears at the same energy as the former one but the couplings are different as can be seen in Table 8.
As it can be seen in Table 5, the last two states in the row of q max = 650 MeV are less stable with the changes of q max and, hence, more uncertain.

Pseudoscalar-baryon (1/2 + ) states, mixed antisymmetric sector
We follow here the same pattern as in the former subsection with the mixed antisymmetric states. In Table 9 we show the channels and the threshold masses. In Table 10 we show the D i j coefficients. In Table 11 we show the states (poles) found in this sector. We find three clear states, one with a large width of almost 200 MeV and two more states with very narrow widths. The first state is actually the same one that we found before, because it couples mostly to bc π as shown in Table 12, and this state appears with M S and M A spins. However, the second state couples mostly to b D and b D s as shown in Table 13, and is then a different state. The state at 7462 MeV couples mostly to cB and cBs , as shown in Table 14.

Vector-baryon (1/2 + ) states, mixed symmetric sector
We show in Table 15 the channels and the thresholds that in this case give rise to degenerate states in sector J P = 1/2 − , 3/2 − . The D i j coefficients are shown in Table 16 and in Table 17 we find four states, three of them with zero width and the last one with a large width. In Tables 18, 19 20 and 21 we show the couplings and wave functions of these states. The state at 7418 MeV could be regarded as a cB * bound state, but it also couples strongly to cB * s . The state at 7501 MeV could be identified with a b D * state, the one at 7595 MeV as a bc ρ bound state, and the state at 7837 MeV and 180 MeV corresponds to a bc K * bound state which decays strongly in bc ω.

Vector-baryon (1/2 + ), mixed antisymmetric sector
In Table 22 we show the coupled channels and the thresholds and in Table 23 we show the D i j coefficients. The states are shown in Table 24. The couplings of the states to the different channels are shown in Tables 25, 26

Pseudoscalar-baryon (3/2) + states
In this case the spin wave function, χ S , is full symmetric and the states generated are in the sector J P = 3/2 − . In Table 29 we show the channels and the thresholds and in Table 30 the D i j coefficients. In Table 31 we show the states obtained for different values of the cut off. We observe three states, one with a width of about 190 MeV and the other two narrow. Inspecting Tables 32, 33, 34 we can see that the first state couples strongly to * bc π , which is open, and this is the reason for the large width. The second state couples mostly to * b D and the third one to * cB .

Wave functions
It is interesting to see how are the wave functions that we have generated. They are a bit different than ordinary wave functions with local potentials decreasing very rapidly as r → ∞. In order to see that, we have to go back to the work of Ref. [94], where we find that the use of Eq. (28) with a G function regularized with a cut off q max is equivalent to solving the Lippmann-Schwinger equation with a potential which is a non-local potential. The wave function in momentum space is particularly easy (see Eqs. (34), (47) of Ref. [94] and Eq. (105) of Ref. [95] generalizing to relativistic energies) where w i ( q ) = m 2 i + q 2 and g is the coupling of the wave function to the channel with particles 1 and 2.
The wave function in coordinate space is given by It is interesting to note that unlike in ordinary potentials, which could be simulated with q max → ∞ and w 1 ( q ) + w 2 ( q ) and sin(qr) providing convergence in the q integration, in our case q max = 650 MeV corresponds to a value where q 2 /m 2 i is very small and the w 1 ( q ) + w 2 ( q ) term does not help in the convergence of the integral, which is determined by q max . One can then see that the shape of the wave function does not depend much on E, which diverts from ordinary wave functions with a rapidly decreasing local potential where the size is roughly given by r = 1/ √ 2μB, with μ the reduced mass and B the binding. We show this in an example of a very bound component, the 7372 MeV state of Table 7, which couples mostly to b D.
In Fig. 2 we show the wave function corresponding to this channel. While naively we would expect a size of around 0.2 fm according to the intuitive formula, we find that the wave function extends much further and even at r 1 fm is not     negligible. This can be better appreciated in Fig. 3 where we plot the wave function squared times r 2 . We see that it peaks around 0.7 fm and still has a sizable strength around 1 fm and beyond.

Conclusions
We have studied the meson-baryon interaction in the sector corresponding to bc quantum numbers. We take the cou-    pled channels that can lead to these quantum numbers and study their interaction in s-wave. The model used for the interaction is based on an extrapolation of the local hidden gauge approach, which uses vector exchange as the source of interaction. The dominant terms come from the exchange of light vector mesons which render the heavy quarks as spectators and the approach automatically satisfies the heavy quark symmetry rules. The interaction is properly unitarized in coupled channels and by looking at poles in the second Riemann sheet we look for the states of the ( ) × 10 5 system. We consider the interaction of pseudoscalars with baryons of J P = 1/2 + , 3/2 + and of vectors with baryons of J P = 1/2 + and distinguish spin mixed symmetric and mixed antisymmetric states in analogy to the and states. We find several states which correspond to bound meson-baryon states with zero or small widths and a few that have a large width, Since the input used to generate these states is the same one used to study c and states of hidden charm, which produced results in excellent agreement with experiment, we have confidence that the predictions made are realistic and encourage the experimental search of such states. In some cases the main decay channels have been identified and this can be useful in the planning of experimental proposals. 11847317 and 11565007. This work is partly supported by the Spanish Ministerio de Economia y Competitividad and European FEDER funds under Contracts No. FIS2017-84038-C2-1-P B and No. FIS2017-84038-C2-2-P B, and the Generalitat Valenciana in the program Prometeo II-2014/068, and the project Severo Ochoa of IFIC, SEV-2014-0398.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: All data generated during this study are contained in this published article.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .