Investigating tetraquarks composed of $us\bar{d}\bar{b}$ and $ud\bar{s}\bar{b}$

In the framework of the quark delocalization color screening model, we investigate tetraquarks composed of $us\bar{d}\bar{b}$ and $ud\bar{s}\bar{b}$ in two structures: meson-meson structure and diquark-antidiquark structure. Neither bound state nor resonance state is found in the system composed of $us\bar{d}\bar{b}$. The reported $X(5568)$ cannot be explained as a molecular state or a diquark-antidiquark resonance of $us\bar{d}\bar{b}$ in present calculation. However, two bound states of the diquark-antidiquark structure are obtained in the tetraquarks system composed of $ud\bar{s}\bar{b}$: an $IJ=00$ state with the mass of $5701$ MeV, and an $IJ=01$ state with the mass of $5756$ MeV, which maybe the partners of $X(5568)$ states. Our results indicate that the diquark-antidiquark configuration would be a good choice for the tetraquarks $ud\bar{s}\bar{b}$ with quantum numbers $IJ=00$ and $IJ=01$. The tetraquarks composed of $ud\bar{s}\bar{b}$ is more possible to form bound states than the one composed of $us\bar{d}\bar{b}$. These bound states are worth investigating in future experiments.


I. INTRODUCTION
In the past few decades, the discovery of numbers of exotic states stimulated extensive interest in understanding the structures of the multiquark hadrons. So far, most tetraquark and pentaquark candidates are composed of hidden charm or bottom quarks. However, the new state X(5568) observed by the D0 collaboration in 2016 [1] was an exception. The X(5568) has a mass m = 5567.8 ± 2.9(stat) +0.9 −1.9 (syst) MeV and width Γ = 21.9 ± 6.4(stat) +5.0 −2.5 (syst) MeV [1]. The decay mode is X(5568) → B 0 s π ± , which indicates that the quark component of the X(5568) should be four different flavors: u, d, s, b. Therefore, the claimed X(5568), if confirmed, would differ from any of the previous observations, as it must be a tetraquark state with usdb or dsūb and their charge-conjugated ones. Unfortunately, this state was not confirmed by other collaborations. The LHCb collaboration [2], the CMS collaboration of LHC [3], the CDF collaboration of Fermilab [4] and the ATLAS Collaboration of LHC [5] all claimed that no evidence for this state was found. Nevertheless, the D0 collaboration's new result still insists on the existence of this tetraquark X(5568) [6]. Clearly, more other measurements are needed.
To search for the tetraquark states with four different flavors, a better state is udsb (or its charge-conjugated one) with replacing the d/s in X(5568) byd/s [23]. Obviously, such state is a partner of X(5568) under the SU (3) flavor symmetry, and their masses are close to each other. But the threshold of udsb is BK, 270 MeV higher than the threshold B s π of X(5568) with usdb. So there is large mass region for this udsb state below threshold and being stable. Besides, Ref. [23] pointed out that if the lowestlying udsb state exists below threshold, it can be definitely observed via the weak decay mode J/ψK − K − π + , with the expectation of hundreds of events in the current LHCb data sample but rejecting backgrounds due to its long lifetime. Therefore, the udsb state would be a more promising detectable tetraquark state. Ref. [24] investigated such state composed of udsb within the chiral quark model, and found the bound state with IJ P = 00 + was possible. Liu et al. also proposed several partner states of X(5568) and estimated the mass difference of these partner states based on the color-magnetic interaction [25], which can provide valuable information on the future experimental search of these states.
It is generally known that quantum chromodynamics (QCD) is the fundamental theory of the strong interaction. Understanding the low-energy behavior of QCD and the nature of the strong interacting matter, however, remains a challenge due to the complexity of QCD. Lattice QCD has provided numerical results describing quark confinement between two static colorful quarks, a preliminary picture of the QCD vacuum and the internal structure of hadrons in addition to a phase transition of strongly interacting matter. But a satisfying description of multiquark system is out of reach of the present calculation. The QCD-inspired models, incorporating the properties of low-energy QCD: color confinement and chiral symmetry breaking, are also power-ful tools to obtain physical insights into many phenomena of the hadronic world. Among many phenomenological models, the quark delocalization color screening model (QDCSM), which was developed in the 1990s with the aim of explaining the similarities between nuclear (hadronic clusters of quarks) and molecular forces [26], has been quite successful in reproducing the energies of the baryon ground states, the properties of deuteron, the nucleon-nucleon (N N ) and the hyperon-nucleon (Y N ) interactions [27]. Recently, this model has been used to study the pentaquarks with hidden-strange [28], hiddencharm and hidden-bottom [29]. Therefore, it is interesting to extend this model to the tetraquark system. In present work, the tetraquark state X(5568) with quark contents usdb and its partner state with udsb are investigated. Besides, two structures, meson-meson and diquark-antidiquark, are considered in this work.
The structure of this paper is as follows. A brief introduction of the quark model and wave functions is given in section II. Section III is devoted to the numerical results and discussions. The summary is shown in the last section.

