Charmed baryon–nucleon interaction

We present a comparative study of the charmed baryon–nucleon interaction based on different theoretical approaches. For this purpose, we make use of (i) a constituent quark model tuned in the light-flavor baryon–baryon interaction and the hadron spectra, (ii) existing results in the literature based both on hadronic and quark-level descriptions, (iii) (2+1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2+1$$\end{document})-flavor lattice QCD results of the HAL QCD Collaboration at unphysical pion masses and their effective field theory extrapolation to the physical pion mass. There is a general qualitative agreement among the different available approaches to the charmed baryon–nucleon interaction. Different from hadronic models based on one-boson exchange potentials, quark-model based results point to soft interactions without two-body bound states. They also support a negligible channel coupling, due either to tensor forces or to transitions between different physical channels, ΛcN\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda _c N $$\end{document}–ΣcN\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \Sigma _c N$$\end{document}. Short-range gluon and quark-exchange dynamics generate a slightly larger repulsion in the 1S0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^1S_0$$\end{document} than in the 3S1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^3S_1$$\end{document}ΛcN\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda _c N$$\end{document} partial wave. A similar asymmetry between the attraction in the two S waves of the ΛcN\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda _c N$$\end{document} interaction also appears in hadronic approaches. A comparative detailed study of Pauli suppressed partial waves, as the 1S0(I=1/2)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^1S_0 (I=1/2)$$\end{document} and 3S1(I=3/2)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^3S_1 (I=3/2)$$\end{document}ΣcN\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Sigma _c N$$\end{document} channels, would help to disentangle the short-range dynamics of two-baryon systems containing heavy flavors. The possible existence of charmed hypernuclei is discussed.


Introduction
There has been an impressive experimental progress in the spectroscopy of heavy hadrons, mainly in the charm sector. The theoretical analysis of hidden and open heavy flavor hadrons has revealed how interesting is the interaction of heavy hadrons, with presumably a long-range part of Yukawa type, and a short-range part mediated by quark-quark and a e-mail: humberto@esfm.ipn.mx b e-mail: valcarce@usal.es c e-mail: carames@usal.es quark-antiquark forces. Some of the recently reported states might appear as bound states or resonances in the scattering of two hadrons with heavy flavor content. See Refs. [1][2][3][4][5][6] for recent overviews and discussions. Thus, the understanding of the baryon-baryon interaction in the heavy flavor sector is a key ingredient in our quest to describing the properties of hadronic matter.
The research programs at various facilities are expected to improve our knowledge on the hadron-hadron interactions involving heavy flavors, particularly in the charm sector. Thus, the LHCb Collaboration at the Large Hadron Collider (LHC) is engaged in an extensive program aimed at the analysis of charmed hadrons produced in the environment of high-energy proton-proton collisions [7]. The observation of five new narrow excited c states has already been reported [8], some of which are suggested as molecules containing a charmed hadron [1][2][3][4][5][6]. The planned installation of a 50 GeV high-intensity proton beam at Japan Proton Accelerator Research Complex (J-PARC) [9,10] intends to produce charmed hypernuclei, in which a Y c baryon ( c or c ) is bound to a nucleus. There are also planned experiments by the PANDA Collaboration at the Facility for Antiproton Ion Research (FAIR) [11,12] to produce charmed hadrons by annihilating antiprotons on nuclei.
In addition to the recent interest in the hadron-hadron interaction involving heavy flavors, there is a long history of speculations as regards bound nuclear systems with a charmed baryon. The observation of events that could be interpreted in terms of the decay of a charmed nucleus [13,14], fostered conjectures about the possible existence of charm analogs of strange hypernuclei [15][16][17]. This resulted in several theoretical estimates about the binding energy and the potential-well depth of charmed hypernuclei based on one-boson-exchange potentials for the charmed baryon-nucleon interaction [18][19][20][21][22]. The current experimental prospects have reinvigorated studies of the lowenergy Y c N interactions [23][24][25][26][27][28][29][30][31][32]. See also the recent reviews [33,34].
As pointed out by Bjorken [35] one should strive to study systems with heavy flavors because due to their size the quark-gluon coupling constant is small and therefore the leading term in the perturbative expansion is enough to describe the system. However, our ability of making firstprinciples analytical calculations of nonperturbative QCD phenomena is very limited. When combined with the lack of experimental information on the elementary Y c N interactions there is room for some degree of speculation in the study of processes involving charmed hadrons. Thus, the situation can be ameliorated with the use of well-constrained models based as much as possible on symmetry principles and analogies with other similar processes, which is still a valid alternative for making progress.
Within such a perspective, in this work we present the first comparative study of the charmed baryon-nucleon interaction based on different theoretical approaches. We employ a widely used constituent quark model (CQM) [36,37] providing a good description of the low-lying spectrum of light and charmed hadrons [38,39] as well as the nucleon-nucleon interaction [36,40]. In addition, we consider different scattered results available in the literature. In particular, we compare to the hadronic description based on one-bosonexchange potentials of Ref.
[26]; the quark-level approach relying on the quark delocalization color screening model (QDCSM) of Ref. [28]; the hybrid model of Ref. [30] based on one-boson-exchange potentials supplemented by a global short-range repulsion of quark origin; and the recent charmed baryon-nucleon potential based on a SU(4) extension of the meson-exchange hyperon-nucleon potential Ã of the Jülich group [41] of Ref. [32]. We will also consider the recent lattice QCD simulations of the Y c N interactions by the HAL QCD Collaboration [42][43][44][45]. However, the lattice QCD simulations are still obtained with unphysical pion masses. They have been extrapolated to the physical pion mass using a chiral effective field theory (EFT) [46]. The paper is organized as follows. In Sect. 2 we outline the basic ingredients of the CQM used to derive the Y c N interactions. We also describe the integral equations of the coupled c N − c N system. In Sect. 3 we present and discuss the results for the c N and c N interactions. We show the results of the CQM in comparison to the available results from other theoretical approaches in the literature. We analyze the consequences of the different approaches for the possible existence of charmed hypernuclei. Finally, in Sect. 4 we summarize the main conclusions of our work.

