Masses and magnetic moments of exotic fully heavy pentaquarks

Inspired by the observation of a resonant state X(6600) of fully charm tetraquark by the CMS experiment of LHCb Collaboration in double J/ψ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J/\psi $$\end{document} decay channel, we perform a systematical study of all configurations of fully heavy pentaquarks PQ1Q2Q3Q4Q5¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{Q_{1}Q_{2}Q_{3}Q_{{4}}\bar{Q_{5}}}$$\end{document} (Qi=c,b\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Q_{i}=c,b$$\end{document}, i=1,2,3,4,5\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,3,4,5$$\end{document}) in their ground states in unified framework of MIT bag model. The color-spin wavefunctions of pentaquarks, classified via Young tableau and presented in terms of the Young–Yamanouchi bases, are used to compute masses and magnetic moments of fully heavy pentaquarks via numerical variational method, predicting a set of masses ranging from 8.229 GeV for the Pccccc¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{cccc\bar{c}}$$\end{document} to 24.770 GeV for the Pbbbbb¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{bbbb\bar{b}}$$\end{document}. Combining with computed masses of fully heavy mesons and baryons, we find that masses of fully heavy hadrons(mesons, baryons, tetraquarks and pentaquarks) with identical flavor rise almost linearly with the number of valence quarks in hadrons, being consistent with the heavy quark symmetry in the heavy-quark limit.


I. INTRODUCTION
Possible existence of multiquark hadrons like tetraquarks (q 2 q2 ) and pentaquarks (q 4 q) was suggested earlier at the birth of quark model [1,2].Later in the 1970s, multiquark states are studied by Jaffe via the MIT bag model [3,4].Since observation of the X(3872) [5] in 2003 by the Belle, many candidates of tetraquarks have been observed, some of which, such as the Z c (3900) [6] and the T cc (3875), are undoubtedly exotic.In 2020, a first candidate of fully charm tetraquark, the X(6900)/the X(6600), has been observed by LHCb in the di-J/ψ invariant mass spectrum and later confirmed by CMS and ATLAS of the LHC at CERN [7][8][9].For experimental searches of light-flavor pentaquark hadrons, it has been a long and nontrivial history, with no undisputed candidates found in over 50 years.Convincing evidence for pentaquarklike structures P c (4450) + and P c (4380) + with a minimal quark content of uudcc was reported by LHCb in a study of Λ 0 b → J/ψpK − (J/ψ → µ + µ − ) decays [10] in 2015, for which the former peak(state) P c (4450) + was resolved further into two states P c (4440) + and P c (4457) + (over 5.4σ) in 2019.More recently, LHCb collaboration reported evidences of two new charmonium pentaquarks with strangeness, P cs (4459) [11], in the J/ψΛ distribution(in Ξ − b → J/ψΛK − decays) in 2021 and P Λ ψs (4338) [12], in flavour-untagged B − → J/ψΛ p decays in 2022.
The purpose of this work is to systematically study all possible ground-states of fully heavy pentaquarks, 3,4,5) in unified framework of MIT bag model.For this, we construct all color-spin wavefunctions of pentaquarks based on S U(2) s × S U(3) c group, classi-fied via Young tableau and expressed in terms of the Young-Yamanouchi bases, to compute masses and magnetic moments of fully heavy pentaquarks in ground states and thereby predict their masses.With the help of numerical variation upon the bag radius, the chromomagnetic mixing among the color-spin configurations are explored in details.Computing further spin-independent masses of the mesons, baryons and tetraquarks all made of identical heavy flavors, we find a linear dependence of the hadron mass upon the heavy quark number N in hadrons.This linear rise of hadron masses is notable for the bottom sector and is demonstrated analytically via a variational method of bag model, being consistent with heavy quark effective theory (HQET) in heavy quark limit.
After introduction, the wavefunctions of the fully heavy pentaquarks with various configurations are classified in Sect.II by the Young tableau and expressed in terms of the Young-Yamanouchi bases.In Sect.III, we briefly review basic relations of MIT bag model including magnetic moments for each spin state, and describe how to apply the variational method to solve the bag model for the heavy hadrons.In Sect.IV, we perform numerical analysis of mass and magnetic moment for fully heavy pentaquarks.Similar computation are given to fully heavy mesons and baryons to show a linear relation for these fully heavy hadrons in Sect.V.The paper ends with summary and conclusions in Sect.VI.

