Regge-like relation and a universal description of heavy–light systems

Using the Regge-like formula (M-mQ)2=πσL\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(M-m_Q)^2=\pi \sigma L$$\end{document} between hadron mass M and angular momentum L with a heavy quark mass mQ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_Q$$\end{document} and a string tension σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document}, we analyze all the heavy–light systems, i.e., D/Ds/B/Bs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D/D_s/B/B_s$$\end{document} mesons and charmed and bottom baryons. Numerical plots are obtained for all the heavy–light mesons of experimental data whose slope becomes nearly equal to 1/2 of that for light hadrons. Assuming that charmed and bottom baryons consist of one heavy quark and one light cluster of two light quarks (diquark), we apply the formula to all the heavy–light baryons including the recently discovered Ωc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Omega _c$$\end{document} and find that these baryons experimentally measured satisfy the above formula. We predict the average mass values of B, Bs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_s$$\end{document}, Λb\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda _b$$\end{document}, Σc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Sigma _c$$\end{document}, Ξc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Xi _c$$\end{document}, and Ωc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Omega _c$$\end{document} with L=2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L=2$$\end{document} to be 6.01, 6.13, 6.15, 3.05, 3.07, and 3.34 GeV, respectively. Our results on baryons suggest that these baryons can be safely regarded as heavy quark–light cluster configuration. We also find a universal description for all the heavy–light mesons as well as baryons, i.e., one unique line is enough to describe both of charmed and bottom heavy–light systems. Our results suggest that instead of mass itself, gluon flux energy is essential to obtain a linear trajectory. Our method gives a straight line for Bc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_c$$\end{document} although the curved parent Regge trajectory was suggested before.


Introduction
Nature has chosen the quantum number 2S+1 L J to classify light and heavy hadrons including light u/d/s and heavy c/b quarks, respectively. This is true for light hadrons which can be treated nonrelativistically but this also holds for heavylight mesons, as analytically derived in Ref. [1] using our semi-relativistic potential model [2,3]. Actually this has been noticed and pointed out in a couple of different contexts, e.g., the string picture/flux tube model [4,5], the quantum mechanical derivation of Regge trajectories [6], the suppression of L S coupling [7,8], the empirical rule of degeneracy among states with the same L [9][10][11][12], etc.
In Ref. [1], we have pointed out that a careful observation of the experimental spectra of heavy-light mesons tells us that heavy-light mesons with the same angular momentum L are almost degenerate. In other words, we have observed that mass differences within a heavy quark spin doublet and between doublets with the same L are very small compared with a mass gap between different multiplets with different L, which is nearly equal to the value of QCD ∼ 300 MeV. This fact is analytically explained by our semi-relativistic potential model [2] which is proposed to describe heavylight mesons, spectra and wave functions. In Refs. [4,5], the authors took a simplest string configuration or flux tube picture of a qq meson based on Nambu's idea [13] of a gluon flux tube for a string. They derived a relation between mass squared and angular momentum. In Ref. [6], the authors also derived a similar relation using quantum mechanics with a Cornell potential model for a meson. In Refs. [7,8], the authors noticed that suppression of L S coupling should occur in the heavy quark symmetry which can be applied to heavylight mesons. In Refs. [9][10][11][12], carefully checking experimental data, they noticed that light hadrons and vector mesons can be well classified by the angular momentum and also find that one unique line is enough to describe vector mesons including φ, ω, and ψ mesons; this they call a universal description.
Low lying mesons are classified in terms of 2S+1 L J , which satisfy the Regge-Chew-Frautch formula M 2 = an +b J +c with constants a, b, and c, principal quantum number n, and total angular momentum J . Or they satisfy the nonrelativistic Regge trajectory for mesons, with the Regge slope α and constants α 0 , β, and β 0 . By using the experimental data, the generalized Regge or the Chew-Frautch formula has been studied in detail in Ref. [9], which gives the values of constants, a, b, and c. Regge trajectories are still effective in determining the spin and parity of newly discovered hadrons. We can have this type of relation which also holds for heavy-light mesons when one notices that mass gaps between states with different L are nearly equal to QC D (see Table 1 in [1]). In this paper, instead of using Eq. (1), as a powerful tool to analyze all the heavy-light systems, we use the formula which was originally derived in Refs. [4,6]. In the next section, following Nambu's picture for hadrons in which quarks are connected by a gluon flux tube [13], we will give a much simpler derivation of this relation between the heavy-light meson mass and angular momentum. Then we apply this formula to heavy-light mesons (D/D s /B/B s ) to obtain linear trajectories in the plane (M − m Q ) 2 vs. L using experimental data as well as theoretical values. Next, we will extend this formula to heavy baryons ( c , b , c , b , c , b , c , and b ), regarding a diquark as a3 color state to see whether it works or not. If this is successful, we can safely say that the heavylight baryons can be well described by a picture in which one heavy quark couples with a diquark. This is one of the motivations of this paper, to see whether diquark picture holds or not since there are some questions to use this concept for baryons. This observation might be checked by lattice gauge theory without fermion.
We also try to check whether a universal description holds, i.e., whether Regge-like lines for X c and X b do overlap or not with X Q being heavy-light systems. The final section is devoted to our conclusions and a discussion, especially of the meaning of the gluon flux energy, M − m Q , and hadron mass M and on the nonlinearity of the parent Regge trajectories for heavy quark systems [14].