The quark-quark interaction
The two-body Y c N interactions are obtained from the chiral constituent quark model of Ref. [36]. The model was pro-posed in the early 1990s in an attempt to obtain a simultaneous description of the light baryon spectrum and the nucleon-nucleon interaction. It was later on generalized to all flavor sectors [37]. In this model, hadrons are described as clusters of three interacting massive (constituent) quarks. The masses of the quarks are generated by the dynamical breaking of the original SU(2) L ⊗ SU(2) R chiral symmetry of the QCD Lagrangian at a momentum scale of the order of CSB = 4π f π ∼ 1 GeV, where f π is the pion electroweak decay constant. For momenta typically below that scale, when using the linear realization of chiral symmetry, light quarks interact through potentials generated by the exchange of pseudoscalar Goldstone bosons (π ) and their chiral partner (σ ): where g 2 ch /4π is the chiral coupling constant, m i are the masses of the constituent quarks, Perturbative QCD effects are taken into account through the one-gluon-exchange (OGE) potential [47]: where λ c are the SU(3) color matrices, r 0 =r 0 /ν is a flavordependent regularization scaling with the reduced mass ν of the interacting pair, and α s is the scale-dependent strong coupling constant given by [37], where α 0 = 2.118, μ 0 = 36.976 MeV and γ 0 = 0.113 fm −1 . This equation gives rise to α s ∼ 0.54 for the lightquark sector, α s ∼ 0.43 for uc pairs, and α s ∼ 0.29 for cc pairs. Finally, any model imitating QCD should incorporate confinement. Although it is a very important term from the spectroscopic point of view, it is negligible for the hadron-hadron interaction. Lattice QCD calculations suggest a screening effect on the potential when increasing the interquark distance [48] which is modeled here by where a c and μ c are the strength and range parameters. Once perturbative (one-gluon exchange) and nonperturbative (confinement and dynamical chiral symmetry breaking) aspects of QCD have been incorporated, one ends up with a quarkquark interaction of the form where n stands for the light quarks u and d. Notice that for the particular case of heavy quarks (c or b) chiral symmetry is explicitly broken and therefore boson exchanges associated to the dynamical breaking of chiral symmetry do not contribute. The parameters of the model are the ones used for the study of the light one-and two-hadron systems [36][37][38][39][40], and for completeness they are quoted in Table 1.
In order to derive the B n B m → B k B l interaction from the basic qq interaction defined above, we use a Born- Oppenheimer approximation where the quark coordinates are integrated out keeping R fixed, the resulting interaction being a function of the two-baryon relative distance. A thorough discussion of the model can be found elsewhere [36,37,49]. We show in Fig. 1 the different diagrams contributing to the charmed baryon-nucleon interaction. While diagrams (a) and (b) are considered in a hadronic description, diagrams (c) and (d) correspond to short-range effects due to quark exchanges that are not mapped in a hadronic description. Diagrams (c) and (d) contain one-gluon exchange contributions that are also missed in hadronic models. To illustrate the capability of the model let us just mention how the obtained N N potentials perfectly describe the S wave phase shifts [40].
In the limit where the two baryons Y c N overlap, the Pauli principle may impose antisymmetry requirements not present in a hadronic description. Such effects, if any, will be prominent for relative S waves, L = 0. The S wave normalization kernel of the two-baryon wave function can be written in the overlapping region (R → 0) as [49] N L=0SI where C(S, I ) is a spin-isospin coefficient and b and b c are the Gaussian parameters for the wave function of the light and charmed quarks, respectively, assumed to be different for the sake of generality. The closer the value of C(S, I ) to 1/3 the larger the suppression of the normalization of the wave function at short distances, generating Pauli repulsion [49,50]. Similarly to Pauli blocked channels, corresponding to C(S, I ) = 1/3, there might exist Pauli suppressed channels, those where C(S, I ) is close to 1/3. This is the case for the c Interaction between two light quarks together with the exchange of identical light quarks. d Interaction between the heavy and a light quark together with the exchange of identical light quarks repulsion. As we will discuss below, this repulsion will be reflected in the phase shifts. Let us finally note that, although we will discuss the dependence of the results on different values of b c , we take a reference value of b c = 0.5 fm.

