# Nucleon-nucleon interaction in a chiral SU(3) quark model revisited

###### Abstract

A dynamical investigation of the nucleon-nucleon () interaction by using the resonating group method (RGM) in a chiral SU(3) quark model has been revisited. The considered quark-quark interaction includes, besides the one-gluon exchange (OGE) and the phenomenological confinement potential, the nonet scalar and pseudo-scalar meson exchanges derived from the spontaneous SU(3) chiral symmetry breaking. The physical consistency requirement that the wave functions of single baryons satisfy the minimums of the Hamiltonian has been strictly imposed in determination of the model parameters. The calculated masses of the octet and decuplet baryon ground states, the binding energy of the deuteron, and the scattering phase shifts up to a total angular momentum are in satisfactory agreement with the experiments.

###### pacs:

13.75.Cs, 12.39.Jh, 14.20.-c, 12.39.Fe## I Introduction

Hadrons are composed of quarks and gluons. It is thus an exciting challenge to understand the phenomena of hadron physics directly from these fundamental degrees of freedom. Despite the progress made in understanding the consequence of quantum chromodynamics (QCD), the theory of strong interactions, the complexity of this theory in its non-perturbative region forces us to employ QCD inspired models in study of the hadron structures and the hadron-hadron interactions. Among these models, the constituent quark model has shown to be quite successful in describing the single baryon properties and the nucleon-nucleon () and hyperon-baryon () interactions Oka:1980 (); Harvey:1981 (); Faessler:1982 (); Shimizu:1984 (); Straub:1988 (); Baruer:1990 (); Valcarce:1994 (); Zhang:1994 (); Zhang:1997 (); Entem:2000 (); Fujiwara:2001 (); Mota:2002 (); Dai:2003 (); Ping:2009 (). In Refs. Huang:2004 (); Huang:2004-2 (); Huang:2005 (); Huang:2007 (); Huang:2008 (), progress has also been made in understanding the kaon-nucleon () and anti-kaon-nucleon () interactions in a chiral constituent quark model. On hadron level, the interaction has been well described in effective field theory Epelbaum:2009 (); Entem:2015 ().

In the constituent quark model study of interaction, the one-gluon exchange (OGE) is found to be one of the most important sources of the short-range repulsion Oka:1980 (); Harvey:1981 (); Faessler:1982 (); Shimizu:1984 (); Straub:1988 (); Baruer:1990 (); Valcarce:1994 (); Zhang:1994 (); Zhang:1997 (); Entem:2000 (); Fujiwara:2001 (); Mota:2002 (). Therefore, to get a proper understanding of the short-range interaction mechanism on a quark level, one needs a credibly determination of the coupling strengths of OGE. In earlier quark model investigations, the OGE coupling constants are usually determined by the mass differences of and , where the masses of single baryons are calculated as the averaged values of the Hamiltonian with the spacial wave functions of constituent quarks described by Gaussian wave functions. The assumption behind this strategy of the OGE couplings determination is that the harmonic oscillator size parameters in the Gaussian wave functions are the same for all single baryons. Under such an assumption, the matrix elements of the kinematic energy and the confinement potential for will be the same as those for . A similar situation applies to and . Consequently, the mass differences of and only come from the one-boson exchanges (OBEs) generated by the quark and chiral-field coupling and the OGE. One can then fix the OGE couplings by the mass differences of and with the parameters in OBEs being predetermined.

The problem for choosing the same size parameter in Gaussian wave functions for all single baryons is that the masses of baryon ground states are not guaranteed to be the minimums of the Hamiltonian, contradictory to the variational principle. In other words, the wave functions chosen for single baryons are not consistent with the model Hamiltonian. Physically, it is hard to understand why different baryons, e.g. and , or and , have exactly the same sizes although their Hamiltonians are different due to their different quantum numbers. In earlier quark model calculations, usually the nucleon is set to be the minimum of the Hamiltonian by a particular choice of the values for other model parameters, e.g. the parameters in the confinement potential. One then needs to be very careful when extends the model from the study of interaction to other baryon-baryon () systems. There may be cases that one needs to introduce additional channels to lower the energy of the considered system. These channels might not be physical ones, but are partially needed to change the internal wave functions of the single baryons. Due cautions should be taken in explaining the configuration structure of any bound states obtained in such cases. In Refs. Ohta:1982 (); Liu:1982 (), Ohta et al. and Liu found that a stability condition of the nucleon should be satisfied to make a meaningful discussion of the interaction. It is natural to expect that such an observation also holds for other systems.