II. MODEL AND WAVE FUNCTIONS
QDCSM has been described in detail in the literatures [26,27]. Here, we just present the salient features of the model. The Hamiltonian for the tetraquark states is shown below: Where S ij is quark tensor operator; Y (x) and H(x) are standard Yukawa functions; T c is the kinetic energy of the center of mass; α s is the quark-gluon coupling constant; g ch is the coupling constant for chiral field, which is determined from the N N π coupling constant through The other symbols in the above expressions have their usual meanings. All model parameters are determined by fitting the meson spectrum we used in this work and shown in Table I. The calculated masses of the mesons in comparison with experimental values are shown in Table II. Besides, a phenomenological color screening confinement potential is used here, and µ ij is the color screening parameter, which is determined by fitting the deuteron properties, N N scattering phase shifts, N Λ and N Σ scattering phase shifts, respectively, with µ uu = 0.45 fm −2 , µ us = 0.19 fm −2 and µ ss = 0.08 fm −2 , satisfying the relation, µ 2 us = µ uu µ ss [30]. When extending to the heavy bottom quark case, there is no experimental data available, so we take it as a adjustable parameter µ bb = 0.001 ∼ 0.0001 fm −2 . We find the results are insensitive to the value of µ bb . So in the present work, we take µ bb = 0.001 fm −2 .
The quark delocalization in QDCSM is realized by   specifying the single particle orbital wave function of QD-CSM as a linear combination of left and right Gaussians, the single particle orbital wave functions used in the ordinary quark cluster model, Here s i , i = 1, 2, ..., n are the generating coordinates, which are introduced to expand the relative motion wavefunction [27]. The mixing parameter ǫ(s i ) is not an adjusted one but determined variationally by the dynamics of the multi-quark system itself. In this way, the multiquark system chooses its favorable configuration in the interacting process. This mechanism has been used to explain the cross-over transition between hadron phase and quark-gluon plasma phase [32]. In this work, the resonating group method (RGM) [33], a well-established method for studying a bound-state or a scattering problem, is used to calculate the energy of all these states. The wave function of the four-quark system is of the form where ψ L , ψ σ , ψ f , and ψ c are the orbital, spin, flavor and color wave functions, respectively, which are given below. The symbol A is the anti-symmetrization operator. For the meson-meson structure, A is defined as where 1 and 3 stand for the quarks in two meson clusters respectively; for the diquark-antidiquark structure, A = 1.
The orbital wave function is in the form of where R 1 and R 2 are the internal coordinates for the cluster 1 and cluster 2. R = R 1 − R 2 is the relative coordinate between the two clusters 1 and 2. The ψ 1 and ψ 2 are the internal cluster orbital wave functions of the clusters 1 and 2, and χ L (R) is the relative motion wave function between two clusters, which is expanded by gaussian bases where s i is called the generate coordinate, n is the number of the gaussian bases, which is determined by the stability of the results. By doing this, the integro-differential equation of RGM can be reduced to an algebraic equation, generalized eigen-equation. Then the energy of the system can be obtained by solving this generalized eigenequation. The details of solving the RGM equation can be found in Ref. [33]. In our calculation, the maximum generating coordinate s n is fixed by the stability of the results. The calculated results are stable when the distance between the two clusters is larger than 6 fm. To keep the dimensions of matrix manageably small, the two clusters' separation is taken to be less than 6 fm. The flavor, spin, and color wave functions are constructed in two steps. First constructing the wave functions for clusters 1 and 2, then coupling the two wave functions of two clusters to form the wave function for tetraquark system. For the meson-meson structure, as the first step, we give the wave functions of the meson cluster. The flavor wave functions of the meson cluster are shown below.
where the superscript of the χ is the index of the flavor wave function for a meson, and the subscript stands for the isospin I and the third component I z . The spin wave functions of the meson cluster are: . (15) and the color wave function of a meson is: Then, the wave functions for the four-quark system with the meson-meson structure can be obtained by coupling the wave functions of two meson clusters. Every part of wave functions are shown below. The flavor wave functions are: The spin wave functions are: The color wave function is: Finally, multiplying the wave functions ψ L , ψ σ , ψ f , and ψ c according to the definite quantum number of the system, we can acquire the total wave functions of the system.
For the diquark-antidiquark structure, the orbital and the spin wave functions are the same with those of the meson-meson structure. For the flavor wave functions, we give the functions of the diquark and antidiquark clusters firstly. =db, Then, the color wave functions of the diquark clusters and the color wave functions of the antidiquark clusters are: After that, the wave functions for the four-quark system with the diquark-antidiquark structure can be obtained by coupling the wave functions of two clusters. Every part of wave functions are shown below. The flavor wave functions are: The color wave functions are: Finally, we can acquire the total wave functions by substituting the wave functions of the orbital, the spin, the flavor and the color parts into the Eq. (10) according to the given quantum number of the system.