The coupled c Nc N system
If we consider the system of two baryons Y c and N in a relative S state interacting through a potential V that contains a tensor force, then there is a coupling to the Y c N D wave so that the Lippmann-Schwinger equation of the system is of the form where t is the two-body amplitude, J , I , and E are the total angular momentum, isospin and energy of the system, and s, s , s are the initial, intermediate, and final orbital angular momentum and spin. p and μ are, respectively, the relative momentum and reduced mass of the two-body system. More precisely, Eq. (8) is only valid for the c N system with isospin 3/2. For this case, the coupled channels of orbital angular momentum and spin that contribute to a given state with total angular momentum J are found in the first two rows of Table 2.
In the case of isospin 1/2, the c N states are coupled to c N states. Thus, if we denote the c N system as channel c and the c N system as channel c , instead of Eq. (8) the Lippmann-Schwinger equation for c N − c N scattering with isospin 1/2 becomes  (2,1) where t c c ;J I is the c N → c N scattering amplitude, t c c ;J I is the c N → c N scattering amplitude, and t c c ;J I is the c N → c N scattering amplitude. The propagators G c (E; p c ) and G c (E; p c ) in Eq. (9) are given by with where the on-shell momenta k c and k c are related by We give in Table 2 the channels ( c , s c ) and ( c , s c ), corresponding to the c N and c N systems, which are coupled in a given state of total angular momentum J for the case of isospin 1/2.