Relation between mass and angular momentum
We take Nambu's picture [13] of a hadronic string, which consists only of gluons and at both ends of which quarks are attached in the case of mesons. The authors in Refs. [4,9] took the same picture to derive the Regge formula: where M is a meson mass, L is an angular momentum and σ is a string tension. Their way of derivation is as follows. See also Refs. [10][11][12] for other applications. Assuming the simple string configuration, we obtain the mass M and angular momentum L given by the following equations: where is a length of a string which connects two light quarks at the ends, and v(r ) = 2r/ is the speed of the flux tube at the distance r from the center of rotation. Equations (4) and (5) are obtained by assuming the simplest configuration of a string which connects two quarks at both ends of a string with the speed of light c and rotates around the center of the mass system. Combining these two equations, we arrive at Eq. (3). Let us apply the above idea to a heavy-light meson. In the heavy quark limit, one considers the situation that a heavy quark is fixed at one end and a light quark rotates around a heavy quark with the speed of light c as described in Fig. 1. Then one obtains the following relation: The right hand side of this equation is just 1/2 of Eq. (3) and our numerical plot does not fit with this equation. Because the Nambu's string consists only of gluons, we need to get rid of the effects of quark masses. In our case, in the limit of heavy quark effective theory, we should subtract only the heavy quark mass, m Q with Q = c, b, from M.
Hence the final expression we should adopt is given by Q q v r l Fig. 1 The schematic diagram for depicting the string connecting heavy and light quarks at the ends and rotating around the heavy quark. A heavy quark is fixed at one point The same form of this equation was derived in Ref. [6] by using the simplified potential model with a couple of intuitive approximations.
Modified and elaborated forms of Eq. (7) have been proposed in Refs. [5,15], among which Ref. [5] gives where σ = 2πσ and κ = a parameter determined by a computer simulation. Reference [15] modifies Eq. (8) so that it does not have a singularity at L = 0 and gives where v 2 is a velocity of the heavy quark in a heavy-light hadron system. From our point of view in this paper, the last terms of Eqs. (8) and (9) are not necessary to analyze the experimental data of heavy-light systems in order to compare them with theoretical models since data or theoretical values with the same L are averaged over isospin and angular momentum L. However, we allow for a constant term on the right hand side of Eqs. (8) or (9) as shown in Eq. (12). Reference [15] further studied detailed mass spectra of Q and Q by including the L S coupling.

Numerical plots for heavy-light systems and universal description
According to Eq. (7), we plot the figures for heavy-light mesons, D/B/D s /B s , as well as charmed and bottom baryons, Q / Q / Q / Q / Q with Q = c, b, for the experimental data listed in PDG [16] and some theoretical models Refs. [17][18][19][20][21]. In the following, when plotting experimental data, we adopt the quark masses used in Ref. [17,18] in Eq. (7) as This is because all the figures with experimental data are compared with a model calculation given by Refs. [17,18].

