Pentaquark interpretation of (cid:2) c states in the quark model

Motivated by the recent observation of the (cid:2) c ( 2910 ) + state in the (cid:3) c ( 2455 ) 0 , ++ π ± spectrum by the Belle Collaboration, we investigate the explanation of (cid:2) c states within the pentaquark framework using the quark delocalization color screening model (QDCSM). To check for bound states and resonance states, we utilize the real-scaling method. Additionally, we calculate the root mean square of cluster spacing to study the structure of the states and further estimate whether a state is a resonance state or not. Our numerical results show that (cid:2) c ( 2910 ) cannot be interpreted as a molecular state, and (cid:3) c ( 2800 ) cannot be explained as the N D molecular state with J P = 1 / 2 − . However, we ﬁnd that (cid:2) c ( 2595 ) can be interpreted as a molecular state with J P = 12 − , where the main component is (cid:3) c π . Similarly, we interpret (cid:2) c ( 2625 ) as a molecular state with J P = 32 − , where the main component is (cid:3) ∗ c π . We also suggest that (cid:2) c ( 2940 ) is likely to be interpreted as a molecular state with J P = 3 / 2 − , where the main


Introduction
Since the discovery of + c by Fermilab in 1976 [1], the charmed baryon family has been enriched through a series of experimental collaborations . These observations have stimulated broad interest in understanding the structures of these states among theoretical groups , particularly  ), and c (2940). As these charmed baryons contain both light and heavy quarks, they provide a transitional link between light and heavy baryons, making them an important subject of study. Furthermore, the analysis of the properties of these states can deepen our understanding of the nonperturbative behavior of quantum chromodynamics (QCD).
In the last two decades, the focus of the various theoretical work and the controversial point is whether these states are excitations or multi-quark states. We can briefly review some of the theoretical work that concerning the single charmed baryons mentioned earlier. These states are studied by assuming that they are traditional three-quark excitations in the constituent quark model [28][29][30], the relativistic flux tube model [31,32], the non-relativistic quark model [33], the relativistic quark-diquark model [34][35][36][37][38], the chiral quark model [39,40], the QCD sum rules [41][42][43], the effective field theory [44,45], the effective Lagrangian method [46], the Faddeev method in momentum space [47], and so on.
In addition to the studies of masses and structures, the studies of the decay of these states were carried out in the effective meson Lagrangian [65][66][67], the 3 P 0 strong decay model [68][69][70][71], and the constituent quark model [72]. Moreover, the view of these states as a mixture of three-quark and five-quark was considered in an unquenched picture [73].
For c (2595), it was investigated as multi-quark state in Ref. [44], where c (2595) was predicted to have a predominant molecular structure with the help the effective field 2 − , 2 ρ in the framework of the chiral quark model. In Ref. [89], an unquenched picture was used to study c (2910) by considering the S−wave D * N channel coupled with the bare udc core ( c (2P)). Their results suggest that c (2910) contains a significant D * N component, and the bare state can cause the D * N binding to be more compact.
An alternative approach to study hadron-hadron interaction and the multi-quark states is the quark delocalization color screening model (QDCSM), which was developed in the 1990s with the aim of explaining the similarities between nuclear and molecular forces [90]. The model gives a good description of N N and Y N interactions and the properties of deuteron [91][92][93][94]. It is also employed to calculate the baryon-baryon and baryon-meson scattering phase shifts in the framework of the resonating group method (RGM), and the exotic hadronic states are also studied in this model. Studies also show that the N N intermediate-range attraction mechanism in the QDCSM, quark delocalization, and color screening, is equivalent to the σ -meson exchange in the chiral quark model, and the color screening is an effective description of the hidden-color channel coupling [95,96]. So it is feasible and meaningful to extend this model to investigate the charmed baryons.
In this work, we explore the pentaquark interpretation of charmed baryons in QDCSM. Both qqq −qc and qqc −qq structures, as well as the coupling of these two structures are taken into account. Our purpose is to investigate whether c (2910) could be explained as a pentaquark state and if the molecular structure is possible. In addition, we also want to see if any other bound or resonance pentaquarks exist or not, which can be used to interpret other charmed baryons.
This paper is organized as follows. After introduction, we briefly introduce the quark model and methods in Sect. 2. Then, the numerical results and discussions are presented in Sect. 3. Finally, the paper ends with summary in Sect. 4.

