Pentaquark states with the $QQQq\bar{q}$ configuration in a simple model

We discuss the mass splittings for the $S$-wave triply heavy pentaquark states with the $QQQq\bar{q}$ $(Q=b,c;q=u,d,s)$ configuration which is a mirror structure of $Q\bar{Q}qqq$. The latter configuration is related with the nature of $P_c(4380)$ observed by the LHCb Collaboration. The considered pentaquark masses are roughly estimated with a simple method. One finds that such states are probably not narrow even if they do exist. This leaves room for molecule interpretation for a state around the low-lying threshold of a doubly heavy baryon and a heavy-light meson, e.g. $\Xi_{cc}D$, if it were observed. As a by product, we conjecture that upper limits for the masses of the conventional triply heavy baryons can be determined by the masses of the conventional doubly heavy baryons.

We discuss the mass splittings for the S-wave triply heavy pentaquark states with the QQQqq (Q = b, c; q = u, d, s) configuration which is a mirror structure of QQqqq. The latter configuration is related with the nature of Pc(4380) observed by the LHCb Collaboration. The considered pentaquark masses are roughly estimated with a simple method. One finds that such states are probably not narrow even if they do exist. This leaves room for molecule interpretation for a state around the low-lying threshold of a doubly heavy baryon and a heavy-light meson, e.g. ΞccD, if it were observed. As a by product, we conjecture that upper limits for the masses of the conventional triply heavy baryons can be determined by the masses of the conventional doubly heavy baryons.

