Heavy quark dominance in orbital excitation of singly and doubly heavy baryons

A mechanism of the heavy quark dominance in the orbital excitation is proposed in this paper which is testified to be reasonable for singly and doubly heavy baryons. In the relativistic quark model, an analysis of the Hamiltonian figures out the mechanism that the excitation mode with lower energy levels is always associated with the heavy quark(s), and the splitting of the energy levels is suppressed by the heavy quark(s). So, the heavy quarks dominate the orbital excitation of singly and doubly heavy baryons. Furthermore, a physical understanding of this mechanism is given in a semi-classical way. Accordingly, the predicted mass spectra of singly and doubly heavy baryons confirm the rationality of this mechanism. In addition, an interesting consequence of this mechanism is that a heavy-light meson is more likely to be produced in the strong decay of the high-orbital excited states, which is supported by experiments. This mechanism is rooted in the breakdown of the mass symmetry. Therefore, it may be also valid for other multi-quark systems, such as the tetraquarks Qqqq and QQqq, or the pentaquarks Qqqqq and QQqqq.


I. INTRODUCTION
Quantum Chromodynamics (QCD), the theory of quarks and gluons, has been developed for 50 years [1].An important prediction of QCD is the existence of heavy baryons.The study of heavy baryon spectroscopy has always been one of the hot topics.It contributes to a deeper understanding of the properties of heavy quarks, such as the heavy quark symmetry, chiral dynamics, feature of strong interactions, and relevant models inspired by QCD.Until now, there have been various methods to analyze baryon spectra (see [1,2] and Refs.therein).Among them, quark potential models are currently an indispensable tool for guiding our complete understanding of the multitude of observed baryons and their properties [3].
Copley, Isgur and Karl first studied the charmed baryons in a quark potential model with hyperfine interactions [4].While the first prediction of the full spectrum of baryons including heavy-quark baryons was presented by Capstick and Isgur in 1986, where the relativistic quark model was developed [5,6].Twenty years later, Roberts and Pervin analyzed the heavy baryons systematically and exhaustively, using a non-relativistic quark model [7].All of the above studies tried to give a unified theoretical description of all baryons including heavy-quark (Q) baryons, by solving the three-quark system explicitly.These works continue to serve as guidelines for experimental and theoretical researches.However, their calculations predict more heavy-quark baryon excited states than those presented in the spectra extracted from data.It is actually a general problem with baryon spectroscopy, which first appeared in the light-quark (q) baryons and was called the 'missing resonance' problem [1][2][3].
A possible solution to this problem is to partially freeze the dynamical degrees of freedom in a three-quark system.For example, if a baryon consists of a quark and a diquark, the reduction of the number of internal degrees of freedom would lead to a more sparsely populated spectrum [3].Based on this idea, Ebert, Faustov and Galkin calculated the spectra of the singly heavy baryons in the heavy quark-light diquark picture in the framework of the QCD-motivated relativistic quark model, where the excitations between the two light quarks are frozen [8,9].It was found that the available experimental data can be well described and the model predicts significantly fewer states than the genuine three-body picture [9].Later, some theoretical works were carried out to investigate the singly heavy baryons with similar methods [10][11][12][13][14].These methods have also been applied to the spectroscopy of doubly heavy baryons, where the light quark-heavy diquark structure is assumed [15][16][17][18][19].However, the reliability of the extension of this method has yet to be tested.'It is telling that this simple diagnostic is difficult to apply since so little is known of the excited baryon spectrum.'[1] For a three-quark system, quark models feature the orbital excitations in two relative coordinates.As shown in Fig. 1, they are often taken to be the Jacobi coordinates ρ and λ, which is consistent with the flavor sectors of the light quark pair (q 1 q 2 ) in a singly heavy baryon or the heavy quark pair (Q 1 Q 2 ) in a doubly heavy baryon (see Sect.II A).In 2016, Yoshida, Hiyama, Hosaka, Oka and Sadato dealt rigorously with the three-quark systems, and analyzed the low-lying excitations of singly and doubly heavy baryons, respectively [20].They concluded that the excitation of the λmode becomes dominant in low-lying states of singly heavy baryons.While the ρ-mode becomes dominant in low-lying states of doubly heavy baryons.In fact, these features have been mentioned earlier [4,7], and in particular, they were expressed systematically in Ref. [7].Later, Chen, Luo, Liu and Matsuki noticed that almost all observed single heavy baryons could be explained as the λ-mode excited states (in the diquark picture) in theory [21].Inspired by the above theoretical achievements, in 2022, the spectra of singly heavy baryons were studied systematically, in which all orbital excitations were assumed in the λmode [22,23] and the ρ-mode was used alone to study the excitations of doubly heavy baryons [24][25][26] in the framework of the relativistic quark model [5,6].The calculated results are satisfactory, especially, the predicted spectra of singly heavy baryons can reproduce almost all experimental data nicely.Then, the rationality of this approach and the physical nature behind it need to be carefully considered.
As shown in Fig. 1, the determination of the Jacobi coordinates λ in a singly heavy baryon and ρ in a doubly heavy baryon is associated with the heavy quark(s).The same is true of the corresponding orbital excitation modes.Hence, it implies that the orbital excitation may be dominated by the heavy quarks.In this paper, we will comprehensively analyze the heavy quark dominance (HQD) mechanism in orbital excitations of singly and doubly heavy baryons, and discuss the reasonability of the HQD approximation.
The remainder of this paper is organized as follows.In Sect.II, the HQD mechanism is analyzed with the Hamiltonian of the relativistic quark model and a semi-classical approach.In Sect.III, the HQD approximation is tested in singly heavy baryons.Sect.IV discusses the HQD approximation in doubly heavy baryons.The signature of the HQD mechanism in the strong decay is analyzed in Sect.V.And Sect.VI is reserved for the conclusions.