Theoretical framework
Herein, QDCSM is employed to investigate the properties of qqqqc systems, and the channel coupling effect is considered. In this sector, we will introduce this model and the way of constructing wave functions.

Quark delocalization color screening model (QDCSM)
The QDCSM is an extension of the native quark cluster model [97][98][99][100]. It has been developed to address multiquark systems. The detail of QDCSM can be found in Refs. [90][91][92][93][94][95]101,102]. Here, we mainly present the salient features of the model. The general form of the pentaquark Hamiltonian is given by where m i is the quark mass, p i is the momentum of the quark, and T C M is the center-of-mass kinetic energy. The dynamics of the pentaquark system is driven by a two-body potential The most relevant features of QCD at its low energy regime: color confinement (V C O N ), perturbative one-gluon exchange interaction (V OG E ), and dynamical chiral symmetry breaking (V χ ) have been taken into consideration.
Here, a phenomenological color screening confinement potential (V C O N ) is used as where a c , V 0 and μ q i q j are model parameters, and λ c stands for the SU(3) color Gell-Mann matrices. Among them, the color screening parameter μ q i q j is determined by fitting the deuteron properties, nucleon-nucleon scattering phase shifts, and hyperon-nucleon scattering phase shifts, respectively, with μ qq = 0.45 fm −2 , μ qs = 0.19 fm −2 and μ ss = 0.08 fm −2 , satisfying the relation, μ 2 qs = μ qq μ ss [103]. Besides, we found that the heavier the quark, the smaller this parameter μ q i q j . When extending to the heavy quark system, the hidden-charm pentaquark system, we took μ cc as a adjustable parameter from 0.01 fm −2 to 0.001 fm −2 , and found that the results were insensitive to the value of μ cc [104]. Moreover, the P c states were well predicted in the work of Refs. [104,105]. So here we take μ cc = 0.01 fm −2 and μ qc = 0.067 fm −2 , also satisfy the relation μ 2 qc = μ qq μ qc .
In the present work, we mainly focus on the low-lying negative parity qqqqc pentaquark states of S-wave, so the spin-orbit and tensor interactions are not included. The onegluon exchange potential (V OG E ), which includes coulomb and color-magnetic interactions, is written as where σ is the Pauli matrices and α s is the quark-gluon coupling constant. However, the quark-gluon coupling constant between quark and anti-quark, which offers a consistent description of mesons from light to heavy-quark sector, is determined by the mass differences between pseudoscalar mesons (spin-parity J P = 0 − ) and vector (spin-parity J P = 1 − ), respectively. For example, from the model Hamiltonian, the mass difference between D and D * is determined by the color-magnetic interaction in Eq. (5), so the parameter α s (qc) is determined by fitting the mass difference between D and D * .
The dynamical breaking of chiral symmetry results in the SU(3) Goldstone boson exchange interactions appear between constituent light quarks u, d and s. Hence, the chiral interaction is expressed as Among them where Y (x) = e −x /x is the standard Yukawa function. The physical η meson is considered by introducing the angle θ p instead of the octet one. The λ a are the SU(3) flavor Gell-Mann matrices. The values of m π , m K and m η are the masses of the SU(3) Goldstone bosons, which adopt the experimental values [106]. The chiral coupling constant g ch , is determined from the π N N coupling constant through Since there exists no strange quark in the qqqqc systems, the K -meson exchange potential will not influence the energies of the pentaquark systems. The other symbols in the above expressions have their usual meanings.
For the purpose of testing the stability of the calculated results against parameter changes, we used two sets of parameters in this work. All the parameters shown in Table 1 are fixed by masses of the ground baryons and mesons. Table 2 shows the masses of the baryons and mesons used in this work.

