No Arabic abstract
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.
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.
The LIGO/Virgo detection of gravitational waves originating from a neutron-star merger, GW170817, has recently provided new stringent limits on the tidal deformabilities of the stars involved in the collision. Combining this measurement with the existence of two-solar-mass stars, we generate a generic family of neutron-star-matter Equations of State (EoSs) that interpolate between state-of-the-art theoretical results at low and high baryon density. Comparing the results to ones obtained without the tidal-deformability constraint, we witness a dramatic reduction in the family of allowed EoSs. Based on our analysis, we conclude that the maximal radius of a 1.4-solar-mass neutron star is 13.6 km, and that smallest allowed tidal deformability of a similar-mass star is $Lambda(1.4 M_odot) = 120$.
We present a proof-of-concept study, based on numerical-relativity simulations, of how gravitational waves (GWs) from neutron star merger remnants can probe the nature of matter at extreme densities. Phase transitions and extra degrees of freedom can emerge at densities beyond those reached during the inspiral, and typically result in a softening of the equation of state (EOS). We show that such physical effects change the qualitative dynamics of the remnant evolution, but they are not identifiable as a signature in the GW frequency, with the exception of possible black-hole formation effects. The EOS softening is, instead, encoded in the GW luminosity and phase and is in principle detectable up to distances of the order of several Mpcs with advanced detectors and up to hundreds of Mpcs with third generation detectors. Probing extreme-density matter will require going beyond the current paradigm and developing a more holistic strategy for modeling and analyzing postmerger GW signals.
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.
Neutron-star (NS) merger simulations are conducted for 38 representative microphysical descriptions of high-density matter in order to explore the equation-of-state dependence of the postmerger ring-down phase. The formation of a deformed, oscillating, differentially rotating very massive NS is the typical outcome of the coalescence of two stars with 1.35 $M_{odot}$ for most candidate EoSs. The oscillations of this object imprint a pronounced peak in the gravitational-wave (GW) spectra, which is used to characterize the emission for a given model. The peak frequency of this postmerger GW signal correlates very well with the radii of nonrotating NSs, and thus allows to constrain the high-density EoS by a GW detection. In the case of 1.35-1.35 $M_{odot}$ mergers the peak frequency scales particularly well with the radius of a NS with 1.6 $M_{odot}$, where the maximum deviation from this correlation is only 60 meters for fully microphysical EoSs which are compatible with NS observations. Combined with the uncertainty in the determination of the peak frequency it appears likely that a GW detection can measure the radius of a 1.6 $M_{odot}$ NS with an accuracy of about 100 to 200 meters. We also uncover relations of the peak frequency with the radii of nonrotating NSs with 1.35 $M_{odot}$ or 1.8 $M_{odot}$, with the radius or the central energy density of the maximum-mass Tolman-Oppenheimer-Volkoff configuration, and with the pressure or sound speed at a fiducial rest-mass density of about twice nuclear saturation density. Furthermore, it is found that a determination of the dominant postmerger GW frequency can provide an upper limit for the maximum mass of nonrotating NSs. The prospects for a detection of the postmerger GW signal and a determination of the dominant GW frequency are estimated to be in the range of 0.015 to 1.2 events per year with the upcoming Advanced LIGO detector.