c N interaction
We show in Fig. 2a the phase shifts for the c N 1 S 0 partial wave as a function of the center of mass (c.m.) kinetic energy. The latest (2 + 1)-flavor lattice QCD simulations by the HAL QCD Collaboration [44] 1 for a pion mass of 570 (410) MeV are denoted by the blue (red) filled circles with their corresponding errors shown by the vertical lines. The black solid line stands for the results of the CQM described in Sect. 2.1. The black dashed line corresponds to the results of the QDCSM of Ref. [28] for a color screening parameter μ = 0.1. In Fig. 2b we present the phase shifts for the c N 3 S 1 partial wave-note that in this case the results of the QDCSM model of Ref. [28] are not available. As can be seen there is a tendency that the attraction obtained by the latest lattice QCD simulations for both c N S waves becomes stronger as the pion mass decreases, moving towards the predictions of the CQM and QDCSM models.
In Fig. 3 we show the dependence of the scattering lengths of the spin-singlet and spin-triplet c N partial waves The solid blue and red lines are just a guide to the eye reported by the HAL QCD Collaboration as a function of the pion mass. The purple diamonds at the physical pion mass stand for the results of the CQM. The repulsive or attractive character of the interaction for the different Y c N partial waves in the CQM is reflected in the scattering lengths and effective range parameters summarized in Table 3.
As can be seen in Figs. 2 and 3, the phase shifts and scattering lengths of the c N 1 S 0 and 3 S 1 partial waves derived by the HAL QCD Collaboration are qualitatively and quantitatively rather similar. Indeed, it was noted in Ref. [44] that the corresponding 1 S 0 and 3 S 1 potentials are almost identical at (a) (b) Fig. 3 a Dependence of the HAL QCD c N 1 S 0 scattering length on the pion mass [44]. The vertical bars include statistical and systematic errors. The purple diamond represents the prediction of the CQM for the physical pion mass. b Same as (a) for the c N 3 S 1 partial wave 410 MeV pion mass and at 570 MeV. These potentials show that the c N interaction is attractive but not strong enough to form two-body bound states. The results of the CQM are slightly different: both partial waves are attractive but without developing two-body bound states. However, the 3 S 1 partial wave is more attractive than the 1 S 0 . This result is due to the short-range dynamics discussed in Sect. 2.1, a consequence of gluon and quark exchanges. It has been outlined long ago in the literature for the N system [51]. If no meson exchanges were considered, the S wave phase shifts of the c N system are very similar to the correspond-ing N N scattering [52]. In both partial waves one obtains typical hard-core phase shifts due to the short-range gluon and quark-exchange dynamics. However, the hard-core radius in the spin-singlet state is larger than in the spin-triplet one [51] leading to a more attractive interaction in the spin-triplet partial wave due to a lower short-range repulsion [53]. In fact, the hard cores caused by the color-magnetic part of the OGE potential have been calculated in Ref. [51], obtaining 0.35 fm for the spin-triplet state and 0.44 fm for the spin-singlet one. If the short-range dynamics is properly considered, this effect has to be transferred to the phase shifts, as concluded by the CQM. This difference stems from the different expectation value in the spin-singlet and spin-triplet c N partial waves of the color-magnetic operator appearing in Eq. (3), σ i · σ j λ c i · λ c j . The matrix elements of this operator are only different from zero when there are quark-exchange effects, as depicted in diagrams (c) and (d) of Fig. 1, 2 giving rise to a genuine quark substructure effect not mapped at the hadronic level.
Reference [44] discusses the qualitative difference between the N and c N interactions due to the absence of K -meson exchanges. The origin of the small spin dependence of the c N interaction is attributed to the heavy D meson mass and the large separation between the c N and c N masses. However, no discussion is found of the role of the short-range dynamics that may contribute to the different behavior of the spin-singlet and spin-triplet c N phase shifts. As will be discussed below, the short-range dynamics also generates a major impact in the c N charmed baryon-nucleon interaction. This is due to additional Pauli suppression, as discussed in Sect. 2.1, in the 1 S 0 (I = 1/2) and 3 S 1 (I = 3/2) c N partial waves, resulting in a strong repulsion.
Recently, Ref. [32] has presented a charmed baryonnucleon potential based on a SU(4) extension of the mesonexchange hyperon-nucleon potential Ã of the Jülich group [41]. Three different models of the interaction were considered, which differ only on the values of the couplings of the scalar σ meson with the charmed baryons. In particular, in a first model the couplings of the σ meson with the charmed baryons are assumed to be equal to those of the and hyperons, and their values are taken from the original YN potential Ã of the Jülich group. In the other two models these couplings are reduced by 15% and 20%, respectively. The c N phase shifts obtained with these models are in qualitative agreement with the CQM results. They predict a higher overall attraction for the 3 S 1 than for the 1 S 0 c N partial wave, unlike the HAL QCD results, predicting similar phase shifts for both partial waves.
There are other studies of the Y c N interactions based on one-boson-exchange potentials at hadronic level [26,27]. Although they do not report explicitly phase shifts or scatter-ing lengths, binding energies of the Y c N two-body systems as a function of the boson-exchange cutoff π are calculated. As can be seen in Tables III and IX of Ref. [26] the J P = 0 + and J P = 1 + states are bound for any value of π . The binding energies of the J P = 1 + state are always a little bit larger than those of the J P = 0 + state. This is due to the similar contribution of the boson-exchange potentials in both partial waves, the difference coming from the channel coupling that enhances the D wave probability. Thus, while for π = 1.2 GeV the probability of the 1 S 0 c N channel in the J P = 0 + state is 98.2%, that of the 3 S 1 c N channel in the J P = 1 + state is 97.6%, with a D wave probability of 1.8%. The small difference between the 1 S 0 and 3 S 1 probabilities in the J P = 0 + and 1 + states, remains almost constant for any value of π . For example, for π = 1.6 GeV they are 80.1% and 79.6%, respectively. However, the D wave probability in the J P = 1 + state augments from 1.8 to 10.1%. Table IV of Ref. [26] reports binding solutions for the individual channels in the J P = 0 + state. As can be seen, the uncoupled 1 S 0 c N state is bound for any value of the cutoff. Unfortunately, binding solutions for the uncoupled c N channel in the J P = 1 + state are not reported. A simplest guess-by-analogy estimation tells us that the results would be the same in both J P states if channel coupling was not considered, as happens for the CQM if the short-range dynamics is neglected.
In a later work [30], the hadron level one-boson-exchange potential was supplemented by an overall short-range repulsion arising from color-magnetic effects evaluated in the heavy quark limit [54][55][56][57]. In general, the results are similar to their previous study, both states J P = 0 + and J P = 1 + being bound or at the edge of binding and obtaining larger binding energies in the 1 + state for the same parametrization. Hence, in both cases [26, 27,30] one expects phase shifts close to 180 degrees at zero energy, being larger for the spintriplet partial wave.
The phase shifts for the 1 S 0 c N interaction reported by the QDCSM model of Ref. [28], dashed line in Fig. 2a, are more attractive than those of the CQM model, although they do not show a bound state. A major difference between the quark model and hadron level approaches has to do with the strength of the channel coupling. The c Nc N transition is rather weak both in the quark-model description of Ref. [28] and the hadronic or hybrid descriptions of Refs. [26,27,30,32]. However, the tensor effects arising from the pseudoscalar or vector meson exchanges become important at hadronic level, 3 while they are negligible in the QDCSM study of Ref. [28]. We have calculated the 3 S 1 c N phase shifts with the CQM just by considering the diagonal inter-action. The results are plotted by the dashed-dotted line in Fig. 2b, where the small contribution of the channel coupling can be seen, in agreement with the QDCSM results of Ref. [28]. It is worth to note that the cc conversion is less important than in the similar system in the strange sector, mainly due to the larger mass difference, namely 168 MeV as compared to 73 MeV in the strange sector. Besides, it comes reduced as compared to the strange sector due to the absence of K -meson exchanges [20], generating a smaller c N − c N transition potential. The small contribution of the channel coupling obtained by the quark-model descriptions, CQM and QDCSM, to the charmed baryon-nucleon interaction is in agreement with the observations of the HAL QCD Collaboration, leading to the conclusion that the c N tensor potential is negligibly weak [44] and that the coupling between c N and c N channels is also weak [45]. Similar conclusions were obtained in Ref. [32].
Reference [46] has extrapolated the results of the HAL QCD Collaboration to the physical pion mass using EFT. The near-identity of the lattice QCD potentials extracted for the 1 S 0 and 3 S 1 c N partial waves [44] persists in the extrapolation to the physical point. As the 3 S 1 -3 D 1 tensor coupling induced by the tensor forces is taken into account in the EFT analysis of Ref. [46], it corroborates its smallness in the spintriplet partial wave, as also derived from the CQM results of Fig. 2b. The EFT extrapolation to the physical pion mass obtains a maximum for the 1 S 0 c N phase shift of around 17-21 degrees. This result is compatible with the predictions of the CQM, as seen in Fig. 4a, where we have calculated the 1 S 0 c N phase shifts for standard quark-model values of b c ∈ [0.2, 0.8] fm. In Fig. 4b we have calculated the scattering length for the same interval of values of b c and we compare with the result of the EFT extrapolation of Ref. [46] at the physical pion mass, the orange vertical line, getting also compatible results. The CQM predicts a slightly larger attraction for the 3 S 1 c N partial wave. This result, which agrees with the conclusions of Ref. [32], is not expected to coincide with the EFT extrapolation of the HAL QCD 1 S 0 and 3 S 1 c N phase shifts, a t ∈ [−0.81, −0.98] fm, due to their identity at unphysical pion masses together with the already mentioned smallness of the tensor force in the spin-triplet partial wave.