I. INTRODUCTION
Recently, the LHCb Collaboration confirmed the Ξ cc baryon [1] which was predicted in the quark model (QM) decades ago [2] and first measured by the SELEX Collaboration in 2002 [3]. Excited conventional heavy hadrons are also observed in recent years [4]. We are still on the way to confirm the conventional structure with three heavy quarks. However, the most exciting thing in hadron physics is that more and more unexpected hadrons are detected [5][6][7][8][9][10][11][12][13][14][15][16] and thus hadrons beyond the conventional quark model should exist. This indicates that we may identify the existence of multiquark states before the confirmation of all the QM states. For example, the P c (4380) and P c (4450) states observed in the J/ψN channel by the LHCb Collaboration [17] can be interpreted as pentaquark states since their high masses are not natural for conventional 3q baryons [18]. If a pentaquark can be confirmed, the existence of other multiquark states should also be possible. It is interesting to investigate theoretically which multiquark systems allow bound states.
Since the pentaquark-like states P c (4380) and P c (4450) are much higher than the J/ψN threshold and the J/ψN interaction is not strong (a lattice simulation gives a value around 0.7 fm for the scattering length [19]), the interpretation with J/ψN scattering states seem not to be appropriate. In the literature, various interpretations, such as Σ ( * ) cD * molecules, diquark-diquark-antiquark pentaquarks or diquark-triquark pentaquarks, have been used to understand their nature [9][10][11][12][13][14][15][16]. Assigning them as qqqcc compact pentaquarks is another possible interpretation, where the cc is a color-octet state. In Refs. [20][21][22], investigations along this line were performed. It was found that interpreting the P c (4380) state as such a pentaquark can be accepted. Further investigations on this type configuration of pentaquarks showed that a stable Λ-type hidden-charm state is also possible [21,23]. In the present study, we are going to consider a mirror-type structure by changing heavy (light) quarks to light (heavy) quarks, pentquarks with the QQQqq configuration where QQQ is always a color-octet triply heavy state. In the following discussions, we use QQQ to denote ccc, bbb, ccb, or bbc. When one needs to distinguish the quark contents, we will also use QQQ to denote ccc or bbb and QQQ to denote ccb or bbc.
In the study of the udsQQ states with a colored uds, one found that the quark-quark and quark-antiquark colorspin interactions can all be attractive [21]. Mixing effects among different color-spin structures further lowered the pentaquark masses and finally resulted in a Λ-type state which has a rather low mass. The mass is not far from the value obtained in Ref. [18] where molecule configuration was adopted and hidden-charm pentaquarks were proposed first. In the present work, we try to understand whether colored QQQ is also helpful to the formation of pentqaurks or not. If the answer is yes, a problem to distinguish the pentaquarks from conventional QQQ baryons may arise. If no, one expects that conventional triply heavy baryons would be found experimentally first.
Up to now, theoretical studies have given the masses of the conventional QQQ baryons in various approaches, although there is still no experimental evidence about them. According to these calculations, the mass of the ground Ω ccc baryon, for example, is in the range of 4.6∼5.0 GeV [4,[24][25][26][27]. If an additional quark-antiquark pair surrounds the heavy color source, there is a possibility that the system has a little higher mass and it looks like an excited Ω ccc state. This possibility would be excluded if the QQQqq state has a much higher mass. We would like to extract some information about this problem from the following investigations.
In this article, we consider the mass splittings of the compact triply heavy pentaquarks QQQqq (Q = b, c; q = u, d, s) 1 and estimate roughly their masses in the framework of a simple quark model. In the following Sec. II, the wave functions will be constructed. In Sec. III, we give the chromomagnetic interaction (CMI) matrices for different types of pentaquarks. In Sec. IV, we show our choice of relevant parameters and analyse the numerical results for the spectrum. The last section is a short summary.

II. CONSTRUCTION OF THE WAVE FUNCTIONS
Before the calculation of CMI matrices for the considered ground pentaquark states, we need to construct all the spin-color wave functions. The configuration we consider contains three heavy quarks (ccc, ccb, bbc, or bbb) and a light quark-antiquark pair (nn, ns, sn, or ss, n = u or d), which is simply denoted as 12345. Table I lists the flavor contents of these 16 systems.  For example, C nn = 1/16(m ∆ − m N ) and C nn = 3/64(m ρ − m π ) may be determined with Eq. (17). However, several effective coupling constants, C ss , C cc , C bb , C bc , and C bc , need to be assigned by models or assumptions. Although the Ξ cc state has been observed, the extraction of C cc needs the mass difference between the J = 1/2 and J = 3/2 ground states. We here simply adopt the assumption C cc =C cc for our evaluations. Similarly, we use C ss =C ss , C bb =C bb , and C bc =C bc , where C bc is obtained with the mass difference between B * c and B c [30]. A variation of the values related with heavy quarks does not induce significant differences [31][32][33]. Table II lists all the coupling strengths we will use. At present, we further assume that these parameters can be applied to different systems so that we may estimate the masses of the studied pentaquark states.
According to Eq. (11), the mass of a pentquark state in the chromomagnetic model is where M and H CM I are the mass of a pentaquark and the corresponding eigenvalue of the choromomagnetic interaction, respectively. By introducing a reference system, the mass of the pentaquark state can be written as Here, M ref and H CM I ref are the mass of the reference system and the corresponding chromomagnetic interaction, respectively. Because none of QQQqq pentaquark states is observed, we choose the threshold of a baryon-meson state with the same quark content as M ref . If the simple model could give correct masses for all the hadron states, the above two formulas should be equivalent. In fact, the model does not involve dynamics and the two approaches result in different multiquark masses. One may consult Refs. [21,34] for some discussions on the difference.  [34,35]. The multiquark masses are also generally overestimated [21,33,34,36,37]. Thus we may treat the masses obtained in this method as theoretical upper limits. This fact indicates that attractions inside hadrons cannot be taken into account sufficiently in the − 8026 ± 7 ± 11 8005 ± 6 ± 11 11211 ± 6 ± 12 11194 ± 5 ± 12 − simple chromomagnetic model. A more reasonable approach is to adopt Eq. (19), where the necessary attractions for conventional hadrons have been implicitly included in their physical masses. In this approach, probably there are two choices on the reference threshold. In this case, we present results with both thresholds. To calculate the baryon-meson thresholds, we use the hadron masses shown in table III. The extraction of the above effective quark masses and coupling parameters also relies on these numbers. Here, we assume that the spin of Ξ cc is 1/2, which is presumed in Ref. [40]. If the spin of Ξ cc is 3/2 (consistent with the prediction in Ref. [39]), the multiquark masses estimated with the threshold relating to Ξ cc would be shifted downward by 64 MeV. When determining thresholds relating to other doubly heavy baryons, we adopt the values obtained in Ref. [39]. However, the mass of Ω * cc in that reference was obtained by adding ∼100 MeV to the mass of Ξ * cc . Consider the fact that the quantum numbers of the observed Ξ cc by LHCb have not been determined, we prefer to use another value obtained in Ref. [27]. A different mass of the doubly heavy baryon (see Refs. [26,41] for a collection) will lead to a different mass of the pentaquark. One may compare the adopted mass with the value in table III to estimate the mass difference. For example, if we got a pentaquark mass M 1 with M Ω bc = 6920 MeV while one wants to adopt M Ω * bc = 7065.7 MeV [27], then one would obtain the mass of that pentaquark by adding (7065.7 − 6983.4) MeV to M 1 , where 6983.4 MeV is the mass of Ω * bc in table III. In order for further discussions, we also present a summary for the theoretical investigations in the literature on the masses of triply heavy conventional baryons in Table IV.

B. Numerical results and global features
Based on the content of heavy quarks, one may classify the pentaquarks into four groups: cccqq, bbbqq, ccbqq, and bbcqq or two types: QQQqq and QQQ qq. With all the parameters given above, we get the spectra of these pentaquark systems which are shown in tables V− VIII. In these tables, the columns labeled with the title "Mass" list masses estimated with Eq. (18). The last one or two columns show results estimated with relevant baryon-meson thresholds. For the QQQ qq states, one may use (QQq)-(Q q) or (QQ q)-(Qq) type threshold for our purpose. The resulting masses may have an uncertainty about 100 MeV, which is not a very large value if compared to the mass around 10 GeV. We present results with both type thresholds.
For the QQQnn states, we do not present the isospin indices in the results because the isoscalar and isovector cases   [21,34,35], we have seen that the multiquark masses estimated with such thresholds can be treated as theoretical lower limits. This feature should be related with the attractions inside conventional hadrons that cannot be taken into account in the present model. Since no evidence for triply heavy baryons is reported, with this feature, we may conversely set upper limits for the masses of the conventional triply heavy baryons. The formula is where  limit. Similarly, from our parameters, we have and further M Ω ccb ≤ 8029 MeV and M Ω bbc ≤ 11261 MeV. These limits are much lower than the baryon masses obtained with Eq. (18). Most of the masses obtained in the literature (table IV) are consistent with such constraints. Of course, these upper limits rely on the input masses of the doubly heavy baryons, but are irrelevant with the pentaquark masses.
In this sense, the masses of conventional triply heavy baryons are constrained by those of conventional doubly heavy baryons. On the other hand, if the experimentally observed triply heavy baryons have masses larger than such values, one may conclude that they probably are not ground QQQ states. In Fig. 1, we show the rough positions of the studied pentaquark states and the thresholds of relevant (QQq)−(Qq), (QQq) − (Q q), or (QQ q) − (Qq) type decay patterns. Each diagram corresponds to a group of pentaquarks. We have   adopted the masses listed in the seventh columns of tables V-VIII. From a study of dibaryon states [66], the additional kinetic energy may be a possible effect in understanding why the predition of the H-dibaryon [67] is inconsistent with experimental results. If quark dynamics are considered, probably reasonable pentaquarks have higher masses than those obtained here. That is why we do not use the masses in the eighth columns of tables VII and VIII for the QQQ qq systems when plotting the diagrams. In Fig. 1, the spectrum relies on the input masses of the doubly heavy baryons, but the relative positions comparing to the thresholds relating to the mass estimation is fixed in the present method. For example, the masses of the J = 5/2 ccbnn states (I = 1 and I = 0) may be changed if one uses a different input mass of Ξ bb but the distance between the threshold of Ξ bbB and them is fixed by the CMI matrices. Of course, the mass splittings between the pentaquark states are also fixed only by the CMI matrices. From tables V− VIII and Fig. 1, one may understand some features for the spectra we study: i). There is a hierarchy around 3200 MeV between the four groups of pentaquarks, which is from the mass difference between the bottom and charm quarks. Within each group, the hierarchy caused by the mass difference m s − m n is about 100 MeV; ii). For each system, the states with highest and lowest masses usually have the lowest spin J = 1/2; iii). The QQQsn states generally have higher masses than the QQQns states, which is a result of the difference in coupling strengths, C Qn = C Qs and C Qn = C Qs .
If the considered compact pentaquarks exist, they will decay into lower hadron states. The easiest decay mode,    (QQq) + (Qq), (QQq) + (Q q), or (QQ q) + (Qq), would be due to quark rearrangements. Fig. 1 shows relevant thresholds for such a mode. When 2J (J is the spin of a pentaquark) is equal to a subscript of the label for the baryon-meson channel, the decay into that channel through S-wave is allowed. The other decay mode is (QQQ) + (qq) or (QQQ ) + (qq). Now, the color of the three heavy quarks is changed by emitting a gluon. Because of the possible constraints from the Pauli principle and/or angular momentum, the spin-flip of QQQ may happen and this decay mode is suppressed for some cases where the spin-flip of QQQ is forbidden in the heavy quark limit. We check this feature in the following analyses.

