Hidden-bottom molecular states from Σb(∗)B(∗)-ΛbB(∗)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varSigma ^{(*)}_bB^{(*)}-\varLambda _bB^{(*)}$$\end{document} interaction

In this work, we study possible hidden-bottom molecular pentaquarks Pb\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_b$$\end{document} from coupled-channel Σb(∗)B(∗)-ΛbB(∗)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varSigma ^{(*)}_bB^{(*)}-\varLambda _bB^{(*)}$$\end{document} interaction in the quasipotential Bethe-Salpeter equation approach. In isodoublet sector with I=1/2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$I=1/2$$\end{document}, with the same reasonable parameters the interaction produces seven molecular states, a state near ΣbB\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \varSigma _bB$$\end{document} threshold with spin parity JP=1/2-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J^P=1/2^-$$\end{document}, a state near Σb∗B\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varSigma ^*_bB$$\end{document} threshold with 3/2-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3/2^-$$\end{document}, two states near ΣbB∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varSigma _bB^*$$\end{document} threshold with 1/2-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1/2^-$$\end{document} and 3/2-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3/2^-$$\end{document}, and three states near Σb∗B∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varSigma _b^*B^*$$\end{document} threshold with 1/2-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1/2^-$$\end{document}, 3/2-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3/2^-$$\end{document}, and 5/2-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$5/2^-$$\end{document}. The results suggest that three states near Σb∗B∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varSigma _b^* B^*$$\end{document} threshold and two states near ΣbB∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varSigma _b B^*$$\end{document} threshold are very close, respectively, which may be difficult to distinguish in experiment without partial wave analysis. Compared with the hidden-charm pentaquark, the Pb\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_b$$\end{document} states are relatively narrow with widths at an order of magnitude of 1 MeV or smaller. The importance of each channel considered is also discussed, and it is found that the ΛbB∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varLambda _b B^*$$\end{document} channel provides important contribution for the widths of those states. In isoquartet sector with I=3/2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$I=3/2$$\end{document}, cutoff should be considerably enlarged to achieve bound states from the interaction, which makes the existence of such states unreliable. The results in the current work are helpful for searching for hidden-bottom molecular pentaquarks in future experiments, such as the COMPASS, J-PARC, and the Electron Ion Collider in China (EicC).


