We employ modern two- and three-nucleon interactions derived in chiral perturbation theory (ChPT) at next to-next-to leading order (N2LO), to calculate the energy per particle of symmetric nuclear matter and pure neutron matter in the framework of the microscopic Brueckner-Hartree-Fock approach. In particular, we present results concerning two optimized versions at N2LO of chiral potentials ( $ \mathrm{N2LO}_{opt}$ ), fitted to properties of light nuclei. We also employ the recently developed $\mathrm{N2LO}_{sat}$ interaction which has been calculated at the same order of ChPT but fitted in a very different way compared with the $\mathrm{N2LO}_{opt}$ interactions. We find that using these potentials, in general, it is not possible to reproduce all the saturation properties of nuclear matter. In particular the behaviour of the symmetry energy ( $E_{sym}$ predicted by the interactions considered is quite soft. This is shown comparing our results with the empirical constraints on $E_{sym}$ obtained from the data analysis of the excitation energies of isobaric analog states in nuclei and from experimental data of the neutron skin thickness of heavy nuclei. We finally confront our results with similar calculations performed by other research groups using nuclear chiral interactions at various order of ChPT and employing different many-body methods.
Read full abstract