skip to main content

Title: Inferring Kilonova Population Properties with a Hierarchical Bayesian Framework. I. Nondetection Methodology and Single-event Analyses

We presentnimbus: a hierarchical Bayesian framework to infer the intrinsic luminosity parameters of kilonovae (KNe) associated with gravitational-wave (GW) events, based purely on nondetections. This framework makes use of GW 3D distance information and electromagnetic upper limits from multiple surveys for multiple events and self-consistently accounts for the finite sky coverage and probability of astrophysical origin. The framework is agnostic to the brightness evolution assumed and can account for multiple electromagnetic passbands simultaneously. Our analyses highlight the importance of accounting for model selection effects, especially in the context of nondetections. We show our methodology using a simple, two-parameter linear brightness model, taking the follow-up of GW190425 with the Zwicky Transient Facility as a single-event test case for two different prior choices of model parameters: (i) uniform/uninformative priors and (ii) astrophysical priors based on surrogate models of Monte Carlo radiative-transfer simulations of KNe. We present results under the assumption that the KN is within the searched region to demonstrate functionality and the importance of prior choice. Our results show consistency withsimsurvey—an astronomical survey simulation tool used previously in the literature to constrain the population of KNe. While our results based on uniform priors strongly constrain the parameter space, those based more » on astrophysical priors are largely uninformative, highlighting the need for deeper constraints. Future studies with multiple events having electromagnetic follow-up from multiple surveys should make it possible to constrain the KN population further.