Resonating group method and wave functions
The resonating group method (RGM) [107,108] and generating coordinates method [109,110] are used to carry out a dynamical calculation. The main feature of the RGM for two-cluster systems is that it assumes two clusters are frozen inside, and only considers the relative motion between the two clusters. So, the conventional ansatz for the two-cluster wave functions is where the symbol A is the anti-symmetrization operator, and A = 1 − P 14 − P 24 − P 34 .
[σ ] = [222] gives the total color symmetry and all other symbols have their usual meanings. φ B and φ M are the q 3 andqq cluster wave functions, respectively. From the variational principle, after variation with respect to the relative motion wave function χ(R) = L χ L (R), one obtains the RGM equation: where H (R, R ) and N (R, R ) are Hamiltonian and norm kernels. By solving the RGM equation, we can get the energies E and the wave functions. In fact, it is not convenient to work with the RGM expressions. Then, we expand the relative motion wave function χ(R) by using a set of gaussians with different centers where L is the orbital angular momentum between two clusters, and S i , i = 1, 2, . . . , n are the generator coordinates, which are introduced to expand the relative motion wave function. By including the center of mass motion: the ansatz Eq. (11) can be rewritten as where χ I 1 S 1 and χ I 2 S 2 are the product of the flavor and spin wave functions, and χ c is the color wave function. These will be shown in detail later. φ α (S i ) and φ β (−S i ) are the single-particle orbital wave functions with different reference centers: With the reformulated ansatz Eq. (15), the RGM Eq. (12) becomes an algebraic eigenvalue equation: where H i, j and N i, j are the Hamiltonian matrix elements and overlaps, respectively. By solving the generalized eigen problem, we can obtain the energy and the corresponding wave functions of the pentaquark systems.
In QDCSM, in order to take into account the mutual distortion or the internal excitations of nucleons in the course of interaction, quark delocalization was introduced to enlarge the model variational space. It is realized by specifying the single particle orbital wave function of QDCSM as a linear combination of left and right Gaussians. The single particle orbital wave functions used in the ordinary quark cluster model are It is important to note that, the mixing parameter is not an adjustable parameter, but is instead determined by the dynamics of the multi-quark system itself. In this way, the multi-quark system could choose its favorable configuration in the interacting process. The intermediate-range attraction and the structure of the formed state can be revealed by studying the variation of delocalization parameter . If the value of the parameter is still large at large distance between two cluster, which indicates that the quarks are willing to run between two clusters, then the state is inappropriate to be explained as a molecular state. Conversely, if the value of the parameter is close to zero at large distance between two cluster, then the formed state is likely to be a molecular state. This mechanism has been used to explain the cross-over transition between hadron phase and quark-gluon plasma phase [111]. For the spin wave function, we first construct the spin wave functions of the q 3 andqq clusters with SU(2) algebra, and then the total spin wave function of the pentaquark system is obtained by coupling the spin wave functions of two clusters together. The spin wave functions of the q 3 andqq clusters are Eq. (19) and Eq. (20), respectively For pentaquark system, the total spin quantum number can be 1/2, 3/2 or 5/2. Considering that the Hamiltonian does not contain an interaction that can distinguish the third component of the spin quantum number, so the wave function of each spin quantum number can be written as follows Similar to constructing spin wave functions, we first write down the flavor wave functions of the q 3 clusters, which are Eur. Phys. J. C (2023) 83:524 Then, the flavor wave functions ofqq clusters are As for the flavor degree of freedom, the isospin of pentaquark systems we investigated in this work is I = 0. The flavor wave functions of pentaquark systems can be expressed as For the color-singlet channel (two clusters are colorsinglet), the color wave function can be obtained by 1 ⊗ 1 Fig. 1 The shape of the resonance state in real-scaling method Finally, we can acquire the total wave functions by combining the wave functions of the orbital, spin, flavor and color parts together according to the quantum numbers of the pentaquark systems.

