It is shown that recently proposed RMF model with $\sigma$ and $\delta$ meson interaction agrees with the observational data and presents an interesting structure with phase transition in the outer part of neutron star core.

The stability of periodically placed slabs occurring in neutron stars (lasagna phase) is examined by exact geometrical methods for the first time. It appears that the slabs are stable against any shape perturbation modes for the whole range of volume fraction occupied by the slab. The calculations are done in the framework of the liquid drop model and obtained results are universal - they do not depend on model parameters like surface tension or charge density. The results shows that the transition to other pasta shapes requires crossing the finite energy barrier

In this paper the analysis of the result of numerical simulations of pasta phases using algebraic topology methods is presented. These considerations suggest that some phases can be further split into (sub)phases and therefore should be more refined in numerical simulations. The results presented in the paper can also be used to relate the Euler characteristic from numerical simulations to the geometry of the phases. The Betti numbers are used as they provide finer characterization of the phases. It is also shown that different boundary conditions give different outcomes.

The rigorous treatment of proton localization phenomenon in asymmetric nuclear matter is presented. The solution of proton wave function and neutron background distribution is found by the use of the extended Thomas-Fermi approach. The minimum of energy is obtained in the Wigner- Seitz approximation of spherically symmetric cell. The analysis of three different nuclear models suggests that the proton localization is likely to take place in the interior of neutron star.