Unified equations of state for cold nonaccreting neutron stars with Brussels-Montreal functionals. V. Improved parametrization of the nucleon density distributions

We previously studied the inner crust and the pasta mantle of a neutron star within the 4th-order extended Thomas-Fermi (ETF) approach with consistent proton shell corrections added perturbatively via the Strutinsky integral (SI) theorem together with the contribution due to pairing. To speed up the computations and avoid numerical problems, we adopted parametrized nucleon density distributions. However, the errors incurred by the choice of the parametrization are expected to become more significant as the mean baryon number density is increased, especially in the pasta mantle where the differences in the energy per nucleon of the different phases are very small, typically a few keV. To improve the description of these exotic structures, we discuss the important features that a nuclear profile should fulfill and introduce two new parametrizations. Performing calculations using the BSk24 functional, we find that these parametrizations lead to lower ETF energy solutions for all pasta phases than the parametrization we adopted before and more accurately reproduce the exact equilibrium nucleon density distributions obtained from unconstrained variational calculations. Within the ETFSI method, all parametrizations predict the same composition in the region with quasi-spherical clusters. However, the two new parametrizations lead to a different mantle structure at mean baryon densities above about 0.07 fm^-3, at which point lasagna is energetically favored. Interestingly, spherical clusters reappear in the pasta region. The inverted pasta phases such as bucatini and Swiss cheese are still found in the densest region above the core in all cases.