The existence of neutron stars with $2M_\odot$ requires the strong stiffness of the equation of state (EoS) of neutron-star matter. We introduce a multi-pomeron exchange potential (MPP) working universally among 3- and 4-baryons to stiffen the EoS. Its strength is restricted by analyzing the nucleus-nucleus scattering with the G-matrix folding model. The EoSs are derived using the Brueckner-Hartree-Fock (BHF) and the cluster variational method (CVM) with the nuclear interactions ESC and AV18. The mass-radius relations are derived by solving the Tolmann-Oppenheimer-Volkoff (TOV) equation, where the maximum masses over $2M_\odot$ are obtained on the basis of the terrestrial data. Neutron-star radii $R$ at a typical mass $1.5M_\odot$ are predicted to be $12.3\!\sim\!13.0$ km. The uncertainty of calculated radii is mainly from the ratio of 3- and 4-pomeron coupling constants, which cannot be fixed by any terrestrial experiment. Though values of $R(1.5M_\odot)$ are not influenced by hyperon-mixing effects, finely-observed values for them indicate degrees of EoS softening by hyperon mixing in the region of $M\!\sim\!2M_\odot$. If $R(1.5M_\odot)$ is less than about 12.4 km, the softening of EoS by hyperon mixing has to be weak. Useful information can be expected by the space mission NICER offering precise measurements for neutron-star radii within $\pm 5\%$.

Neutron stars with large masses $\sim 2M_{\odot}$ require the hard stiffness of equation of state (EoS) of neutron-star matter. On the other hand, hyperon mixing brings about remarkable softening of EoS. In order to solve this problem, a multi-pomeron exchange potential (MPP) is introduced as a model for the universal many-body repulsion in baryonic systems on the basis of the Extended Soft Core (ESC) baryon-baryon interaction. The strength of MPP is determined by analyzing the nucleus-nucleus scattering with the G-matrix folding model. The interactions in $\Lambda\!N$, $\Sigma\!N$ and $\Xi\!N$ channels are shown to be consistent with experimental indications. The EoS in neutron-star matter with hyperon mixing is obtained from ESC in addition of MPP, and mass-radius relations of neutron stars are derived. The maximum mass is shown to reach $2M_{\odot}$ even in the case of including hyperon mixing on the basis of model-parameters determined by terrestrial experiments.

A multi-pomeron exchange potential (MPP) is proposed as a model for the universal many-body repulsion in baryonic systems on the basis of the Extended Soft Core (ESC) bryon-baryon interaction. The strength of MPP is determined by analyzing the nucleus-nucleus scattering with the G-matrix folding model. The interaction in $\Lambda N$ channels is shown to reproduce well the experimental $\Lambda$ binding energies. The equation of state (EoS) in neutron matter with hyperon mixing is obtained including the MPP contribution, and mass-radius relations of neutron stars are derived. It is shown that the maximum mass can be larger than the observed one $2M_{\odot}$ even in the case of including hyperon mixing on the basis of model-parameters determined by terrestrial experiments.

A multi-pomeron exchange potential (MPP) is proposed as a model for the three-body repulsion indicated in neutron-star matter, which works universally among three- and four-baryons. Its strength is determined by analyzing the nucleus-nucleus scattering with the G-matrix folding model. The EoS in neutron matter is obtained including the MPP contribution. The neutron-star mass is calculated by solving the TOV equation. The maximum mass is obtained to be larger than the observed one $1.97 M_{solar}$ on the basis of the experimental data.