C. Stability of various states
We explore the possible compact structure QQQqq in this article, where QQQ is always a color octet state. Since the mass estimation for the pentaquark states depends on effective quark masses or the masses of the unobserved doubly heavy baryons, what we obtain here are only the rough positions of such states. More accurate values need further dynamical calculations. Before those studies, it is beneficial to discuss preliminarily the stability of such pentaquark states and to find out possible interesting exotic states.
The basic idea is to consider the allowed two-body strong decays and the size of the phase space for decay. From the above results, all the pentaquark states seem to have rearrangement decays and are probably not stable. Since hadrons are composed of quarks, the masses and decay properties are finally determined by the quark-quark interactions. It is helpful to understand the hadron-level features from the inner interactions. If heavy quark spin symmetry plays a role in the decay processes [68][69][70], probably narrow pentaquarks are still possible. We will check this possibility. In the study of the QQQQ systems [33], we checked the effective interaction between two heavy quarks in the mixing case by varying the effective coupling strengths, through which we may roughly guess whether the tetraquark states are stable or not. Here, we also perform a similar study.

The cccqq and bbbqq states
For these states, the spin of ccc or bbb (QQQ) is always 1/2 from the symmetry consideration, irrespective of the total spin of the system, 1/2 or 3/2. Since the spin of the color-singlet Ω ccc or Ω bbb is 3/2, the decay of such pentaquarks into a triply heavy baryon plus a light meson involves an emission of a chromomagnetic gluon and the spin-flip of QQQ. In the heavy quark limit, such a decay process is suppressed. In principle, no symmetry principle suppresses the decay into a doubly heavy baryon plus a heavy-light meson and one expects that narrow QQQqq states would not exist if the pentaquark masses are high enough. To check this point, one may calculate overlapping factor between the initial state wave function and the final state wave function. From the recoupling formula in the spin space, for the cccnn case (other QQQqq cases are similar), one has where the superscripts indicate the spins. Obviously, there is no suppressed coupling from the comparison between these wave functions and those in Eqs. (5) and (6). One may roughly get relative coupling amplitudes between different channels with these wave functions, but could not derive the absolute partial widths without knowing the coupling strength. Therefore, the possible nature that the QQQqq states have generally broad widths cannot be excluded from this coupling analysis. In a multiquark state, the interactions between each pair of quark components together affect the final mass which determines the phase space of decay. After the complicated coupling among different color-spin structures, the property of the interaction between two quark components becomes unclear. It is helpful to understand whether the effective interaction is attractive or repulsive by introducing a measure. To find it, we may study the effects induced by the artificial change of coupling strengths in the Hamiltonian. When a coupling strength is decreased, the resulting mass may be increased or decreased. If the effective interaction between the considered components is attractive (repulsive), the mass would be shifted upward (downward) and vice versa. For the purpose to see the effects, we may define a dimensionless variable where ∆C ij is the variation of a coupling strength in Eq. (11) and ∆M = ∆ Ĥ CM I is the corresponding variation of a pentaquark mass. Then the positive (negative) K ij corresponds to a repulsive (attractive) effective interaction between the ith and jth components. When ∆C ij is small enough, K ij tends to be a constant (∂M )/(∂C ij ). With these constants, we may write an equivalent expression for the pentaquark mass as where M 0 is a constant for a given pentaquark system and its value can be obtained based on the estimation approach, Eq. (18) or (19). One should note that this formula does not mean that M is linearly related with C ij . Now we apply this definition of K ij to the cccnn states. By reducing one relevant coupling strength to 99.99% of its original value, we get a K ij . Each time we change only one C ij while others remain unchanged. The results are collected in table IX where the order of states is the same as that in table V. A positive (negative) number indicates that the corresponding color-spin interaction is effectively repulsive (attractive). In the bcbc or bbbb case, for example, one may understand that a relatively stable tetraquark is favored if the effective bc or bb interaction is attractive. In the present case, not only one quark-quark interaction exists. Usually, it is hard to guess whether a pentaquark state is stable or not just from the signs of K ij . The sum of their contributions gives the final mass shift. From table IX, it is obvious that the effective color-spin interactions between each pair of components in the high mass J = 1/2 cccnn states (I = 1 and I = 0) are repulsive and the states should not be stable. By checking the CMI matrix in Eq. (12) where the interaction between quark components can be attractive, one understands that the mixing effect may change the nature of each interaction. In addition to judge the stability of a pentaquark with Fig.  1, one may equivalently use Eq. (24) to judge whether the decay into a (ccn) − (cn) channel happens or not. Since K cc C cc + K cn C cn > 8/3C cc + 16/3C cn (CMI of Ξ * cc ) and K cn C cn > 16/3C cn (CMI of D * ), the decay into Ξ * cc D * (and thus Ξ cc D * and Ξ cc D) through S-wave is allowed. If this state exists, its width should not be narrow. For a J = 3/2 state (I = 1 or I = 0), the nn and cn interactions are both attractive, which is obvious from Eq. (13) or (24), and the width of this state should be relatively narrower. From masses or effective interactions, S-wave Ξ cc D * and Ξ * cc D decay channels are both opened. For the remaining J = 1/2 states, the attraction appears mainly between a charm quark and a light quark. Now, the allowed S-wave decay channels are Ξ cc D and possibly Ξ cc D * . The same K cc for different states reflects the fact that the spin of ccc is always 3/2. From these discussions, it seems that all these cccnn pentaquarks are not narrow states even if they do exist.
For the other cccqq states, only one isospin is allowed, either 1/2 or 0. The higher J = 1/2 cccns state has three S-wave channels Ξ * cc D * s , Ξ cc D * s , and Ξ cc D s and the J = 3/2 state has two S-wave channels Ξ cc D * s and Ξ * cc D s . The allowed S-wave decay channels for the remaining J = 1/2 cccns state are Ξ cc D s and possibly Ξ cc D * s . Effective contributions from each pair of quark components can also be obtained from the values shown in table IX. In principle, these cccns states should not be narrow. The cccsn states have similar decay properties, but D s (D * s ) is replaced by D (D * ) and Ξ cc (Ξ * cc ) is replaced by Ω cc (Ω * cc ). If one replaces only Ξ cc (Ξ * cc ) by Ω cc (Ω * cc ) in the decay products of the cccns states, the decay patterns of the cccss states are obtained. From the above discussions, the widths of all the cccqq should not be narrow.
Similarly, one gets rough decay properties of the bbbqq states. It is easy to understand that their widths should not be narrow either, since all of them have S-wave rearrangement decay patterns and no suppressed channels are found.
If the interactions between light quark components have large contributions to the pentaquark mass, the state may be low enough and probably stable, because the coupling constant in the CMI model is proportional to 1/(m i m j ). Two examples are the udscc state proposed in Refs. [21,23] and the QQqqq states studied in Ref. [34]. In the present QQQqq case, it is easy to read out contributions from each pair of interaction from table IX. Obviously, the qq effective interactions are not attractive enough (even repulsive). The effective interactions for heavy-heavy and heavy-light components do not provide enough attraction, either. Thus, the masses of QQQqq states are not low enough, which results in the conclusion that such pentaquark states are probably not narrow even if they exist. Of course, further conclusion needs accurate hadron masses and the determination of decay amplitudes.

