The Skyrme model is a low energy effective field theory of strong interactions where nuclei and baryons appear as collective excitations of pionic degrees of freedom. In the last years, there has been a revival of Skyrme's ideas and new related models, some of them with BPS bounds (topological lower energy bounds), have been proposed. It is the aim of this paper to review how they can be applied to the study of neutron stars allowing for a description by means of topological solitons. We will focus on different aspects as the equation of state or the mass-radius relation, where we find that high maximal masses are supported.

Using a solitonic model of nuclear matter, the BPS Skyrme model, we compare neutron stars obtained in the full field theory, where gravitational back reaction is completely taken into account, with calculations in a mean-field approximation using the Tolman-Oppenheimer-Volkoff approach. In the latter case, a mean-field-theory equation of state is derived from the original BPS field theory. We show that in the full field theory, where the energy density is non-constant even at equilibrium, there is no universal and coordinate independent equation of state of nuclear matter, in contrast to the mean-field approximation. We also study how neutron star properties are modified by going beyond mean field theory, and find that the differences between mean field theory and exact results can be considerable. Further, we compare both exact and mean-field results with some theoretical and phenomenological constraints on neutron star properties, demonstrating thus the relevance of our model even in its most simple version.

The BPS Skyrme model has been demonstrated already to provide a physically intriguing and quantitatively reliable description of nuclear matter. Indeed, the model has both the symmetries and the energy-momentum tensor of a perfect fluid, and thus represents a field theoretic realization of the "liquid droplet" model of nuclear matter. In addition, the classical soliton solutions together with some obvious corrections (spin-isospin quantization, Coulomb energy, proton-neutron mass difference) provide an accurate modeling of nuclear binding energies for heavier nuclei. These results lead to the rather natural proposal to try to describe also neutron stars by the BPS Skyrme model coupled to gravity. We find that the resulting self-gravitating BPS Skyrmions provide excellent results as well as some new perspectives for the description of bulk properties of neutron stars when the parameter values of the model are extracted from nuclear physics. Specifically, the maximum possible mass of a neutron star before black-hole formation sets in is a few solar masses, the precise value depending on the precise values of the model parameters, and the resulting neutron star radius is of the order of 10 km.