skip to main content

Title: The impact of the crust equation of state on the analysis of GW170817

The detection of GW170817, the first neutron star-neutron star merger observed by Advanced LIGO and Virgo, and its following analyses represent the first contributions of gravitational wave data to understanding dense matter. Parameterizing the high density section of the equation of state of both neutron stars through spectral decomposition, and imposing a lower limit on the maximum mass value, led to an estimate of the stars’ radii ofkm andkm (Abbottet al2018Phys. Rev. Lett.121161101). These values do not, however, take into account any uncertainty owed to the choice of the crust low-density equation of state, which was fixed to reproduce the SLy equation of state model (Douchin and Haensel 2001Astron. Astrophys.380151). We here re-analyze GW170817 data and establish that different crust models do not strongly impact the mass or tidal deformability of a neutron star—it is impossible to distinguish between low-density models with gravitational wave analysis. However, the crust does have an effect on inferred radius. We predict the systematic error due to this effect using neutron star structure equations, and compare the prediction to results from full parameter estimation runs. For GW170817, this systematic error affects the radius estimate by 0.3 km, approximatelyof the neutron stars’ radii.

more » « less
Award ID(s):
Author(s) / Creator(s):
; ;
Publisher / Repository:
IOP Publishing
Date Published:
Journal Name:
Classical and Quantum Gravity
Page Range / eLocation ID:
Article No. 025008
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We demonstrate how to quantify the frequency-domain amplitude and phase accuracy of waveform models,δAandδφ, in a form that could be marginalized over in gravitational-wave inference using techniques currently applied for quantifying calibration uncertainty. For concreteness, waveform uncertainties affecting neutron-star inspiral measurements are considered, and post-hoc error estimates from a variety of waveform models are made by comparing time-domain and frequency-domain analytic models with multiple-resolution numerical simulations. These waveform uncertainty estimates can be compared to GW170817 calibration envelopes or to Advanced LIGO and Virgo calibration goals. Signal-specific calibration and waveform uncertainties are compared to statistical fluctuations in gravitational-wave observatories, giving frequency-dependent modeling requirements for detectors such as Advanced LIGO Plus, Cosmic Explorer, or Einstein Telescope. Finally, the distribution of waveform error for the GW170817 posterior is computed from tidal models and compared to the constraints onδφorδAfrom GWTC-1 by Edelmanet al.In general,δφandδAcan also be interpreted in terms of unmodeled astrophysical energy transfer within or from the source system.

    more » « less
  2. An advanced LIGO and Virgo’s third observing run brought another binary neutron star merger (BNS) and the first neutron-star black hole mergers. While no confirmed kilonovae were identified in conjunction with any of these events, continued improvements of analyses surrounding GW170817 allow us to project constraints on the Hubble Constant (H0), the Galactic enrichment fromr-process nucleosynthesis, and ultra-dense matter possible from forthcoming events. Here, we describe the expected constraints based on the latest expected event rates from the international gravitational-wave network and analyses of GW170817. We show the expected detection rate of gravitational waves and their counterparts, as well as how sensitive potential constraints are to the observed numbers of counterparts. We intend this analysis as support for the community when creating scientifically driven electromagnetic follow-up proposals. During the next observing run O4, we predict an annual detection rate of electromagnetic counterparts from BNS of0.430.26+0.58(1.971.2+2.68) for the Zwicky Transient Facility (Rubin Observatory).

    more » « less
  3. Abstract

    The multi-messenger detection of the gravitational-wave signal GW170817, the corresponding kilonova AT2017gfo and the short gamma-ray burst GRB170817A, as well as the observed afterglow has delivered a scientific breakthrough. For an accurate interpretation of all these different messengers, one requires robust theoretical models that describe the emitted gravitational-wave, the electromagnetic emission, and dense matter reliably. In addition, one needs efficient and accurate computational tools to ensure a correct cross-correlation between the models and the observational data. For this purpose, we have developed the Nuclear-physics and Multi-Messenger Astrophysics framework NMMA. The code allows incorporation of nuclear-physics constraints at low densities as well as X-ray and radio observations of isolated neutron stars. In previous works, the NMMA code has allowed us to constrain the equation of state of supranuclear dense matter, to measure the Hubble constant, and to compare dense-matter physics probed in neutron-star mergers and in heavy-ion collisions, and to classify electromagnetic observations and perform model selection. Here, we show an extension of the NMMA code as a first attempt of analyzing the gravitational-wave signal, the kilonova, and the gamma-ray burst afterglow simultaneously. Incorporating all available information, we estimate the radius of a 1.4Mneutron star to be$$R=11.9{8}_{-0.40}^{+0.35}$$R=11.980.40+0.35km.

    more » « less
  4. Abstract

    We measure the molecular-to-atomic gas ratio,Rmol, and the star formation rate (SFR) per unit molecular gas mass, SFEmol, in 38 nearby galaxies selected from the Virgo Environment Traced in CO (VERTICO) survey. We stack ALMA12CO (J= 2−1) spectra coherently using Hivelocities from the VIVA survey to detect faint CO emission out to galactocentric radiirgal∼ 1.2r25. We determine the scale lengths for the molecular and stellar components, finding a ∼3:5 relation compared to ∼1:1 in field galaxies, indicating that the CO emission is more centrally concentrated than the stars. We computeRmolas a function of different physical quantities. While the spatially resolvedRmolon average decreases with increasing radius, we find that the mean molecular-to-atomic gas ratio within the stellar effective radiusRe,Rmol(r<Re), shows a systematic increase with the level of Hi, truncation and/or asymmetry (HIperturbation). Analysis of the molecular- and the atomic-to-stellar mass ratios withinRe,Rmol(r<Re)andRatom(r<Re), shows that VERTICO galaxies have increasingly lowerRatom(r<Re)for larger levels of HIperturbation (compared to field galaxies matched in stellar mass), but no significant change inRmol(r<Re). We also measure a clear systematic decrease of the SFEmolwithinRe, SFEmol(r<Re), with increasingly perturbed Hi. Therefore, compared to field galaxies from the field, VERTICO galaxies are more compact in CO emission in relation to their stellar distribution, but increasingly perturbed atomic gas increases theirRmoland decreases the efficiency with which their molecular gas forms stars.

    more » « less
  5. Abstract

    We present the second data release of gravitational waveforms from binary neutron star (BNS) merger simulations performed by the Computational Relativity (CoRe) collaboration. The current database consists of 254 different BNS configurations and a total of 590 individual numerical-relativity simulations using various grid resolutions. The released waveform data contain the strain and the Weyl curvature multipoles up to=m=4. They span a significant portion of the mass, mass-ratio, spin and eccentricity parameter space and include targeted configurations to the events GW170817 and GW190425.CoResimulations are performed with 18 different equations of state, seven of which are finite temperature models, and three of which account for non-hadronic degrees of freedom. About half of the released data are computed with high-order hydrodynamics schemes for tens of orbits to merger; the other half is computed with advanced microphysics. We showcase a standard waveform error analysis and discuss the accuracy of the database in terms of faithfulness. We present ready-to-use fitting formulas for equation of state-insensitive relations at merger (e.g. merger frequency), luminosity peak, and post-merger spectrum.

    more » « less