In this work, we re-investigate the interaction in the framework of resonating group method (RGM) within a chiral SU(3) quark model. The quark-quark interaction includes the OGE, the phenomenological confinement potential, and the nonet scalar and pseudo-scalar meson exchanges generated from the spontaneous SU(3) chiral symmetry breaking. The major difference between the present work with earlier quark model calculations is that the quark-quark interaction employed in the present work is constrained to describe the energies of single baryons, the binding energy of deuteron, and the scattering phase shifts in a rather consistent manner without introducing any additional parameters. Specifically, in the present work, the harmonic oscillator size parameters for constituent quarks are not treated as predetermined parameters and taken to be the same for all single baryons. Instead, they are determined by variational method in calculation of the energies of single baryon ground states, which ensures that all single baryons are minimums of the Hamiltonian. The model parameters in Hamiltonian are then adjusted to simultaneously match the calculated energies of single baryons, the binding energy of the deuteron, and the predicated scattering phase shifts with their experimental values. We mention that this is the first chiral quark model investigation where the octet and decuplet baryon ground states and the interactions are handled in a consistent manner. It is expected that the results will be more reliable when one extends the present model from the study of interaction to other systems, especially to the search for dibaryons, which is expected to be done in our next step work.

The paper is organized as follows. In the next section, we review the main aspects of the chiral SU(3) quark model, the description of the baryon ground states, the formulation of RGM, and the determination of the model parameters. The results for masses of octet and decuplet baryon ground states, the binding energy of deuteron, and the scattering phase shifts are shown and discussed in Sec. III. Finally, the summary and conclusions are drawn in Sec. IV.

## Ii Formulation

### ii.1 The chiral SU(3) quark model

The idea and details of the chiral SU(3) quark model can be found in Refs. Zhang:1997 (); Huang:2004 (); Huang:2004-2 (). Here we just present the main features of this model.

The quark and chiral field interaction Lagrangian in the flavor SU(3) case can be obtained by a linear generalization of the SU(2) linear -model, which gives

(1) |

Here is the quark field, and are nonet pseudo-scalar and scalar fields, is the Gell-Mann matrix of the flavor SU(3) group, and the quark and chiral-field coupling constant. Clearly, this Lagrangian is invariant under the infinitesimal chiral SU(3)SU(3) transformation. In practice, a form factor will be inserted into the vertices of quark and chiral field coupling to describe the chiral-field structure,

(2) |

with the cutoff mass indicating the chiral symmetry breaking scale amk91 (); abu91 (); emh91 (). Equations (1) and (2) result in the following chiral field induced potentials between the -th and -th quarks:

(3) |

with

(4) | ||||

(5) |

and

(6) | ||||

(7) |

where

(8) | ||||

(9) | ||||

(10) | ||||

(11) | ||||

(12) | ||||

(13) | ||||

(14) | ||||

(15) | ||||

(16) | ||||

(17) |

and and are masses of the -th and -th constituent quarks. and are related to the mass of scalar meson, , and the mass of pseudo-scalar meson, , by

(18) | ||||

(19) |

The relations for in Eqs. (18) and (19) come from an explicit treatment of the mass difference of and quarks in the derivation of potentials of Eqs. (4) and (5) for - and -exchange from the quark and chiral field interaction Lagrangian of Eq. (1).

For pseudo-scalar meson exchanges, the fact that the physical and ’ are mixing states of and is considered in the calculation:

(20) |

with the mixing angle taken to be the empirical value .

Apart from the potentials from OBEs given in Eq. (3), to study the hadron structure and hadron-hadron dynamics, one still needs the potential from OGE,