« less
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
The Astrophysical Journal
Page Range or eLocation-ID:
Article No. 58
DOI PREFIX: 10.3847
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT Joint multimessenger observations with gravitational waves and electromagnetic (EM) data offer new insights into the astrophysical studies of compact objects. The third Advanced LIGO and Advanced Virgo observing run began on 2019 April 1; during the 11 months of observation, there have been 14 compact binary systems candidates for which at least one component is potentially a neutron star. Although intensive follow-up campaigns involving tens of ground and space-based observatories searched for counterparts, no EM counterpart has been detected. Following on a previous study of the first six months of the campaign, we present in this paper the next five months of the campaign from 2019 October to 2020 March. We highlight two neutron star–black hole candidates (S191205ah and S200105ae), two binary neutron star candidates (S191213g and S200213t), and a binary merger with a possible neutron star and a ‘MassGap’ component, S200115j. Assuming that the gravitational-wave (GW) candidates are of astrophysical origin and their location was covered by optical telescopes, we derive possible constraints on the matter ejected during the events based on the non-detection of counterparts. We find that the follow-up observations during the second half of the third observing run did not meet the necessary sensitivity tomore »constrain the source properties of the potential GW candidate. Consequently, we suggest that different strategies have to be used to allow a better usage of the available telescope time. We examine different choices for follow-up surveys to optimize sky localization coverage versus observational depth to understand the likelihood of counterpart detection.« less

    High-fidelity radio interferometric data calibration that minimizes spurious spectral structure in the calibrated data is essential in astrophysical applications, such as 21 cm cosmology, which rely on knowledge of the relative spectral smoothness of distinct astrophysical emission components to extract the signal of interest. Existing approaches to radio interferometric calibration have been shown to impart spurious spectral structure to the calibrated data if the sky model used to calibrate the data is incomplete. In this paper, we introduce BayesCal: a novel solution to the sky-model incompleteness problem in interferometric calibration, designed to enable high-fidelity data calibration. The BayesCal data model supplements the a priori known component of the forward model of the sky with a statistical model for the missing and uncertain flux contribution to the data, constrained by a prior on the power in the model. We demonstrate how the parameters of this model can be marginalized out analytically, reducing the dimensionality of the parameter space to be sampled from and allowing one to sample directly from the posterior probability distribution of the calibration parameters. Additionally, we show how physically motivated priors derived from theoretical and measurement-based constraints on the spectral smoothness of the instrumental gains can be usedmore »to constrain the calibration solutions. In a companion paper, we apply this algorithm to simulated observations with a HERA-like array and demonstrate that it enables up to four orders of magnitude suppression of power in spurious spectral fluctuations relative to standard calibration approaches.

    « less
  3. Abstract Pulsar timing array (PTA) experiments are becoming increasingly sensitive to gravitational waves (GWs) in the nanohertz frequency range, where the main astrophysical sources are supermassive black hole binaries (SMBHBs), which are expected to form following galaxy mergers. Some of these individual SMBHBs may power active galactic nuclei, and thus their binary parameters could be obtained electromagnetically, which makes it possible to apply electromagnetic (EM) information to aid the search for a GW signal in PTA data. In this work, we investigate the effects of such an EM-informed search on binary detection and parameter estimation by performing mock data analyses on simulated PTA data sets. We find that by applying EM priors, the Bayes factor of some injected signals with originally marginal or sub-threshold detectability (i.e., Bayes factor ∼1) can increase by a factor of a few to an order of magnitude, and thus an EM-informed targeted search is able to find hints of a signal when an uninformed search fails to find any. Additionally, by combining EM and GW data, one can achieve an overall improvement in parameter estimation, regardless of the source’s sky location or GW frequency. We discuss the implications for the multi-messenger studies of SMBHBs withmore »PTAs.« less
  4. Abstract

    The waveform of a compact binary coalescence is predicted by general relativity. It is therefore possible to directly constrain the response of a gravitational-wave (GW) detector by analyzing a signal’s observed amplitude and phase evolution as a function of frequency. GW signals alone constrain the relative amplitude and phase between different frequencies within the same detector and between different detectors. Furthermore, if the source’s distance and inclination can be determined independently, for example from an electromagnetic (EM) counterpart, one can calibrate the absolute amplitude response of the detector network. We analyze GW170817’s ability to calibrate the LIGO/Virgo detectors, finding a relative amplitude calibration precision of approximately20% and relative phase precision of(uncertainty) between the LIGO Hanford and Livingston detectors. Incorporating additional information about the distance and inclination of the source from EM observations, the relative amplitude of the LIGO detectors can be tightened to  ∼%. Including EM observations also constrains the absolute amplitude precision to similar levels. We investigate the ability of future events to improve astronomical calibration. By simulating the cumulative uncertainties from an ensemble of detections, we find that with several hundred events with EM counterparts, or several thousandmore »events without counterparts, we reach percent-level astronomical calibration. This corresponds to  ∼5–10 years of operation at advanced LIGO and Virgo design sensitivity. It is to be emphasized that directin situmeasurements of detector calibration provide significantly higher precision than astronomical sources, and already constrain the calibration to a few percent in amplitude and a few degrees in phase. In this sense, our astronomical calibrators only corroborate existing calibration measurements. Nonetheless, it is remarkable that we are able to use an astronomical GW source to characterize properties of a terrestrial GW instrument, and astrophysical calibration may become an important corroboration of existing calibration methods, providing a completely independent constraint of potential systematics.

    « less

    The association of GRB170817A with a binary neutron star (BNS) merger has revealed that BNSs produce at least a fraction of short gamma-ray bursts (SGRBs). As gravitational wave (GW) detectors push their horizons, it is important to assess coupled electromagnetic (EM)/GW probabilities and maximize observational prospects. Here, we perform BNS population synthesis calculations with the code mobse, seeding the binaries in galaxies at three representative redshifts, $z$ = 0.01, 0.1, and 1 of the Illustris TNG50 simulation. The binaries are evolved and their locations numerically tracked in the host galactic potentials until merger. Adopting the microphysics parameters of GRB170817A, we numerically compute the broad-band light curves of jets from BNS mergers, with the afterglow brightness dependent on the local medium density at the merger site. We perform Monte Carlo simulations of the resulting EM population assuming either a random viewing angle with respect to the jet, or a jet aligned with the orbital angular momentum of the binary, which biases the viewing angle probability for GW-triggered events. We find a gamma-ray detection probability of $\sim\!2{{\rm per\ cent}},10{{\rm per\ cent}},\mathrm{and}\ 40{{\rm per\ cent}}$ for BNSs at $z$ = 1, 0.1, and 0.01, respectively, for the random case, rising to $\sim\!75{{\rm per\ cent}}$more »for the $z$ = 0.01, GW-triggered aligned case. Afterglow detection probabilities of GW-triggered BNS mergers vary in the range of $\sim \! 0.3 \!-\! 0.5{{\rm per\ cent}}$, with higher values for aligned jets, and are comparable across the high- and low-energy bands, unlike gamma-ray-triggered events (cosmological SGRBs) which are significantly brighter at higher energies. We further quantify observational biases with respect to host galaxy masses.

    « less