Tetraquark states in the bottom sector and the status of the $Y_b$(10890) state

We have done the exploratory study of bottom tetraquarks ($[bq\bar b \bar q];{q\in u,d}$) in the diquark-antidiquark framework with the inclusion of spin hyperfine, spin-orbit and tensor components of the one gluon exchange interaction. Our focus here is on the $Y_b$(10890) and other exotic states in the bottom sector. We have predicted some of the bottom counterparts to the charm tetraquark candidates. Our present study shows that if $Z_b(10610)$ and $Z_b(10650)$ are diquark-diantiquark states then they have to be first radial excitations only and we have predicted $Z_b(10650)$ state as first radial excitation of tetraquark state $X_b$ (10.143-10.230). We have identified $X_b$ state with $J^{PC}= 1^{+-}/0^{++}$ as being the analogue of $Z_c(3900)$. An observation of the $X_b$ will provide a deeper insight into the exotic hadron spectroscopy and is helpful to unravel the nature of the states connected by the heavy quark symmetry. We particularly focus on the lowest P wave $[bq][\bar b\bar q]$ states with $J^{PC}=1^{--}$ by computing their leptonic, hadronic and radiative decay widths to predict the status of still controversial $Y_b$(10890) state. Apart from this, we have also shown here the possibility of mixing of P wave states. In the case of mixing of $1^{--}$ state with different spin multiplicities, we found that predicted masses of the mixed P states differ from $Y_b$(10890) state only by $\pm20$ MeV energy difference which can be helpful to resolve further the structure of $Y_b$(10890).