Stabilization method
To provide the necessary information for experiments to search for exotic hadron states, the coupling calculation between the bound channels and open channels is indispensable. As we know that some channels are bound because of the strong attractions of the system. However, these states will decay to the corresponding open channels by coupling to open channels and become resonance states. Besides, some states will become scattering state by the effect of coupling to both the open and closed channels. The stabilization method, which is also called the real-scaling method, is one of the effective ways to look for the genuine resonance states. This method was proved to be a reliable and mature approach for estimating the energies of long-lived metastable states of electron-atom, electron-molecule, and atom-diatom complexes [112]. It was firstly applied to quark model by Emiko Hiyama [114] to search for P c states.
In this approach, the distance between baryon and meson clusters is labeled as S i , and the largest one is S m . We scale the finite volume for expansion by increasing S m continuously. The continuum states (scattering states) will fall down to the corresponding thresholds, while the genuine resonance states will survive after coupling to the scattering states and keep stable with increasing S m . In this way, there will be two cases of the genuine resonance states, which are shown in Fig. 1. (a) If the energy of a scattering state is far away from that of the resonance state which means there is a weak coupling (or no coupling) between the resonance states and the scattering states, the shape of the resonance state is a stable straight line. (b) If the energy of a scattering state is close to the resonance state, there will be a strong coupling between the resonance state and the scattering state, which would appear as an avoid-crossing structure between two declining lines. The avoid-crossing structure will reproduce stably with the continual increase of S m . In this case, one can obtain the possible resonance state by the correspond- Table 3 The energies of the qqqqc pentaquark system with quantum numbers J P = 1 2 − (unit: MeV) Parameter Structure ing avoid-crossing structure. Furthermore, by analyzing the components of avoid-crossing structures, we can determine the formation mechanism of the resonance state. The main component of the avoid-crossing structure is usually bound in single-channel calculation and the avoid-crossing structure is below the threshold of the main component, which can account for the mechanism of the resonance state. Otherwise, if the main component of the avoid-crossing structure is unbound in the single-channel calculation, then it is usually a false resonance state. By this mean, we can distinguish resonance states from scattering states. More details can be found in Refs. [113][114][115][116].

The results and discussions
In this work, we investigate the S-wave qqqqc pentaquark systems in the framework of QDCSM using two sets of parameters. The quantum numbers of these systems are I = 0, J P = 1/2 − , 3/2 − and 5/2 − . Two structures qqq−qc and qqc −qq, as well as the coupling of these two structures are taken into consideration. To find out if there exists any bound state, we carry out a dynamic bound-state calculation. The real-scaling method is employed to investigate the genuine resonance states. Moreover, the calculation of root mean square (RMS) is helped to explore the structure of the bound states or resonance states on the one hand, and to further estimate whether the observed states are resonance states or scattering states on the other hand. The numerical results of different systems are listed in Tables 3, 5 and 7, respectively. The first column headed with Parameter represents the results of two sets of parameters. The second column headed with Structure includes qqq −qc and qqc −qq two kinds. The third and forth columns, headed with χ f i and χ σ j , denote the way how wave functions constructed, which can be seen in Sect. 2.2. The fifth column headed with Channel gives the physical channels involved in the present work. The sixth column headed with E th (T heo.) refers to the theoretical value of non-interacting baryon-meson threshold. The seventh column headed with E sc shows the energy of each single channel. The values of binding energies (E B = E sc − E th (T heo.)) are listed in the eighth column only if E B < 0 MeV. E sc is the corrected energy after the mass correction, which will be introduced in detail later. E ccs is the energy of the channel coupling of single spatial structure, while E cct takes two spatial structures into account in channel coupling. E ccs and E cct are the corrected energies, corresponding to different coupling modes. In addition, the proportion of each channel in the channel coupling calculation is shown at the bottom of Tables 3, 5 and 7. With the help of the proportion of each channel, we can further study the influence of the channel coupling.
To reduce the theoretical errors, we can shift the mass to E sc = M 1 (E x p.) + M 2 (E x p.) + E B , where the experimental values of a baryon M 1 (E x p.) and a meson M 2 (E x p.) are used. The above formula is used to handle single channel mass correction. When we deal with the mass correction of the coupled channels, a modified formula will be used.
where p i is the proportion of various physical channels.