Heavy-Light Mesons
In this subsection, we plot figures in (M − m Q ) 2 vs. L for D, B, D s , and B s mesons taken from experimental data [16] as well as the model calculations of Ref. [17]. D/B mesons: Using Tables 1 and 2, the results are given in Figs. 2 and 3 for D and B mesons, separately.
To compare Eqs. (3) with (7), we give the numerical value of the coefficient of L in Eq. (3) obtained by Afonin,   which is taken from Table 4 of Ref. [9] with radial quantum number n. The calculations of the EFG model [17] for D and B mesons are given by Fig. 3.
Looking at linear equations written on Fig. 2, we obtain As you can see, the coefficients of L are nearly equal to 1/2 of that of light hadrons in Eq. (11). Hence we can conclude that D and B mesons satisfy Eq. (7), support an approximate rotational symmetry of heavy-light mesons claimed in Ref. [1], and the string picture for heavy-light mesons works well. Two lines in Fig. 2 can be written in one figure in the upper row of Fig. 4, which nicely shows the universal description of these mesons. This means that two lines almost overlap, irrespective of the heavy quark flavors. In the same way, we draw a figure for the EFG model [17] which is also given in the upper right of Fig. 4. From now on, we plot two lines for heavy-light systems X c and X b in one figure, which makes it easy to compare the two lines and we can draw conclusions on whether their slopes are close to 1/2 and whether they overlap to confirm a universal description. Finally, we predict the average mass of B with L = 2 to be 6.009 GeV using Eq. (12). D s /B s mesons: In the lower row of Fig. 4, we plot mass squared vs. L for D s and B s mesons using Tables 1 and 2, which presents a similar behavior to that of D and B mesons; it is obvious that they also satisfy Eq. (7). The values of the slope for D s and B s are close to each other, which means D s and B s satisfy a universal description. To show model calculations for D s /B s mesons, we plot the figure for these mesons of the GI model in the lower right of Fig. 4, which shows the GI model well satisfies the linear equation (7) of (M − m Q ) 2 vs. L and two lines almost overlap, i.e., a universal description is confirmed. As one can see, the experimental data is much better than the model calculation in regard to a universal description, which may be due to the fact that the EFG model does not explicitly respect heavy quark symmetry. We also predict the average mass of B s with L = 2 as 6.129 GeV using M B s − m b 2 = 0.650L + 0.261 written on the lower left of Fig. 4.