Introduction
A plethora of new kind of states which have been observed recently has inspired extensive interest in revealing the underlying structure of these newly observed states. Exploration of these states will improve our understanding of non-perturbative QCD. In recent years a significant experimental progress has been achieved regarding discoveries of bottomonium-like and charmoniumlike charged manifestly exotic resonances Z b (10610), Z b (10650) [1,2,3,4,5,6,7], Z c (3900) [8,9,10,11,12] and Z c (4020/4025) [13,14,15,16].Their production mechanism and decay rates are not compatible with a standard quarkonium interpretation. A huge effort in understanding the nature of these new states and in building a new spectroscopy is forthcoming.
In the recent years strong experimental evidence from B and charm factories has been accumulating for the existence of exotic new quarkonia states, narrow resonances called X, Y, Z particles which do not seem to have a simple qq structure. Their masses and decay modes show that they contain a heavy quark-antiquark pair, but their quantum numbers are such that they must also contain a light quark-antiquark pair [17]. The theoretical challenge has been to determine the nature of these resonances. Their production mechanism, masses, decay widths, spin-parity assignments and decay modes have been revisited recently [18,19,20]. The term exotica labels states which have an identical number of quarks and antiquarks but defy an ordinary meson classification. Many exotic states in the charm sector with cc content have been discovered by Belle and others [8,21]. While there are most likely many more which are yet unknown and many of them should also reflect in the bb sector according to heavy quark symmetry. The non-discovery of the respective bb partners of the charmonium-like exotica would be even more enigmatic. Belle collaboration has extended the study of the XYZ exotic state family to the bottomonium sector by claiming the observations of two exotica states in Υ (5S) decays [3]. The CMS experiment also searched for the bottomonium partner of X(3872) at hadron colliders [22] in the Υ (1S)ππ decay mode and found no evidence for the X b state while the ratio of the cross section X b to Υ (2S) shows upper limit in the range of (0.9-5.4)% at 95% confidence level for X b masses between 10-11 GeV. Those are the first upper limits on the production of a possible X b state at a hadron collider. Currently there are pending, unanswered questions concerning the exotic spectroscopy in the heavy quark sectors especially in the bottom sector. To promote the endeavor of understanding the heavy exotic states, the exploration of the bottom sector is important. Motivated by the BaBar's discovery of large Y (4260) → π + π + J/ψ signal discovered in the charmonium mass region, Belle experiment have searched for similar state in the bottomonium sector [23]. They observed partial decay widths Υ (5S) → π + π + Υ (nS) (n = 1, 2, 3) associated with the peak in the π + π + Υ (nS) cross section hundreds of times larger than the theoretical predictions [1] and the corresponding measured rates for the Υ (4S) [24]. This observation suggests the presence of a new, non-conventional hadronic state in the bottom sector equivalent of the Y (4260) of the charm sector with mass around 10.890 GeV [26] which is referred as Y b (10890) state. Indeed, there exist three candidates up to date, namely the states labeled Y b (10890), Z b (10610) and Z b (10650), observed by Belle [3]. Not only new states are waiting to be discovered but also the existence of Y b (10890) needs to be established or refuted. The Y b (10890) is a potential exotic state still remains to be confirmed since its observation first reported by the Belle collaboration [1,27]. Looking into the interest in this case, present study is particularly focus on the negative parity 1 −− exotic states. Apart from its spin parity, study of its di-leptonic, hadronic and radiative decay widths also help us to solve the puzzling features of this state. The interpretation of the hidden bottom four quark state as a tetraquark exotic states has been advanced and has been studied in considerable detail [25,26,28,29,30,31,32,33,34,35,36]. The experimental search for tetraquark states is a very difficult problem, since exotic candidates are nothing but the resonances immersed in the excited hadron spectra and moreover they usually decay to several hadrons. Their mass and decay products put them in the category of quarkonia-like resonances but their masses do not fit into the conventional quark model spectrum of quark-antiquark mesons [37,38]. However, to confirm a new resonance it is necessary to study all its properties with high level of accuracy including its mass and width. In this work, we develop phenomenology to study some of the theoretical problems of multiquarks and predict multiquark bound states and resonances. In particular, as a benchmark, we study in detail the heavy-light antilight-antiheavy systems who are expected to produce tetraquarks. Despite the intense experimental attempts, these resonances are still mysterious and complicated and we still lack of a comprehensive theoretical framework. In particular, the most popular phenomenological models proposed to explain the internal structure of these particles are the compact tetraquark in the constituent diquark-antidiquark picture and the loosely bound di-meson molecular picture. Following Gell-Mann ′ s suggestion of the possibility of diquark stucture [39], various authors have introduced effective degrees of freedom of diquarks in order to describe tetraquarks as composed of a constituent diquark and diantiquark using QCD sum rules [40,41]. This concept of diquark was even used to account for some experimental phenomena [42]. The authors of refs [43,44] studied the tetraquark systems in the diquarkantidiquark picture using the chromo-magnetic interactions. In the same way Maiani et al. [45,46,47] also studied tetraquarks and pentaquarks systems by considering this concept of diquark. In their study they have included the spin spin interactions. On the other hand Ebert et. al. [48,25,49] employed the relativistic quark model based on the quasi-potential approach in order to find the mass spectra of hidden heavy tetraquark systems. Unlike Maiani et al., they ignored the spin−spin interactions inside the diquark and anti-diquark. The presence of a coherent diquark structure within tetraquarks helps us to treat the problem of four-body to that of two two-body interactions. In the present case, we employ the diquark and anti-diquark picture in the beauty sector and compute the mass spectra of the diquarkantidiquarks [bqbq]; q ∈ u, d in the ground and orbitally excited states with the inclusion of both S = 0 and S = 1 diquarks. We present the formalism of the study of hidden bottom tetraquark states in section 2. In section 3, we discuss the Y b states and their decay properties. We conclude and discuss our findings in section 4.

Theoretical framework
In this paper we shall take a different path and investigate different ways in which the experimental data can be reproduced. We have treated the four particle system as two-two body systems interacting through effective potential of the same form of the two body interaction 3 potential of Eq.(1). The existence of exotic hadrons of the diquarks-diantiquarks pair called tetraquarks or diquakonia is a problem which was foremost raised about 20 years ago and was used to describe scalar mesons below 1GeV in 1977 by R. Jaffe [50,51]. He suggested the idea of strongly correlated two-quarks-two-antiquarks states to baryon-antibaryon channels where the MIT bag model used to predict the quantum numbers and the masses of prominent states. There are two types of diquarks one is S = 0 good (scalar) diquarks and another one is S = 1 bad (vector) diquarks. We have available lattice results which favors the evidence of an attractive diquark(antidiquark) channel for the good diquarks (color antitriplet, flavour antisymmetric) with spin S = 0 in accordance with Jaffe's proposal. On the other hand there is no lattice results available for an attractive channel for the bad diquarks i.e. with spin S = 1. Here, we use the fact the effective QCDlagrangian is independent of spin in the heavy quark limit and we incorporate the diquark with S = 1 also in computing the mass spectra. There are many methods to estimate the mass of a hadron, among which phenomenological potential model is a fairly reliable one especially for heavy hadrons [52,53,54,55]. In the present study, the non-relativistic interaction potential we have used is the Cornell potential which consists of a central term V (r) which is just the sum of the Coulomb (vector) and linear confining (scalar) parts given by The value of the α s , the running coupling constant is determined by [56] α s (µ 2 ) = 4π where µ = 2m a m b /(m a + m b ), Λ = 0.413 GeV, M B is the background mass and n f is number of flavours [56]. The model parameters we have used in the present study are same as in refs [56,57]. The constituent quark masses employed here are: m u = m d = 0.33GeV and m b = 4.88GeV . The degeneracy of these exotic states are removed by including the spin-dependent part of the usual one gluon exchange potential [58,59,60,61]. The potential description extended to spin-dependent interactions results in three types of interaction terms such as the spin-spin, the spin-orbit and the tensor part. Accordingly, the spin dependent part V SD is given by The coefficient of these spin dependent terms of Eq.(4) can be written in terms of the vector (V V ) and scalar (V S ) parts of the static potential described in Eq.(1) as Where M i , M j correspond to the masses of the respective constituting two-body systems. The Schrödinger equation with the potential given by Eq. (1) is numerically solved using the Mathematica notebook of the Runge-Kutta method [62] to obtain the energy eigen values and the corresponding wave functions.

