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.

We study a scalar field theory with a non-standard kinetic term minimally coupled to gravity. We establish the existence of compact boson stars, that is, static solutions with compact support of the full system with self-gravitation taken into account. Concretely, there exist two types of solutions, namely compact balls on the one hand, and compact shells on the other hand. The compact balls have a naked singularity at the center. The inner boundary of the compact shells is singular, as well, but it is, at the same time, a Killing horizon. These singular, compact shells therefore resemble black holes.