(21) |

with

(22) | ||||

(23) |

and a phenomenological confinement potential of which the frequently used linear and quadratic types are considered in the present work,

(24) |

The total Hamiltonian for a multiquark system can then be written as

(25) |

with being the three-momentum of the -th quark, and being the number of quarks for the system considered.

### ii.2 Baryon ground states

The anti-symmetrized wave function for decuplet or octet ground state baryon in spin-flavor-color-orbit space can be written as

(26) |

where , and are wave functions of baryon in orbit-flavor, spin and color spaces, respectively, with the superscripts S, MX denoting symmetric and mixed-symmetric under interchange of any pair of quarks in the corresponding space. Note that the baryon wave functions in orbit and flavor spaces are coupled together, since the mass of strange quark is different from that of up (down) quark (),

(27) | ||||

(28) |

where is the flavor of the -th quark, the SU(3) Clebsch-Gordon coefficient in flavor space, and the orbit wave function for three quarks with flavor contents . The orbit wave function is determined by the variational method with the trial wave function taken to be a product of Gaussian functions,

(29) |

With the assumption that the harmonic oscillator frequency is the same for , and quarks which ensures that the three-quark center-of-mass motion is irrelevant to the calculation, the size parameter for strange quark, , is related to that for up quark, , by

(30) |

In earlier quark model calculations, is treated as a pre-determined parameter and taken as the same for all single baryons. In the present work, for each baryon will be obtained by minimize the corresponding baryon mass from the model calculation, i.e.

(31) |

which ensures that each single baryon is the solution of the Hamiltonian of Eq. (25).

### ii.3 RGM for system

The following Jacobi coordinates are defined to construct the total wave function of the system:

(32) | ||||

(33) | ||||

(34) | ||||

(35) | ||||

(36) |

with () being the internal coordinates for two nucleon clusters, and bing the relative coordinate of two clusters.

Following the cluster model calculations mka77 (); mok81 (), the RGM wave function is written as

(37) |

with being the trial wave function of relative motion of two nucleon clusters, and being the total spin and isospin of the system. is the internal wave function of each single cluster, taken to be the same as that from Eq. (26). The wave function of total center-of-mass motion is stripped off as it is irrelevant to the calculation. The symbol is the anti-symmetrizing operator for interchange any pair of constituent quarks between two clusters,

(38) |

where is the permutation operator of the -th and -th quarks that are from different clusters. Substituting into the following projection equation

(39) |

with and being the inner energy of nucleon and the relative energy between two nucleon clusters, respectively, one gets the RGM equation for the unknown relative motion wave function

(40) |

with

(41) |

where the Hamiltonian kernel and normalization kernel can, respectively, be calculated by

(42) |

In practical calculation, the unknown wave function for two-cluster relative motion is projected into partial waves,

(43) |

For a bound state problem, the -wave relative wave function is expanded as

(44) |

with

(45) |

where is called the generate coordinate, and the -th modified spherical Bessel function. The information about the unknown wave function of two-cluster relative motion is now exhibited by the coefficients ’s with properly chosen values of the generate coordinates. Performing a variational procedure, one deduces the -th partial-wave equation for a bound-state problem,

(46) |

with

(47) | |||

(48) |

Solving Eq. , one gets the binding energy and the corresponding wave function of the two-cluster system.

For a scattering problem, the radial part of the relative motion wave function is expanded as

(49) |

with

(50) |

where is the -th spherical Hankel functions, is the momentum of the two-cluster relative motion, and is a cutoff radius beyond which all the strong interactions can be disregarded. The complex parameters and are determined by the smoothness condition at and ’s satisfy . Performing a variational procedure, the -th partial-wave equation for the scattering problem can be deduced as

(51) |

with

(52) | ||||

(53) |

and

(54) |

where the RGM kernel is defined in Eq. . Solving Eq. , the -matrix element and the phase shifts can be obtained by

(55) |

### ii.4 Model parameters

Model I | Model II | |
---|---|---|

( conf.) | ( conf.) | |