J
The energies of the systems with J P = 1/2 − are listed in Table 3 that none of these channels is bound. The energies of c ρ and * c ρ channels are below the corresponding theoretical threshold, with the binding energies −3 to −6 MeV and −13 to −15 MeV, respectively. However, the channels of the system are influenced by each other, so it is unavoidable to take into account the channel coupling effect.
To give a better understanding of channel coupling, we couple the channels with the same spatial structure, as well as all the channels with two spatial structures. From Table 3, we can see that for the qqq −qc structure, a bound state with the corrected mass 2809-2810 MeV is obtained, which is close to the charmed baryon c (2800). Similar conclusions can be found in Refs. [52,53,57,58,65]. However, this state can decay to the c π channel, which has a lower energy. For the qqc −qq structure, the lowest energy of this system is 2614-2615 MeV, which is 8-9 MeV lower than the lowest energy channel c π . So a bound state with the main component c π is obtained here. However, it is unclear whether the bound single channels c ρ and * c ρ become resonance states or not. So the channel coupling of all channels is needed and the real-scaling method is employed to explore the resonance states. By solving the Schrodinger equation with channel coupling, we can obtain a series of eigenvalues theoretically. All these eigenvalues are used to draw for real-scaling method, although some extreme high energies are not visible in the images that finally presented. On the other hand, only the lowest energy is presented in Table 3, because whether the system can form a bound state depends on whether the lowest energy is below the lowest threshold.
We calculate the energy eigenvalues of the qqqqc system by taking the value of S m from 6.0 to 12.0 fm, to see if there is any stable state. The results of the two sets of parameters are consistent with each other. To save space, we only show the results of PARAM I here. The stabilization plots of the energies of the qqqqc system with quantum numbers J P = 1/2 − are shown in Fig. 2. The continuum states fall off towards their respective thresholds, which are marked with red lines. When the distance between baryon and meson clusters is small, the threshold structure is not obvious enough at close distance. This phenomenon is caused by a lack of computing space and this situation will improve with the increase of S m . And for genuine resonance states, which appear as avoid-crossing structure will be marked with blue lines. Bound states below the lowest threshold of the systems are also marked with blue lines. In addition, the main component of the avoid-crossing structures is usually bound in single-channel calculations, which can account for the mechanism of the resonance state.
The RMS can also be used to further estimate whether the observed states are resonance states or scattering states. It is worth noting that, the scattering states have no real RMS since the relative motion wave functions of the scattered states are non-integrable in the infinite space. If we calculate the RMS of scattering states in a limited space, we can only obtain a value that increases with the expansion of computing space. Although the wave function of a resonance state is also nonintegrable, we can calculate the RMS of the main component of the resonance state, which is usually bound in singlechannel calculation and its wave function is integrable. In this way, we can calculate the RMS of various states to identify the nature of these states by keep expanding the computing space. Besides, the structure of a multi-quark system can also be estimated by calculating the RMS. The results of RMS of the single channel and channel-coupling are listed in the Table 4.
After taking into account the full channel coupling, the lowest energy of the J P = 1/2 − system is pushed down to 2598-2599 MeV, and the corrected mass is 2575 MeV, as shown in Table 3. After calculating the composition, we find that the channel of c π is the main component with the proportion of 90.0%, while the channel of N D and the rest channels are the minor components, accounting for 8.1% and 1.9%, respectively. From Fig. 2 we can see that the energy of this state is very stable with the increase of S m , which confirms that it is a bound state. Besides, the RMS of this state is stable with the increase of the computing space, which further confirms the bound state conclusion. The value of the RMS of this state is 1.4 fm, indicating that the two clusters  are not too close to each other. As discussed in Sect. 2.2, the delocalization parameter can be used to reveal the internal structure of multi-quark state. The variation of delocalization parameter of different states are shown is Fig. 3. To be clear, the variation of delocalization parameter for c π with J P = 1/2 − and * c π with J P = 3/2 − are similar and that causes these two lines to overlap. The delocalization parameter of c π with J P = 1/2 − is quite small around 1.4 fm, which means that the quarks are not likely to run between baryon and meson clusters. All these properties show that the bound state with J P = 1/2 − is inclined to be a molecular state, and its mass is close the c (2595). Since the c (2595) is located very close to the c π threshold, this observation leads us naturally to consider a predominant baryon-meson structure of this lowest-lying odd parity charmed baryon. Here, we prefer to interpreted c (2595) as the molecular state with the main component of c π , and the quantum number is J P = 1/2 − . The similar conclusions can be seen in Refs. [44,52,62].
As mentioned above, a quasi-bound state is obtained by coupling N D and N D * channels. From Fig. 2 we can see that the avoid-crossing structures appear at the place around the threshold of N D. However, it is difficult to estimate if there is a resonance state, because the avoid-crossing structure is too close to the threshold of N D. We calculate the RMS of the single N D and N D * channel, as well as the case by channel-coupling (labeled as E ccs (2776)). We find that the values of RMS are very large and they are not stable with the increase of the computing space, which indicates that they are scattering states. So, we conclude that the avoidcrossing structure around the threshold of N D or N D * is not a resonance state. It is because that the decay rate of different channels to open channels is different, which leads to the different slope of the dots representing energies in the real-scaling figure and forms the avoid-crossing structure. So, we cannot explain c (2800) as the N D molecular state with J P = 1/2 − in present work.
It is worth noting that in the middle of Fig. 2, there exists an avoid-crossing structure around 2905 MeV, the mass of which is close to the newly observed c (2910). We're interested in whether it a resonance state and can be used to explain the c (2910). At the bottom of Table 4, we present the proportion of each channel of this E(2905) and its main components are N D * and N D. Although E(2905) is above the thresholds of N D * and N D, it is still possible to be a color-structure resonance state because the effect of the hidden-color channel-coupling is included in the QDCSM. However, the value of RMS of this state is very large and it is variational with the increase of the computing space, which indicates that it is a scattering state. So it cannot be used to explain the c (2910) in this work.
For the c ρ and * c ρ, the single channel calculation shows that both of them are bound states, and the value of RMS of each channel is also consistent with this conclusion. However, after the full channel-coupling, we cannot find any avoid-crossing structure below the threshold of c ρ or * c ρ. It is reasonable. There are several channels below c ρ and * c ρ, which will push the energy of these two states above the thresholds. So, these two bound states disappear after taking into account the effect of channel-coupling.
The energies of qqqqc pentaquark system with quantum numbers J P = 3 2 − are listed in the Table 5. For the qqq −qc spatial structure, the N D * channel is bound in the single channel calculation. At the same time, for the qqc−qq struc- ture, the single channel calculation shows that both the c ρ and * c ρ are bound states, with the binding energy of −77 to −79 MeV and −14 to −16 MeV, respectively, while the c ω and * c π channels are unbound. After coupling all possible channels, a bound state is obtained, whose energy is 2624-2628 MeV (2635-2636 MeV after mass correction). The proportion of each channel is * c π : 95.6%, c ω : 2.9% and the rests: 1.5%, which means that the main component of this bound state is * c π . Since the results of the two sets of parameters are consistent, the following discussion is mainly based on PARAM I. Figure 4 and Table 6 show the stabilization plots of the energies and the RMS of the qqqqc system with J P = 3/2 − , respectively. It is obvious in Fig. 4 that there is a stable state under the lowest threshold, which is marked by the blue line. The value of RMS of this state is stable with the increase of the computing space and it is 1.4 fm, indicating that the two clusters are not too close to each other. Meanwhile, from Fig. 3, the delocalization parameter of * c π with J P = 3/2 − is quite small around it's RMS. All these properties show that the bound state with J P = 3/2 − tends to be a molecular state, and its mass is close the c (2625). So it is possible to interpreted the c (2625) as a molecular state with J P = 3/2 − dominated by * c π channel. The similar explanation could be found in Refs. [48,49,52].
In Fig. 4, there are five red lines, which represent the thresholds for each single channel. A little below the threshold line of N D * , we can see that stable avoid-crossing structures repeated periodically there. After calculating the composition, we find that the main component is the N D * channel with the proportion of 66.5%, while the proportion of the * c π channel is 25.1% and the one of rest channels is 8.4%. The theoretical energy of this structure is 2849 MeV, lower than the threshold of N D * . So these avoid-crossing structures may represent a resonance state. By using the proportion of each channel, the corrected mass 2933 MeV is obtained for this state. Besides, the calculation of the RMS of this state shows that it is stable with the increase of the computing space and the RMS is 1.9 fm, which indicates that it is a resonance state with the molecular structure. The delocalization parameter of N D * from Fig. 3 also confirms the molecular structure of N D * state. Clearly, the corrected mass of this resonance state is close to the c (2940). So the c (2940) is likely to be interpreted as a molecular state with J P = 3/2 − , and the main component is N D * . This conclusion is consistent with the work of Refs. [50,51,[53][54][55][56][57]60,66,67].
Particularly, there is another repeated avoid-crossing structure below the threshold of c ρ. The calculated mass of this state is 3160 MeV, and the main component is the c ρ channel with the proportion of about 72%. The corrected mass is 3140 MeV, and the RMS of this state is 1.4 fm. All these properties show that it is also a resonance state, which is worth searching in future work. Besides, the delocalization parameter of this state around 1.4 fm is not very small, as shown in Fig. 3. This indicates that quarks are likely to run between c and ρ clusters, so c ρ with J P = 3/2 − is improper to be a molecular state.
As for PARAM II, the results of real-scaling method and RMS are consistent with that of PARAM I, and two resonance states are obtained. the corrected mass of the N D * and the c ρ resonance states are 2933 MeV and 3142 MeV, respectively.    to interpret this state directly as a molecular state. Around 1.4 fm, delocalization parameter of this state is still large, as shown in Fig. 3. This indicates a strong attraction between * c and ρ despite the distance between two hadrons. On the basis of the above discussion, the * c ρ state is not suitable to be a molecular state. Although it can decay to some D-wave channels, like N D, N D * , c ω, c ρ, and so on, it is still possible to be a resonance, which is worthy of experimental search and research.