The four-quark state in diquark-antidiquark picture
In this section, we calculate the mass spectra of tetraquarks with hidden bottom as the bound states of two clusters (Qq andQq), (Q = b; q = u, d). We think of the diquarks as two correlated quarks with no internal spatial excitation. Because a pair of quarks cannot be a color singlet, the diquark can only be found confined into hadrons and used as effective degree of freedom. Heavy light diquarks can be the building blocks of a rich spectrum of exotic states which can not be fitted in the conventional quarkonium assignment. Maiani et al [45] in the framework of the phenomenological constituent quark model considered the masses of hidden/open charm diquark-antidiquark states in terms of the constituent diquark masses with their spin-spin interactions included. We discuss the spectra in the framework of a non-relativistic hamiltonian including chromo-magnetic spin-spin interactions between the quarks (antiquarks) within a diquark(antidiquark. Masses of diquark (antidiquark) states are obtained by numerically solving the Schrödinger equation with the respective two body potential given by Eq.(1) and incorporating the respective spin interactions described by Eq.(4) Table 1 Mass spectra of four quark states in the diquark-antidiquark picture (For L 1 = 0, L 2 = 0)(in GeV)) . perturbatively.
In the diquark-antidiquark structure, the masses of the diquark/diantiquark system are given by: Further, the same procedure is adopted to compute the binding energy of the diquark-antidiquark bound system as Where Q and q represents the heavy quark and light quark respectively. In the present paper, d andd represents diquark and antidiquark respectively. While E d , Ed, E dd are the energy eigen values of the diquark, antidiquark and diquark-antidiquark system respectively. The spin-dependent potential (V SD ) part of the hamiltonian described by Eq.(4) has been treated perturbatively. Details of the computed results are listed in Table 1, 2, 3, 4 and 5 for the low lying positive parity and negative parity states respectively.