The predetermined model parameters are: (1) the masses of mesons, which are taken to be their experimental values except for meson, whose mass is treated as a parameter to be fitted by the binding energy of deuteron and the scattering phase shifts; (2) the cutoff mass in the form factor of Eq. (2), which is taken to be 1100 MeV as usual Dai:2003 (); Huang:2004 (); Huang:2004-2 (); Huang:2005 (); Huang:2007 (); Huang:2008 (); (3) the constituent , , quark masses, which are chosen to be MeV, MeV Zhang:1997 (); Dai:2003 (); Huang:2004 (); Huang:2004-2 (); Huang:2005 (); Huang:2007 (); Huang:2008 (); (4) the quark and chiral field coupling constant, which is fixed by the relation

(56) |

with taken to be the empirical value .

The left parameters are the coupling constants of OGE, the parameters in confinement potential, and the mass of meson. They are adjusted to match the baryon masses calculated at the minimums of the Hamiltonian, the binding energy of deuteron, and the scattering phase shifts with their experimental values. The fitted values of those parameters are listed in Table 1, where model I and model II refer to models with linear and quadratic confinement as presented in Eq. (24), respectively.

One sees from Table 1 that the OGE coupling constants fixed in the present work are a little bit larger than those in earlier quark model calculations, but they are still qualitatively consistent with those usually expected in QCD at the squared momentum transfer GeV, namely .

Note that in earlier quark model calculations, the harmonic oscillator size parameter in the Gaussian wave functions is also treated as a predetermined parameter, and the masses of all the octet and decuplet baryons are then calculated by using the same as that for nucleon. As mentioned in Sec. I, the problem in such calculations is that the calculated energies of singe baryon ground states other than nucleon are not minimums of the Hamiltonian, contradictory to the variational principle. Therefore, cautions should be taken when one uses the same set of parameters to study the systems other than . This issue will be further discussed in the next section.

## Iii Results and discussions

As shown in Table 1, we have nine adjustable model parameters, namely the sigma meson mass , the OGE coupling constants and , the confinement strengths , , and zero point energies , , . With the parameter values listed in Table 1, we get the masses of baryon ground states by variational method of Eq. (31), the deuteron binding energy by solving Eq. (46), the scattering phase shifts by solving Eq. (51) and calculating Eq. (55).

Expt. [MeV] | ||||||||
---|---|---|---|---|---|---|---|---|

Theo. [MeV] | ||||||||

[fm] ( conf.) | ||||||||

( conf.) |

Table 2 shows our results for energies and corresponding size parameters of octet and decuplet baryon ground states. One sees that our theoretical masses calculated for all baryon ground states are consistent with their experimental values in both linear confinement and quadratic confinement models. One also observes that the size parameters in the Gaussian wave functions of Eq. (29) are quite different for various baryons, and they are also different in models with linear or quadratic confinement.

Model I | Model II | Expt. |

( conf.) | ( conf.) | |

In earlier quark model calculations, the size parameter is treated as a pre-determined parameter and it is the same for all , and other single baryons. In our opinion, once the Hamiltonian is given, the baryon wave functions should be obtained in principle by solving a Schrödinger equation. Of course this is impractical due to the complexity of the quark-quark interacting potential [cf. Eq. (25)]. An approximated solution could then be obtained by using Gaussians as trial wave functions and determining the baryon masses and sizes by variational method. This has, for the first time, been strictly imposed in the present work without introducing any additional parameters compared with earlier quark model calculations Zhang:1997 (); Dai:2003 ().

