The possible presence of amorphous and heterogeneous phases in the inner crust of a neutron star is expected to reduce the electrical conductivity of the crust, with potentially important consequences on the magneto-thermal evolution of the star. In cooling simulations, the disorder is quantified by an impurity parameter which is often taken as a free parameter. We aim to give a quantitative prediction of the impurity parameter as a function of the density in the crust,performing microscopic calculations including up-to-date microphysics of the crust. A multi-component approach is developed at finite temperature using a compressible liquid drop description of the ions with an improved energy functional based on recent microscopic nuclear models and optimized on extended Thomas-Fermi calculations. Thermodynamic consistency is ensured by adding a rearrangement term and deviations from the linear mixing rule are included in the liquid phase. The impurity parameter is consistently calculated at the crystallization temperature as determined in the one-component plasma approximation for the different functionals. Our calculations show that at the crystallization temperature the composition of the inner crust is dominated by nuclei with charge number around $Z \approx 40$, while the range of the $Z$ distribution varies from about 20 near the neutron drip to about 40 closer to the crust-core transition. This reflects on the behavior of the impurity parameter that monotonically increases with density up to around 40 in the deeper regions of the inner crust. Our study shows that the contribution of impurities is non-negligible, thus potentially having an impact on the transport properties in the neutron-star crust. The obtained values of the impurity parameter represent a lower limit; larger values are expected in the presence of non-spherical geometries and/or fast cooling dynamics.

Context. In the cooling process of a non-accreting neutron star, the composition and properties of the crust are thought to be fixed at the finite temperature where nuclear reactions fall out of equilibrium. A lower estimation for this temperature is given by the crystallization temperature, which can be as high as $\approx 7\times 10^9$ K in the inner crust, potentially leading to sizeable differences with respect to the simplifying cold-catalyzed matter hypothesis. Aims. We extend the recent work by Fantina et al. (2019) on the outer crust, to the study of the crystallization of the inner crust and the associated composition in the one-component plasma approximation. Methods. The finite temperature variational equations for non-uniform matter in both the liquid and the solid phases are solved using a compressible liquid-drop approach with parameters optimized on four different microscopic models which cover the present uncertainties in nuclear modeling. Results. We consider separately the effect of the different nuclear ingredients with their associated uncertainties, namely the nuclear equation of state, the surface properties in the presence of a uniform gas of dripped neutrons, and the proton shell effects arising from the ion single-particle structure. Our results suggest that the highest source of model dependence comes from the smooth part of the nuclear functional. Conclusions. We show that shell effects play an important role at the lowest densities close to the outer crust, but the most important physical ingredient to be settled for a quantitative prediction of the inner crust properties is the surface tension at extreme isospin values.

The neutron star crustal EoS and transition point properties are computed within a unified meta-modeling approach. A Bayesian approach is employed including two types of filters: bulk nuclear properties are controlled from low density effective field theory (EFT) predictions as well as the present knowledge from nuclear experiments, while the surface energy is adjusted on experimental nuclear masses. Considering these constraints, a quantitative prediction of crustal properties can be reached with controlled confidence intervals and increased precision with respect to previous calculations: {$\approx 11\%$} dispersion on the crustal width and {$\approx 27\%$} dispersion on the fractional moment of inertia. The crust moment of inertia is also evaluated as a function of the neutron star mass, and predictions for mass and radii are given for different pulsars. The possible crustal origin of Vela pulsar glitches is discussed within the present estimations of crustal entrainment, disfavoring a large entrainment phenomenon if the Vela mass is above $1.4M_\odot$. Further refinement of the present predictions requires a better estimation of the high order isovector empirical parameters, e.g. $K_{sym}$ and $Q_{sym}$, and a better control of the surface properties of extremely neutron rich nuclei.

The influence of the uncertainties of the equation of state empirical parameters on the neutron stars crust-core phase transition is explored within a meta-modeling approach, in which the energy per particle is expanded as a Taylor series in density and asymmetry around the saturation point. The phase transition point is estimated from the intersection of the spinodal instability region for dynamical fluctuations with the chemical equilibrium curve. Special attention is paid to the inclusion of high-order parameters of the Taylor series and their influence on the transition point. An uncorrelated prior distribution is considered for the empirical parameters, with bulk properties constrained through effective field theory predictions, while the surface parameters are controlled from a fit of nuclear masses using the extended Thomas Fermi approximation. The results show that the isovector compressibility $K_{sym}$ and skewness $Q_{sym}$ have the most significant correlations with the transition point, along with the previously observed influence of the $L_{sym}$ parameter. The estimated density and pressure of the crust-core transition are $n_t = (0.071 \pm 0.011) fm^{-3}$ and $P_t = (0.294 \pm 0.102) MeV fm^{-3}$.