An official website of the United States government
ABSTRACT Observations of gravitational waves emitted by merging compact binaries have provided tantalizing hints about stellar astrophysics, cosmology, and fundamental physics. However, the physical parameters describing the systems (mass, spin, distance) used to extract these inferences about the Universe are subject to large uncertainties. The most widely used method of performing these analyses requires performing many Monte Carlo integrals to marginalize over the uncertainty in the properties of the individual binaries and the survey selection bias. These Monte Carlo integrals are subject to fundamental statistical uncertainties. Previous treatments of this statistical uncertainty have focused on ensuring that the precision of the inferred inference is unaffected; however, these works have neglected the question of whether sufficient accuracy can also be achieved. In this work, we provide a practical exploration of the impact of uncertainty in our analyses and provide a suggested framework for verifying that astrophysical inferences made with the gravitational-wave transient catalogue are accurate. Applying our framework to models used by the LIGO–Virgo–KAGRA collaboration and in the wider literature, we find that Monte Carlo uncertainty in estimating the survey selection bias is the limiting factor in our ability to probe narrow population models and this will rapidly grow more problematic as the size of the observed population increases.
more »
« less
Flexible and Accurate Evaluation of Gravitational-wave Malmquist Bias with Machine Learning
Abstract Many astronomical surveys are limited by the brightness of the sources, and gravitational-wave searches are no exception. The detectability of gravitational waves from merging binaries is affected by the mass and spin of the constituent compact objects. To perform unbiased inference on the distribution of compact binaries, it is necessary to account for this selection effect, which is known as Malmquist bias. Since systematic error from selection effects grows with the number of events, it will be increasingly important over the coming years to accurately estimate the observational selection function for gravitational-wave astronomy. We employ density estimation methods to accurately and efficiently compute the compact binary coalescence selection function. We introduce a simple pre-processing method, which significantly reduces the complexity of the required machine-learning models. We demonstrate that our method has smaller statistical errors at comparable computational cost than the method currently most widely used allowing us to probe narrower distributions of spin magnitudes. The currently used method leaves 10%–50% of the interesting black hole spin models inaccessible; our new method can probe >99% of the models and has a lower uncertainty for >80% of the models.
more »
« less
- Award ID(s):
- 1912594
- PAR ID:
- 10363838
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 927
- Issue:
- 1
- ISSN:
- 0004-637X
- Format(s):
- Medium: X Size: Article No. 76
- Size(s):
- Article No. 76
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Standardizing the definition of eccentricity is necessary for unambiguous inference of the orbital eccentricity of compact binaries from gravitational wave observations. In previous works, we proposed a definition of eccentricity for systems without spin-precession that relies solely on the gravitational waveform, is applicable to any waveform model, and has the correct Newtonian limit. In this work, we extend this definition to spin-precessing systems. This simple yet effective extension relies on first transforming the waveform from the inertial frame to the coprecessing frame, and then adopting an amplitude and a phase with reduced spin-induced effects. Our method includes a robust procedure for filtering out spin-induced modulations, which become non-negligible in the small eccentricity and large spin-precession regime. Finally, we apply our method to a set of Numerical Relativity and Effective One Body waveforms to showcase its robustness for generic eccentric spin-precessing binaries. We make our method public via Python implementation ingw_eccentricity.more » « less
-
null (Ed.)ABSTRACT Compact white dwarf (WD) binaries are important sources for space-based gravitational-wave (GW) observatories, and an increasing number of them are being identified by surveys like Extremely Low Mass (ELM) and Zwicky Transient Facility (ZTF). We study the effects of non-linear dynamical tides in such binaries. We focus on the global three-mode parametric instability and show that it has a much lower threshold energy than the local wave-breaking condition studied previously. By integrating networks of coupled modes, we calculate the tidal dissipation rate as a function of orbital period. We construct phenomenological models that match these numerical results and use them to evaluate the spin and luminosity evolution of a WD binary. While in linear theory the WD’s spin frequency can lock to the orbital frequency, we find that such a lock cannot be maintained when non-linear effects are taken into account. Instead, as the orbit decays, the spin and orbit go in and out of synchronization. Each time they go out of synchronization, there is a brief but significant dip in the tidal heating rate. While most WDs in compact binaries should have luminosities that are similar to previous traveling-wave estimates, a few per cent should be about 10 times dimmer because they reside in heating rate dips. This offers a potential explanation for the low luminosity of the CO WD in J0651. Lastly, we consider the impact of tides on the GW signal and show that the Laser Interferometer Space Antenna (LISA) and TianGO can constrain the WD’s moment of inertia to better than $$1{{\ \rm per\ cent}}$$ for centi-Hz systems.more » « less
-
Abstract LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmological distances; from the early inspirals of stellar-mass black holes that will ultimately venture into the ground-based detectors’ view to the death spiral of compact objects into massive black holes, and many sources in between. Central to realising LISA’s discovery potential are waveform models, the theoretical and phenomenological predictions of the pattern of gravitational waves that these sources emit. This White Paper is presented on behalf of the Waveform Working Group for the LISA Consortium. It provides a review of the current state of waveform models for LISA sources, and describes the significant challenges that must yet be overcome.more » « less
-
Abstract The delay time distribution of neutron star mergers provides critical insights into binary evolution processes and the merger rate evolution of compact object binaries. However, current observational constraints on this delay time distribution rely on the small sample of Galactic double neutron stars (with uncertain selection effects), a single multimessenger gravitational wave event, and indirect evidence of neutron star mergers based on r -process enrichment. We use a sample of 68 host galaxies of short gamma-ray bursts to place novel constraints on the delay time distribution and leverage this result to infer the merger rate evolution of compact object binaries containing neutron stars. We recover a power-law slope of α = − 1.83 − 0.39 + 0.35 (median and 90% credible interval) with α < −1.31 at 99% credibility, a minimum delay time of t min = 184 − 79 + 67 Myr with t min > 72 Myr at 99% credibility, and a maximum delay time constrained to t max > 7.95 Gyr at 99% credibility. We find these constraints to be broadly consistent with theoretical expectations, although our recovered power-law slope is substantially steeper than the conventional value of α = −1, and our minimum delay time is larger than the typically assumed value of 10 Myr. Pairing this cosmological probe of the fate of compact object binary systems with the Galactic population of double neutron stars will be crucial for understanding the unique selection effects governing both of these populations. In addition to probing a significantly larger redshift regime of neutron star mergers than possible with current gravitational wave detectors, complementing our results with future multimessenger gravitational wave events will also help determine if short gamma-ray bursts ubiquitously result from compact object binary mergers.more » « less
