Low Density Neutron Star Matter with Quantum Molecular Dynamics: The Role of Vector Interactions

The effect of isospin-dependent nuclear forces on the inner crust of neutron stars is modeled within the framework of Quantum Molecular Dynamics (QMD). To successfully control the density dependence of the symmetry energy of neutron-star matter below nuclear saturation density, a coupling potential between the $\omega$ and $\rho$ meson fields is introduced. This approach is inspired by the baryon density and isospin density-dependent repulsive Skyrme force of asymmetric nuclear matter. In isospin-asymmetric nuclear matter, the system shows nucleation, as nucleons are arranged into shapes resembling nuclear pasta. The dependence of clusterization in the system on the isospin properties is also explored by calculating two-point correlation functions. We show that, as compared to previous results that did not involve the $\omega$-$\rho$ potential, the energy symmetry slope $L$ is successfully controlled by varying the $\omega$-$\rho$ coupling strength. Nevertheless, the effect of changing the slope of the nuclear symmetry energy $L$ on the crust-core transition density does not seem significant. To the knowledge of the authors, This is the first implementation of a $\omega$-$\rho$ coupling in a QMD model for the inner crust of neutron stars.