Mixing of P wave states
In the limit of heavy quark, the spin of the light and heavy degrees of freedom are separately conserved by the strong interaction. So hadrons containing a heavy quark can be simultaneously assigned the quantum numbers S Qq , m Qq ,Qq, mQ q . Since dynamics depends only on spin of the light degrees of freedom, the hadron will appear in degenerate multiplets of total spin S that can be formed from diquark and antidiquark and accordingly we can classify the states in the convenient way. In the present study, we find that the masses of orbitally excited state with relative angular momentum L = 1 and total spin S = 0, 1, 2 corresponding to 1 P 1 , 3 P 1 , 5 P 1 and 3 P 0 are close to each other in the mass region around 10.850-11.201GeV. The importance of the - Ali [26] Brink [32] Semay [33] Mass (GeV) J PC linear combination of scalar and axial vector states was noted by Rosner [63] and he emphasized in the context of the constituent quark model the individual conservation of heavy and light degree of freedom in the heavy quark systems. In the mass spectra shown in Table 2, two 1 P 1 states with masses 10.931 GeV and 10.882 GeV and another 5 P 1 state with mass 10.853 GeV are there. Similarly, in the mass spectra shown in Table 3, there are two 1 P 1 states with masses 10.201 GeV and 10.145 GeV respectively, 3 P 1 state with mass 11.163 GeV and 5 P 1 state with mass 11.117 GeV. Generally mixing is done through Where, U −1 is given by Table 2 Mass spectra of four quark states in the diquark-antidiquark picture(L 1 = 1, L 2 = 0)(in GeV)) .  Table 3 Mass spectra of four quark states in the diquark-antidiquark picture(L 1 = 1, L 2 = 1)(in GeV)) .   Table 4 1 st radially excited mass spectra of four quark states in the diquark-antidiquark picture(For L 1 = 0, L 2 = 0)(in GeV)) . . So we have shown here the possibility that these states might be getting mixed up with each other to yield mixed states according to [63,64,65] Where, |α and |β are same parity states. The |P J ′ > and |P J > are the lower and higher eigen states respectively as given in Ref. [66]. For a finite mixing angle (or mixing probability 1) the masses of the |P J ′ > and |P J > states will lie only between the masses of the |α > and |β > states. Accordingly, we get mixed states at 10.914 GeV and 10.898 GeV for mixing of two states 1 P 1 (10.931 GeV) and 1 P 1 (10.882 GeV). Similarly for the mixing of 1 P 1 (10.931) and 5 P 1 (10.853), we obtained states at 10.905 GeV and 10.879 GeV and for that for the mixing of 1 P 1 (10.882) and 5 P 1 (10.853) states, we obtained mixed states at 10.871GeV and 10.862 GeV. In the same way we obtained mixed states for other combinations also. These mixed states are listed in Table 6. The masses of 1 −− mixed states lie very close to the 10.890 resonance. We have also computed the leptonic and hadronic and radiative decay widths for these mixed states.

Y b (10890) state and its decay properties
The prominent exotic state Y b (10890) with J P C = 1 −− was first observed by the Belle collaboration [1,4] and to date, it remains to be confirmed by independent experiments. The anomalously large production cross sections for e + e − → Υ (1S; 2S; 3S)π + π − measured at Υ (5S) was not in good agreement with the line shape and production rates for the conventional bb Υ (5S)) state. An important issue is whether the puzzling events seen by Belle stem from the decays of the Υ (5S) or from another particle Y b having a mass close enough to the mass of the Υ (5S). This results motivated theorists to resolve the puzzling features of this peak which lies approximately at mass 10.890 GeV. Currently, there are two competing theoretical explanations: the tetraquark interpretation on the one hand [26,29,30,31] and the re-scattering model [67] on the other. The tetraquark model can explain the enhancement and the resonant structure via Zweig allowed decay processes and coupling to intermediate states while the re-scattering model is based on the decay Υ (5S)) → B * B * and a subsequent recombination of the B mesons. A detailed study on this state is available in literature [26,29,28]. The state Υ (10890) is usually referred as the Υ (5S), since its mass is close to the mass of the 5S state predicted by potential models. However, a different proposal has been put forwarded by the authors of Ref. [26], in which they call this state as Y b and this state is being a P-wave tetraquark analogous to the Y (4260), though the current experimental situation regarding the peak about 10.890 GeV is still debatable. However, here we found three vector states with J P C = 1 −− whose mass is around 10.890 GeV i.e. 10.882GeV, 10.853 GeV and 10.931 GeV. We also found 1 −− mixed 1 P 1 states lie at 10.914 and 10.898 GeV. To resolve the Y b further we have calculated the di-electronic, hadronic and radiative decay widths of these states.

