Low-lying charmed and charmed-strange baryon states

In this work, we systematically study the mass spectra and strong decays of 1P and 2S charmed and charmed-strange baryons in the framework of non-relativistic constituent quark models. With the light quark cluster–heavy quark picture, the masses are simply calculated by a potential model. The strong decays are studied by the Eichten–Hill–Quigg decay formula. Masses and decay properties of the well-established 1S and 1P states can be reproduced by our method. Σc(2800)0,+,++\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Sigma _c(2800)^{0,+,++}$$\end{document} can be assigned as a Σc2(3/2-)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Sigma _{c2}(3/2^-)$$\end{document} or Σc2(5/2-)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Sigma _{c2}(5/2^-)$$\end{document} state. We prefer to interpret the signal Σc(2850)0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Sigma _c(2850)^0$$\end{document} as a 2S(1/2+)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2S(1/2^+)$$\end{document} state although at present we cannot thoroughly exclude the possibility that this is the same state as Σc(2800)0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Sigma _c(2800)^0$$\end{document}. Λc(2765)+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda _c(2765)^+$$\end{document} or Σc(2765)+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Sigma _c(2765)^+$$\end{document} could be explained as the Λc+(2S)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda _c^+(2S)$$\end{document} state or Σc1+(1/2-)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Sigma ^+_{c1}(1/2^-)$$\end{document} state, respectively. We propose to measure the branching ratio of B(Σc(2455)π)/B(Σc(2520)π)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {B}(\Sigma _c(2455)\pi )/\mathcal {B}(\Sigma _c(2520)\pi )$$\end{document} in the future, which may disentangle the puzzle of this state. Our results support Ξc(2980)0,+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Xi _c(2980)^{0,+}$$\end{document} as the first radial excited state of Ξc(2470)0,+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Xi _c(2470)^{0,+}$$\end{document} with JP=1/2+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J^P=1/2^+$$\end{document}. The assignment of Ξc(2930)0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Xi _c(2930)^0$$\end{document} is analogous to Σc(2800)0,+,++\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Sigma _c(2800)^{0,+,++}$$\end{document}, i.e., a Ξc2′(3/2-)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Xi ^\prime _{c2}(3/2^-)$$\end{document} or Ξc2′(5/2-)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Xi ^\prime _{c2}(5/2^-)$$\end{document} state. In addition, we predict some typical ratios among partial decay widths, which are valuable for experimental search for these missing charmed and charmed-strange baryons.