In Fig. 1 we show our calculated baryon mass as a function of the variational size parameter for each baryon ground state. One sees that the minimums of the masses of various baryons are located at different values of the size parameter . For octet baryons, the size parameters resulted from the linear confinement model are close to those from the quadratic confinement model. While for decuplet baryons, the size parameters resulted from the linear confinement model are much bigger than those from the quadratic confinement model. This indicates that the decuplet baryons are much more sensitive than the octet baryons to the choice of the type of the phenomenological confinement potential. One also observes that in both models, the size parameters for decuplet baryons are always much bigger than those for octet baryons, which is distinct from the earlier quark model calculations in literature, where the size parameters for all octet and decuplet baryons are taken to be the same. Under such an assumption, the OGE coupling constants are claimed to be determined by the mass differences of - and -, with the masses of all , , and baryons calculated at the same size parameter . In these calculations, one may adjust the parameters in the confinement potential to make the masses of and stable against the parameter , but then the , and other baryons are obviously not the minimums of the Hamiltonian, as indicated by Table 2 and Fig. 1. When one extends the model from the study of interaction to other systems especially when the decuplet baryons are involved, there might be cases that one needs to introduce additional channels to lower the energy of the considered system by changing the internal wave functions of the single baryons. Due caution should be exercised in explaining the configuration structure of any bound states obtained in such kind of calculations.

Table 3 shows the binding energy of deuteron calculated in the present work with a comparison with the experimental value. One sees that our calculated values in both linear confinement model and quadratic confinement model are in good agreement with the data. This can be easily understood, since in the parameters fitting procedure, we found that this quantity is rather sensitive to the mass of meson, and one thus can fine tune to make the theoretical binding energy of deuteron close to the experimental data.

Figures 2-8 show the scattering phase shifts and Fig. 9 shows the mixing parameters for coupled partial waves up to a total angular momentum calculated in the present work, where the red solid and blue dashed lines represent the results calculated in models with linear confinement and quadratic confinement, respectively. The scattered symbols are results form SAID’s partial wave analysis Workman:2016 (), with the full circles representing their single energy analysis and the empty diamonds representing their energy dependent solution (SM16). One sees that the overall agreement of our results with SAID’s partial wave analysis is satisfactory. At least they are not worse than the earlier quark model calculations, while in the present work we have put much more strict constraints on the model parameters.

Figure 2 shows the -wave phase shifts. For partial wave, the theoretical phase shifts are in perfect agreement with the data. This channel corresponds to the deuteron quantum numbers, which means that the coupling to the -partial wave caused by the tensor force resulted from OGE and one-pion exchange is important. Given the size parameter for nucleon determined by variational method of Eq. (31), and the mass of sigma meson fine tuned to get the deuteron binding energy, the phase shifts for this partial wave are found to be automatically in agreement with the data in the present work. The differences of the phase shifts from the linear confinement model and quadratic confinement model are found to be negligible. For partial wave, the calculated phase shifts represented by the red dotted and blue dash-dotted lines of the left subfigure indicate that there is a lack of attraction in this partial wave. It is known that the required attraction in this partial wave is supplied by the coupling to the partial wave Zhang:1997 (), caused by the tensor force offered by OGE and one-pion exchange. In the present work, since and have different size parameters, it is rather complicated to perform a rigorous calculation of the coupling between and partial waves in the framework of RGM. Here for the sake of simplicity, we estimate the effect from this coupling by calculating the off-diagonal transition matrix elements at nucleon’s size parameter. The resulted phase shifts are shown as red solid line (linear confinement model) and blue dashed line (quadratic confinement model) in left subfigure of Fig. 2. It is clearly seen that the coupling to the partial wave improves the phase shifts significantly. Nevertheless, there are still noticeable discrepancies in the energy region of MeV, indicating a redundant attraction obtained from the coupling to the partial wave. A more complete analysis of this partial wave requires a strict calculation of the coupling of and partial waves in RGM with the difference of the size parameters for and been properly taken into account. Such an investigation requires a new development of the RGM and is postponed until our future work.

The -wave phase shifts are shown in Fig. 3. One sees that the phase shifts for -wave and -wave are described satisfactorily. The phase shifts are too attractive, while the phase shifts are a little bit less attractive in the energy region of MeV, indicating the need for a moderate spin-orbit force, a topic deserving special treatment outside this general work.

Figure 4 shows the -wave phase shifts. Apart from the partial wave, which seems too attractive in the energy region of MeV, the phase shifts for other partial waves are described quite well.

