skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: AT2017gfo: Bayesian inference and model selection of multicomponent kilonovae and constraints on the neutron star equation of state
ABSTRACT The joint detection of the gravitational wave GW170817, of the short γ-ray burst GRB170817A and of the kilonova AT2017gfo, generated by the the binary neutron star (NS) merger observed on 2017 August 17, is a milestone in multimessenger astronomy and provides new constraints on the NS equation of state. We perform Bayesian inference and model selection on AT2017gfo using semi-analytical, multicomponents models that also account for non-spherical ejecta. Observational data favour anisotropic geometries to spherically symmetric profiles, with a log-Bayes’ factor of ∼104, and favour multicomponent models against single-component ones. The best-fitting model is an anisotropic three-component composed of dynamical ejecta plus neutrino and viscous winds. Using the dynamical ejecta parameters inferred from the best-fitting model and numerical–relativity relations connecting the ejecta properties to the binary properties, we constrain the binary mass ratio to q < 1.54 and the reduced tidal parameter to $$120\lt \tilde{\Lambda }\lt 1110$$. Finally, we combine the predictions from AT2017gfo with those from GW170817, constraining the radius of a NS of 1.4 M⊙ to 12.2 ± 0.5 km (1σ level). This prediction could be further strengthened by improving kilonova models with numerical-relativity information.  more » « less
Award ID(s):
2011725 2020275 1630782
PAR ID:
10250779
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
505
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
1661 to 1677
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract GW170817 is the first binary neutron star (NS) merger detected in gravitational waves (GWs) and photons, and so far remains the only GW event of its class with a definitive electromagnetic counterpart. Radio emission from the structured jet associated with GW170817 has faded below the sensitivity achievable via deep radio observations with the most sensitive radio arrays currently in operation. Hence, we now have the opportunity to probe the radio re-brightening that some models predict, which should emerge at late times from the interaction of the dynamically stripped merger ejecta with the interstellar medium. Here we present the latest results from our deep radio observations of the GW170817 field with the Karl G. Jansky Very Large Array (VLA), 4.5 yr after the merger. Our new data at 3 GHz do not show any compelling evidence for emission in excess to the tail of the jet afterglow (<3.3 μ Jy), confirming our previous results. We thus set new constraints on the dynamical ejecta afterglow models. These constraints favor single-speed ejecta with energies ≲10 50 erg (for an ejecta speed of β 0 = 0.5), or steeper energy–speed distributions of the kilonova ejecta. Our results also suggest larger values of the cold, nonrotating maximum NS mass in equal-mass scenarios. However, without a detection of the dynamical ejecta afterglow, obtaining precise constraints on the NS equation of state remains challenging. 
    more » « less
  2. Abstract We study the production of very light elements (Z< 20) in the dynamical and spiral-wave wind ejecta of binary neutron star mergers by combining detailed nucleosynthesis calculations with the outcome of numerical relativity merger simulations. All our models are targeted to GW170817 and include neutrino radiation. We explore different finite-temperature, composition-dependent nuclear equations of state, and binary mass ratios, and find that hydrogen and helium are the most abundant light elements. For both elements, the decay of free neutrons is the driving nuclear reaction. In particular, ∼0.5–2 × 10−6Mof hydrogen are produced in the fast expanding tail of the dynamical ejecta, while ∼1.5–11 × 10−6Mof helium are synthesized in the bulk of the dynamical ejecta, usually in association with heavyr-process elements. By computing synthetic spectra, we find that the possibility of detecting hydrogen and helium features in kilonova spectra is very unlikely for fiducial masses and luminosities, even when including nonlocal thermodynamic equilibrium effects. The latter could be crucial to observe helium lines a few days after merger for faint kilonovae or for luminous kilonovae ejecting large masses of helium. Finally, we compute the amount of strontium synthesized in the dynamical and spiral-wave wind ejecta, and find that it is consistent with (or even larger than, in the case of a long-lived remnant) the one required to explain early spectral features in the kilonova of GW170817. 
    more » « less
  3. Kilonovae, one source of electromagnetic emission associated with neutron star mergers, are powered by the decay of radioactive isotopes in the neutron-rich merger ejecta. Models for kilonova emission consistent with the electromagnetic counterpart to GW170817 predict characteristic abundance patterns, determined by the relative balance of different types of material in the outflow. Assuming that the observed source is prototypical, this inferred abundance pattern in turn must matchr-process abundances deduced by other means, such as what is observed in the solar system. We report on analysis comparing the input mass-weighted elemental compositions adopted in our radiative transfer simulations to the mass fractions of elements in the Sun, as a practical prototype for the potentially universal abundance signature from neutron star mergers. We characterize the extent to which our parameter inference results depend on our assumed composition for the dynamical and wind ejecta and examine how the new results compare to previous work. We find that a dynamical ejecta composition calculated using the FRDM2012 nuclear mass and FRLDM fission models with extremely neutron-rich ejecta (Ye= 0.035) along with moderately neutron-rich (Ye= 0.27) wind ejecta composition yields a wind-to-dynamical mass ratio ofMw/Md= 0.47, which best matches the observed AT2017gfo kilonova light curves while also producing the best-matching abundance of neutron capture elements in the solar system, though, allowing for systematics, the ratio may be as high as of order unity. 
    more » « less
  4. null (Ed.)
    ABSTRACT The detection of the optical transient AT2017gfo proved that binary neutron star mergers are progenitors of kilonovae (KNe). Using a combination of numerical-relativity and radiative-transfer simulations, the community has developed sophisticated models for these transients for a wide portion of the expected parameter space. Using these simulations and surrogate models made from them, it has been possible to perform Bayesian inference of the observed signals to infer properties of the ejected matter. It has been pointed out that combining inclination constraints derived from the KN with gravitational-wave measurements increases the accuracy with which binary parameters can be estimated, in particular breaking the distance-inclination degeneracy from gravitational wave inference. To avoid bias from the unknown ejecta geometry, constraints on the inclination angle for AT2017gfo should be insensitive to the employed models. In this work, we compare different assumptions about the ejecta and radiative reprocesses used by the community and we investigate their impact on the parameter inference. While most inferred parameters agree, we find disagreement between posteriors for the inclination angle for different geometries that have been used in the current literature. According to our study, the inclusion of reprocessing of the photons between different ejecta types improves the modeling fits to AT2017gfo and, in some cases, affects the inferred constraints. Our study motivates the inclusion of large ∼ 1-mag uncertainties in the KN models employed for Bayesian analysis to capture yet unknown systematics, especially when inferring inclination angles, although smaller uncertainties seem appropriate to capture model systematics for other intrinsic parameters. We can use this method to impose soft constraints on the ejecta geometry of the KN AT2017gfo. 
    more » « less
  5. 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. 
    more » « less