The ccbqq and bbcqq states
Without the constraint from the Pauli principle, the spin of ccb or bbc (QQQ ) can be both 3/2 and 1/2. The resulting spectrum is more complicated. When the angular momentum of the pentaquark is 5/2, the spin of QQQ must be 3/2. The decay of the pentaquark into Ω QQQ plus a light meson is forbidden or suppressed, but the decay into Ω * QQQ plus a light vector meson not. The latter decay is forbidden only by kinematics. When the angular momentum of the pentaquark is 1/2 or 3/2, the spin of the inside QQQ can be both 1/2 and 3/2. From the eigenvectors in tables VII and VIII, the spin of QQQ in some states is dominantly 1/2. The decays of the 1/2 or 3/2 pentaquarks into Ω QQQ and a qq meson are not suppressed, either. For the (QQq) + (Q q) decay mode, it is easy to understand that each overlapping factor does not vanish from Eqs. (7), (8), and (22). (The nonvanishing overlapping factor for the decay of a J = 5/2 pentaquark is obvious.) For the (QQ q) + (Qq) decay mode, the nonvanishing overlapping factor can be understood similarly since the spin of QQ in a physical (QQ q) state can be both 1 and 0. Even if the spin of QQ in (QQ q) is only zero, the factor still does not vanish. Therefore, no symmetry requirement can suppress the decay of the QQQ qq states. Whether such a state can decay or not depends only on kinematics.
To understand the effective quark-quark or quark-antiquark color-spin interactions, we list K ij of each pair of interactions in table X where ∆C ij = 0.00001C ij . Contributions to the pentaquark mass from each pair of components are easy to obtain. For example, the effective cc or bb interaction is always repulsive. From the table, the interactions in all the highest states (J = 1/2) are repulsive and such states are certainly not stable. For the lowest states (J = 1/2) in each system, the chromomagnetic interactions except for cc or bb are all attractive. To understand whether the interactions result in low enough and thus possible narrow pentaquarks, what we need is to check whether the lowest states with various spins can decay.
First, we focus on the lowest ccbnn states. The thresholds of the decay mode containing a (bcn) baryon are generally lower than those of the mode containing a (ccn) baryon. For the case J = 5/2, the states (I = 1 and 0) probably do not have the S-wave decay pattern Ξ * cc B * , but they may decay into Ξ * bc D * and Ω * ccb ρ or Ω * ccb ω through S-wave. Their widths should not be narrow. For the lowest J = 3/2 states (I = 1 and 0), the channels Ξ bc D * , Ξ * bc D, and Ξ bc D * are all opened. They can also decay into Ω * ccb π, Ω ccb ρ, or Ω ccb ω through S-wave. Their widths should not be narrow, either. For the lowest J = 1/2 states, the decay patterns Ξ bc D * , Ξ bc D, Ξ bc D, Ω ccb π, Ω ccb ω, etc. are all allowed, which should result in broad pentaquark states if they exist. Therefore, narrow ccbnn pentaquarks are not expected. However, this conclusion depends on the pentaquark masses. If the masses estimated with the Ξ bc D are more reasonable, one may conclude that relatively narrow J = 5/2 or J = 3/2 pentaquark states are still possible.
Next, one may analyse similarly the cases ccbns, ccbsn, and ccbss. The basic conclusion is that relatively narrow J = 5/2 pentaquarks are possible only when their masses are low enough.
Finally, the bbcqq systems have similar spectra to the ccbqq systems, but the thresholds of the decay mode containing a bbq baryon are lower than those of the mode containing a bcq baryon. If pentaquark masses shown in Fig. 1 are reasonable, no narrow or stable bbcqq states should exist.
From the above results, one basically finds that all the QQQ qq pentaquark states are probably not stable, although relatively narrow states cannot be excluded completely. The color-spin interaction itself does not lead to enough attraction. Further theoretical and experimental investigations on doubly heavy and triply heavy states may provide more information.