III. THE RESULTS AND DISCUSSIONS
In present work, we investigate tetraquarks with quark components: usdb and udsb in two structures, mesonmeson and diquark-antidiquark. The quantum numbers of the tetraquarks we study here are I = 0, 1, J = 0, 1 and the parity is P = +. The orbital angular momenta are set to zero because we are interested in the ground states. To check whether or not there is any bound state in such tetraquark system, we do a dynamic boundstate calculation. Both the single-channel and channelcoupling calculations are carried out in this work. All the general features of the calculated results are as follows.

A. Tetraquarks composed of usdb
For tetraquarks composed of usdb, the isospin is I = 1. The energies of the states with J = 0, 1 are calculated and the results are listed in Table III and IV. In the tables, the second column gives the index of the wave functions of each channel. The columns headed with E sc and E cc represent the energies of the single-channel and channel-coupling calculation respectively. For mesonmeson structure, there are two additional columns, the column headed with "Channel" denotes the physical contents of the channel and the coulmn headed with E th denotes the theoretical threshold of the channel. From the Table III, we can see that the energies of every single channel approach to the corresponding theoretical threshold. The channel-coupling cannot help too much. Energies are still above the threshold of the lowest channel (B 0 s π + for IJ = 10 and B * 0 s π + for IJ = 11), which indicates that no bound usdb state with meson-meson structure is formed in our quark model calculation.
With regard to the diquark-antidiquark structure, the energies are listed in Table IV. The channels with different flavor-spin-color configurations have different energies and the coupling of them is rather stronger than that of the meson-meson structure. However, the energy of the IJ = 10 state is still higher than the theoreti-cal threshold of the lowest channel B 0 s π + , 5506.9 MeV. Similarly, the energy of the IJ = 11 state is higher than the theoretical threshold of the lowest channel B * 0 s π + , 5532.9 MeV. Thus, there is no bound state with diquarkantidiquark structure in the present calculation.
Nevertheless, the colorful subclusters diquark (qq) and antidiquark (qq) cannot fall apart because of the color confinement, so there may be a resonance state with diquark-antidiquark structure. To check the possibility, we perform an adiabatic calculation of energy for both the IJ = 10 and IJ = 11 states. The results are shown in Fig. 1, where the horizontal axis S is the distance between two subclusters and the vertical axis stands for the energy of the system at the corresponding distance S. It is obvious in Fig. 1 that the energy of both the IJ = 10 and IJ = 11 states is increasing when the two subclusters fall apart, which indicates that the two subclusters tend to clump together. In other words, the odds are the same for the states being meson-meson structure, diquark-antidiquark structure or other structures. As mentioned above, the energy of the state is higher than the theoretical threshold of the lowest channel, so neither the state of IJ = 10 nor the state of IJ = 11 is a resonance state in QDCSM.
Therefore, the X(5568) cannot be explained as a molecular state or a diquark-antidiquark resonance of usdb in the present calculation. Our results are consistent with the analysis of Ref. [19] and Ref. [20]. In Ref [19], the four-quark system usdb with both mesonmeson structure and diquark-antidiquark structure was studied in the framework of the chiral quark model by using the Gaussian expansion method, and no candidate of X(5568) was found. In Ref. [20], Burns and Swanson explored a lot of possible explanations of the X(5568) signal, a tetraquark, a hadronic molecule or a threshold effect and found that none of them can be a candidate of the observed state.