Charmed and bottom baryons
Regarding two light quarks inside a heavy baryon as a light cluster (diquark), we can apply the formula Eq. (7) to charmed and bottom baryons, in which only either the c or the b quark is included. For the baryons, we take models of Refs. [18][19][20][21] to compare with experiments [16], which includes the effect of heavy quark symmetry. Here we should mention that there are other models for calculating the mass spectrum of heavy baryons; see also the pioneering work [23][24][25]. We also have to take care of the total spin of a diquark, S qq . We call this system a good type when S qq = 0, and a bad type when S qq = 1 according to Ref. [22].  Table 3, we list the present experimental data and EFG model calculated results for Q baryons. Their Regge-like lines are given by the upper row of Fig. 5 for experimental and theoretical values. We take Table 3 The experimental [16] and EFG model values [18] for Q and Q baryons with a good diquark (S qq = 0). Here, the constituent quark masses [17] we adopted to analyze the experimental data as well as EFG calculations are given by Eq. ( theoretical values from Ref. [18]. From these figures, we can conclude that Eq. (7) is satisfied, the slopes are close to 1/2, and a unified description holds. We can predict the average values of some states which have not yet been observed. For instance, by using the line for b , we predict the average mass of b (3/2 + , 5/2 + ) with L = 2 to be 6.145 GeV.
c / b baryons (I = 1, S qq = 1): Q baryons have S qq = 1, i.e., a bad type. We present the experimental data [16] and EFG model results [18] for Q baryons in Table 4. From this table, we can see that only two experimental data for b with L = 0 were observed, i.e., only one point in the L vs. (M − m b ) 2 plane, and experimental data for c with L = 0, 1 have been measured. Hence, we plot figures only for c of experimental and c,b of theoretical values [18] in the lower row of Fig. 5. From these figures, we can conclude again that Eq. (7) is satisfied for c , the slope is close to 1/2, a unified description holds, and the model calculations obey the same rules as the experimental data. We predict the average mass of c with L = 2 to be 3.053 GeV using the linear equation for c written on the lower left of Fig. 5.  Table 3. Because of the lack of experimental data, we can only plot a figure for experimental values of c with L = 0, 1, 2. The theoretical values of Ref. [18] are taken. Their figures are given by the left two of Fig. 6, respectively. From these figures, we can conclude again that Eq. (7) is satisfied for c and the slope is close to 1/2. The model calculations obey the same rules as the experimental data, including a unified description. We predict the average mass of c with L = 2 as 3.068 GeV using the linear equation for c written on the leftmost of Fig. 6. c / b baryons (I = 1/2, S qq = 1): Q baryons have S qq = 1, i.e., a bad type. We also present the experimental data and model calculation for Q baryons in Table 4. There is no definite experimental data for Q with Q = c, b, and hence we can only plot figures for Q of the model calculations given by Ref. [18]. The figure is given in the rightmost one of Fig. 6, which shows that Eq. (7) is satisfied, the slope is close to 1/2, and a unified description holds for the theoretical data. c / b baryons (I = 0, S qq = 1): Q baryons have S qq = 1, i.e., a bad type. Recently the LHCb Collaboration observed six 0 c in the + c K − invariant mass spectrum [26]. However, there is only one b state with L = 0 listed in [16]. Hence, we only list the present experimental data and some of the theoretical results in Table 5. We also plot the left of Fig. 7 for c of experimental data and the right of Fig. 7 for c and b of the theoretical calculation in Ref. [18]. Besides, since there might be higher radial excited states in the recent discovery by LHCb, we also plot figures with principal quantum number n = 1, 2 for theoretical values [18][19][20][21] in Fig. 8. We predict the average mass of c with L = 2 as 3.338 GeV using the linear equation for c written on the left of Fig. 7.
Since experimental data are lacking, the slope is not reliable enough to determine which model is most preferable. Reference [18] satisfies Eq. (7) (EFG), the slope is close to 1/2 and a universal description holds for theoretical data. As for two values of n, there are two figures obtained by the same group [20,21], whose figures, referred to STRV and STRV2 Table 4 The experimental [16] and EFG model values [18] for Q and Q baryons with a bad diquark (S qq = 1). Here, the constituent quark masses [17] we adopted to analyze experimental data as well as EFG calculations are given by Eq. ( in Fig. 8, respectively, are slightly different from each other. The second has a larger slope than the first one and hence, it can be rejected from our point of view. The other two diagrams plotted with the results from CL [19] and EFG [18] models are similar to the plot obtained from Ref. [20] and it may be possible that these three could be close to experiments, which we expect to have in the future.

Conclusions and discussion
In this article, using the Regge-like formula (M − m Q ) 2 = πσ L of Eq. (7)  We have adopted only the confirmed values cited in Ref. [16] to plot figures. Light quarks, u, d, and s, form the so-called chiral particles, e.g., π , K , etc. These quarks have the current quark masses, i.e., very tiny masses. When these quarks are dressed with gluon clouds, they become constituent quarks. In our paper, we treat light quarks as current ones in the string picture so that their masses vanish in the chiral limit. Hence, we should subtract only the heavy quark mass from the hadron. Light current quark masses should not be subtracted from the hadron because a part of a light current quark mass comes from the gluon energy. We have numerically checked the cases in which a light constituent quark mass or a light For heavy-light mesons with the same L, we have used the average mass value of the L-wave states for a plot, i.e., the average value of 2S+1 L J = 1 S 0 and 3 S 1 states for Swave, the average value of 1 P 1 , 3 P 1 , 3 P 0 and 3 P 2 states for P-wave, etc. For heavy-light baryons with the same L, we only have a singlet J = 1/2 for the S-wave, we have used the average value of the J = 1/2 and 3/2 states for the P-wave, the average value of 3/2 and 5/2 states for the D-wave with S qq = 0 (a good type), etc. If there are isospin nonsinglet states, we, of course, have averaged over the isospin states.
Numerical plots have been obtained for all the heavy-light mesons of the experimental data whose slopes become nearly equal to 1/2 of that for the light mesons as expected. A universal description also holds for all the heavy-light mesons, i.e., one unique line is enough to describe both the charmed and bottom heavy-light systems. Surprisingly enough it has been found that both of D/B and D s /B s have had almost the same unique lines as can be seen from Fig. 4. We have also checked the theoretical model of Ref. [17] to see whether this model also obeys the above rules which are satisfied by the experimental data and have found that the EFG model also supports their slopes being 1/2 and a universal description. We have predicted the averaged mass values of B and B s with L = 2 to be 6.01 and 6.13 GeV, respectively.  Regarding that charmed and bottom baryons consist of one heavy quark and one light cluster of two light quarks (diquark), we have applied the formula Eq. (7) to all the heavy-light baryons including the recently discovered c and have found that experimental values of all the heavylight baryons, Q / c / Q / Q / Q with Q = c, b, if they exist, have well satisfied the formula, (M − m Q ) 2 = πσ L and the coefficient of L is close to 1/2, to be compared with that for light hadrons. Since there exist experimental data both for c and b , we have found a universal description only for Q . We have also checked the model calculations of Refs. [18][19][20][21] to see whether they satisfy Eq. (7), have a slope close to 1/2, and a universal description holds and we have found that they really do satisfy all these rules except for Ref. [20]. Because of the unknown assignments of five Fig. 8 Plots of values calculated in EFG [18], CL [19], STRV [20], STRV2 [21] models for Q baryons with principal quantum number n = 1, 2. The best fit lines are given by an equation c , we have also provided Fig. 8 with n = 1, 2, which can be used for future analysis.
When looking at figures for a universal description, one notices that the slope for the bottom heavy-light system has a smaller value than that for the charmed one. This is understandable because the heavy-light system with a b quark is dominated by m b compared with that with m c . The slope for the bottom system is much closer to 1/2 than that for charmed one.
We have also predicted the average mass values of b , c , c , and c with L = 2 as 6.15, 3.05, 3.07, and 3.34 GeV, respectively, by using the straight lines for these baryons. These baryons should be averaged over two spin states (3/2 + , 5/2 + ) for b and c and over six spin states (1/2 + , 3/2 + , 3/2 + , 5/2 + , 5/2 + , 7/2 + ) for c and c . Our results of all the heavy-light baryons suggest that heavylight baryons can be safely regarded and treated as heavy quark-light cluster configuration.
Finally, we would like to comment on the reason why we could analyze all the heavy-light systems. This can be done because we have ignored L S and SS couplings which cause small mass splittings among states with the same L. Otherwise we would have immediately faced serious problems, e.g., other than Q and Q with S qq = 0 (good type), the heavy-light baryons with S qq = 1 (bad type) have cumbersome interactions as has been pointed out in Ref. [15]. Hence, our way of analysis of heavy-light systems is an important and powerful tool to analyze experimental data as well as model calculations since by analyzing data, we can judge whether the experimental data observed are reliable and which model should be adopted or is reliable.
Nonlinearity of the parent Regge trajectoris is observed in model calculations of Ref. [14]. Our analysis suggests that instead of the mass, the gluon flux energy should be used to relate it to the angular momentum for heavy quark systems. Actually, nonlinearity of the parent Regge trajectory for B c obtained in Ref. [14] can be remedied by adopting Eq.
(2), which gives a linear trajectory. That is, for heavy quark systems, the heavy quark mass dominates and determines the curves in the ordinary Regge trajectories, while from our point of view, the gluon flux energy determines the behavior of heavy quark systems.
Future measurements of higher orbitally and radially excited states and their masses of heavy quark systems by LHCb and forthcoming BelleII are waited for to test our observation.