D. Discussions
In obtaining the rough masses of the studied pentaquark states, we have several uncertainties. First, the values of coupling constants C ij in the model Hamiltonian are determined from the conventional hadrons. Whether they can be used to multiquark states is still an open question, as noticed in Refs. [71,72]. Secondly, the mass shifts from CMI energies for a system and thus the pentaquark masses are not well determined. If one uses the effective quark masses determined from conventional hadrons, one obtains hadron masses with overestimated values. To reduce the uncertainty, we phenomenologically adopt a reference hadron-hadron system for our purpose. Whether the selected system is appropriate or not is also a question to be answered. At present, we tend to use a system with high threshold in order to include possible additional effects like the additional kinetic energy [66]. The uncertainty for the obtained pentaquark masses in this method will be checked by the future experimental measurements. Thirdly, the adopted model does not involve dynamics and much information is hidden in constant parameters. The parameters should be different from system to system. However, as a preliminary work for multiquark system where the few-body problem is difficult to deal with, the present investigation on basic features of multiquark spectra is still helpful for further theoretical and experimental studies on multiquark properties. We want to emphasize that the estimated masses do not indicate that all of these states should be bound. Whether such states exist or not needs the experimental measurements to confirm. On the other side, if one compact pentaquark state could be observed, the masses of its partner states may be predicted with our results. If the pentaquark masses shown in Fig. 1 are all reasonable, a question arises to distinguish a compact multiquark state from a molecule state. In understanding the nature of the P c (4380) state, both the (cqq)(cq) molecule configuration and the (cc) 8c (qqq) 8c pentaquark configuration are not contradicted with experimental measurements. The masses of the proposed udscc in these two configurations are also consistent [18,21]. For the present QQQqq case, the situation is different. Our study does not favor low mass cccqq or bbbqq pentaquarks while the investigations at the hadron level [73,74] indicate that the molecule states such as Ξ cc D and Ξ cc D * are possible. If a low-lying pentaquark-like state around such thresholds were observed, the molecule interpretation would be preferred. However, for the possible states around the thresholds of Ξ ccB and Ξ ccB * , one needs to use other information like branching ratios to understand whether it is a compact pentaquark or a molecule, which can be studied in future works. Now we move on to the question how to distinguish an orbitally excited QQQ state from an isoscalar pentaquarklike state if a high mass QQQ-like state was observed. The conventional P -wave states and the pentaquark states may have similar decay patterns if their mass difference is not large. We may check the gaps between the upper limits for the masses of ground QQQ baryons and the lowest (QQn)(Qn) thresholds above which pentaquarks probably exist. They are around 590 MeV for the ccc case, 730 MeV for the bbb case, 600 MeV for the ccb case, and 650 MeV for the bbc case. Usually, an orbital excitation needs an energy around 300 MeV, which is smaller than these gaps. Therefore, one does not need to worry about the exotic nature of the observed QQQ baryon once it is below relevant (QQn)(Qn) thresholds.