The phase shifts for the -wave are shown in Fig. 5. The partial wave is well described. The and partial waves are perfectly described in the MeV energy region, and they are a little bit too attractive when MeV. The partial wave is a little bit too attractive in the MeV energy region.

Figures 6-8 show the results for , , and partial waves up to a total angular momentum . The phase shifts for all the partial waves are described quite well, except that there are minor deviations in a few partial waves.

Figure 9 shows the mixing parameters for the relevant coupled partial waves. The coupling of these spin-triplet partial waves are due to the tensor forces stemming from the OGE and pseudo-scalar meson exchanges. One sees that the mixing parameters for all the coupled partial waves considered in the present work are in good agreement with the experimental data, except that for the coupling of -, which agrees with data only at very low energies. Note that the phase shifts for the and partial waves are not quite well described in the present work, as shown in Figs. 3 and 5.

In a word, the scattering phase shifts for , , , , , , and partial waves up to a total angular momentum and mixing parameters for the relevant coupled partial waves have been satisfactorily described in the present work. Note that the model has only 9 adjustable parameters as tabulated in Table 1. It is known that the confinement potential itself results in negligible contributions between two color-singlet clusters Huang:2004 (); Huang:2004-2 (); Huang:2005 (); Huang:2007 (). Therefore, among these 9 adjustable parameters, only , and are relevant to the scattering observables. Since and have already been determined in fitting the masses of octet and decuplet single baryons, apart from the fact that is fine tuned to reproduce the binding energy of the deuteron, the calculation of the scattering phase shifts for 26 partial waves and mixing parameters for the relevant coupled partial waves in the present work is almost parameter-free. With this in mind, it is fair to say that the description of the scattering data in the present work is promising.

In earlier quark model calculations, the couplings of , and the mass of sigma meson are independent of the choice of confinement type (linear or quadratic) which is assumed to be proportional to the color-color operator as shown in Eq. (24). There, the masses of all single baryons are calculated at the same size parameter . Of course by doing so the single baryons other than and are not minimums of the model Hamiltonian. In spite of that, the and get the same contributions from the kinematic energy and the confinement potential, so do and . Then the parameters and can be uniquely determined by the mass differences of and once the quark-chiral field coupling is fixed by Eq. (56) and the cutoff mass is selected to be around the chiral symmetry breaking scale. The parameters in the confinement potential are just used to make the and stable against the size parameter and to adjust the matrix elements of Hamiltonian calculated at close to the empirical masses of all single baryon ground states. Since the and are independent of the choice of confinement type, and the confinement potential itself does not considerably contribute to the interaction, the sigma meson mass determined by fitting the scattering data will then also be independent of the choice of confinement type. When one extends the model to study other systems composed of two color-singlet clusters, no sizable differences will be observed when different types of confinement potential are chosen.

The situation is quite different in the present work. One sees that the parameters , and are different in linear confinement model and quadratic confinement model. The reason for that is the following. We want to describe the baryon ground states, the binding energy of deuteron, and the scattering processes in a consistent manner. This requires the size parameter in the wave functions of each single baryon to be determined by minimization of the Hamiltonian. As can be seen in Table 2 and Fig. 1, the resulted size parameters are different for different baryons in both the linear confinement model and the quadratic confinement model. Therefor, unlike in the earlier quark model calculations, the kinematic energies and the confinement potential now also contribute to the mass differences of and due to the unequal size parameters of these baryons. The mass differences remained for contributions of OGE are thus usually different for linear confinement model and quadratic confinement model, resulting in different OGE coupling constants. The difference of the parameter is then easy to be understood as it is partially needed to compensate the difference of the OGE interaction.