c N interaction
In Fig. 5 we show the I = 1/2 c N phase shifts. Figure 5a presents the prediction of the CQM model for the 1 S 0 partial wave. There are no data available to compare with. The strong repulsion observed in the c N 1 S 0 (I = 1/2) interaction is a consequence of Pauli suppression effects arising in spin-isospin saturated channels [50], as discussed in Sect. 2.1. Results of other theoretical approaches for this partial wave would help to disentangle the role of the short-  Figure 5b shows the phase shifts for the 3 S 1 partial wave. The black solid line stands for the results of the CQM. The latest (2 + 1)-flavor lattice QCD simulations by the HAL QCD Collaboration [44] for a pion mass of 570 (410) MeV are shown by the blue (red) filled circles with their corresponding errors. As in the c N interaction, the tendency can be seen that the attraction becomes stronger as the pion mass decreases, the phase shifts moving towards the results of the CQM. One observes that lattice QCD simulations predict the attraction in the c N 3 S 1 (I = 1/2) channel to be stronger than in the equivalent c N channel. This conclusion also  Table 3. Note that the scattering lengths of the I = 1/2 c N system are complex because the lower c N channel is always open. Figure 6a shows the inelasticity for the c N 3 S 1 (I = 1/2) partial wave derived by the HAL QCD Collaboration [45] for a pion mass of 570 (410) MeV by blue (red) filled circles with their corresponding errors. Figure 6b shows the inelasticity obtained with the CQM for the c N 3 S 1 and 1 S 0 I = 1/2 partial waves. Although the coupling between the c N and c N channels in the 3 S 1 partial wave is small, see Fig. 2b, the inelasticity predicted by the CQM is larger than the HAL QCD simulation. . 6 a Inelasticity for the c N 1 S 0 (I = 1/2) and c N 3 S 1 (I = 1/2) partial waves predicted by the CQM as a function of the c.m. kinetic energy. a Inelasticity for the c N 3 S 1 (I = 1/2) partial wave of the HAL QCD Collaboration [45] at m π = 570 MeV (blue filled circles) and m π = 410 MeV (red filled circles) as a function of the c.m. kinetic energy. The vertical line at each point represents the statistical error of the lattice QCD simulations. The solid blue and red lines are just a guide to the eye In Fig. 7 we show the I = 3/2 c N phase shifts. The 1 S 0 c N channel presents an attraction comparable to the 1 S 0 c N system. The scattering length is still far from the standard values of the N system, in the order of −2.9 to −2.6 fm, which may allow for the existence of three-body bound states as we will discuss below. The 3 S 1 (I = 3/2) c N channel presents a strong repulsion, a consequence again of quark-Pauli effects arising in spin-isospin saturated channels. As mentioned above for the c N 1 S 0 (I = 1/2) state, it would be convenient to have results of other theoretical approaches for the phase shifts of the c N 3 S 1 (I = 3/2) partial wave, to scrutinize the short-range dynamics.
Finally, in Fig. 8 we present the CQM results for the total cross section for the c N , c N (I = 1/2), and c N (I = 3/2) scattering. One of the most interesting applications of the charmed baryon-nucleon interaction is the study of the possible existence of charmed hypernuclei. The binding energy of c hypernuclei has been analyzed in Ref. [44] using the HAL QCD c N interaction for m π = 410 MeV, where it has been noted that for nuclei with A = 12−58 the Coulomb repulsion is not much stronger than the strong binding energy, which leads to the possible existence of c hypernuclei in light or medium-heavy nuclei. On the contrary, Refs. [29,30] concluded to the existence of light c hypernuclei. Moreover, Ref. [32] concluded the existence of c hypernuclei for all nuclei studied, from 5 He to 209 Pb.
Regarding the possibility of a J = 1/2 charmed hypertriton with the HAL QCD Y c N interactions there is a delicate balance. On the one hand, it would be favored, bearing in mind the tendency that the c N attraction becomes stronger as the pion mass decreases. On the other hand, since the average c N ( N ) potential that it is relevant for the charmed hypertriton (hypertriton) is dominated by the spinsinglet channel [58], the considerably smaller 1 S 0 c N scattering length compared to the N system goes against its existence. The balance could be tilted if the spin dependence of the c N interaction induced by the short-range dynamics would slightly enhance the attraction in the spin-triplet partial wave as compared with the spin-singlet one. Then the existence of J = 3/2 c hypernuclei might be considered seriously. The isoscalar J = 3/2 state is dominated by the more attractive spin-triplet interaction [59], which together with the reduction of the kinetic energy associated with the c induced by its larger mass as compared to the , could lead to a slightly bound J = 3/2 c hypernucleus [29]. In this regard, it is important to keep in mind that the isoscalar J = 3/2 N N state is close to threshold; see Table V and Fig. 2 of Ref. [60].
The recent few-body calculation of Ref. [30], employing the strongly attractive one-boson-exchange interactions discussed above leading already to c N bound states, leads to several c N N bound states with binding energies of the order of 20 MeV. As has been discussed above, one made use of a slightly more attractive interaction for the c N spintriplet partial wave than for the spin-singlet partial wave. This generates an isoscalar J = 3/2 c N N ground state instead of J = 1/2, see Fig. 11 of Ref. [30].
The order of the isoscalar c N N J = 1/2 and J = 3/2 channels is also reversed with respect to the strange sector in the CQM model, the J = 3/2 being the most attractive one. This difference can easily be associated with the importance of the -conversion in the strange sector [59]. When the N -N potential is disconnected, the J = 3/2 channel is almost not modified, while the J = 1/2 loses great part of its attraction. Thus, the ordering between the J = 1/2 and J = 3/2 channels is reversed in such a way that the hypertriton would not be bound (see Fig. 6(a) of Ref. [60]). As we have already discussed, the cc conversion is less important than in the strange sector, giving rise to a softer c Nc N transition potential. Thus, the calculation of Ref. [29], making use of the CQM phase shifts presented in Figs. 2 and 5, i.e., without two-body bound states, obtained an isoscalar J = 3/2 charmed hypernucleus with a binding energy of 0.27 MeV. After correcting exactly by the Coulomb potential the final binding energy obtained was 0.14 MeV. Different from the hadron level calculation of Ref. [30], in the CQM model the J = 1/2 c N N is unbound. Let us finally note that the hard-core radius of the c N interaction, relevant for the study of charmed hypernuclei [20], in the CQM is fixed by the short-range dynamics [51].
There are not few-body calculations with the QDCSM Y c N interactions of Ref. [28]. However, a simple reasoning hints towards the possible existence of a J = 1/2 charmed hypertriton in this model. As one can see in Fig. 2a, the c N 1 S 0 phase shifts predicted by the QDCSM are similar to those of the c N 3 S 1 partial wave obtained with the CQM, see Fig. 2b. As the channel coupling is negligible in both cases, with the QDCSM one would obtain a scattering length for the c N 1 S 0 state of about −2.31 fm, see Table 1. This scattering length is within the order of that of the 1 S 0 N system, between −2.9 and −2.6 fm, which is a key ingredient for the existence of the hypertriton. The possible existence of a J = 1/2 charmed hypertriton in the QDCSM would be reinforced by the reduced kinetic energy contribution of the c baryon. It might be at an disadvantage by the lack of the N − N coupling that, as seen in Fig. 6(a) of Ref. [60], is of basic importance to get the hypertriton in quark-model based descriptions.
Reference [32] has also studied the possible existence of bound states of the c in different nuclei. One makes use of the c self-energy as an effective c -nucleus mean-field potential in a Schrödinger equation to get the bound state energies. c hypernuclei from 5 c He to 209 c Pb are studied. Even the less attractive model for the Y c N interaction of those discussed in Sect. 3.1, where the couplings of the σ meson with the charmed baryons are reduced 20% as compared to the original YN potential Ã of the Jülich group, is able to bind the c in all the nuclei considered. This is in contrast with the HAL QCD Collaboration results [44], which suggest that only light-or medium-mass c nuclei could really exist. The conclusions of this work come to reinforce the results obtained with the CQM in Ref. [29]. On the one hand they arrive at the same conclusion as regards the negligible contribution of the c N − c N coupling, and on the other hand they support the possible existence of light charmed hypernuclei.