Leptonic decay width of J
In the conventional bb systems, the decay widths are determined by the wave functions at the origin for ground state while for the P waves the derivation of these wave functions at the origin are used. We have used the same Van-Royen-Weisskopf formula but with a slight modification. Since tetraquark size is larger than that of quarkonia,to take into account larger size of tetraquark, we have modified wavefunctions by including a quantity σ, size parameter whose value varies from σ ∈ [ 1 2 , . These tetraquarks wave functions will affect the decay amplitudes and thereby influencing the decay rates. The partial electronic decay widths Γ ee[bu] and Γ ee[bd] of the tetraquark states Y bu and Y bd made up of diqaurks and antidiquarks (for up quark and down quark respectively) are given by the well known Van Royen Weiss-kopf formula for P waves [69] Here, α is the fine structure coupling constant and σ < 1 and < e Q > is the effective charge of Qq diquark system given by [70] For computing the leptonic decay width, we have employed the numerically obtained radial solutions while the authors [69,26] have used value calculated by using QQ-onia package [71] giving |R ′ 11 (0)| 2 = 2.067GeV 5 . Our calculated results for leptonic decay widths for Y bu and Y bd are shown in Table 7 with available theoretical data. Since all the vector 1 −− states are P-waves, the value of R ′ (0) will not change as the masses of the 8 diquarks remain the same. Hence, value of leptonic decay width does not change significantly as it only varies with the mass. However in the case of mixed states the contributions from the radial wavefnctions will be noticeable.

Hadronic decay width of
In this section, we have studied the hadronic decay of the 1 −− P wave Y b (10890) state. We discuss the two body hadronic decays i.e. Y b (q) → B * q (k)B * q (l).These are zweig allowed processes and involve essentially the quark rearrangements. For calculating domiant two body hadronic decay widths of the 1 −− Y b (10890) state, the vertices are given as [72] and the corresponding decay widths are given by Here | − → k | is the center of mass momentum given by Where, M is the mass of the decaying particle and M k , M l are the masses of the decay products. The decay constant F is the non-perturbative quantity and to evaluate it is the beyond the scope in our approximation. We adopted the same approach used in [26,73]and estimate them using the known two body decays of Υ (5S) which are described by the same vertices as given in [72]. To extract the value of F and | − → k |, we have used the values of decay widths for the decays Υ (5S) → B q (k)B q (l), B q (k)B * q (l), B * q (k)B * q (l) from Particle data group [74]. The extracted value of the F and | − → k | are shown in the Table 8 along with the decay width results. To take into account different hadronic size of tetraquark we have included a quantity σ that already discussed earlier. The results for hadronic decay widths are shown in table 5 which differ from the corresponding PDG [74] values of Υ (5S). Out of these three P wave states, computed value of hadronic decay width for Y b (10853) is of the order of 50 MeV as against the PDG value of 110 ± 13 MeV and consistent with the BELLE measurements. For other two states we are getting more higher values than they actually should have. So out of these three states, we predict only the state with mass 10.853 GeV as a Y b (10890) state.

Radiative decay width of
We study the radiative decays of these states using the idea of Vector Meson Dominance (VMD) which describe interactions between photons and hadronic matter [75] and we hope that this will increase an insight about these tetraquark states. The transition matrix element for radiative decay of Y b → χ b + γ is given with the use of VMD and the decay width is given by Where, λ is the center of mass momentum and f ρ = 0.152GeV 2 [76]. Similarly, we have computed radiative decay Y b → η b + γ. The present results are shown in the Table 9 with available theoretical data. There is no experimental data available for the radiative decay of Y b (10890) and we look forward to see the experimental support in favour of our predictions.

Results and Discussions
We have computed the mass spectra of hidden bottom four quark states in diquark-antidiquark picture which are listed in table 1, 2, 3, 4 and 5. We have taken various combinations of the orbital and spin excitations to compute the mass spectra. The computed mass spectra are compared with other available theoretical results in figure 1. Apart from this we mainly have paid attention to Y b (10890) state and have computed leptonic, hadronic and radiative decay width of Y b which are listed in tables 7, 8 and 9 respectively. Apart from this, we have also done mixing of 1 −− P waves which are also listed in the respective Tables. The core of the present study is that the color diquark is handled as a constituent building block. We predicted some of the bottom tetraquark states as counterpart in the charm sector. It is necessary to highlight that the observation of the bottom counterparts to the new anomalous charmonium-like states Table 6 Mixed P wave states(in GeV)) .   [26] 0.08 ± 0.03 [26] 0.12 [73] is very important since it will allow to distinguish between different theoretical descriptions of these states. In this viewpoint, it would be also valuable to look for the analogue in the bottom sector as states related by heavy quark symmetry may have universal behaviours. The predicted bottom counterparts are shown in the Fig. 2 for better understanding. In the present study, we have noticed that mass difference between predicted X b (10233) and χ b1 (9892) which of the same order of magnitude of the mass difference between X(3872) and χ c1 (3510) of the charm sector This kind of similarity between charm and bottom sector is very interesting. We found that mass difference between X b (10143) and its first radially excited state X b (10650) states is ∼ 510M eV as similar to charmonia which is about 590 MeV. In the same way, we found that mass difference between X b (10233) and its first radially excited state X b (10683) states is ∼ 450M eV . So by taking the the evidence from these results, we can say that 4-quark state in the bottom sector analogous to charm sector should exist. We have predicted some of the radially excited states which are listed in Table  10. Accordingly, we predicted Z b (10650) state as the first radial excitation of either X b (10143) (0 ++ ) state or X b (10233)(1 +− ) state. The authors of Ref. [77] studied the masses of the S-wave [bq][bq] tetraquark states with the inclusion of chromomagnetic interaction and  [41]. These predictions of X b state and its production rates in hadron-hadron collisions have indicated a promising prospect to find the X b at hadron collider in particular the LHC and we suggest our experimental colleagues to perform an analysis. Such attempt will likely lead to the discovery of the Xb and thus enrich the list of exotic hadron states in the heavy bottom sector. An observation of the X b will provide a deeper insight into the exotic hadron spectroscopy and is helpful to unravel the nature of the states connected by the heavy quark symmetry. Similarly, there exist other radial excited states in the region 11.095-11.151 GeV corresponding to 2P states. We look forward to see experimental search for these states. The authors in Refs [79,80] have predicted Z b (10650) state as di-mesonic molecular state in the ground state. From our present study, we suggest that if Z b states are diquark -diantiquark states then they are not the ground state of bottomonium-like four quark state but the first radially excited state of its ground state which lies in the 10.100 − 10.300GeV which is in agreement with the results reported by the authors of Ref [81]. The same presumption was made by authors of Table 9 Radiative decay widths (in keV)) . Ref [46,54] to explain Z(4430) state as an excitation of state Z c (3900)/Z c (3885) in the charm sector. So in conjecture with this, our prediction regarding Z b (10650) state is just a straightforward extension to beauty sector and we observe that the Z b (10650) is also a radially excited state of still unmeasured X b state just like that of authors of Ref [81] who predicted Z b (10610) state as the the radial excitation of X b (10100) such that the mass difference is M Z b (10650) − M X b (10143) ∼ 510M eV which is very close to mass difference between Υ (2S) − Υ (1S) = 560M eV . To have a clear-cut picture about the discussion made regarding the bottom exotic states, the above discussed exotic states are displayed in the Fig. 2 with analogous states at the charm sector.
The comparison between the bottom tetraquark states and charm tetraquark states accentuates the resemblance between presumptions made in the present study, namely the existence of a X b (10143) as a ground state of Z b (10650) and presumption related to existence of ground state of Z(4430) made in Refs [46,54]. The pre-sumption of bottom tetraquark states analogous to charm spectra should stimulate searches for these states in both the beauty and the charm sector within the mass range around 10100 − 10300M eV and 3500 − 3870M eV respectively. The searching of these states would be not only able to find unobserved state shown in Fig. 1 but also be able to detect many more prominent states in these mass range. As Y b (10890) state with quantum number 1 −− is of our keen of interest, in this study we have predicted three P wave 1 −− states in the mass region around 10.850-10.931 GeV. We have observe that P wave state with mass 10.853 GeV as the Y b state. The calculated partial electronic decay widths for P wave Y b is about 0.03-0.12 keV which is in agreement with the available experiment data [6] and other theoretical predictions [26,73]. Our present calculation show that the leptonic width of Y b is much lower than that of the width of conventional state Υ (5S)(0.31 ± 0.07keV ) [74]. From this we can say that Υ (10890) peak is different from the Υ (5S) and possibly may be the Y b (10890) only.
We have also computed the two body hadronic decays of Y b . The total hadronic decay width is of the order of 50 MeV which is lower than the total decay width of Υ (5S) = 110M eV state. So this narrow width state Y b (10890) can be tetraqurk state only rather than being the conventional bb state. We have also computed the radiative decay widths of Y b , but due to lack of experimental results we can not make any concrete conclusion here. These results can be guidelines for future studies. In the absence of experimental data, we can't make any conclusion regarding mixing of P wave states but we expect our results could be helpful to understand the structure of these states. Our computed masses of 1 −− mixed states i.e. 1 P 1 and 5 P 1 states lie very close to Y (10890) state by at most an order of ±20 MeV. So we look forward to see the experimental search for these states with very high precision as these states are very closely spaced. The experiments should have in principle the sensitivity to detect and also to explore the nature of such near-lying states. The present study of mixing is an attempt to signify its importance to further resolve mystery of Y b (10890). If the status of Y b (10890) is confirmed then it will be a major step in the direction of testing the models and provide theorists with vital input to present a credible explanation of this new form of hadrons.