All though the parameters , and are different in linear confinement model and quadratic confinement model, one sees from Figs. 2-8 that the phase shifts from these two models are very close. In Fig. 10 we show the diagonal matrix elements of OGE and the total Hamiltonian in the generator coordinate method (GCM) calculation for and partial waves, which can be qualitatively regarded as the effective OGE and Hamiltonian between two clusters. In this figure, the red solid and blue dashed lines represent the matrix elements for Hamiltonian in models with linear and quadratic confinement, respectively. The red dotted and blue dash-dotted lines correspond to the results for OGE in models with linear and quadratic confinement, respectively. The horizontal axis denotes the generator coordinate which can qualitatively describe the distance between two clusters. One sees that the interaction of OGE at short range is more repulsive in linear confinement model than that in quadratic confinement model, so is the total Hamiltonian. Nevertheless, these two models result in quite similar phase shifts, implying that the calculated phase shifts in the energy region considered are not that sensitive to the short range interaction. For deeply bound states, one may see difference of the results from these two models. The OGE and Hamiltonian matrix elements in the present work are also different from those in earlier quark model calculations. This to some extent also indicates that the scattering data alone cannot uniquely determine the parameters in a constituent quark model. To get more reliable quark-quark interactions for the extension to study other systems in a parameter-free way, one should constrain the model by fitting as many data as possible. In this regard, the present work, which describes the energies of single baryon ground states, the binding energy of the deuteron, and the scattering phase shifts up to a total angular momentum in a rather consistent manner, may be considered as a step in the right direction.

Using the present model to get predictions for other systems is beyond the scope of the present work. But based on the parameters listed in Table 1 and the baryon size parameters shown in Table 2 and Fig. 1, it is still possible to make some interesting arguments without any detailed calculations. As we all know, the sigma meson exchange contributes attraction to all baryon-baryon systems. Its mass, , determined in the present work is 608 MeV in linear confinement model and 625 MeV in quadratic confinement model, both higher than that used in Ref. Dai:2003 (), 595 MeV. For system, the attraction from the meson exchange is not reduced in the present work since nucleon has a smaller size parameter, fm, compared with fm used in Ref. Dai:2003 (). However, for decuplet baryons, the size parameters determined by the variational method in the present work are much bigger than fm. It is then expected that a much smaller attraction from the meson exchange will be obtained in the present model when study the interactions of decuplet baryons, e.g. or . Detailed investigations of these systems in the present model are in progress, and the results will be reported subsequently.

In the present work, the short-range quark-quark interaction is dominated by OGE and quark exchange effects. In literature, there are authors studied the short-range interaction as stemming from vector-meson exchanges on quark level Dai:2003 (); Shimizu:2000 (). As mentioned in Ref. Huang:2005 (), it is still a controversial and challenging problem whether OGE or vector-meson exchange is the right mechanism or both of them are important for describing the short-range quark-quark interaction. We postpone the study of this issue by incorporating the vector-meson exchanges in the quark-quark interaction Hamiltonian to our next step work.

## Iv Summary and conclusions

The aim of the present work is to perform a consistent description of the masses of single baryons and the scattering data by using the same Hamiltonian in a chiral SU(3) quark model. The considered quark-quark interactions includes, besides the OGE potential and the phenomenological confinement potential, the potentials from the nonet scalar meson exchanges and the nonet pseudo-scalar meson exchanges derived from the spontaneous SU(3) chiral symmetry breaking. The masses of single baryons are calculated by using the Gaussian trial wave functions with the size parameters determined by the variational method, which ensures that all single baryons are minimums of the model Hamiltonian. The interaction is dynamically investigated by using the resonating group method. It is found that the calculated masses of the octet and decuplet baryon ground states, the binding energy of the deuteron, and the scattering phase shifts up to a total angular momentum and mixing parameters for the relevant coupled partial waves are in satisfactory agreement with the experiments. The present model may serve as a good starting point to achieve a consistent and unified description of the single baryon properties and the baryon-baryon dynamics. Investigations of baryon spectroscopy and the interaction dynamics for systems other than in a completely parameter-free way within the present model are planned for our future work.

###### Acknowledgements.

The authors thank the discussions with Profs. Z. Y. Zhang, P. N. Shen, Y. B. Dong and J. X. Wang. This work is partially supported by the National Natural Science Foundation of China under Grants No. 11475181, No. 11635009 and No. 11475015, the Youth Innovation Promotion Association of CAS under Grant No. 2015358, and the Key Research Program of Frontier Sciences of CAS under Grant No. Y7292610K1.## References