B. Tetraquarks composed of udsb
For tetraquarks composed of udsb, four states with the quantum numbers IJ = 00, 01, 10 and 11 are studied. The energies of the meson-meson structure and the diquark-antidiquark structure are listed in Tables V and VI, respectively. For the meson-meson structure, the results are similar to that of the tetraquarks of usdb. Table V shows that the energies of every single channel are above the corresponding theoretical threshold. The effect of channel-coupling is very small except for the IJ = 01 state. For the states with IJ = 00, IJ = 10, and IJ = 11, all energies are above the threshold of the lowest channel (B 0 K + for IJ = 00, B 0 K + for IJ = 10, and B * 0 K + for IJ = 11) even by the channel-coupling calculation. While for the state with IJ = 01, the energy is about 1.0 MeV lower than the threshold of the lowest channel B * 0 K + after channel-coupling. However, the binding energy is not very large, so there maybe a weak molecular bound state of udsb with quantum numbers of IJ = 01, and the mass of this state is about 5813 MeV.
For the diquark-antidiquark structure, all the possible channels are shown in Table VI. One can see that the en-ergy of each single channel is higher than the theoretical threshold of the corresponding channel, which are shown in Table V. Although the effect of the channel-coupling is much stronger than that of the meson-meson structure, the energy of the IJ = 10 and IJ = 11 states are still above the theoretical threshold of the corresponding channel. So there is no any bound states for the IJ = 10 or IJ = 11 state. In order to check if there is any resonance state, we also perform the adiabatic calculation of energy for both the IJ = 10 and IJ = 11 states. The results are shown in Fig. 2. The case is similar to the tetraquarks composed of usdb. The energy of both the IJ = 10 and IJ = 11 states is increasing when the two subclusters diquark (qq) and antidiquark (qq) fall apart, which indicates that the two subclusters tend to clump together. So there is no resonance state with quantum numbers IJ = 10 and IJ = 11.
However, things are different for the IJ = 00 state and the IJ = 01 state. The energy of the IJ = 00 state is about 5701 MeV, 74 MeV lower than the theoretical threshold of the B 0 K + , which indicates that the IJ = 00 state of the diquark-antidiquark structure can be a bound state. Ref. [24] also found that the bound state with IJ = 00 was possible. Meanwhile, the energy of the IJ = 01 state is 58 MeV lower than the theoretical threshold of the B * 0 K + , so the IJ = 01 state is also bound here. Thus, both the IJ = 00 state and the IJ = 01 state of diquark-antidiquark structure can form bound states.
By comparing with the results of the meson-meson structure, we note that the energy of the IJ = 01 state of the diquark-antidiquark structure is about 5756 MeV, which is much lower than that of the meson-meson structure shown in Table V. This shows that the IJ = 01 state prefers to be a bound state of the diquark-antidiquark structure. Moreover, the IJ = 00 state of the diquarkantidiquark structure is easier to form the bound state than the one of the meson-meson structure. All these indicate that the diquark-antidiquark configuration maybe a good choice for some tetraquarks. Some work have been done to explain the exotic XY Z states depending on the diquark-antidiquark configuration. Ref. [34] proposed the hypothesis that the diquarks and antidiquarks in tetraquarks were separated by a potential barrier to explain the properties of exotic resonances such as X and Z. Ref. [35] presented a dynamical picture to explain the nature of some exotic XY Z states based on a diquark-antidiquark open-string configuration. The picture combined the advantages of diquark-based models, which can accommodate much of the known XY Z spectrum, with the experimental fact that such states are both relatively narrow and are produced promptly. Thus both the IJ = 00 and IJ = 01 states of the diquarkantidiquark structure we obtain here maybe the narrow resonance states. The study of the decay width of these states is our further work.
Contrasting with the tetraquarks composed of usdb, we find the tetraquarks composed of udsb is more likely to form bound state. The reasons are as follows. First, the diquark pair of two light quarks (ud) or two heavier quarks (sb) is usually more stable than the one of two quarks with larger mass difference like us or db pair. Our results show that the tetraquarks composed of udsb of the diquark-antidiquark structure is most possible to form bound states, which just supports this point. Secondly, the lowest threshold of udsb is BK, 270 MeV higher than the threshold B s π of usdb. So there is large mass region for this udsb state below threshold and being stable. This conclusion also supports the assumption of Ref. [23], which proposed such particle with the quark component of udsb (or its charge-conjugated one) as a partner of X(5568) of usdb under the SU (3) flavor symmetry.

