Constraints on the nuclear symmetry energy from asymmetric-matter calculations with chiral NN and 3N interactions

The nuclear symmetry energy is a key quantity in nuclear (astro)physics. It describes the isospin dependence of the nuclear equation of state, which is commonly assumed to be almost quadratic. In this work, we confront this standard quadratic expansion of the equation of state with explicit asymmetric nuclear-matter calculations based on a set of commonly used Hamiltonians including two- and three-nucleon forces derived from chiral effective-field theory. We study, in particular, the importance of nonquadratic contributions to the symmetry energy, including the nonanalytic logarithmic term introduced by Kaiser [Phys. Rev. C 91, 065201 (2015)]. Our results suggest that the nonquadratic contribution to the symmetry energy can be systematically determined from the various Hamiltonians employed, and we obtain 0.74-+0.11 0.08MeV (or -1.02-(0.08)+0.11 MeV for the potential term with the effective-mass contribution) at nuclear saturation density, while the logarithmic contribution to the symmetry energy is relatively small and model-dependent. We also employ the meta-model approach to study the impact of the higher-order contributions on the neutron-star crust-core transition density, and find a 5% correction.