Fully heavy pentaquarks
for which the corresponding parts of the wavefunctions have to be antisymmetric due to Pauli principle and colorlessness of . Consider the pentaquark P cccb b, for instance.Since the component ccc in P cccb b are identical in flavor, the antisymmetry requirement for a fermionic hadron (color singlet) implies that the wavefunction of P cccb b has to be antisymmetric in orbital, color and spin under exchange of each cc pair among the component ccc.To respect the overall symmetry of the hadronic wavefunctions mentioned above, we utilize the Young tableau to represent the irreducible representations (bases) of the permutation group and thereby classify the pentaquark configurations with certain symmetry.For a given classification, we use the Young-Yamanouchi basis, which corresponds to the Young tableau, to construct explicitly the configuration describing the pentaquark we consider.
First of all, we consider color wavefunctions of the pentaquark , which are singlets in color space.Utilizing direct product of the fundamental (color) representations [3] c and [ 3] c , one can classify, in the language of group theory, the color wavefunctions of a pentaquark as below: in which the allowed color singlet is {3 c ⊗ 3c } only.
One is then forced to consider the color triplet in the direct product, [ [3] . ( Here, the subscript labels the irreducible representation of the color group S U(3) c .This yields three color-singlet configurations of if one combines each of color-triplets (3 c ) listed in Eq. ( 2) with the antitriplet ( 3c ) of the remaining antiquark Q5 .We express the Young tableau representing the obtained colorless configurations as and correspond them to the following color wavefunctions explicitly: In the spin space, the direct product of five fermions represented in terms of Young tableau can be written as For the spin J = 5/2, 3/2 and 1/2 of the pentaquarks, the configurations can be similarly represented in terms of Young tableau [5], [4,1], and [3,2] of one, four, and five dimensions.For the pentaquark with J = 5/2, we write the spin wavefunction as, In the case of the J = 3/2, the spin wavefunctions become We see from the above analysis that to describe fully heavy pentaquarks with one antiquark Q5 having a given flavor and color, the antiquark can be temporarily omitted from the spin degrees of freedom.Thus, we shall identify the spin states shown in Eqs. ( 8) through (10) via the Young-Yamanouchi bases associated with the Young tableau [4], [3,1], and [2,2].Given that the color and spin states represented in terms of the Young tableau respect certain symmetry, one can construct the combined color-spin state of a pentaquark which is fully antisymmetric under the exchange of any pair among the heavy quarks 1 = Q 1 , 2 = Q 2 , 3 = Q 3 , and 4 = Q 4 (for short).
We start from the color singlets in Eq. ( 3), and combine them with spin states by the outer product of the permutation group, S 4 , resulting in the color ⊗ spin states for the quarks 1, 2, 3 and 4.After isolating spin of the antiquark, one can deduce the outer product between Young tableaus [2,1,1] of the color singlets and that of the spin states [4], [3,1], [2,2].For the state expression in terms of the first four quarks (1234) mentioned above, one can obtain the Young tableau reps. of color ⊗ spin states and write them as that in Ref. [51] (see Eq. ( 17)) According to the method described in Ref. [52], all the possible Young-Yamanouchi bases (i.e.ψ = ϕ P 2 χ P 1 ) for the pentaquarks we shall address in this work can be obtained from the couplings(⊗) of the color and spin degrees of freedom.Collectively, we express these color-spin coupled bases, again, in terms of the Young tableau: : in which (for J = 5/2) and (for J = 3/2) In the case of J = 1/2 , the color ⊗ spin bases in Eq. ( 11) are explicitly, Based on the Young tableau in Eq. ( 11) and Pauli principle for each pair of two fermions, we can construct all allowed color-spin wavefunctions and list them in Table I for pentaquarks consisting of four heavy quarks and a heavy antiquark.Based on similar principle, one can construct the states involved in chromomagnetic interactions (CMI) (18) with the help of color and spin factors (19,20) via the color bases (3) and spin bases (8)(9)(10), respectively.Following Ref. [53], we shall use the scaling ratios of color factors to evaluate the matrices of the binding energy, as discussed in section III.TABLE I: Color-spin bases of pentaquarks q 1 q 2 q 3 q 4 q5 with J P quantum number.

III. MIT BAG MODEL AND CMI
MIT bag model describes hadron as a sphere bag containing confined valence quarks.The model also includes perturbative inter-quark interactions via considering the lowest-order gluon exchange among quarks, which is referred as chromomagnetic interaction [54,55].Mass of hadrons for a bag of radius R is [54,55], where the first term is total sum of the kinematic energy of relativistic quark i with mass m i , the second is bag volume energy with B the bag constant, the third is zero point energy with coefficient Z 0 , the forth is binding energy (M BD ) between heavy quarks or between heavy and strange quarks [56,57], and the fifth is chromomagnetic interaction.In Eq. ( 15), the bag radius R is bag radius to be determined variationally, and the quark momentum x i in unit of R −1 satisfies a boundary condition on bag surface, which is to be solved iteratively in this work.The interaction in this work enters among quarks via two parts, ∆M = M BD + M CMI .The first part, the binding energy M BD = M<N B MN (M, N = s, c and b), proposed in Refs.[56,57], rises mainly from short-range chromoelectric interactions between heavy quarks or within heavystrange pairs as they are (relatively) massive and move nonrelativistically.This part was supported by our previous work [53], which reconciles bag dynamics of the light and heavy hadrons.Five binding energies B cs , B cc , B bs , B bc and B bb in the color 3c rep., extracted from the bag mass corrections to the baryons when heavy pairs (cs, cc, bs, bc and bb) involved, can be scaled to the binding energies B MN in 6 c reps. via scaling factor ratios between two reps.[53,56].As there are two reps.( 3c and 6 c ) for color wavefunctions of quark pairs, as discussed in section II, the binding energies could be matrices in color-spin space given in Table I, where spin bases are orthogonal.
The second part of interaction in Eq. ( 15) is chromomagnetic interaction M CMI , which can be due to the short-range gluon exchange [58].Similar to the interaction of magnetic moments of the quark spins, this part has the form of where i and j denote the quark (anti-quark) indices, λ the Gell-Mann matrices, σ the Pauli matrices, and C i j the coupling parameters of the CMI.
To compute the color and spin factor in Eq. ( 18), we employ the following formula of the matrix element: where n, m and x, y indicate components of the basis vectors of color and spin wavefunctions of a hadron, respectively.The symbols c and χ denote bases of color and spin vectors associated with a quark, respectively.The matrices of two factors can be calculated with the help of Eqs. ( 19)-( 20) once colorspin wavefunctions ψ = ϕ P χ P were determined.
In the MIT bag model, the CMI parameter C i j is known analytically, given by [54] where the reduced magnetic moment μi and the (running) strong coupling α s are given by μi and Here, α i = ω i R and λ i = m i R [54].The function of α s (R) in Eq. ( 23) is discussed in Ref. [53].The function F(x i , x j ) in Eq. ( 24) is rational in terms of parameters x i and x j , where From Eq. ( 22), one can write the magnetic moment µ i of the quark (anti-quark) i as, in bag model, with electric charge Q i .As a result, the magnetic moment of a hadron with color-spin wavefunction ψ is, where g i = 2, and S iz is the third component of spin of the individual quark i [59].We shall calculate the magnetic moments of the pentaquarks in unit of the magnetic moment of the proton µ p and further transform the obtained results into that in unit of µ N with the help of measured proton magnetic moment µ p = 2.79285µ N [60,61].Note that Eq. ( 28) also holds true in the case of chromomagnetic mixing, for which one can expand ψ in terms of the color-spin bases to find magnetic moment by Eq. (28).For the spin wavefunctions in section II, one can derive each diagonal element of magnetic moment.The results are listed in Table II in detail.
For the model parameters, we use the values in our previous work [53] which reconcile dynamics of the light and heavy hadrons.These values are Z 0 , the constant B and masses of the quarks (the nonstrange quarks n = u, d, the strange quark s, the charm quark c and the bottom quark b) in this work, TABLE II: Sum rule of magnetic moments for pentaquarks q 1 q 2 q 3 q 4 q5 .

Spin basis µ χ
Meanwhile, using the binding energies B QQ ′ for the quark pair QQ ′ in color antitriplet rep.[53], 3c : one can obtain the corresponding values of the binding energies for color-6 c rep., via multiplying the color-factor ratios (scaling factors).The results of the binding energies between quark pairs (denoted in subscripts) are: Given the parameter inputs in Eqs. ( 29), (30) and (31), one can apply variational method to Eq. ( 15) to determine the bag radius R and the respective momentum x i via Eq.( 17) for a given hadron with color-spin wavefunction ψ, given in Table I.It is then straightforward to calculate masses and magnetic moments of fully heavy pentaquarks we address in this work, as detailed in section IV.

IV. MASSES AND MAGNETIC MOMENTS OF FULLY HEAVY HADRONS
Before discussing fully heavy pentaquarks, we first consider triply (fully) heavy baryons, the fully heavy binding systems of three quarks, which include the five baryons Ω ccc , Ω ccb , Ω * ccb , Ω cbb , Ω * cbb and Ω bbb , as listed in Table III.For these baryons, we apply the same framework of MIT bag model (15) with the inputs in Eqs. ( 29), ( 30) and ( 31) to compute the mass M and the bag radii R 0 of the fully heavy baryons, with other works cited for comparison.Given the color-spin bases [53] (see Eq. ( 15)), defined as ϕ B χ B 1 and ϕ B χ B 2 for J P = 3/2 + and 1/2 + , respectively, the CMI matrices and binding energies are well determined.When pentaquarks included, the computed results could help us further understand the general features of fully heavy systems, as we shall discuss in Sect.V.
In Table III, the second and third rows present bag radii R 0 and masses of triplet heavy baryons, and the data followed rows present the masses (in unit of MeV) cited from various works.One sees that the predicted masses of the Ω ccc , for instance, ranges from 4520 MeV to 5130 MeV for which various methods including QCD sum rules [62,63,71,78] are employed.Up to 610 MeV, our prediction is close to the prediction 4810 ± 100 MeV in Ref. [71].The lattice QCD was also applied to study the Ω ccc in Ref. [67], with the prediction of 4796 ± 8 ± 18 MeV.Comparing with the bag model computation of the doubly heavy baryons and tetraquarks [53] and the fully charm tetraquarks cccc in Ref. [80], our results for fully heavy baryons and pentaquarks are reasonable, considering that the inputs of the parameters are fitted only from light baryons and heavy mesons.
Using the method described in sections II and III, one can construct a total mass formula (15) for the fully heavy hadrons containing the CMI and binding energies and calculate masses, magnetic moments and eigenvectors for all fully heavy systems containing the fully heavy mesons, baryons and pentaquarks, especially of the ten pentaquarks P bbbb b, P bbbbc , P cccc b, P ccccc , P bbbc b, P bbbcc , P cccb b, P cccbc , P ccbb b and P ccbbc , whose color-spin bases have been listed in Table I.Among these heavy hadrons, more attention is given to fully bottom and charm systems, which attracts extensive attention [81][82][83][84][85][86].
In Table IV, we present our computed results for the masses and magnetic moments of the fully bottom system P bbbb b and the fully charm system P ccccc , compared to other works via the CMI model [81] and constituent quark model [82], which adopt the same color-spin wavefunctions, and via the lattice-QCD inspired model [83] and the chiral quark model [84].The comparison is also given with the predictions via the QCD sum rules [85,86].Agreement of our predictions with that by the constituent quark model [82] is achieved.Given the masses of triply heavy baryons computed in Table III and that of the measured heavy mesons [60], one finds that our predicted masses of fully heavy pentaquarks are all above the thresholds of the heavy baryons and mesons, as shown in FIG. 1.
Further more, we examine the respective mass splittings ∆M between spin multiplets due to the lowest two quan-  tum numbers J for the fully heavy mesons, tetraquarks and pentaquarks and show them collectively in Table V.Here, J = {0, 1} or J = {1/2, 3/2}.For the fully heavy systems with identical quarks, it is found that ∆M are suppressed significantly with the heavy quark number N increasing and hadron becoming heavier.We see that the mass splittings are narrower for the bottom sector compared to the charm sector.
For the fully heavy pentaquarks with both bottom and charm flavors, we list the corresponding results computed via the MIT bag model in Tables VI, VII, VIII, and IX, in comparison with other similar works cited [81,82].The obtained masses are also compared to the respective thresholds and plotted in FIG. 2. All hadrons computed are above the thresholds associated with the baryons and mesons, meaning that these fully heavy pentaquarks are unstable against strong decays to two-hadron final states.
The instability of pentaquarks can be investigated by analyzing the eigenvectors obtained in this study.In order to examine the color-spin wavefunction ψ of pentaquark state, we utilize the corresponding eigenvectors and bases defined in Table I, and express it in the form as shown below: Here, c 1 represents the overlap between the wavefunctions of the pentaquark and the specific baryon ⊗ meson component, which corresponds to a scattering state.If the pentaquark state couples strongly to a scattering state, indicating that the probability |c 1 | 2 of the collapse of ψ into a color-singlet (6) approaches 1, we can infer that it is unstable against strong decay with a very broad width.
For the pentaquark states presented in Table IV and VI, the values of |c 1 | 2 are determined to be 1/3.Similarly, in Table IX, all the calculated values are approximately 1/3.Notably, the pentaquark states denoted by asterisks in Table VII and VIII have dominant components of scattering states with |c 1 | 2 values greater than 0.83, while the remaining states range from 0.34 to 0.80.Based on the above discussions, we conclude that the pentaquark states P bbbb Q, P cccc Q, and P ccbb Q (Q = c, b) are likely to possess a compact structure with relatively narrow widths.
On the other hand, the bag radius R 0 ranging from 3.48 GeV −1 = 0.7 fm to 5.08 GeV −1 = 1.0 fm, which are in the order of the typical hadrons in size (about 1 fm), implies TABLE VII: Predicted spectra of pentaquarks P bbbc b and P bbbcc with comparisons from various works.Bag radius R 0 is in GeV −1 .Masses are in GeV.Magnetic moments µ are in unit of µ N .The states denoted by asterisks couples strongly to scattering states.

V. MASS PATTERN AND LINEARITY UPON QUARK NUMBER
To explain the mass pattern obtained in this work, we consider question as to if the bag radius R leads to the larger volume energy M V = 4πR 3 B/3 for the larger heavy quark number N, so that the hadron masses tend to be of suprathreshold.For this, we ignore the mass splittings due to the CMI and rewrite     Eq. ( 15) in the form that hadrons consist of the identical flavor where M stands for the spin-independent mass, N the total number of valence heavy quarks, and B QQ the binding energy given in Eq. (30).
Owing to the scaling of the color factor, as mentioned in Sect.III, the fully heavy hadrons with N identical quarks happen to have total binding energy of NB QQ .Given Eq. (33), one can use it to examine each part of the energy for the fully heavy systems with identical quark number N = 2 ∼ 5.The numerical results are listed in Table X.One sees that the volume energy M V grows quickly while the suppression of zero point energy M Z = −Z 0 /R becomes weaker significantly with more valence quarks involved.This implies that fully heavy hadrons tend to be heavier when N becomes larger and they are all above thresholds of the decaying final states.TABLE X: Numerical results for the mass components (all in GeV), the kinematic energy ω Q (in GeV), the volume energy M V , zero point energy M Z and spin-independent mass M of fully heavy systems (mesons, baryons and pentaquarks) in detail.The parameters of the bag radius R 0 (GeV −1 ) and x Q are also shown for the heavy quark Q.We further examine this mass pattern for fully heavy systems of the mesons, baryons, tetraquarks and pentaquarks within the same framework of MIT bag model.Beside the masses of baryons and pentaquarks which are computed via the MIT bag model, we cite the results for the mesons and tetraquarks in Refs.[53] and [80], respectively.The obtained results are plotted in FIG. 3.For the fully bottom and charm hadrons one sees that hadron masses rise almost linearly with N ranging from 2 to 5. It is notable that our mass predictions of the mesons and tetraquark cccc are in good agreements with experimental data, especially, the computed mass 6572 MeV of the fully charm tetraquark cccc agrees well with the measured mass 6552 ± 10 ± 12 MeV of the newly-discovered resonance X(6600) [8].

System
Inspired by Eq. ( 33) and the numerical results above, we promote the heavy quark number N to be as large as 20 mathematically, despite that a hadron with 20 valence quarks, if exist, may result in a very large and thereby unphysical bag radius R. We apply this procedure to N = 2, 3, 4 and 5 to test if the linear dependence of the hadron mass holds true as far as this work is involved.Similar to FIG. 3, we plot the numer-ical results in FIG.4(a) for the spin-independent masses v.s.quark number N for fully bottom and charm systems respectively, with the data from N = 2 ∼ 20.From these results, we find the following numerical fits for the fully bottom (denoted by subscript b) and charm (denoted by subscript c) systems (in GeV, 2 ≤ N ≤ 5): In FIG.4(a), one observes the linear lines numerically, corresponding to Eq. ( 34), but the marching of it with the data of the spin-independent masses in Table X is approximated.To see this, we apply Eq. ( 33) to the fully light hadron systems with binding energy ignored, B QQ = 0, and obtain the results plotted in FIG.4(b).The linear dependence breaks slightly when N is small, as it should be for the light hadrons.Notice that the (solid) line for the bottom sector is nearly linear in contrast with that for the charm sector, this implies that linearity of mass dependence holds true at heavy quark limit.
In the following, we show such a mass linearity upon N in heavy quark limit analytically.We expand the mass formula (33) in term of the inverse heavy-quark mass m and minimize the mean mass Eq. ( 35) to find the bag radius, The solution of Eq. ( 36) in the large-N limit is R 0 ∼ N 1/5 .This gives rise to the following mass formula for fully heavy hadrons with ΛN = NB QQ + (4πB/3) N 3/5 − Z 0 N −1/5 , by analogy with the mass relation m H Q = m Q + Λ + ∆m 2 2m Q in heavy quark effective theory (HQET).This explains the linear mass dependence of fully heavy hadrons with heavy quark flavor N, as demonstrated in FIG. 3.

VI. CONCLUSIONS AND REMARKS
In this work, we employ a unified framework of MIT bag model to perform a systematical study of all ground-state configurations of fully heavy pentaquarks and other fully heavy hadrons.With the help of the color-spin wavefunctions of pentaquarks, which are classified via Young tableaux which respect symmetry of the S U(2) s ⊗ S U(3) c group and present in terms of the Young-Yamanouchi bases, we have computed  [53], and that of tetraquarks are from Ref. [80].
the masses and magnetic moments of the fully charm, bottom and bottom-charm pentaquarks in ground states via numerical variational method that applies to the MIT bag model.Our computation predicts a set of masses of fully heavy pentaquarks, ranging from 8.229 GeV for the P ccccc to 24.770 GeV for the P bbbb b.Combining with the similar bag-model computation of masses of fully heavy systems, we find that the masses of fully heavy hadrons (mesons, baryons, tetraquarks and pentaquarks) with identical flavors depend almost linearly upon the number N of the valence quarks in hadrons, and we demonstrate this linear behavior of the hadron masses using the MIT bag model, being consistent with the heavy quark symmetry at the heavy-quark limit.Our results also indicate that the heavier the fully heavy hadrons, the narrower the mass splittings of their lowest two multiplets (in J).
We remark that with the increase of the heavy quark num- ber N, the bag radius R 0 of the fully heavy hadrons can rise up as high as 1 fm, and it results in the larger volume energy M V and the less suppression of zero point energy M Z , which makes them suprathreshold and unstable against strong decay.On the other hand, the shorter bag radius of the compact pentaquarks implies some larger possibilities for their formation via hadronization in experiments, compared to the formation of the hadrons containing light degrees of freedom.We hope that our predictions can help the future LHCb experiments to search the fully heavy pentaquarks or baryons addressed in this work.

FIG. 3 :
FIG.3: Masses of fully heavy systems ranging from mesons (N = 2) to pentaquarks (N = 5) with N the number of quarks and antiquarks, calculated in the unified framework of MIT bag model.The results of mesons are quoted from Ref.[53], and that of tetraquarks are from Ref.[80].

FIG. 4 :
FIG. 4:The calculated spin-independent masses of fully heavy and light systems by Eq.(33).Quark number N is from 2 up to 20 for mathematical extrapolation.

TABLE III :
Computed masses (MeV) of triply heavy baryons, compared to other works cited.Bag radius R 0 is in GeV −1 .

TABLE IV :
Predicted spectra of pentaquarks P bbbb b and P ccccc with comparisons from various works.Bag radius R 0 is in GeV −1 .Masses are in GeV.Magnetic moments µ are in unit of µ N .

TABLE V :
[80]ulated masses M (in GeV) of fully heavy mesons[53], tetraquarks[80]and pentaquarks, with their mass splittings ∆M (in GeV) due to the lowest two J P quantum numbers.

TABLE VI :
Predicted spectra of pentaquarks P bbbbc and P cccc b with comparisons from various works.Bag radius R 0 is in GeV −1 .Masses are in GeV.Magnetic moments µ are in unit of µ N .

TABLE VIII :
Predicted spectra of pentaquarks P cccb b and P cccbc with comparisons from various works.Bag radius R 0 is in GeV −1 .Masses are in GeV.Magnetic moments µ are in unit of µ N .The states denoted by asterisks couples strongly to scattering states.

TABLE IX :
Predicted spectra of pentaquarks P ccbb b and P ccbbc with comparisons from various works.Bag radius R 0 is in GeV −1 .Masses are in GeV.Magnetic moments µ are in unit of µ N .