Introduction
At present, the particle data group (PDG) lists nine charmed and ten charmed-strange baryons [1]. They are c (2286) + ,  [2] and Belle [3]. Very recently, the Belle Collaboration updated the measurements of c (2580) 0,+ , c (2645) 0,+ , c (2790) 0,+ , c (2815) 0,+ , and c (2980) 0,+ [4]. On the other hand, new decay modes for the higher excited charmed baryon states have been found by experiment. For instance, the decay channel of D + was first found for c (3055) + and c (3080) + , and the following ratios of branching fractions were first reported by Belle several months ago [5]: where the uncertainties are statistical and systematic. Obviously, these new measurements are very useful to understand the nature of these excited charmed baryon states.
Theoretically, the charmed baryons which contain one heavy quark and two light quarks occupy a particular position in the baryon physics. Since the chiral symmetry and heavy quark symmetry (HQS) can provide some qualitative insight into the dynamics of charmed baryons, the investiga-tion of charmed baryons should be more helpful for improving our understanding of the confinement mechanism. The spectroscopy of charmed baryons has been investigated in various models. So far, the several kinds of quark potential models [6][7][8][9][10], the relativistic flux tube (RFT) model [11,12], the coupled channel model [13], the QCD sum rule [14][15][16], and the Regge phenomenology [17] have been applied to study the mass spectra of excited charmed baryons, and so did lattice QCD [18,19]. The strong decay behaviors of charmed baryons have been studied by several methods, such as the heavy hadron chiral perturbation theory (HHChPT) [20,21], the chiral quark model [22,23], the 3 P 0 model [24], and a non-relativistic quark model [25]. The decays of 1P c and c baryons have also been investigated by a light front quark model [26,27], a relativistic three-quark model [28], and the QCD sum rule [29].
Although many experimental and theoretical efforts have been made for the research of charmed baryons, most of the 1P and 2S charmed baryons are not yet established. Several candidates, including c (2765) + , c (2800) 0,+,++ , c (2930) 0 , and c (2980) 0,+ are still in controversy. c (2765) + was first observed by the CLEO Collaboration in the decay channel of c (2765) + → + c π + π − [30], and confirmed by Belle in the c (2455)π mode [31,32]. Because both + c and + c excitations can decay through + c π + π − and c (2455)π , we even do not know whether the observed charmed baryon signal around 2765 MeV is the + c or + c state, or their overlapping structure [33]. In the e + e − annihilation process, an isotriplet state, c (2800) 0,+,++ , was observed by Belle in the channel of + c π , and it was tentatively identified as the c2 state with J P = 3/2 − [34]. Interestingly, another neutral resonance was later found by BaBar in the process of B − → * 0 cp → + c π −p with the mass 2846 ± 8 ± 10 MeV and decay width 86 +33 −22 MeV [35]. The higher mass and the weak evidence of J = 1/2 indicate that the signal observed by BaBar might be different from the Belle observation. In this paper, we denote the signal discovered by BaBar by c (2850) 0 . c (2930) 0 , which was only seen by BaBar in the decay mode + c K − [36], still needs more confirmation. c (2980) 0,+ was first reported by Belle in the channels + c K − π + and + c K 0 S π − [37], and was later confirmed by Belle [4,38] and BaBar [39] in the channels c (2580)π , c (2645)π and c (2455)K , respectively. However, the decay widths reported by Refs. [4,[37][38][39] were quite different from each other. More experimental information as regards the charmed baryons can be found in Refs. [40][41][42][43].
Obviously, a systematic study of masses and decays is required for these unestablished charmed baryons. More importantly, most of 2S and 1P charmed baryons have not yet been detected by any experiments. Such a research can also help the future experiments to find them. In the present work, we will focus on both the mass spectra and the strong decays of low-lying 1P and 2S charmed baryons. We pay attention to only the charmed baryons inside of which degrees of freedom of two light quarks are frozen. It means that two light quarks are not considered to be excited, neither radially nor orbitally. As illustrated in Ref. [44], these kinds of charmed excitations carry lower-excited energies, which means these excited charmed baryons may more likely be detected by experiments. Fortunately, our results indicate that most of the observed charmed baryons can be accommodated in this way.
The paper is arranged as follows. In Sect. 2, the masses of low-excited charmed baryons are calculated by the nonrelativistic quark potential model. In Sect. 3, the Eichten-Hill-Quigg (EHQ) decay formula, which is employed to study the strong decays of excited charmed baryons, is introduced. The properties of low-lying charmed baryon states are fully discussed in Sect. 4. Finally, the paper ends with a conclusion and an outlook. Some detailed calculations and definitions are collected in the appendices.

Treating charmed baryon system as a two-body problem
To study the baryon dynamics, one crucial question which should be answered is "What are the relevant degrees of freedom in a baryon?" [45]. In constituent quark models, a baryon system consists of three confined quarks. Thus the dynamics of a baryon resonance is surely more complex than a meson. Due to the HQS, however, the dynamics of charmed baryons could be greatly simplified. The HQS suggests that the couplings between a c quark and two light quarks are weak [46]. Therefore, two light quarks in a charmed baryon could first couple with each other to form a light quark cluster. 1 Then the light quark cluster couples with a charm quark, and a charmed baryon resonance forms. With this assumption, two light quarks have the same status to a c quark, including the average distances to a c quark. In the light cluster-heavy quark picture, the dynamics of heavy baryon can be simplified. In the non-relativistic constituent quark model, the spin-independent part of the Hamiltonian is where the Cornell potential [47] is used as a phenomenological confining term. The Jacobi coordinates are usually taken to deal with the three-body problem for baryons. ρ, λ, and R are related to the quark positions by where the indices, 1, 2, and Q, are for two light quarks and a heavy quark, respectively. The momenta p ρ , p λ and p R , which are conjugate to the Jacobi coordinates above, can be defined easily. Now the spin-independent Hamiltonian becomes where m ρ = m 1 m 2 /(m 1 + m 2 ), m λ = (m 1 + m 2 )m Q /M, and M = m 1 +m 2 +m Q . According to the definitions above, the relative motion between two light quarks is usually called the ρ-mode, while the one between the center of mass of the two light quarks and the heavy quark is called the λ-mode.
As emphasized above, the average distances to c quark should be equal for two light quarks in a cluster, i.e., τ ≡ r 1Q = r 2Q (see Fig. 1). In practice, we should solve the following Schrödinger equation for the mass of a heavy baryon: Since the excited mode between two light quarks is not considered in this paper, the light quark cluster can be treated as a block with the antitriplet color structure and peculiar size. Specifically, a color singlet baryon system should be formed as 3 q 1 ⊗ 3 q 2 ⊗ 3 Q 3 cl. ⊗ 3 Q 1 Q−cl. . In this way, a heavy baryon could be treated as a quasi-two-body system. In the light quark cluster-heavy quark picture, Isgur has discussed the similarity of dynamics between heavy baryons and heavy-light mesons [48]. It should be stressed that the scenario of a light cluster-heavy quark picture is not contradictory to the one-gluon exchange interaction. Since the Fig. 1 A single heavy baryon in the light cluster-heavy quark picture. With the SU(3) flavor symmetry, the relation of τ = λ 2 + ρ 2 /4 is obvious color-spin interaction is proportional to the inverse of the quark masses, two light quarks in the heavy-baryon system are expected to strongly couple to each other [49]. Thus, they may develop into a quark cluster. In fact, the existence of a light quark cluster correlations was partly confirmed by the lattice QCD [50] and the Bethe-Salpeter equation [51]. Thus, we expect the light quark cluster to be an effective degree of freedom for a charmed baryon. Since the ρ mode of a cluster is not considered here, the Schrödinger equation (3) is simplified as 2.2 Adopted effective potentials As a whole, the light quark cluster which occupies an antitriplet color structure interacts with the c quark. Then we would like to substitute λ (distance between light cluster and c quark) for τ (distance between light quark and c quark).
To this end, the effective potential [52] H conf describes the interaction between the cluster and c quark, where ν is an adjustable parameter. This approximation can greatly decrease the computational complexity. As shown in Tables 2 and 3, the mass spectra given in this way are reasonable for the low-lying excited charmed baryons. As a two-body problem, we treat the masses of different kinds of clusters as parameters and first fix them before calculating the masses of the low-lying charmed baryons. According to the flavor and spin, the light cluster can be classified into two kinds: one is the "scalar" cluster, and another is the "vector" cluster. Constrained by the Pauli exclusion principle, the total wave function of the light quark cluster should be antisymmetric in the exchange of two quarks. Because the spatial and color parts of this light quark cluster are always symmetric and antisymmetric, respectively, the function |flavor × |spin should be symmetric. Therefore, the scalar light quark cluster [qq] (S = 0) is always flavor antisymmetric, and the axial-vector light quark cluster {qq} (S = 1) is flavor symmetric. In Jaffe's terminology [49], the "scalar" and "vector" quark clusters are named the "good" and "bad" quark clusters, respectively. The masses of the "good" light quark cluster [qq] and [qs] are taken from our previous work where m [qq] and m [qs] were fixed as 450 and 630 MeV by the RFT model [11], respectively. Following Jaffe's method [49], the bad light quark cluster masses can be evaluated by the following relationships: the B-type baryons for convenience. Due to an antitriplet color structure, the spin-dependent interactions between light cluster and c quark are expected to be the same as the meson systems. In a constituent quark model [53], the spin-dependent interaction is written as The color contact interaction H cont Q−cl. is usually given by the following form: where S Q and S cl. refer to the heavy quark and light cluster spins. A Gaussian-smeared function (σ/ √ π) 3 e −σ 2 λ 2 is normally used forδ σ (λ) [54]. If the SU(3) flavor symmetry is kept for the charmed baryons, we may modify the color contact interaction as where the mass of a light quark cluster is just replaced by a unit. This assumption is supported by the mass differences of the 1S B-type charmed baryons, The mass differences shown above are mainly due to the color contact interaction in the quark potential model. Clearly, these values are almost independent of the light quark cluster masses. The color tensor interaction in Eq. (6) is Finally, H SO Q−cl. denotes the spin-orbit interaction which contains two terms. One is the color magnetic interaction which arises from one-gluon exchange, where S denotes the spin of a baryon, S = S Q + S cl. . Another spin-orbit interaction is the Thomas-precession term, . (11) To reflect the importance of the heavy quark symmetry, we rewrite the spin-dependent interactions as The degrees of freedom of the light quark cluster are characterized by its total angular momentum j cl. , i.e., j cl. = S cl. + L. Obviously, the orbital angular momentum L of a charmed baryon in the present picture is defined by the angular momentum between light quark cluster and c quark, i.e., L = L λ . The tensor operator is defined asŜ 12 With the confining term of Eq. (5), the coefficients V ss , V 1 , V 2 , and V t in Eq. (12) are defined by

Getting masses of charmed baryons
In our calculation, the following Schrödinger equation is solved for the nS state: The confining and contact terms have been given by Eqs. (5) and (8). For the orbital excitations, all spin-dependent interactions are treated as the leading-order perturbations. Our calculation indicates that the color contact interaction can be ignored for the orbital excitations. Two bases are employed to extract the mass matrix elements. One is defined by the eigenstates |S cl. , L , j cl. , S Q , J ( j j coupling scheme) and the other is by |S cl. , S Q , S, L , J (L S coupling scheme). The relation between these two bases is Here j l = L + j cl. and subindices and h of the last column mean low and high states in mass after includingṼ ss andṼ t interactions Due to S cl. = 0, only V 2 S Q · j cl. contributes to the masses of G-type charmed baryons. With a bad light quark cluster, however, B-type charmed baryons have more complicated splitting structures. Within the framework of the heavy quark effective theory, the spin of an axial-vector light quark cluster, S cl. , first couples with the orbital angular momentum L. As illustrated in Fig. 2, in the heavy quark limit m c → ∞, there are only three states which are characterized by j cl. for 1P charmed baryons. When the heavy quark spin S Q couples with j cl. , the degeneracy is resolved and the five states appear. These are two J P = 1/2 − , two J P = 3/2 − , and one J P = 5/2 − states. Lastly, the states with the same J P mix with each other by the interactions of V ss S Q · S cl. and V tŜ12 , and physical states are formed.
We now turn to a calculation of the mass matrix in the j j coupling scheme. For 1P states with J P = 1/2 − , the mass matrix is given by Similarly, for two states with J P = 3/2 − , For the J P = 5/2 − state, In the following, we denote |S cl. , L , j cl. , S Q , J as | j cl. , J P for brevity. Then the notations | 1/2 and | 3/2 appearing above are defined by Table 1 Values of the parameters of the non-relativistic quark potential model. The unit of b is GeV ν+1 which varies depending on each value of ν The mass matrix of 1D states can also be obtained by the similar procedure. As shown above, there are seven parameters in the non-relativistic quark potential model, which are m Q , m cl. , b, α, γ , ν, and C Qqq . All values of parameters are listed in Table 1. If the SU(3) flavor symmetry is taken into account for the charmed and charmed-strange baryons, the dynamics of + c states should be like c . The case of c and c is alike. Accordingly, the same value of γ is selected for the G-type charmed baryons, as well as the case of B-type.
We have adopted the typical values in the quark potential models for m c , b, α, and ν (see Table 1). It is an effective method to investigate charmed baryons in the heavy quarklight quark cluster picture. We do not expect the values of ν to be the same for G-type and B-type baryons. Here, ν of + c / c is slightly larger than c / c . The predicted masses of the low-excited charmed baryons are collected in Tables  2 and 3.
As mentioned earlier, the nonzero off-diagonal elements in mass matrices of 1/2 | H S | 1/2 and 3/2 | H S | 3/2 cause the mixing between two states with the same J P but different j cl. . However, the mechanism of mixing effects in hadron physics is still unclear. In principle, a physical hadron state with a specific J P comprises all possible Fock states with the same total spin and parity. As the most famous member of the XYZ family, X (3872) may be explained as a mixture of charmonium and the molecular state with J PC = 1 ++ [55]. Here we take the | j cl. , J P basis to describe the mixing for the B-type baryons. Then two physical states characterized by different masses can be denoted   Here we have denoted the physical states by their masses (see Table 3). The mixing angles for other states in Table 3 with the same J P are listed in Table 4.
The uncertainty may exist in the mixing angles. Firstly, the loop corrections to the spin-dependent one-gluon-exchange potential may be important for the heavy-light hadrons. As Table 4 The mixing angles for the 1P, 2P, and 1D c / c states  [53] can be well explained by the corrected spin-dependent potential [59,60]. If we use this type of potential in our calculation, of course, the mixing angle will change. Secondly, the mixing angles depend on the parameters. Thirdly, there are other mechanisms, e.g., hadron loop effects [61], which may contribute to the mixing phenomenon in hadron physics. Anyway, we expect that the mixing angles in Table 4 reflect the main features of the mixing states. Due to the uncertainties of the mixing angles, however, we ignore the mixing effects as the first step to the study of the decays of charmed excitations in the next subsection. Obviously, it is a good approximation only when the mixing effects are not large. Fortunately, this crude procedure is partially supported by the former analysis of charmed mesons [62][63][64]. If the decay properties obtained in this way describe the principal characteristics of the mixing states, the angles obtained by the potential model may be overestimated.

Simple harmonic oscillator (SHO) β values
In the next section, the Okubo-Zweig-Iizuka (OZI) allowed decays will be calculated for the 1P and 2S charmed baryons where the SHO wave functions are used to evaluate the transition factors via the 3 P 0 model. We will also discuss the mixing effects for the decays of the relevant states. Following the method of Ref. [65], all values of the SHO wave function scale, denoted β in the following, are calculated (see Table 5). The values of β reflect the distances between the light quark cluster and the c quark. In our calculation of strong decays, we will consider the structures of light diquarks. What is more important is that the possible final states of an excited charmed baryon may contain a light flavor meson, a charmed meson, a light flavor baryon, e.g., π , K, D, p, and . For the β of these hadrons, the following potential will be used: where Here, the parameters α s and b are taken as 0.45 and 0.145 GeV ν+1 as in Table 1, respectively. To reproduce the masses of the light quark clusters in Table 1, the masses of u/d, s are fixed as 0.195 and 0.380 GeV. While σ and C are treated as adjustable parameters, the masses of the π/ρ, K /K * , D/D * , p/ , and families are fitted with experimental data. Meanwhile, the values of β for the corresponding states are also obtained; they are collected in Table 6.
Before ending this section, we briefly summarize the complicated deduction presented here. Firstly, the dynamics of heavy baryon is simplified as a two-body system when the symmetric configuration is considered. Secondly, the mass matrices were calculated in the j j coupling scheme. By solving the Schrödinger equation, we obtain the mass spectra and mixing angles for the relevant states. For estimating the twobody strong decays in the next section, finally, we also present the values of the SHO wave function scale for all initial and final states.

Strong decays
In this section, we will use the formula provided by Eichten, Hill, and Quigg (EHQ) [66] to extract the decay widths of excited charmed baryons. Since the dynamical behavior of the heavy-light hadrons is governed by the light degrees of freedom in the limit of heavy quark symmetry, a doublet formed by two states with the same j cl. but different J shall have similar decay properties. More specifically, the transitions between two doublets should be determined by a single amplitude which is proportional to the products of four Clebsch-Gordan coefficients [46]. Some typical ratios of excited charmed baryons with negative parity were pre- dicted by this law [46]. Later, a more concise formula (the EHQ formula) was proposed for the widths of heavy-light mesons [66]. The EHQ formula has been applied systematically to the decays of excited open-charm mesons [62][63][64].
Recently, the EHQ formula has been extended to the study of the decay properties of 1D c and c states [11].
For the charmed baryons, the EHQ formula can be written as where ξ is the flavor factor given in Table 13 in Appendix B. q = | q| denotes the three-momentum of a final state in the rest frame of an initial state. A and B represent the initial and final heavy-light hadrons, respectively. C denotes the light flavor hadron (see Fig. 3). The explicit expression ofβ is given in Eq. (A.11) in Appendix A. In addition, C s Q , j B ,J B j C , j A ,J A is a normalized coefficient given by the following equation: where j C ≡ s C + . The symbols s C and represent the spin of the light hadron C and the orbital angular momentum relative to B, respectively. The transition factors M j A , j B j C , (q) involved in the concrete dynamics can only be calculated by various phenomenological models. For the decays of heavylight mesons, transition factors have been calculated by the relativistic chiral quark model [67] and the 3 P 0 model [62,64,68]. In our work, we will employ the 3 P 0 model [69][70][71] to obtain the transition factors. More details of an estimate of the transition factors are given in Appendix A.

Experimentally well-established 1S and 1P states
At present, all the ground states and 1P G-type charmed states have been experimentally established [1]. These states have been observed, at least, by two different collaborations, and their properties including masses and decays have been well determined. With good precision, the strong decays of these states provide a crucial test of our method.
Among the 1S charmed baryons, the measurements of c (2455) and c (2520) have been largely improved [2,3] (see Table 7). In our calculation, the mass and decay width of c (2520) ++ measured by CDF will be taken as input data to fix the constant γ peculiar to the 3 P 0 model. With the transition factor for the process c (2520) → c (2286) + π [see Eq. (A.12) in Appendix A], the value of γ is fixed as 1.296. 2 As shown in Tables 8 and 9, the predicted widths of other ground and 1P G-type charmed baryons are well consistent with experiments. Our results of mass spectra and decay widths indicate that c (2595) + , c (2625) + , c (2790) 0,+ , and c (2815) 0,+ can be accommodated with the 1P Gtype charmed baryons. c (2595) + and c (2790) 0,+ can be classified into the 1/2 − states while c (2625) + and c (2815) 0,+ into the 3/2 − states. The predicted mass splittings between the 1P 1/2 − and 3/2 − states are 25 MeV and 27 MeV for the c and c baryons, respectively, which are also consistent with the experiments. The assignments of c (2595) + , c (2625) + , c (2790) 0,+ , and c (2815) 0,+ are also supported by other work [9][10][11][12] in which the light quark cluster scenario was also employed. In addition, the mass spectra obtained by different types of the quark potential models in the three-body picture also support these assignments [7,8,[56][57][58]. However, the investigations by QCD sum rules indicate that these 1P candidates may have more complicated structures [14][15][16]. Especially, the work by Chen et al. suggested that c (2595) + and c (2625) + form the heavy doublet˜ c1 (1/2 − , 3/2 − ) (the same assignments as the case of c (2790) 0,+ and c (2815) 0,+ ) [16], which is    Table 10 is about 63.52 MeV, which is also in agreement with the measurements [30][31][32]. Furthermore, the signal of c (2765) + has been observed in the c (2455)π intermediate state, while there is no clear evidence for the decay of c (2765) + through c (2520)π [31,32]. This is also consistent with our results of the |1P, 1/2 − h . Based on the combined analysis of the mass spectrum and strong decays, we therefore conclude that c (2765) + could be regarded as a good candidate of c1 (1/2 − ). Considering the uncertainties of the quark potential models, the masses obtained by Refs. [9,10,56] are not contradictory to our assignment to c (2765) + .
According to the predicted masses in Table 3, c (2800) 0,+,++ could be assigned to either the |1P, 3/2 − l , or the |1P, 3/2 − h , or the |1P, 5/2 − states. When we consider the decay properties of these three states (see Table  10), the possibility of assignment to the |1P, 3/2 − l state can be excluded since the Belle Collaboration observed this state in the + c π mode. 3 At present, the Belle Collaboration tentatively identified c (2800) 0,+,++ as members of the c2 (3/2 − ) isospin triplet, which agrees with our results of both mass spectrum and strong decays. When the measured mass of c (2800) 0 (2806 MeV) is used for the c2 (3/2 − ) state, the predicted width is about 40.1 MeV, which is comparable with the experiment [34]. However, we notice that the quantum number J P of c (2800) 0,+,++ has not yet been measured. Then the possibility of this state as the c2 (5/2 − ) candidate cannot be excluded by our results since the decay mode of + c π is dominant for this state. In addition, the predicted mass and total width of the c2 (5/2 − ) state are also compatible with experimental data of the c (2800) 0,+,++ baryon. Therefore, we would like to point out that the signal of c (2800) 0,+,++ found by Belle might be their overlapping structure. We hope the future experiments measure the following branching ratios to disentangle this state.
For the c2 (3/2 − ) state, For the c2 (5/2 − ) state, As mentioned earlier, the signal c (2850) 0 discovered by the BaBar Collaboration may be a J = 1/2 state. If c (2850) 0 is the 1/2 + (2S) state, the corresponding ratios (see Sect. 4.5) are different from Eqs. (22)- (25). So the measurements of  Table 10). Since the decay mode of c (2520) π is the largest for the c1 (3/2 − ) state, we suggest to search this channel for this state in future experiments. In the heavy quark limit, the branching ratio for the  Table 11 favor c (2930) 0 as the c2 (3/2 − ) or c2 (5/2 − ) state. Then c (2930) 0 might be regarded as the strange partner of c (2800) 0,+,++ by our results. Interestingly, the mass difference between c (2930) 0 and c (2800) 0,+,++ is about 130 MeV, which is comparable with the mass differences among sextet states of the ground charmed baryons [21]. With a chiral quark model, Liu et al. also analyzed the c (2930) 0 by the two-body strong decays [23]. Their results support c (2930) 0 as the | 2 c P λ , 1/2 − or | 4 c P λ , 1/2 − state. Since the heavy quark symmetry was not considered in Ref. [23], the notations of charm-strange baryons in Ref [23] are different from our c0 (1/2 − ) and c1 (1/2 − ). Although the results in Table 11 indicate that the + c K decay mode dominates the decay of c0 (1/2 − ) state, the mass of this state is predicted to be about 2840 MeV, which is much smaller than 0 c (2930). In addition, the + c K decay mode is forbidden for the c1 (1/2 − ) state. Thus, according to our results, 0 c (2930) is unlikely to be a 1P state with J P = 1/2 − .
Another charm-strange baryon, c (2980) 0,+ , is slightly higher than the predicted mass range of 1P c states. This state has been observed in c (2455)K , c (2580)π , c (2645)π , and nonresonant + cK π decay channels. However, it was not seen in the decay modes of + cK and c π [37-39]. Comparing the mass and decay properties of c (2980) 0,+ with our results, the possibility of a 1P c state might be excluded. As shown in the next subsection, c (2980) 0,+ could be a good 2S c candidate. Based on our results on strong decays, we find that the c1 (1/2 − ) and c1 (3/2 − ) are quite narrow (see Table 11).   [35] 4.4 2S + c and 0,+ c states According to the mass spectrum (see Table 2), c / c (2765) + can also be regarded as the first radial (2S) excitation of the c (2286) + with J P = 1/2 + . Interestingly, the results of strong decays in Table 12 do not contradict this assignment. Our calculation indicates that the decay channel c (2455)π is a dominant decay channel for the + c (2S) state. This is in line with the observations by Belle [31,32]. At present, both 1/2 + (2S) + c and 1/2 − (1P) + c are possible for the assignment of c / c (2765) + . However, there is the very important feature for experiments to distinguish these two assignments in future. Specifically, we suggest to search c / c (2765) + in the channel of c (2520)π . As shown in Table 12, the channel c (2520)π is large enough to find the + c (2S) state. On the other hand, this mode seems to be too small to be detected for the c1 (1/2 − ) (see Table 10). Explaining the criteria concretely, we give the following branching ratios for these two states.
For the c (2S) state, For the c1 (1/2 − ) state, The branching ratio of B( c (2520)π )/B( c (2455)π ) for the c1 (1/2 − ) state is roughly an order of magnitude smaller than c (2S). If c (2765) is the 2S excitation, c (2980) could be a good candidate as its charm-strange analog [21] as seen in Table 12 is predicted for c (2980), which can be tested by future experiments. Recently, the ratio of branching fractions ≈ 75% (30) has been estimated by the Belle Collaboration [4]. Combining this with the predicted partial widths of c (2815) and c (2645) in Tables 8 and 9, the branching fraction B( c (2980) + → c (2580) 0 π + ) is evaluated to be about 40%, which is consistent with our direct result of 41.8%.  Table 3, masses of the 2S c (1/2 + , 3/2 + ) states are predicted as 2850 and 2876 MeV, respectively. The neutral c (2850) 0 found by the BaBar Collaboration in the decay channel B − → c (2850) 0p → + c π −p [35] can be regarded as the 2S c state with J P = 1/2 + . The mass and width of the neutral c (2800) 0 and c (2850) 0 are collected here: For lack of experimental information, at present, PDG treated c (2850) 0 and c (2800) 0,+,++ as the same state [1]. As pointed out by the BaBar Collaboration [35], however, there are indications that these two signals detected by Belle [34] and BaBar [35] are two different * c states. The main reasons are listed as follows: 1. Although the widths of c (2800) 0,+,++ and c (2850) 0 are consistent with each other, their masses are 3σ apart.
Our results also indicate that c (2800) 0,+,++ and c (2850) 0 are the different c excited states. One notices that the predicted mass of 1/2 + (2S) c state in this work and in Ref. [9] are around 2850 MeV. Even the results in Refs. [10,56] are only about 50 MeV larger than the measurements. Due to the intrinsic uncertainties of the quark potential model, it is appropriate to assign c (2850) 0 as a 2S 1/2 + state. More importantly, the predicted decay width of c (1/2 + , 2S) state is 118.36 MeV, which is comparable with the measurement by BaBar [35]. The partial width of c π is 35.11 MeV, which can explain why c (2850) 0 was first found in this channel. We find that the decay modes of c (2455)π and c (2520)π are also large. Finally, we give the following branching ratios: which can be tested by future experiments. If c (2850) 0 is the 1/2 + (2S) state, the mass of its doublet partner in the heavy quark effective theory is predicted as 2876 MeV (denoted c (2880)). According to the predicted decay widths in Table 12, this state might also be broad. + c π , c (2455)π , and c (2520)π are also dominant for the decay of c (2880). The ratio of ( c (2455)π )/ ( c (2520)π ) for c (2880) is different from c (2850), whose numerical value is given by Even though the strange partners of c (2850) and c (2880) have not been found by any experiments, their decay properties are calculated and presented in Table 12. Our results indicate that + c K , c (2580)π , and c (2455)K are the dominant decay modes of the c (3000) state with J P = 1/2 + , while + c K and c (2645)π are those of the c (3030). 4 Besides the masses and decay widths, the following branching ratios may also be valuable for future experiments:

Summary and outlook
In principle, both ρ and λ modes can be excited in a baryon system. For charmed baryons, the excitation energies of the ρ and λ modes are different due to the heavier mass of a c quark. For the ordinary confining potential, such as the linear or harmonic form, the excited energy of the ρ mode is larger than the λ mode [44]. Hence the low-excited charmed baryons may be dominated by the λ mode excitations. Recently, the investigation by Yoshida et al. confirmed this point [75]. Furthermore, they find that the ρ and λ modes are well separated for the charmed and bottom baryons, which means the component of the ρ mode can be ignored for the low-excited charmed baryons. Interestingly, Refs. [9][10][11] have also shown that the masses of existing charmed baryons can be explained by the λ mode. Hence, our study of strong decays of the lowexcited charmed baryons is an important complement to this work [9][10][11]75].
Up to now, several candidates of the 1P and 2S charm and charm-strange baryons have been found by experiments, and some of them are still open to debate. To better understand these low-excited charmed baryons, in this paper, we carry out a systematical study of the mass spectra and strong decays for the 1P and 2S charmed baryon states in the framework of the non-relativistic constituent quark model. The masses have been calculated in the potential model where the charmed baryons are simply treated as a quasi-two-body system in a light quark cluster picture. The strong decays are computed by the EHQ decay formula where the transition factors are determined by the 3 P 0 model. When calculating the decays, the inner structure of a light quark cluster has also been considered. Except for the unique parameter γ of the QPC model, the parameters in the potential model and in the EHQ decay formula have the same values.
The well-established ground and 1P G-type charmed baryons provide a good test to our method. The experimental properties including both masses and widths for these states can be well explained by our results. This success has made us more confident of our predictions for other 1P and 2S states. Our main conclusions are as follows.
The analysis of the mass and decay properties supports that c (2980) 0,+ is the 2S excitation (the first radial excited state of c (2470)). The existence of c (2930) 0 is still in dispute. If it exists, the assignments of c2 (3/2 − ) and c2 (5/2 − ) are possible. In other words, it could be regarded as a strange partner of c (2800) 0,+,++ . Some useful ratios of partial decay widths are also presented for c (2980) 0,+ and c (2930) 0 .
Although both the masses and the strong decays have been explained in the heavy quark-light quark cluster picture for the observed 2S and 1P candidates, it is not the end of the story of the study of the excited charmed baryon states. Investigation of the ρ mode excited states with higher energies are also important to identify the effective degrees of freedom of charmed baryons. However, this topic needs much laborious work and is beyond the scope of the present work. In addition, the quark model employed here neglects the effect of virtual hadronic loops. In the future, a more reasonable scheme for studying the properties of heavy baryons will be obtained by the unquenched quark model. Another topic which is left as a future task is to calculate the sum rules among the branching fractions of charmed baryons by applying the technique found in Ref. [76].
In the following, we will show how to obtain the partial wave amplitudes by the 3 P 0 strong decay model for the decays of excited charmed baryons. As an example, the process c (2520) → c (2280)π will be constructed and the transition factor for the EHQ formula will be extracted.
As pointed out in Sect. 3, there are two possible decay processes for an excited charmed baryon state (see Fig. 3). The final states of the left figure contain a charmed baryon and a light meson. The right one contains a charmed meson and a light baryon. If a baryon decays via the so-called 3 P 0 mechanism, a quark-antiquark pair is created from the vacuum and then regroups two outgoing hadrons by a quark rearrangement process. In the non-relativistic limit, the transition operatorT of the 3 P 0 model is given bŷ )ω (4,5) ϕ (4,5) 0 where ω (4,5) and ϕ (4,5) 0 are the color and flavor wave functions of the q 4q5 pair created from the vacuum. Thus, ω (4,5) = (RR + GḠ + BB)/ √ 3 and ϕ (4,5) 0 = (uū + dd + ss)/ √ 3 are color and flavor singlets. The pair is also assumed to carry the quantum number of 0 ++ , suggesting that they are in a 3 P 0 state. The χ (4,5) 1,−m represents the pair production in a spin triplet state. The solid harmonic polynomial Y m 1 ( k) ≡ | k|Y m 1 (θ k , φ k ) reflects the momentum-space distribution of the q 4q5 . γ is a dimensionless constant which expresses the strength of the quark-antiquark pair created from the vacuum. The value of γ is usually fixed by fitting the well-measured partial decay widths.
When the mock state [77] is adopted to describe the spatial wave function of a meson, the helicity amplitude M j A , j B , j C (q) can easily be constructed in the L S basis [71]. The mock state for an A meson is As for the left decay process in Fig. 3, the wave function of a B baryon can be constructed in the same way. The wave function of a C meson is Here, the symbols of i (i = A, B, and C) represent the Clebsch-Gordan coefficients for the initial and final hadrons, which arise from the couplings among the orbital, spin, and total angular momentum and their projection of l z and s z to j z . More specifically, i (i = A, B, and C) are given by respectively. The helicity amplitude M j A , j B , j C (q) is defined by where q represents the momentum of an outgoing meson in the rest frame of a meson A. For comparison with experiments, one obtains the partial wave amplitudes M L S (q) via the Jacob-Wick formula [78], Then the decay width (A → BC) is derived analytically in terms of the partial wave amplitudes in the A rest frame, where i = A, B, C and j = 1, 2, . . . , 5. The color matrix element ω 235 B ω 14 C |ω 45 0 ω 123 A is a constant which can be absorbed into the parameter γ . The flavor matrix element ξ = ϕ 235 B ϕ 14 C |ϕ 45 0 ϕ 123 A will be presented in the next subsection. To obtain the analytical amplitudes, the SHO wave functions are employed to describe the spatial wave function of a hadron. In the momentum space, the SHO radial wave function, ψ n Lm (q), is given by Based on Eq. (A.7), we obtain the amplitude: ; ε 1 = m 1 m 1 + m 2 ; ε 2 = m 3 m 1 + m 3 + m 5 ; where a phase space factor (2π E B E C /M A ) 1/2 γ is omitted. One notices that the unitary rotation between the L S coupling and j j coupling (Eq. (15)) should be performed to reduce the transition factors of 1P state with the same J P . More details for calculating the decay amplitudes of an excited baryon in the 3 P 0 model can be found in Refs. [24,79].