- (1) M. Oka and K. Yazaki, Phys. Lett. B 90, 41 (1980).
- (2) M. Harvey, Nucl. Phys. A 352, 326 (1981).
- (3) A. Faessler, F. Fernández, G. Lübeck, and K. Shimizu, Phys. Lett. B 112, 201 (1982).
- (4) K. Shimizu, Phys. Lett. B 148, 418 (1984).
- (5) U. Straub, Z. Y. Zhang, K. Braüer, A. Faessler, S. B. Khadkikar, and G. Lübeck, Nucl. Phys. A 483, 686 (1988).
- (6) K. Braüer, A. Faessler, F. Fernández, and K. Shimizu, Nucl. Phys. A 507, 599 (1990).
- (7) A. Valcarce, A. Buchmann, F. Fernández, and A. Faessler, Phys. Rev. C 50, 2246 (1994).
- (8) Z. Y. Zhang, A. Faessler, U. Straub, and L.Ya. Glozman, Nucl. Phys. A 578, 573 (1994).
- (9) Z. Y. Zhang, Y. W. Yu, P. N. Shen, L. R. Dai, A. Faessler, and U. Straub, Nucl. Phys. A 625, 59 (1997).
- (10) D. R. Entem, F. Fernández, and A. Valcarce, Phys. Rev. C 62, 034002 (2000).
- (11) Y. Fujiwara, T. Fujita, M. Kohno, C. Nakamoto, and Y. Suzuki, Phys. Rev. C 65, 014002 (2001).
- (12) R. D. Mota, A. Valcarce, F. Fernández, D. R. Entem, and H. Garcilazo, Phys. Rev. C 65, 034006 (2002).
- (13) L. R. Dai, Z. Y. Zhang, Y. W. Yu, and P. Wang, Nucl. Phys. A 727, 321 (2003).
- (14) J. L. Ping, H. X. Huang, H. R. Pang, F. Wang, and C. W. Wong, Phys. Rev. C 79, 024001 (2009).
- (15) F. Huang, Z. Y. Zhang, and Y. W. Yu, Phys. Rev. C 70, 044004 (2004).
- (16) F. Huang and Z. Y. Zhang, Phys. Rev. C 70, 064004 (2004).
- (17) F. Huang and Z. Y. Zhang, Phys. Rev. C 72, 024003 (2005).
- (18) F. Huang, W. L. Wang, Z. Y. Zhang, and Y. W. Yu, Phys. Rev. C 76, 018201 (2007).
- (19) F. Huang, W. L. Wang, and Z. Y. Zhang, Int. J. Mod. Phys. A 23, 3057 (2008).
- (20) E. Epelbaum, H. W. Hammer, and Ulf-G. Meißner, Rev. Mod. Phys. 81, 1773 (2009).
- (21) D. R. Entem, N. Kaiser, R. Machleidt, and Y. Nosyk, Phys. Rev. C 91, 014002 (2015).
- (22) S. Ohta, M. Oka, A. Arima, and K. Yazaki, Phys. Lett. B 119, 35 (1982).
- (23) K. F. Liu, Phys. Lett. B 114, 22 (1982).
- (24) A. M. Kusainov, V. G. Neudatchin, and I. T. Obukhovsky, Phys. Rev. C 44, 2343 (1991).
- (25) A. Buchmann, E. Fernández, and K. Yazaki, Phys. Lett. B 269, 35 (1991).
- (26) E. M. Henley and G. A. Miller, Phys. Lett. B 251, 453 (1991).
- (27) M. Kamimura, Suppl. Prog. Theor. Phys. 62, 236 (1977).
- (28) M. Oka and K. Yazaki, Prog. Theor. Phys. 62, 556 (1981).
- (29) R. L. Workman, W. J. Briscoe, and I. I. Strakovsky, Phys. Rev. C 94, 065203 (2016).
- (30) K. Shimizu and L. Ya. Glozman, Phys. Lett. B 477, 59 (2000).