skip to main content


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

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):
1806962
PAR ID:
10303213
Author(s) / Creator(s):
; ;
Publisher / Repository:
IOP Publishing
Date Published:
Journal Name:
Classical and Quantum Gravity
Volume:
37
Issue:
2
ISSN:
0264-9381
Page Range / eLocation ID:
Article No. 025008
Format(s):
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

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

    We present measurements of [Fe/H] and [α/Fe] for 128 individual red giant branch stars (RGB) in the stellar halo of M31, including its Giant Stellar Stream (GSS), obtained using spectral synthesis of low- and medium-resolution Keck/DEIMOS spectroscopy (and 6000, respectively). We observed four fields in M31's stellar halo (at projected radii of 9, 18, 23, and 31 kpc), as well as two fields in the GSS (at 33 kpc). In combination with existing literature measurements, we have increased the sample size of [Fe/H] and [α/Fe] measurements from 101 to a total of 229 individual M31 RGB stars. From this sample, we investigate the chemical abundance properties of M31's inner halo, findingand. Between 8 and 34 kpc, the inner halo has a steep [Fe/H] gradient (−0.025 ± 0.002 dex kpc−1) and negligible [α/Fe] gradient, where substructure in the inner halo is systematically more metal-rich than the smooth component of the halo at a given projected distance. Although the chemical abundances of the inner stellar halo are largely inconsistent with that of present-day dwarf spheroidal (dSph) satellite galaxies of M31, we identified 22 RGB stars kinematically associated with the smooth component of the stellar halo that have chemical abundance patterns similar to M31 dSphs. We discuss formation scenarios for M31's halo, concluding that these dSph-like stars may have been accreted from galaxies of similar stellar mass and star formation history, or of higher stellar mass and similar star formation efficiency.

     
    more » « less
  5. Abstract

    We develop a Newtonian model of a deep tidal disruption event (TDE), for which the pericenter distance of the star,rp, is well within the tidal radius of the black hole,rt, i.e., whenβrt/rp≫ 1. We find that shocks form forβ≳ 3, but they are weak (with Mach numbers ∼1) for allβ, and that they reach the center of the star prior to the time of maximum adiabatic compression forβ≳ 10. The maximum density and temperature reached during the TDE follow much shallower relations withβthan the previously predictedρmaxβ3andTmaxβ2scalings. Belowβ≃ 10, this shallower dependence occurs because the pressure gradient is dynamically significant before the pressure is comparable to the ram pressure of the free-falling gas, while aboveβ≃ 10, we find that shocks prematurely halt the compression and yield the scalingsρmaxβ1.62andTmaxβ1.12. We find excellent agreement between our results and high-resolution simulations. Our results demonstrate that, in the Newtonian limit, the compression experienced by the star is completely independent of the mass of the black hole. We discuss our results in the context of existing (affine) models, polytropic versus non-polytropic stars, and general relativistic effects, which become important when the pericenter of the star nears the direct capture radius, atβ∼ 12.5 (2.7) for a solar-like star disrupted by a 106M(107M) supermassive black hole.

     
    more » « less