Summary
In this work, we try to explain some charmed baryons within the pentaquark framework. The S-wave pentaquark systems qqqqc with I = 0, J P = 1 2 − , 3 2 − and 5 2 − are systematically investigated in the QDCSM. The dynamic bound state calculation is carried out to search for any bound state in the qqqqc systems. Both the single channel and the channel coupling calculation are performed to explore the effect of the multi-channel coupling. Meanwhile, the real-scaling method is employed to examine the existence of the resonance states and the bound states. We also calculate the RMS of cluster spacing to study the structure of the states and estimate if the state is resonance state or not.
The numerical results show that the effect of the channel coupling is important for forming a bound state and deepening the bondage to some extent. We can draw the following conclusions: (1) Three bound states are obtained in present work, among which c (2595) can be interpreted as the molecular state with J P = 1  icated, whose mass is 3140-3142 MeV. All these charmed states are worth searching in future work. In describing the multi-quark system, the channel coupling effect has to be taken into account, especially for the resonance state, where the coupling to the open channels will shift the mass of the resonance state, or destroy it. The real-scaling method may be an effective method to pick up the genuine resonance states from the states with discrete energies. Besides, from the above discussion of the charmed baryons, we would like to note that there exist different points of view to the structure of these states. To explore the structure of exotic hadrons, the unquenched quark model may be another critical approach.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: The data have been illustrated in the figures and tables, so they are not necessary to be deposited. Data may be made available upon request.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 . SCOAP 3 supports the goals of the International Year of Basic Sciences for Sustainable Development.