IV. SUMMARY
In summary, we investigate tetraquarks composed of usdb and udsb in the framework of QDCSM. Two structures, meson-meson and diquark-antidiquark, are considered. Our results show that there is no bound state or resonance state composed of usdb. The reported X(5568) cannot be explained as a molecular state or a diquark-antidiquark resonance of usdb in present calculation. In contrast, two bound states are obtained for the tetraquarks system composed of udsb: an IJ = 00 state with the mass of 5701 MeV, and an IJ = 01 state with the mass of 5756 MeV, which maybe the better tetraquark candidates with foure different flavors. These two bound states are of the diquarkantidiquark structure. For the system with IJ = 00 and IJ = 01, it is obvious that the state of the diquarkantidiquark structure is more likely to form bound state than that of the meson-meson structure, which indicates that the diquark-antidiquark configuration would be a good choice for the tetraquarks udsb with IJ = 00 and IJ = 01. During the calculation, we find that the effect of the channel-coupling in the diquark-antidiquark structure is much stronger than that in the meson-meson structure, and the channel-coupling plays an important role in forming bound states in the diquark-antidiquark structure.
Meanwhile, our results also show that the tetraquarks composed of udsb is more possible to form bound states than the one composed of usdb. Thus, if the X(5568) does exist, the tetraquarks composed of udsb must be a more stable state. If the X(5568) is proved to be nonexistent, it is still possible for the existence of such tetraquarks with udsb components. This conclusion is in accordance with the point of Ref. [23], which proposed the state composed of udsb (or its charge-conjugated one) as a partner of X(5568) of usdb. Ref. [23] also pointed out that if the lowest-lying udsb state exists below threshold, it can be definitely observed via the weak decay mode J/ψK − K − π + , with the expectation of hundreds of events in the current LHCb data sample but rejecting backgrounds due to its long lifetime. Therefore, the udsb state would be a promising detectable tetraquark state. We hope that experiments will help to discover these interesting tetraquark states.