No Arabic abstract
The increasing richness of data related to cold dense matter, from laboratory experiments to neutron-star observations, requires a framework for constraining the properties of such matter that makes use of all relevant information. Here, we present a rigorous but practical Bayesian approach that can include diverse evidence, such as nuclear data and the inferred masses, radii, tidal deformabilities, moments of inertia, and gravitational binding energies of neutron stars. We emphasize that the full posterior probability distributions of measurements should be used rather than, as is common, imposing a cut on the maximum mass or other quantities. Our method can be used with any parameterization of the equation of state (EOS). We use both a spectral parameterization and a piecewise polytropic parameterization with variable transition densities to illustrate the implications of current measurements and show how future measurements in many domains could improve our understanding of cold catalyzed matter. We find that different types of measurements will play distinct roles in constraining the EOS in different density ranges. For example, better symmetry energy measurements will have a major influence on our understanding of matter somewhat below nuclear saturation density but little influence above that density. In contrast, precise radius measurements or multiple tidal deformability measurements of the quality of those from GW170817 or better will improve our knowledge of the EOS over a broader density range.
We show how observations of gravitational waves from binary neutron star (BNS) mergers over the next few years can be combined with insights from nuclear physics to obtain useful constraints on the equation of state (EoS) of dense matter, in particular, constraining the neutron-matter EoS to within 20% between one and two times the nuclear saturation density $n_0approx 0.16 {text{fm}^{-3}}$. Using Fisher information methods, we combine observational constraints from simulated BNS merger events drawn from various population models with independent measurements of the neutron star radii expected from x-ray astronomy (the Neutron Star Interior Composition Explorer (NICER) observations in particular) to directly constrain nuclear physics parameters. To parameterize the nuclear EoS, we use a different approach, expanding from pure nuclear matter rather than from symmetric nuclear matter to make use of recent quantum Monte Carlo (QMC) calculations. This method eschews the need to invoke the so-called parabolic approximation to extrapolate from symmetric nuclear matter, allowing us to directly constrain the neutron-matter EoS. Using a principal component analysis, we identify the combination of parameters most tightly constrained by observational data. We discuss sensitivity to various effects such as different component masses through population-model sensitivity, phase transitions in the core EoS, and large deviations from the central parameter values.
Using a Bayesian approach, we combine measurements of neutron star macroscopic observables obtained by astrophysical and gravitational observations, to derive joint constraints on the equation of state (EoS) of matter at supranuclear density. In our analysis we use two sets of data: (i) the masses and tidal deformabilities measured in the binary neutron star event GW170817, detected by LIGO and Virgo; (ii) the masses and stellar radii measured from observations of nuclear bursts in accreting low-mass X-ray binaries. Using two different parametrizations of the equation of state, we compute the posteriorprobability distributions of the EoS parameters, using which we infer the posterior distribution for the radius and the mass of the two neutron stars of GW170817. The constraints we set on the radii are tighter than previous bounds.
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.
X-ray pulse profile modeling of PSR J0740+6620, the most massive known pulsar, with data from the NICER and XMM-Newton observatories recently led to a measurement of its radius. We investigate this measurements implications for the neutron star equation of state (EoS), employing a nonparametric EoS model based on Gaussian processes and combining information from other x-ray, radio and gravitational-wave observations of neutron stars. Our analysis mildly disfavors EoSs that support a disconnected hybrid star branch in the mass-radius relation, a proxy for strong phase transitions, with a Bayes factor of $6.9$. For such EoSs, the transition mass from the hadronic to the hybrid branch is constrained to lie outside ($1,2$) $M_{odot}$. We also find that the conformal sound-speed bound is violated inside neutron star cores, which implies that the core matter is strongly interacting. The squared sound speed reaches a maximum of $0.75^{+0.25}_{-0.24}, c^2$ at $3.60^{+2.25}_{-1.89}$ times nuclear saturation density at 90% credibility. Since all but the gravitational-wave observations prefer a relatively stiff EoS, PSR J0740+6620s central density is only $3.57^{+1.3}_{-1.3}$ times nuclear saturation, limiting the density range probed by observations of cold, nonrotating neutron stars in $beta$-equilibrium.
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.