Outlook
We have performed a comparative study of the charmed baryon-nucleon interaction based on different theoretical approaches. For this purpose, we make use of: (i) a constituent quark model tuned in the light-flavor baryon-baryon interaction and the hadron spectra, (ii) hadronic descriptions based on one-boson-exchange potentials, (iii) a quark delocalization color screening model, (iv) (2 + 1)-flavor lattice QCD results of the HAL QCD Collaboration at unphysical pion masses and their effective field theory extrapolation to the physical pion mass. There is a general qualitative agreement among the different available approaches to the charmed baryon-nucleon interaction. Quark-model based results point to soft interactions without two-body bound states. They also support a negligible channel coupling, due either to tensor forces or transitions between different physical channels, c N − c N . The short-range dynamics of the CQM model, fixing the hard-core radius of the S wave interactions, generates a slightly larger repulsion in the 1 S 0 than in the 3 S 1 c N partial wave. A similar asymmetry between the attraction in the two S waves of the c N interaction also appears in hadronic approaches.
Pauli suppression effects generate a major impact in the c N charmed baryon-nucleon interaction, resulting in a strong repulsion in the 1 S 0 (I = 1/2) and 3 S 1 (I = 3/2) partial waves. A comparative detailed study of Pauli suppressed partial waves, as the 1 S 0 (I = 1/2) and 3 S 1 (I = 3/2) c N channels, would help to disentangle the short-range dynamics of two-baryon systems containing heavy flavors. Quarkmodel approaches predict a small contribution of the channel coupling to the charmed baryon-nucleon interaction, concluding that the c N tensor potential is negligibly weak and that the coupling between c N and c N channels is also weak.
In the light of the results for the Y c N interactions, the possible existence of charmed hypernuclei has been discussed. The order of the isoscalar J = 1/2 and J = 3/2 channels is reversed in the charm with respect to the strange sector. While the existence of an isoscalar J = 1/2 c N N charmed hypernucleus is not likely, that of an isoscalar J = 3/2 state seems more feasible. In any case, the possible existence of c hypernuclei in light or medium-heavy nuclei is a firm prediction of quark-model and hadronic approaches to the Y c N interaction.
The understanding of the baryon-baryon interaction in the heavy flavor sector is a key ingredient in our quest to describing the properties of hadronic matter. The study of unknown two-baryon systems could benefit from wellconstrained models based as much as possible on symmetry principles and analogies with other similar processes. Subsequently, lattice QCD simulations could incorporate firmly established predictions to validate our understanding of lowenergy Quantum Chromodynamics in the multiquark sector.