Introduction
It is one of the most important topic in hadron physics community to search for the hadronic exotic states beyond the conventional quark model. Among the theoretical pictures in the market, molecular state is a competitive one to explain existing candidates of exotic states, such as the XY Z particles and P c states [1]. A molecular state is analogous to a nucleus, especially the deuteron, that is, a loosely bound state of two or more hadrons. It immediately leads to a conclusion a e-mail: junhe@njnu.edu.cn (corresponding author) that a molecular state is close to the threshold of constituent hadrons. In practice, the study of the molecular state also focuses on resonance structures near thresholds. Vise versa, if we can find more structures near thresholds, especially those with corresponding relationship, it will strongly support existence of molecular states. In the current work, we will provide predictions of hidden-charm pentaquarks P b , which are partners of the hidden-charm P c states.
The observation of hidden-charm pentaquarks at LHCb is a great breakthrough of the study of exotic states [2,3]. It is also an important support on the molecular state picture. Three narrow resonance structures were reported at LHCb in an update measurement as P c (4457) and P c (4440) states near Σ cD * threshold and a P c (4312) state near Σ cD threshold [2]. Combined with P c (4380) near Σ * cD threshold suggested in the first observation [3], it exhibits a good pattern of the S-wave molecular states from interactions corresponding to the thresholds. Such observation confirms the prediction of existence of the hidden-charm pentaquark in some models [4][5][6][7][8] A lot of theoretical interpretations of these structures emerged after the experimental observation. Due to the strong correlation between these structures and the thresholds, the molecular state is the most popular picture to explain the P c states [9][10][11][12][13][14][15][16][17][18], though other interpretations can not be excluded [19][20][21][22]. In Ref. [23], authors even proposed existence of seven hidden-charm molecular states as a complete heavy-quark spin symmetry multiplet. In our previous works [24,25], we systematically investigate coupledchannel Σ 1/2 − and 3/2 − were also produced with the same parameters, but the result suggests that their effects on the experimental observable may be small. Besides, the decay pattern was also discussed, and the ΛD * channel is found dominant in the decays of these states. Now that experimentally observed P c states were well interpreted in the molecular state picture, we can predict hidden-bottom states above 11 GeV. Though there are a large amount of works about P c states reported in the literature, the studies about the hidden-bottom pentaquark are still inadequate. There are a few incidental studies about P b states in the works to interprete the P c states in the molecular picture, such as within constituent quark model [17,26], the chiral effective field theory [9,27], and the Bethe-Salpeter Eq. [28]. The width and decay pattern of the P b states were also discussed in Refs. [17,[28][29][30][31]. Due to very large mass of the hidden-bottom pentaquark, it is relatively difficult to search for the P b state in experiment compared with the P c state. In Ref. [32], pion and photon induced productions of hiddenbottom pentaquarks were studied, the calculation suggests that it is possible to search for these states at COMPASS, J-PARC and EicC. Hence, it is interesting to perform a systematical study about the hidden-charm pentaquarks based on experimental information and theoretical analysis about the P c states.
In this work, we will investigate coupled-channel Σ interaction in the quasipotential Bethe-Salpeter equation (qBSE) approach to find possible hidden-bottom molecular states. The interaction was described in the one-boson-exchange model with the help of the effective Lagrangians within the heavy quark symmetry and chiral limit as in Refs. [24,25], where the P c states were interpreted. The masses and widths of molecular states are predicted by finding poles in complex energy plane. The decay channels of predicted states will be discussed also.
This article is organized as follows. After introduction, the details of theoretical frame of coupled-channel Σ ) interactions is presented in Sect. 2. In Sect. 3, the single-channel results of the states with isospin I = 1/2 and I = 3/2 are given first. Then, coupled-channel results are presented, and the importance of the channels considered are discussed. Finally, summary and discussion will be given in Sect. 4.

