GW170817: observation of gravitational waves from a binary neutron star inspiral
TL;DR: The association of GRB 170817A, detected by Fermi-GBM 1.7 s after the coalescence, corroborates the hypothesis of a neutron star merger and provides the first direct evidence of a link between these mergers and short γ-ray bursts.
Abstract: On August 17, 2017 at 12∶41:04 UTC the Advanced LIGO and Advanced Virgo gravitational-wave detectors made their first observation of a binary neutron star inspiral. The signal, GW170817, was detected with a combined signal-to-noise ratio of 32.4 and a false-alarm-rate estimate of less than one per 8.0×10^{4} years. We infer the component masses of the binary to be between 0.86 and 2.26 M_{⊙}, in agreement with masses of known neutron stars. Restricting the component spins to the range inferred in binary neutron stars, we find the component masses to be in the range 1.17-1.60 M_{⊙}, with the total mass of the system 2.74_{-0.01}^{+0.04}M_{⊙}. The source was localized within a sky region of 28 deg^{2} (90% probability) and had a luminosity distance of 40_{-14}^{+8} Mpc, the closest and most precisely localized gravitational-wave signal yet. The association with the γ-ray burst GRB 170817A, detected by Fermi-GBM 1.7 s after the coalescence, corroborates the hypothesis of a neutron star merger and provides the first direct evidence of a link between these mergers and short γ-ray bursts. Subsequent identification of transient counterparts across the electromagnetic spectrum in the same location further supports the interpretation of this event as a neutron star merger. This unprecedented joint gravitational and electromagnetic observation provides insight into astrophysics, dense matter, gravitation, and cosmology.
Citations
More filters
TL;DR: A Bayesian statistical framework is developed and employed that enables us to implement constraints on the equation of state from laboratory measurements of nuclei and state-of-the-art chiral effective field theory methods and it is found that lower bounds on the neutron star tidal deformability will very strongly constrain microscopic models of the dense matter equations of state.
Abstract: We confront observational data from gravitational wave event GW170817 with microscopic modeling of the cold neutron star equation of state. We develop and employ a Bayesian statistical framework that enables us to implement constraints on the equation of state from laboratory measurements of nuclei and state-of-the-art chiral effective field theory methods. The energy density functionals constructed from the posterior probability distributions are then used to compute consistently the neutron star equation of state from the outer crust to the inner core, assuming a composition consisting of protons, neutrons, electrons, and muons. In contrast to previous studies, we find that the 95% credibility range of predicted neutron star tidal deformabilities (136<Λ<519) for a 1.4 solar-mass neutron star is already consistent with the upper bound deduced from observations of the GW170817 event. However, we find that lower bounds on the neutron star tidal deformability will very strongly constrain microscopic models of the dense matter equation of state. We also demonstrate a strong correlation between the neutron star tidal deformability and the pressure of beta-equilibrated matter at twice saturation density.
200 citations
TL;DR: In this paper, a nonparametric representation of the neutron-star equation of state based on Gaussian processes conditioned on nuclear theory models is proposed. But the model is not considered.
Abstract: Observations of neutron stars, whether in binaries or in isolation, provide information about the internal structure of the most extreme material objects in the Universe. In this work, we combine information from recent observations to place joint constraints on the properties of neutron star matter. We use (i) lower limits on the maximum mass of neutron stars obtained through radio observations of heavy pulsars, (ii) constraints on tidal properties inferred through the gravitational waves neutron star binaries emit as they coalesce, and (iii) information about neutron stars' masses and radii obtained through X-ray emission from surface hot spots. In order to combine information from such distinct messengers while avoiding the kind of modeling systematics intrinsic to parametric inference schemes, we employ a nonparametric representation of the neutron-star equation of state based on Gaussian processes conditioned on nuclear theory models. We find that existing astronomical observations imply ${R}_{1.4}=12.3{2}_{\ensuremath{-}1.47}^{+1.09}\text{ }\text{ }\mathrm{km}$ for the radius of a $1.4\text{ }\text{ }{\mathrm{M}}_{\ensuremath{\bigodot}}$ neutron star and $p(2{\ensuremath{\rho}}_{\mathrm{nuc}})=3.{8}_{\ensuremath{-}2.9}^{+2.7}\ifmmode\times\else\texttimes\fi{}1{0}^{34}\text{ }\text{ }\mathrm{dyn}/{\mathrm{cm}}^{2}$ for the pressure at twice nuclear saturation density at the 90% credible level. The upper bounds are driven by the gravitational wave observations, while X-ray and heavy pulsar observations drive the lower bounds. Additionally, we compute expected constraints from potential future astronomical observations and find that they can jointly determine ${R}_{1.4}$ to $\mathcal{O}(1)\text{ }\text{ }\mathrm{km}$ and $p(2{\ensuremath{\rho}}_{\mathrm{nuc}})$ to 80% relative uncertainty in the next five years.
199 citations
TL;DR: In this article, the effects of viscosity, which is likely to arise in the merger remnant due to magnetohydrodynamical turbulence, were considered for a massive neutron star surrounded by a torus, which was a canonical remnant formed after the binary neutron star merger.
Abstract: We perform long-term general relativistic neutrino radiation hydrodynamics simulations (in axisymmetry) for a massive neutron star (MNS) surrounded by a torus, which is a canonical remnant formed after the binary neutron star merger. We take into account the effects of viscosity, which is likely to arise in the merger remnant due to magnetohydrodynamical turbulence. As the initial condition, we employ the azimuthally averaged data of the MNS-torus system derived in a three-dimensional, numerical-relativity simulation for the binary neutron star merger. The viscous effect plays key roles for the remnant evolution and mass ejection from it in two phases of the evolution. In the first $t\lesssim10$ ms, a differential rotation state of the MNS is changed to a rigidly rotating state, and as a result, a sound wave, which subsequently becomes a shock wave, is formed in the vicinity of the MNS due to the variation of the quasi-equilibrium state of the MNS. The shock wave induces significant mass ejection of mass $\sim(0.5-2.0)\times 10^{-2}M_\odot$ for the alpha viscosity parameter of $0.01-0.04$. For the longer-term evolution with $\sim 0.1-10$ s, a significant fraction of the torus material is ejected. The ejecta mass is likely to be of order $10^{-2}M_\odot$, so that the total mass of the viscosity-driven ejecta could dominate that of the dynamical ejecta of mass $\lesssim 10^{-2}M_\odot$. The electron fraction, $Y_e$, of the ejecta is always high enough ($Y_e\gtrsim0.25$) that this post-merger ejecta is lanthanide-poor; hence, the opacity of the ejecta is likely to be $\sim 10-100$ times lower than that of the dynamical ejecta. This indicates that the electromagnetic signal from the ejecta would be rapidly evolving, bright, and blue if it is observed from a small viewing angle ($\lesssim 45^\circ$) for which the effect of the dynamical ejecta is minor.
198 citations
Institut des Hautes Études Scientifiques1, Istituto Nazionale di Fisica Nucleare2, University of Jena3, University of Parma4, University of Pisa5, University of Turin6, Max Planck Society7, University of Paris8, Moscow Institute of Physics and Technology9, IAC10, University of Milan11, University of Milano-Bicocca12, Princeton University13, University of Cambridge14, Cardiff University15, Albert Einstein Institution16, University of Valencia17, Sapienza University of Rome18, Radboud University Nijmegen19
TL;DR: In this paper, an effective one-body (EOB) waveform model for non-precessing (spin-aligned) and tidally interacting compact binaries is presented.
Abstract: We present TEOBResumS, a new effective-one-body (EOB) waveform model for nonprecessing (spin-aligned) and tidally interacting compact binaries. Spin-orbit and spin-spin effects are blended together by making use of the concept of centrifugal EOB radius. The point-mass sector through merger and ringdown is informed by numerical relativity (NR) simulations of binary black holes (BBHs) computed with the SpEC and bam codes. An improved, NR-based phenomenological description of the postmerger waveform is developed. The tidal sector of TEOBResumS describes the dynamics of neutron star binaries up to merger and incorporates a resummed attractive potential motivated by recent advances in the post-Newtonian and gravitational self-force description of relativistic tidal interactions. Equation-of-state-dependent self-spin interactions (monopole-quadrupole effects) are incorporated in the model using leading order post-Newtonian results in a new expression of the centrifugal radius. TEOBResumS is compared to 135 SpEC and 19 bam BBH waveforms. The maximum unfaithfulness to SpEC data ¯ F—at design Advanced LIGO sensitivity and evaluated with total mass M with a variance of 10M⊙≤M≤200M⊙—is always below 2.5×10−3 except for a single outlier that grazes the 7.1×10−3 level. When compared to bam data, ¯ F is smaller than 0.01 except for a single outlier in one of the corners of the NR-covered parameter space that reaches the 0.052 level. TEOBResumS is also compatible, up to merger, to high-end NR waveforms from binary neutron stars with spin effects and reduced initial eccentricity computed with the bam and thc codes. The data quality of binary neutron star waveforms is assessed via rigorous convergence tests from multiple resolution runs and takes into account systematic effects estimated by using the two independent high-order NR codes. The model is designed to generate accurate templates for the analysis of LIGO-Virgo data through merger and ringdown. We demonstrate its use by analyzing the publicly available data for GW150914.
198 citations
California Institute of Technology1, University of Chicago2, Lawrence Berkeley National Laboratory3, University of California, Berkeley4, Northwestern University5, Cardiff University6, Max Planck Society7, University of Copenhagen8, Queen's University Belfast9, Columbia University10, Rochester Institute of Technology11
TL;DR: In this paper, the authors used Gaussian Process Regression to model the light curves and spectra of the electromagnetic transient AT2017gfo and obtained a 90 per cent upper bound on the mass ratio q ≲ 1.38 and a lower bound Λ ≳ 197.
Abstract: The detection of the binary neutron star merger GW170817 together with the observation of electromagnetic counterparts across the entire spectrum inaugurated a new era of multimessenger astronomy. In this study, we incorporate wavelength-dependent opacities and emissivities calculated from atomic-structure data enabling us to model both the measured light curves and spectra of the electromagnetic transient AT2017gfo. Best fits of the observational data are obtained by Gaussian Process Regression, which allows us to present posterior samples for the kilonova and source properties connected to GW170817. Incorporating constraints obtained from the gravitational wave signal measured by the LIGO-Virgo Scientific Collaboration, we present a 90 per cent upper bound on the mass ratio q ≲ 1.38 and a lower bound on the tidal deformability of Λ ≳ 197, which rules out sufficiently soft equations of state. Our analysis is a path-finder for more realistic kilonova models and shows how the combination of gravitational wave and electromagnetic measurements allow for stringent constraints on the source parameters and the supranuclear equation of state.
197 citations
Cites background or methods from "GW170817: observation of gravitatio..."
...These results can be compared to estimates obtained from a reanalysis of GW170817 (De et al. 2018), which incorporates quasi-universal relations for the tidal deformability and obtains 90% lower bounds on the tidal deformability Λ̃ & 117 and 90% upper bounds on the mass ratio q ....
[...]
...The detection of the binary neutron star GW170817 together with the observation of electromagnetic counterparts across the entire spectrum inaugurated a new era of multi-messenger astronomy....
[...]
...Best-fits of the observational data are obtained by Gaussian Process Regression, which allows us to present posterior samples for the kilonova and source properties connected to GW170817....
[...]
...We initially took the photometry from the UV to K−band from (Andreoni et al. 2017; Arcavi et al. 2017; R. Chornock et al. 2017; Cowperthwaite et al. 2017; Drout et al. 2017; Evans et al. 2017; Kasliwal et al. 2017; Tanvir et al. 2017; Pian et al. 2017; Troja et al. 2017; Smartt et al. 2017; Utsumi et al. 2017; Valenti et al. 2017) from phases +0.467d to +25.19d after GW170817 and at each epoch created the broadest spectral energy distribution possible....
[...]
...In this article, we obtained constraints on the GW170817 progenitors mass ratio and tidal deformability, which are more stringent than those obtained purely from gravitational-wave observations....
[...]
References
More filters
TL;DR: In this article, the authors present a cosmological analysis based on full-mission Planck observations of temperature and polarization anisotropies of the cosmic microwave background (CMB) radiation.
Abstract: This paper presents cosmological results based on full-mission Planck observations of temperature and polarization anisotropies of the cosmic microwave background (CMB) radiation. Our results are in very good agreement with the 2013 analysis of the Planck nominal-mission temperature data, but with increased precision. The temperature and polarization power spectra are consistent with the standard spatially-flat 6-parameter ΛCDM cosmology with a power-law spectrum of adiabatic scalar perturbations (denoted “base ΛCDM” in this paper). From the Planck temperature data combined with Planck lensing, for this cosmology we find a Hubble constant, H0 = (67.8 ± 0.9) km s-1Mpc-1, a matter density parameter Ωm = 0.308 ± 0.012, and a tilted scalar spectral index with ns = 0.968 ± 0.006, consistent with the 2013 analysis. Note that in this abstract we quote 68% confidence limits on measured parameters and 95% upper limits on other parameters. We present the first results of polarization measurements with the Low Frequency Instrument at large angular scales. Combined with the Planck temperature and lensing data, these measurements give a reionization optical depth of τ = 0.066 ± 0.016, corresponding to a reionization redshift of . These results are consistent with those from WMAP polarization measurements cleaned for dust emission using 353-GHz polarization maps from the High Frequency Instrument. We find no evidence for any departure from base ΛCDM in the neutrino sector of the theory; for example, combining Planck observations with other astrophysical data we find Neff = 3.15 ± 0.23 for the effective number of relativistic degrees of freedom, consistent with the value Neff = 3.046 of the Standard Model of particle physics. The sum of neutrino masses is constrained to ∑ mν < 0.23 eV. The spatial curvature of our Universe is found to be very close to zero, with | ΩK | < 0.005. Adding a tensor component as a single-parameter extension to base ΛCDM we find an upper limit on the tensor-to-scalar ratio of r0.002< 0.11, consistent with the Planck 2013 results and consistent with the B-mode polarization constraints from a joint analysis of BICEP2, Keck Array, and Planck (BKP) data. Adding the BKP B-mode data to our analysis leads to a tighter constraint of r0.002 < 0.09 and disfavours inflationarymodels with a V(φ) ∝ φ2 potential. The addition of Planck polarization data leads to strong constraints on deviations from a purely adiabatic spectrum of fluctuations. We find no evidence for any contribution from isocurvature perturbations or from cosmic defects. Combining Planck data with other astrophysical data, including Type Ia supernovae, the equation of state of dark energy is constrained to w = −1.006 ± 0.045, consistent with the expected value for a cosmological constant. The standard big bang nucleosynthesis predictions for the helium and deuterium abundances for the best-fit Planck base ΛCDM cosmology are in excellent agreement with observations. We also constraints on annihilating dark matter and on possible deviations from the standard recombination history. In neither case do we find no evidence for new physics. The Planck results for base ΛCDM are in good agreement with baryon acoustic oscillation data and with the JLA sample of Type Ia supernovae. However, as in the 2013 analysis, the amplitude of the fluctuation spectrum is found to be higher than inferred from some analyses of rich cluster counts and weak gravitational lensing. We show that these tensions cannot easily be resolved with simple modifications of the base ΛCDM cosmology. Apart from these tensions, the base ΛCDM cosmology provides an excellent description of the Planck CMB observations and many other astrophysical data sets.
10,728 citations
TL;DR: In this paper, the authors present results based on full-mission Planck observations of temperature and polarization anisotropies of the CMB, which are consistent with the six-parameter inflationary LCDM cosmology.
Abstract: We present results based on full-mission Planck observations of temperature and polarization anisotropies of the CMB. These data are consistent with the six-parameter inflationary LCDM cosmology. From the Planck temperature and lensing data, for this cosmology we find a Hubble constant, H0= (67.8 +/- 0.9) km/s/Mpc, a matter density parameter Omega_m = 0.308 +/- 0.012 and a scalar spectral index with n_s = 0.968 +/- 0.006. (We quote 68% errors on measured parameters and 95% limits on other parameters.) Combined with Planck temperature and lensing data, Planck LFI polarization measurements lead to a reionization optical depth of tau = 0.066 +/- 0.016. Combining Planck with other astrophysical data we find N_ eff = 3.15 +/- 0.23 for the effective number of relativistic degrees of freedom and the sum of neutrino masses is constrained to < 0.23 eV. Spatial curvature is found to be |Omega_K| < 0.005. For LCDM we find a limit on the tensor-to-scalar ratio of r <0.11 consistent with the B-mode constraints from an analysis of BICEP2, Keck Array, and Planck (BKP) data. Adding the BKP data leads to a tighter constraint of r < 0.09. We find no evidence for isocurvature perturbations or cosmic defects. The equation of state of dark energy is constrained to w = -1.006 +/- 0.045. Standard big bang nucleosynthesis predictions for the Planck LCDM cosmology are in excellent agreement with observations. We investigate annihilating dark matter and deviations from standard recombination, finding no evidence for new physics. The Planck results for base LCDM are in agreement with BAO data and with the JLA SNe sample. However the amplitude of the fluctuations is found to be higher than inferred from rich cluster counts and weak gravitational lensing. Apart from these tensions, the base LCDM cosmology provides an excellent description of the Planck CMB observations and many other astrophysical data sets.
9,745 citations
TL;DR: This is the first direct detection of gravitational waves and the first observation of a binary black hole merger, and these observations demonstrate the existence of binary stellar-mass black hole systems.
Abstract: On September 14, 2015 at 09:50:45 UTC the two detectors of the Laser Interferometer Gravitational-Wave Observatory simultaneously observed a transient gravitational-wave signal. The signal sweeps upwards in frequency from 35 to 250 Hz with a peak gravitational-wave strain of $1.0 \times 10^{-21}$. It matches the waveform predicted by general relativity for the inspiral and merger of a pair of black holes and the ringdown of the resulting single black hole. The signal was observed with a matched-filter signal-to-noise ratio of 24 and a false alarm rate estimated to be less than 1 event per 203 000 years, equivalent to a significance greater than 5.1 {\sigma}. The source lies at a luminosity distance of $410^{+160}_{-180}$ Mpc corresponding to a redshift $z = 0.09^{+0.03}_{-0.04}$. In the source frame, the initial black hole masses are $36^{+5}_{-4} M_\odot$ and $29^{+4}_{-4} M_\odot$, and the final black hole mass is $62^{+4}_{-4} M_\odot$, with $3.0^{+0.5}_{-0.5} M_\odot c^2$ radiated in gravitational waves. All uncertainties define 90% credible intervals.These observations demonstrate the existence of binary stellar-mass black hole systems. This is the first direct detection of gravitational waves and the first observation of a binary black hole merger.
9,596 citations
Journal Article•
TL;DR: The first direct detection of gravitational waves and the first observation of a binary black hole merger were reported in this paper, with a false alarm rate estimated to be less than 1 event per 203,000 years, equivalent to a significance greater than 5.1σ.
Abstract: On September 14, 2015 at 09:50:45 UTC the two detectors of the Laser Interferometer Gravitational-Wave Observatory simultaneously observed a transient gravitational-wave signal. The signal sweeps upwards in frequency from 35 to 250 Hz with a peak gravitational-wave strain of 1.0×10(-21). It matches the waveform predicted by general relativity for the inspiral and merger of a pair of black holes and the ringdown of the resulting single black hole. The signal was observed with a matched-filter signal-to-noise ratio of 24 and a false alarm rate estimated to be less than 1 event per 203,000 years, equivalent to a significance greater than 5.1σ. The source lies at a luminosity distance of 410(-180)(+160) Mpc corresponding to a redshift z=0.09(-0.04)(+0.03). In the source frame, the initial black hole masses are 36(-4)(+5)M⊙ and 29(-4)(+4)M⊙, and the final black hole mass is 62(-4)(+4)M⊙, with 3.0(-0.5)(+0.5)M⊙c(2) radiated in gravitational waves. All uncertainties define 90% credible intervals. These observations demonstrate the existence of binary stellar-mass black hole systems. This is the first direct detection of gravitational waves and the first observation of a binary black hole merger.
4,375 citations
TL;DR: This second gravitational-wave observation provides improved constraints on stellar populations and on deviations from general relativity.
Abstract: We report the observation of a gravitational-wave signal produced by the coalescence of two stellar-mass black holes. The signal, GW151226, was observed by the twin detectors of the Laser Interferometer Gravitational-Wave Observatory (LIGO) on December 26, 2015 at 03:38:53 UTC. The signal was initially identified within 70 s by an online matched-filter search targeting binary coalescences. Subsequent off-line analyses recovered GW151226 with a network signal-to-noise ratio of 13 and a significance greater than 5 σ. The signal persisted in the LIGO frequency band for approximately 1 s, increasing in frequency and amplitude over about 55 cycles from 35 to 450 Hz, and reached a peak gravitational strain of 3.4+0.7−0.9×10−22. The inferred source-frame initial black hole masses are 14.2+8.3−3.7M⊙ and 7.5+2.3−2.3M⊙ and the final black hole mass is 20.8+6.1−1.7M⊙. We find that at least one of the component black holes has spin greater than 0.2. This source is located at a luminosity distance of 440+180−190 Mpc corresponding to a redshift 0.09+0.03−0.04. All uncertainties define a 90 % credible interval. This second gravitational-wave observation provides improved constraints on stellar populations and on deviations from general relativity.
3,448 citations