From Finite Nuclei to Neutron Stars: The Essential Role of High-Order Density Dependence in Effective Forces

A unified description of finite nuclei and equation of state of neutron stars presents both a major challenge and also opportunities for understanding nuclear interactions. Inspired by the Lee–Huang–Yang formula of hard-sphere gases, we develop effective nuclear interactions with an additional high-order density dependent term. While the original Skyrme force SLy4 is widely used in studies of neutron stars, there are not satisfactory global descriptions of finite nuclei. The refitted SLy4' force can improve descriptions of finite nuclei but slightly reduces the radius of neutron star of 1.4M_\odot with M_\odot being the solar mass. We find that the extended SLy4 force with a higher-order density dependence can properly describe properties of both finite nuclei and GW170817 binary neutron stars, including the mass-radius relation and the tidal deformability. This demonstrates the essential role of high-order density dependence at ultrahigh densities. Our work provides a unified and predictive model for neutron stars, as well as new insights for the future development of effective interactions.