Theoretical frame
In the qBSE approach, we will use the one-boson-exchange interaction of two bottom hadrons as dynamical kernel. In the current work, we will adopt the Lagrangians with heavy quark limit and chiral symmetry, and the channels with hidden-charm mesons are ignored as in Ref. [24,25] to keep the consistence. The peseudoscalar P, vector V and scalar σ exchanges will be considered, and the effective Lagrangians depicting the couplings of light mesons and bottom mesons or bottom baryons are required and will be presented in the below.
First, we consider the couplings of light mesons to heavylight bottom mesons P = (B 0 , B + , B + s ). The Lagragians were constructed in the literature as [33][34][35][36], where The P and V are the pseudoscalar and vector matrices as The explicit forms of the Lagrangians for the couplings of light mesons to bottom baryons can be written as [37], where B iμ is defined as and the bottomed baryon matrices are defined as In the calculation, the masses of particles are chosen as suggested central values in the Review of Particle Physics (PDG) [38]. The mass of broad σ meson is chosen as 500 MeV. The coupling constants involved was cited from the literature [15,37,39,40], and listed in Table 1, With the vertices obtained from the above Lagrangians, the potential of couple-channel interaction can be constructed. Because six channels are invovled in the current work, it is tedious and fallible to give explicit 36 potential elements and input them into code. Instead, in this work, we input vertices Γ and propagators P into code directly, and the potential can be obtained as The propagators are defined as usual as where the form factor f (q 2 ) is adopted to compensate the offshell effect of exchanged meson as f ( with m e being the m P,V,σ and q being the momentum of the exchanged meson. The cutoff is rewritten as a form of Λ e = m +α e 0.22 GeV. The f I is the flavor factor for certain meson exchange of certain interaction, and the explicit values are listed in Table 2.
With the potential kernel obtained, we use the qBSE to solve the scattering amplitude [41][42][43][44][45]. After partialwave decomposition and spectator quasipotential approximation, the 4-dimensional Bethe-Saltpeter equation in the Minkowski space can be reduced to a 1-dimensional equa-tion with fixed spin-parity J P as [42], where the sum extends only over nonnegative helicity λ .
Here, the reduced propagator with the spectator approximation can be written as with p h,l and m h,l being the momenta and masses of heavy or light constituent particles. The partial wave potential is defined with the potential of interaction obtained in the above in Eq. (6) as where η = P P 1 P 2 (−1) J −J 1 −J 2 with P and J being parity and spin for system, B ( * ) meson or Σ ( * ) b baryon. The initial and final relative momenta are chosen as p = (0, 0, p) and p = (p sin θ, 0, p cos θ). The d J λλ (θ ) is the Wigner dmatrix. we also adopt an exponential regularization by introducing a form factor into the propagator as [42] In the current work, the relation of the cutoff Λ r = m + α r 0.22 GeV with m being the mass of the exchanged meson is also introduced into the regularization form factor to suppress large momentum, i.e., the short-range contribution of the π exchange as warned in Ref. [46].

Numerical results
The 1-dimensional integral equation can be transformed into a matrix equation as M = V + V G 0 M by Gauss discretization. The molecular states can be found by searching for the pole of scattering amplitude M in complex energy plane at |1 − V (z)G(z)| = 0 with z = W + iΓ /2 equaling to system energy W at real axis [42]. In addition, we take two free parameters α e and α r as α for simplification.

Single-channel results
Each experimental observed P c state is close to a threshold, respectively [2,3]. It suggests that each of these states should be mainly from a single-channel interaction in the molecular state picture, which is confirmed by previous study in Table 1 The coupling constants adopted in our calculation. The λ and λ S,I are in the unit of GeV −1 Others are in the unit of 1 Table 2 The flavor factors f I for certain meson exchanges of certain interaction. The values in bracket are for the case of I = 3/2 if the values are different from these of I = 1/2 Ref. [25]. In this work, we present the single-channel results first.
In the current work, we consider all states with spin parities which can be produced from S-wave interaction, The results for isodoublet with I = 1/2 are illustrated in Fig. 1. The Λ b B ( * ) with 1/2 − , 3/2 − and Λ b B with 1/2 − are also calculated, however, large α beyond reasonable limit is required to produce bound states.
The results suggest that bound states can be produced in all seven cases in a range of α from 0 to 3.5. All states appear at α about 0.5, and binding becomes deeper with the increase of α, and gradually reach to a binding energy about 20 MeV or the next threshold at α = 2.5−3.5 . The binding energies will continue to increase, but we no longer present such results. The trends of curves for three states produced from the Σ * b B * interaction with different spin parities are almost the same. Such phenomenon can also be found for two curves for two states from the Σ b B * interaction. Compared with the results for P c states [24,25], one can find that the values of parameter α to produce the hidden-bottom molecular states is relatively smaller.
In Fig 2, we present the results for isoquartet states with I = 3/2. In the calculation, we also consider seven cases as for isodoublet. No bound state is produced form Σ b B interaction with (1/2 − ) and Σ * b B interaction from (3/2 − ) even if α is taken to 9. Except these two states, left five states can be produced from single-channel interaction as shown in 2, but with considerablly large α. The production of bound states, , needs a value of α at least 4, which is larger than the maximum value of α required for binding of isodoublet molecular states with I = 1/2. It indicates that these three states are hardly to be found if the isoboublet states have small binding energies. The rest two states Σ b B * (3/2 − ) and Σ * b B * (5/2 − ) appears at α = 3 and 2.5, respectively. it implies that these two Generally speaking, if we assume that the P b states are also loosely bound states as P c states, the possibility of existence of isoquaret states is vey small. In Ref [26], within the frame of constituent quark model, the molecular states with I = 3/2 were also not found.

Coupled-channel results
In the previous calculation, the bound states from single channel calculation exhibit as poles at real axis of complex energy plane, that is, the widths of these states are zero. In the case of the P c states, we found that the experimental width can be well reproduced with inclusion of coupledchannel effect [25]. In the above single-channel calculation, seven bound states are produced. Those states can be coupled to each other by exchanges of light mesons. Besides, the ΛD * channel is also found important for the width of P c states [25,30,31]. In the following, coupled-channel results for the Σ ( * ) b B ( * ) − Λ b B ( * ) interaction will be given. Here, we first give an example to show a general picture of coupled-channel results. Since there is no experimental data about the P b states, we should choose a parameter to present the results. The only free parameter in our model is α, and in the single-channel calculation the α dependences of the masses of seven bound states exhibit a similar trend. Hence, we choose α as 1.5 to illustrate the poles from coupledchannel Σ Fig. 3. The values of log |1 − V (z)G 0 (z)| with variation of complex energy z is adopted to show the positions of poles of coupled-channel scattering amplitude because M = (1 − V G 0 ) −1 V . And we present the results for spin parities 5/2 − , 3/2 − , and 1/2 − in a range from 11.06 to 11.16 GeV for real part of complex energy Re(z) and -2 to 2 MeV for imaginary part I m(z). 11.08 11.10 11.12 11.14 11.16 11.06 11.08 11.10 11.12 11.14 11.16 One can find that there are still seven poles produced as in the single-channel calculation. It suggests that only states with an S-wave interaction can be produced for three spin parities considered. In the case with J P = 1/2 − , there exist three poles near the Σ b B (S wave), Σ b B * (S and D waves) and Σ * b B * (S, D and G waves) thresholds, respectively. No pole appears near Σ * B (P wave) threshold for 1/2 − . In the case with J P = 3/2 − , we also have three poles near the Σ * b B (S and D waves), Σ b B * (S and D waves) and Σ * b B * (S, D and G waves) thresholds, respectively. And no pole appears near Σ b B (P wave) threshold. For spin parity 5/2 − there is only one pole near the Σ * b B * thresholds, and only this channel can produce a pole with 5/2 − in S wave.
The P c (4457) and P c (4440) is close to each other near Σ cD * threshold, which was even taken as one resonance structure in the first observation of the P c states [2]. In the P b case, our results exhibit a more serious overlapping between two states near Σ b B * as in the single channel calculation. A shown in Fig. 3, these two poles with 1/2 − and 3/2 − have almost the same mass and width. Furthermore, the masses of three poles near the Σ * b B * threshold are also close very much, which are 11149.2 MeV, 11149.7 MeV and 11150.1 MeV corresponding to the two molecular states with 1/2 − , 3/2 − and 5/2 − , respectively. Different from the P c states [25], here, the state with 5/2 − stands out background obviously while other two poles with J P = 1/2 − and 3/2 − are very dimly, and may be difficult to be found at experiment.
As shown in the figure, the poles acquire imaginary parts after coupled-channel effects are included in the calculation. However, these states are generally very narrow, with imaginary parts smaller than 0.5 MeV. By using the rela- Table 3 The masses and widths of molecular states at different values of α. The "CC" means full coupled-channel calculation. The values of the complex position means mass of corresponding threshold subtracted by the position of a pole, M th − z, in the unit of MeV. The two short line "−−" means the coupling does not exist. The imaginary part of some poles are shown as "0.00", which means too small value under the current precision chosen tion Γ = −2 Im(z), it means a small width about 1 MeV or smaller, which is much smaller than the P c states with similar binding energy. It is worth mentioning that at α = 1.5, except for Σ * b B * (5/2 − ), the masses of six molecular states after inclusion of coupled-channel effect are in good agreement with the results obtained under the frame of the constituent quark model in Ref [17].
In the above, we only present the poles of molecular states obtained by the coupled-channels calculation at a value of α = 1.5. Because there is no experimental data, in the following, we will present the results with different vales of α from 0.5 to 2.0 to show the dependence of results on the parameter in second and third columns of Table 3. The twochannel calculation results are also listed in the fourth to eighth columns to show the role of each channel on widths of molecular states. Here, to emphasize the nearest threshold, we replace the real part of pole by z → M th − z with M th being the mass of nearest higher threshold.
In the first column, we list thresholds with certain spin parity, and the result of pole under the corresponding threshold with different α is given in the second and third columns with full coupled-channel Σ ( * ) b B ( * ) −Λ b B ( * ) interaction. One can find, except a small width is acquired, the results are similar to those from the single-channel calculation. Hence, we take such channel as production channel of this pole. All poles appear on threshold at about α = 0.5 and leave the threshold with the increasing of α. If we choose a binding energy about 10 MeV, the widths of most states are very small, about 1 MeV or smaller.
In the fourth to eighth columns, we consider two-channel result with the coupling between the production channel and a channel below it. The imaginary part reflects the strength of couplings between two channels. Since the pole is mainly from the production channel, the effect of a channel on the pole can be also estimated from the two-channel resutls. Because the width with smaller α is very small, in the followings, we focus on the results at larger α, 1.5 and 2.0.
Three states near Σ * b B * threshold, which is the highest threshold of channels considered in the current work, can decay into five channels. Among these decay channels, the Λ b B * channel has strongest couplings to this three states. For the two states near Σ b B * threshold, there are four decay channels, the Λ b B * channel is much stronger than other channels for 3/2 − . For 1/2 state, both Σ * b B and Λ b B * channels couples strongly to the Σ b B * channel. Among three possible decay channels of the state near Σ * b B thresholds, only Λ b B * channel provides large width, and other channels only give very small imaginary part of the position. For the Σ b B(1/2 − ) case, only Λ b B ( * ) channels involves, among which, the Λ b B * channel is still dominant one. Hence, for all seven states, the Λ b B * channel is the most important one in all channels considered, which is consistent with the results of the hidden-charm pentaquarks in Refs. [25,30,31].

Summary and discussion
In this work, the masses and widths of hidden-bottom molecular pentaquarks are predicted from coupled-channel Σ ( * ) b B ( * ) − Λ b B ( * ) interaction in the qBSE approach with the help of effective Lagrangians with heavy quark and chiral symmetries. The results suggest that seven molecular states can be produced from the interactions. All states appear at α about 1, which corresponds to reasonable radii of the constituent hadrons, about 0.5 fm.
Among the seven states, three of them are near the Σ * b B * threshold, and the masses of these three states are very close. The two sates with 1/2 − and 3/2 − are very weak compared with the state with 5/2 − . Hence, these three states should exhibit as one resonance structure without partial wave analysis. Even with partial wave analysis, the states with 1/2 − and 3/2 − are difficult to be distinguished from the one with 5/2 − . The two states near Σ b B * are also mixing together but a partial-wave analysis will be helpful to distinguish them. Hence, the results suggest that four resonance structures may be observed in experiment, though there exist seven molecular states from coupled-channel Σ Compared with hidden-charm P c states, the widths of P b states are much smaller, about 1 MeV or smaller. And the calculation suggests that the Λ b B * channel has strong couplings to the molecular states, which is due to the strong couplings of vertex Λ b Σ ( * ) b π . The small width have both advantage and disadvantage in experimental observation of such states. The small width makes the production possibility small, which needs high luminosity of experimental facility. However, a small width also makes the peak of state stand out obviously from background in experiment. In Ref. [32], we study the possibility to search for such states in pion and photon induced productions. The results suggest that with small widths the measurement of the P b states is promising at the such as the COMPASS J-PARC, especially the Electron Ion Collider (EicC) in China.