A. Hamiltonian and the wave function
In our previous studies on the spectroscopy of singly and doubly heavy baryons, we adopted the relativistic quark model [5,6].In the calculation, the Gaussian expansion method (GEM) and the infinitesimally-shifted Gaussian (ISG) basis functions have been used as well [27].For the convenience of discussion, a brief introduction of the Hamiltonian and the wave function is presented here.In the relativistic quark model, the Hamiltonian for a three-quark system reads, where the interaction terms H con f i j , H so i j and H hyp i j are the confinement, spin-orbit and hyperfine interactions, respectively.The confinement term H con f i j includes the one-gluon exchange potentials G i j and linear confining potentials S i j , in which the smearing transformations are performed.And the interactions have been modified with the momentum-dependent factors.
With the requirement of the flavor S U(3) F subgroups for the light quark pair (q 1 q 2 ), the singly heavy baryons belong to either a sextet (6 F ) of the flavor symmetric states, or an antitriplet ( 3F ) of the flavor antisymmetric states.Thus, Λ Q and Ξ Q belong to the 3F sector; Σ Q , Ξ ′ Q and Ω Q belong to the 6 F sector.Here Q denotes charm (c) quark or bottom (b) quark.Analogously, in the doubly heavy baryons, Ξ bc and Ω bc belong to the 3F sector, Ξ bb(cc) , Ξ ′ bc , Ω bb(cc) and Ω ′ bc belong to the 6 F sector, with the requirement of the heavy quark pair (Q 1 Q 2 ) as shown in Fig. 1.
Based on the above discussion and the heavy quark effective theory (HQET) [7,28,29], the spin and orbital wave function of the state is written as l ρ (l λ ), L and s are the quantum numbers of the relative orbital angular momentum l ρ (l λ ), total orbital angular momentum L, and total spin of the quark pair s, respectively.j denotes the quantum number of the coupled angular momentum of L and s, so that the total angular momentum J = j± 1 2 .And M J is the 3rd component of J.Then, the baryon state is simply labeled with nL(J P ) j , in which n is the quantum number of the radial excitation.
The orbital part of the above wave function (|l ρ m ρ |l λ m λ ) L is expanded in terms of the Gaussian basis functions.Then, the matrix element of each Hamiltonian term can be obtained on these bases with the Jacobi coordinates ρ and λ as shown in Fig. 1 (for details, see [5,22,23]).Thus, the eigenvalues of the Hamiltonian for the orbital excited states can be calculated.Particularly, the contribution of each term to the eigenvalues can be analyzed, which will be done in the next subsection.
Now we analyze the excitation energies of the λ-mode and the ρ-mode with L = 1.Tables I and II in the appendix list the details of each term contribution to the eigenvalues of the Hamiltonian.For comparison, Types I, II and III in each table are presented, which correspond to the light-quark baryon, singly heavy baryon and doubly heavy baryon, respectively.H in the last column is the total energy (the mass of the state).∆ in the brackets is the excitation energy of each term with respect to the ground state.
As shown in Tables I and II, the energy levels of each (l ρ , l λ )nL mode are dominated by the sum of H 0 , G i j and S i j .While H hyp and H so are responsible for the splitting of the levels.This feature can be seen more clearly in Figs. 2 and  3.For the convenience of the following discussion, we define whose eigenvalues stand for the energy levels of the excitation modes (l ρ , l λ ) L .Then, we use ∆H mode to represent the excitation energies (relative energy levels) of the modes with respect to the pure S -wave mode (0, 0) 0 .
Compared with Type I, as shown in Fig. 2, m 3 in Type II is replaced by a heavy quark mass (m b = 4977 MeV), which leads to two results in Type II: (1) The excitation energy of the (0, 1)1P mode is suppressed; (2) The corresponding energy splitting is also reduced.As a whole, the energy level of the λ-mode and the energy splitting are both suppressed due to the mass of the heavy quark.The same is true in Type III, nevertheless, here the excitation energy of the ρ-mode is suppressed which is caused by the heavy quark pair.The similar phenomenon can be seen for the 6 F sector as shown in Fig. 3.The conclusion in this work is the same as that of Refs.[7,20].We can uniformly summarize that the lower orbital excitation energy is caused by the heavy quark(s) in singly and doubly heavy baryons.This is the HQD mechanism.
Type I This mechanism can be understood in a semi-classical way.The orbital excitation generates the rotational kinetic energy E r.k., which increases the kinetic energies of quarks and the radial distance between quarks.The former contributes directly to H 0 .While the latter raises the values of G i j and (For the 6 S i j due to the one-gluon exchange potential G(r) = − 4α s (r) 3r and linear confinement potential Here b ′ and c ′ are constants.α s (r) is the running coupling constant and has been parameterized [6].Therefore, the rotational kinetic energy E r.k.determines ∆H mode finally.While the rotational kinetic energy in a given state E r.k.depends on the moment of inertia of the state I through E r.k.= L 2 2 I .So, a larger I will cause a smaller E r.k., which will eventually lead to a lower ∆H mode .
For a two-body system rotating in a plane, the following relationships exist, m A r A = m B r B and r A + r B = r.Here m A(B) denotes the mass of A(B) particle.r A(B) is the corresponding distance to the center of mass.Then, the moment of inertia of the system For the P-wave orbital excitation of singly heavy baryons, the λ-mode means (l ρ , l λ )=(0,1).If we take the baryon as a two-body system with m A = m Q ≫ m q and m B = m q 1 + m q 2 ∼ 2m q , I λ ∼ 2m q r 2 λ is obtained based on the above discussion.In the ρ-mode with (l ρ , l λ )=(1,0), the rotation occurs in the two light quarks.And I ρ ∼ 1 2 m q r 2 ρ as m A ∼ m B ∼ m q .In this way, for both of the ( 3F ) and (6 F ) sectors, we can obtain the moment of inertia I for Types I, II and III.By using the calculated values of r 2 ρ 1/2 and r 2 λ 1/2 , we can estimate the value of the rotational kinetic energy E r.k. and compare it with the excitation energy ∆H mode .
The used variables and comparison results are listed in Table III.It is shown that for Type II (the singly heavy baryons), the estimated values of E r.k. in the λ-mode are lower than those in the ρ-mode, being the bigger moment of inertia I in the λ-mode due to the influence of the heavy quark.By comparison, it can be seen that the estimated E r.k. is near ∆H mode .For Type III (the doubly heavy baryons), we can conclude that the heavy quarks dominate the lower E r.k. in the ρ-mode, which qualitatively agrees with ∆H mode as shown in Table III.
The above discussion confirms our understanding that the rotational kinetic energies E r.k. of the excitation modes involving heavy quark(s) are significantly lower, which determines the relative energy levels of the excitation modes.There-fore, the HQD mechanism has been understood in the semiclassical way.On the other hand, the suppression of the energy levels splitting for heavy baryons can be explained in the heavy quark limit [7,30,31].We may conclude that the excitation mode with lower energy levels is always associated with the heavy quark(s), and the energy level splitting is also suppressed by the heavy quark(s).So, the heavy quarks dominate the orbital excitation of singly and doubly heavy baryons.

III. THE TEST OF HQD APPROXIMATION IN SINGLY HEAVY BARYONS
A. The orbital excitation of singly heavy baryons For further investigating the HQD approximation in singly heavy baryons, the evolution of the excitation energies with m 3 for different modes is plotted in the cases of the P-, D-and F-wave states.At the same time, the relative energies of the (1, 1)1S ( 12 + ) and (2, 2)1S ( 12 + ) states with respect to the pure S -wave state (0, 0)1S ( 12 + ) are also given for comparison, as shown in Fig. 4.
In the 3F sector, the energies of the (1, 1)1S ( 12 (2, 2)1S ( 12 + ) states are obviously higher than that of the pure S -wave state.For the P-, D-and F-wave states, the excitation energies of most modes are quickly suppressed as m 3 increases, except for the ρ-mode.And the λ-mode is indeed the one with the lowest excitation energy.
It is reasonable to believe that the excitation mode with the lowest energy is most stable and has the greatest probability of being observed experimentally.Then, all excitation modes with higher energies can be approximately ignored.Thus, the 'true' ground state comes from the pure S -wave state, and the 'real' orbital excited state originates from the λ-mode, which is dominated by the heavy quark(s).This is the HQD approximation as mentioned before in this paper.
The situation of the 6 F sector is similar to that of the 3F sector that the λ-mode is dominant.But the phenomenon is not so obvious when m 3 ∼ 1.6 GeV for the P-wave states, where the excitation energy of the (0, 1)1P( 1 2 + ) 1 state is only slightly lower than that of the (1, 0)1P( 1 2 + ) 1 state.It implies that c quark (m c = 1.628GeV) is not heavy enough and hardly treated as a heavy quark, which was also pointed out in some theoretical works [32].Even so, the HQD approximation is still applicable to the 6 F sector.

B. Predicted excited states with HQD approximation
With the HQD approximation, the excited spectra of singly heavy baryons were systemically investigated with a unified set of parameters [22,23].The results show that the predicted masses reproduce the experimental data well.As shown in Table IV, for the well established baryons in experiment [33,34], the calculated masses are very close to the data, apart from Λ c (2940) + [22,35].The maximum difference of the mass values between theory and experiment is generally no more than 20 MeV.
For the observed baryons whose J P values are not yet confirmed [33,[36][37][38][39][40], the calculated results are also helpful for their rational assignments.As shown in Table V, these baryons have been divided into several groups.In particular, some members in the same family were grouped together and labeled with 1P(6 F ).With the help of the mass spectra, one can see them more clearly.As shown in Fig. 5 ) 2 , and increase in energy.Some baryons of these 1P states have been observed as shown in Table V and Fig. 5.However, because their energies are so close to each other, the identification of these baryons becomes very difficult.So they can only be labeled with 1P(6 F ), as shown in Table V.
Fig. 5 shows the predicted spectra agree well with the experimental data.And the 'missing resonance' problem does not seem to exist here.It makes one more confident that the HQD mechanism is reasonable and the HQD approximation is feasible.In addition, Fig. 5 shows an obvious shell structure of the spectrum.It implies that L might be approximated as a good quantum number.If this is true, the γ transitions would occur between the neighbouring orbits.This feature has been partially observed in experiment.For example, Ξ c (2815) 0 → Ξ 0 c γ is the γ transition between the 1P( 3 2 − ) and 1S ( 1 2 + ) states [33].

IV. HQD MECHANISM IN DOUBLY HEAVY BARYONS
The doubly heavy baryons Ξ cc (Ω cc ), Ξ bb (Ω bb ) and Ξ ′ bc (Ω ′ bc ) belong to the 6 F sector, and Ξ bc (Ω bc ) belongs to the 3F sector, as mentioned in Sect.II A. To examine the HQD mechanism in the 6 F sector, we set m 2 = m 1 and change m 1 from 0.3 GeV to 5.1 GeV.For the 3F sector, m 2 is set to be m c m b × m 1 , which stands for Ξ bc (or Ω bc ) baryon when m 1 approaches to 5.0 GeV (see the caption in Fig. 6).
As shown in Fig. 6, the orbital excitation of doubly heavy baryons presents the same feature as that of singly heavy baryons.In this case, nevertheless, the mode with the low-   est energy is the ρ-mode and the lowest orbital excitation energy is still dominated by the heavy quarks.Hence, the HQD mechanism remains valid and the HQD approximation is feasible for doubly heavy baryons.
Based on the HQD approximation, the spectra of doubly heavy baryons were investigated in the ρ-mode [24][25][26].The predicted ground state mass of Ξ cc is 3640 MeV, which is very close to the experimental data 3621.6±0.4MeV [33].The spectrum also presents an obvious shell structure as shown in Fig. 7.It is expected that more doubly heavy baryons will be observed and can be used to further test the HQD approximation.
In addition, some works have noticed the excitation of the heavy quark pair is important and can not be ignored [16,18].In those works, the spectra of the doubly heavy baryons therefore contain the excitation of the heavy quark pair.And the heavy quark pair was not treated as a point-like diquark, which was used in some other works [58][59][60].Unlike both of them, the HQD approximation freezes the orbital excitation of the λ-mode and only considers the ρ-mode, which predicts fewer excited states than the former [16,18], and lowers the excited energy levels than the latter [58][59][60].

V. SIGNATURE OF HQD MECHANISM IN STRONG DECAY
If the HQD mechanism is reasonable, an interesting phenomenon would occur in the strong decay.The calculation shows that for a singly heavy baryon, the root-mean-square value of r λ increases with the orbital excitation [22,23].It means the heavy quark is further away from the light quark pair, as the orbital quantum number increases.For a highorbital excited state, when the strong decay occurs, the heavy quark is therefore likely to capture the light antiquark which comes from the creation of a quark-antiquark pair in the vac-  I, and the open circles are the ones in Table II.The three bigger red circles represent Λ c (2940) + , Ξ c (3123) + , and Ω c (3120) + , respectively, which can not be assigned very well.More excited states can be found in Refs.[22,23].
uum, according to the 3 P 0 model [61,62].So, the heavy-light meson is more likely to be produced.
From the Particle Data Group (PDG) [33], one can see the D 0 p decay mode occurs in Λ c (2860) + and Λ c (2880) + , which are identified as the D-wave states in theory [9,22,63].Similarly, the ΛD + decay modes are also observed experimentally for Ξ c (3055) + and Ξ c (3080).So, they are assigned to the 1D( 3 2 ) doublet states of the 3F sector in Table V.For Λ + c , the ground state of the Λ c family, the exotic mode of pγD has been observed in experiment.But its branching ratio is very small.Additionally, the pD 0 mode is observed in the strong decay of Λ c (2940) + , which means it should be a high-orbital excited state.
In Table V and Fig. 5, Ξ b (6327) 0 and Ξ b (6333) 0 is identified as the D-wave doublet states in this work.In their observed strong decay modes, however, the heavy-light meson is absent [33].From Ref. [23], we find their r 2 λ 1/2 values are 0.743 fm and 0.751 fm, respectively, which are smaller than those (0.825 fm and 0.843 fm) of the D-wave states Ξ c (3055) + and Ξ c (3080).It implies the heavy-light meson is more difficult to be produced in the strong decay of the bottom baryons, as compared to the charmed baryons.
So far, the two-body strong decay mode containing a heavylight meson is mainly observed in the high-orbital excited states of the charmed baryons.The same phenomenon is expected to occur in doubly heavy baryons, which might be helpful in searching for the doubly heavy baryons with the high-orbital excitation in experiments.

VI. CONCLUSIONS
The singly and doubly heavy baryons are special, due to the large mass values of charm and bottom quarks.Since these heavy quarks are much more massive than the light quarks, they break the mass symmetry.So, the S U(3) F flavor symmetry valid for light-quark baryons are not well applied here.On the other hand, the heavy quark limit is not fully applicable here either, because the masses of the heavy quarks are much less than infinity.Inspired by the progress of the heavy baryon spectroscopy in recent decades, we speculate a mechanism of the heavy quark dominance (HQD) in the orbital excitation of  singly and doubly heavy baryons.
By analyzing the Hamiltonian of the relativistic quark model for a three-quark system, we tried to understand the HQD mechanism.It is shown that the excitation mode with lower energy levels is always associated with the heavy quark(s), and the splitting of the energy levels is suppressed by the heavy quark(s) as well.In other words, the heavy quarks dominate the orbital excitation of singly and doubly heavy baryons.In a semi-classical way, we show that the rotational kinetic energy generated by the orbital excitation is positively related to the excitation energy of the mode.And the smaller rotational kinetic energy is determined by the heavy quark(s).So, the heavy quarks eventually dominate the lower orbital excitation for singly and doubly heavy baryons.Thus, we get a physical understanding of this mechanism.
The HQD mechanism is then tested by comparing various orbital excitation modes for singly heavy baryons.It is shown that the excitation mode dominated by the heavy quark(s) has lower energy levels, therefore, it is most stable and has the greatest probability of being observed experimentally, which justifies the HQD approximation.With the HQD approximation, the predicted excited spectra of singly heavy baryons can reproduce the available experimental data nicely.Additionally, the shell structure of the mass spectra implies L can be regarded as a good quantum number approximately, which is supported by the γ transition in experiment.Thus, the HQD mechanism is proved to be reasonable in singly heavy baryons.For doubly heavy baryons, their orbital excitations are also dominated by the heavy quarks.The corresponding mass spectra with the HQD approximation are expected to be tested in future experiments.The HQD mechanism also predicts an interesting phenomenon that the heavy-light meson is more likely to be produced in the strong decay of the highorbital excited states for singly and doubly heavy baryons.This kind of strong decay mode may be the signature of the HQD mechanism as well.
The HQD mechanism is rooted in the breakdown of the mass symmetry in a multi-quark system.Therefore, it should be valid not only for singly and doubly baryons, but also for the multi-quark systems containing one or two heavy quarks, such as the tetraquarks Qqqq and QQqq, or the pentaquarks Qqqqq and QQqqq.cc .More excited states can be found in Refs.[24][25][26].

1 FIG. 1 :
FIG.1: Jacobi coordinates for the three-quark systems.The heavy quarks are represented by bigger balls and the light quarks by smaller balls.They are numbered for ease of use in calculations.

FIG. 2 :
FIG. 2: Relative energy levels ∆H mode and the energy splittings of the ρ-mode and the λ-mode for the 3F sector with Types I, II and III.The energy values are in MeV.
, for Σ c(b) , Ξ ′ c(b) and Ω c(b) , there are five excited states for the 1Pwave, and their energies are almost degenerate.The five states are 1P

2 -FIG. 4 : 2 + 1 2+
FIG. 4:The evolution of excitation energies with m 3 for different modes of singly heavy baryons.Here only the modes with J = L ± 1 2 are selected.m 3 changes from 0.2 GeV to 5.0 GeV.The excitation energies are measured from the pure S -wave states 1S ( 1 2 + ) 0 for the 3F and

FIG. 5 :
FIG.5: Calculated spectra of singly heavy baryons and the relevant experimental data.'++','+', '0' and '-' in the brackets indicate the charge states of baryons.The solid black circles denote the baryons in TableI, and the open circles are the ones in TableII.The three bigger red circles represent Λ c (2940) + , Ξ c (3123) + , and Ω c (3120) + , respectively, which can not be assigned very well.More excited states can be found in Refs.[22,23].

2 -FIG. 6 :
FIG. 6: Evolution of excitation energies with m 1 for different modes of doubly heavy baryons, where m 1 changes from 0.3 GeV to 5.1 GeV.The value of m 3 is equal to 0.22 GeV.m 2 is set as m 2 = m 1 for the 6 F sector.For the 3F sector, m 2 is set as m 2 = mc m b × m 1 , which corresponds to Ξ bc or Ω bc when m 1 = 4.977 GeV.