V. SUMMARY
In this paper, we have investigated the mass splittings for the pentaquark states with the configuration QQQqq in a chromomagnetic model, where QQQ can be color-octet ccc, bbb, ccb, or bbc. The values of their masses cannot be determined accurately since the model does not involve dynamics. We estimate roughly the masses with two methods. Since no symmetry constraints suppress the S-wave rearrangement decays into (QQq) + (Qq), whether the decays happen or not depends only on kinematics. Although the obtained results have uncertainties, it seems that all the studied pentaquarks are above the respective lowest (QQq)-(Qq) type thresholds and thus do not have low masses and narrow widths. On the other hand, the proposed molecules around the thresholds of Ξ cc D, Ξ cc D * , Ξ ccB , and Ξ ccB * are still possible [73]. Based on our results, once a state around the threshold of Ξ cc D were observed, its molecular nature other than compact multiquark nature is favored.
To see the important mixing effects from different color-spin structures, we express the mass of an obtained pentaquark state with effective chromomagnetic interactions by varying coupling constants in the model Hamiltonian. The relevant coefficients K ij 's are tabulated in tables IX and X. From these tables, it is easy to understand the contributions to the pentaquark mass from each pair of quark interactions. When comparing with the CMI matrices in Sec. III, one finds that the mixing may change the interaction strengths significantly (even signs). Actually, the value of K ij can also be obtained with the eigenvectors in tables V-VIII and the CMI matrices in Sec. III. With the explicit expressions, one can see how the mixing effects contribute to a state. For example, the effective cn interaction in the lower J = 1/2 cccnn state (I = 1 or 0) reads which equals to −9.99C cn if the error bars in the eigenvector are included. As a byproduct of the pentaquark study, we have obtained a conjecture for the mass inequalities for conventional doubly heavy baryons and conventional triply heavy baryons. Such relations may be tested in the future experimental results.
We hope the present study may stimulate further investigations about properties of conventional doubly or triply heavy baryons and multiquark states on both the theoretical side and the experimental side. Search for the exotic triply heavy baryons can be performed at LHC.