skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 10:00 PM ET on Friday, December 8 until 2:00 AM ET on Saturday, December 9 due to maintenance. We apologize for the inconvenience.

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. 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
  3. 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
  4. Abstract

    The Neutron Star Interior Composition Explorer (NICER) records exceptional data on pulsars’ energy-dependent X-ray pulse profiles. However, in searching for evidence of pulsations, Guillot et al. (2019) introduce a procedure to select an ordered subset of data that maximizes a detection statistic (the H-test). I show that this can degrade subsequent analyses using an idealized model with stationary expected count rates from both noise and signal. Specifically, the data-selection procedure biases the inferred mean count rate to be too low and the inferred pulsation amplitude to be too high, and the size of these biases scales strongly with the amount of data that is rejected and the true signal amplitude. The procedure also alters the H-test’s null distribution, rendering nominal significance estimates overly optimistic. While the idealized model does not capture all the complexities of real NICER data, it suggests that these biases could be important for NICER’s observations of J0740+6620 and other faint pulsars (observations of J0030+0451 are likely less affected). I estimate that these effects may introduce a bias of(10%)on average in the inferred modulation depth of lightcurves like J0740+6620's, and may be as large as(50%)for fainter pulsars. However, the change for a single data set like J0740+6620 is expected to be a shift between −5% and +20%. This could imply that the lower limit on J0740+6620's radius is slightly larger than it should be, although preliminary investigations suggest the radius constraints shift to larger radii by(1%)with the same overall statistical precision using real J0740+6620 data.

    more » « less
  5. Abstract

    We present the stellar population properties of 69 short gamma-ray burst (GRB) host galaxies, representing the largest uniformly modeled sample to date. Using theProspectorstellar population inference code, we jointly fit photometry and/or spectroscopy of each host galaxy. We find a population median redshift ofz=0.640.32+0.83(68% confidence), including nine photometric redshifts atz≳ 1. We further find a median mass-weighted age oftm=0.80.53+2.71Gyr, stellar mass of log(M*/M) =9.690.65+0.75, star formation rate of SFR =1.441.35+9.37Myr−1, stellar metallicity of log(Z*/Z) =0.380.42+0.44, and dust attenuation ofAV=0.430.36+0.85mag (68% confidence). Overall, the majority of short GRB hosts are star-forming (≈84%), with small fractions that are either transitioning (≈6%) or quiescent (≈10%); however, we observe a much larger fraction (≈40%) of quiescent and transitioning hosts atz≲ 0.25, commensurate with galaxy evolution. We find that short GRB hosts populate the star-forming main sequence of normal field galaxies, but do not include as many high-mass galaxies as the general galaxy population, implying that their binary neutron star (BNS) merger progenitors are dependent on a combination of host star formation and stellar mass. The distribution of ages and redshifts implies a broad delay-time distribution, with a fast-merging channel atz> 1 and a decreased neutron star binary formation efficiency from high to low redshifts. If short GRB hosts are representative of BNS merger hosts within the horizon of current gravitational wave detectors, these results can inform future searches for electromagnetic counterparts. All of the data and modeling products are available on the Broadband Repository for Investigating Gamma-ray burst Host Traits website.

    more » « less