We study the localization of protons in the core of neutron stars for ten realistic nuclear models that share a common behaviour of nuclear symmetry energy which saturates and eventually decreases at high densities. This results in the low proton fraction of beta-stable neutron star matter. Protons form a small admixture in the neutron star core, which is localized at sufficiently high densities. For every model we calculate the density $n_{loc}$ above which the localization effect is present. Our results indicate that localization occurs at densities above $0.5-1.0 fm^{-3}$. The phase with localized protons occupies a spherical shell or a core region inside neutron stars which contains significant fraction of all nucleons. Proton localization is of great importance for astrophysical properties of neutron stars as it strongly affects transport coefficients of neutron star matter and can produce spontaneous magnetization in neutron stars.

Neutron stars are discussed as laboratories of physics of strong gravitational fields. The mass of a neutron star is split into matter energy and gravitational field energy contributions. The energy of the gravitational field of neutron stars is calculated with three different approaches which give the same result. It is found that up to one half of the gravitational mass of maximum mass neutron stars is comprised by the gravitational field energy. Results are shown for a number of realistic equations of state of neutron star matter.

We study properties of the proton component of neutron star matter for a number of realistic nuclear models. Protons which form a few percent admixture tend to be localized in potential wells corresponding to neutron matter inhomogeneities created by the protons in the neutron medium. We calculate the energy of the Wigner-Seitz cell enclosing a single localized proton. The neutron background is treated in the Thomas-Fermi approximation and the localized proton is described by the Gaussian wave function. The neutron density profile is obtained by solving the appropriate variational equation. This approach gives lower energies of localized protons than obtained previously with less sophisticated methods.