No Arabic abstract
We propose a new equation of state for nuclear matter based on a generalized Skyrme model which is consistent with all current constraints on the observed properties of neutron stars. This generalized model depends only on two free parameters related to the ranges of pressure values at which different submodels are dominant, and which can be adjusted so that mass-radius and deformability constraints from astrophysical and gravitational wave measurements can be met. Our results support the Skyrme model and its generalizations as good candidates for a low energy effective field-theoretic description of nuclear matter even at extreme conditions such as those inside neutron stars.
We study the coupling of nuclear matter described by the BPS Skyrme model to generalized gravity. Concretely, we consider the Starobinsky model which provides the leading-order correction to the Einstein-Hilbert action. Static solutions describing neutron stars are found both for the full field theory and for the mean-field approximation. We always consider the full Starobinsky model in the nonperturbative approach, using appropriately generalized shooting methods for the numerical neutron star calculations. Many of our results are similar to previous investigations of neutron stars for the Starobinsky model using other models of nuclear matter, but there are some surprizing discrepancies. The Newtonian mass relevant for the surface redshift, e.g., results larger than the ADM mass in our model, in contrast to other investigations. This difference is related to the particularly high stiffness of nuclear matter described by the BPS Skyrme model and offers an interesting possibility to distinguish different models of nuclear matter within generalized gravity.
The equation of state (EOS) of dense matter is an essential ingredient for numerical simulations of core-collapse supernovae and neutron star mergers. The properties of matter near and above nuclear saturation density are uncertain, which translates into uncertainties in astrophysical simulations and their multi-messenger signatures. Therefore, a wide range of EOSs spanning the allowed range of nuclear interactions are necessary for determining the sensitivity of these astrophysical phenomena and their signatures to variations in input microphysics. We present a new set of finite temperature EOSs based on experimentally allowed Skyrme forces. We employ a liquid drop model of nuclei to capture the non-uniform phase of nuclear matter at sub-saturation density, which is blended into a nuclear statistical equilibrium EOS at lower densities. We also provide a new, open-source code for calculating EOSs for arbitrary Skyrme parametrizations. We then study the effects of different Skyrme parametrizations on thermodynamical properties of dense astrophysical matter, the neutron star mass-radius relationship, and the core collapse of 15 and 40 solar mass stars.
First proposed in 2013 by Yagi and Yunes, the quasi-universal emph{I-Love-Q relations} consist of a set of relations between the moment of inertia, the spin-induced quadrupole moment and the electric quadrupolar tidal deformability of neutron stars which are independent of the Equation of State (EoS) within an accuracy of $sim1%$. In this work, we show that these relations hold for different Skyrme-based nuclear matter EoS and also for the star-like solutions of different Einstein-BPS-Skyrme-models, some of which do not even present a barotropic equation of state. Further, other quasi-universal relations are analyzed, and together with recent GW observations, we use them to select the generalized Skyrme model that better reproduces observations. Our results reaffirm both the universality of the emph{I-Love-Q} relations and the suitability of generalized Skyrme models to describe nuclear matter inside neutron stars.
The increasing number and precision of measurements of neutron star masses, radii, and, in the near future, moments of inertia offer the possibility of precisely determining the neutron star equation of state. One way to facilitate the mapping of observables to the equation of state is through a parametrization of the latter. We present here a generic method for optimizing the parametrization of any physically allowed EoS. We use mock equations of state that incorporate physically diverse and extreme behavior to test how well our parametrization reproduces the global properties of the stars, by minimizing the errors in the observables mass, radius, and the moment of inertia. We find that using piecewise polytropes and sampling the EoS with five fiducial densities between ~1-8 times the nuclear saturation density results in optimal errors for the smallest number of parameters. Specifically, it recreates the radii of the assumed EoS to within less than 0.5 km for the extreme mock equations of state and to within less than 0.12 km for 95% of a sample of 42 proposed, physically-motivated equations of state. Such a parametrization is also able to reproduce the maximum mass to within 0.04 M_sun and the moment of inertia of a 1.338 M_sun neutron star to within less than 10% for 95% of the proposed sample of equations of state.
We present a novel method for revealing the equation of state of high-density neutron star matter through gravitational waves emitted during the postmerger phase of a binary neutron star system. The method relies on a small number of detections of the peak frequency in the postmerger phase for binaries of different (relatively low) masses, in the most likely range of expected detections. From such observations, one can construct the derivative of the peak frequency versus the binary mass, in this mass range. Through a detailed study of binary neutron star mergers for a large sample of equations of state, we show that one can extrapolate the above information to the highest possible mass (the threshold mass for black hole formation in a binary neutron star merger). In turn, this allows for an empirical determination of the maximum mass of cold, nonrotating neutron stars to within 0.1 M_sun, while the corresponding radius is determined to within a few percent. Combining this with the determination of the radius of cold, nonrotating neutron stars of 1.6 M_sun (to within a few percent, as was demonstrated in Bauswein et al., PRD, 86, 063001, 2012), allows for a clear distinction of a particular candidate equation of state among a large set of other candidates. Our method is particularly appealing because it reveals simultaneously the moderate and very high-density parts of the equation of state, enabling the distinction of mass-radius relations even if they are similar at typical neutron star masses. Furthermore, our method also allows to deduce the maximum central energy density and maximum central rest-mass density of cold, nonrotating neutron stars with an accuracy of a few per cent.