-
Real-time gravitational-wave inference for binary neutron stars using machine learning
Authors:
Maximilian Dax,
Stephen R. Green,
Jonathan Gair,
Nihar Gupte,
Michael Pürrer,
Vivien Raymond,
Jonas Wildberger,
Jakob H. Macke,
Alessandra Buonanno,
Bernhard Schölkopf
Abstract:
Mergers of binary neutron stars (BNSs) emit signals in both the gravitational-wave (GW) and electromagnetic (EM) spectra. Famously, the 2017 multi-messenger observation of GW170817 led to scientific discoveries across cosmology, nuclear physics, and gravity. Central to these results were the sky localization and distance obtained from GW data, which, in the case of GW170817, helped to identify the…
▽ More
Mergers of binary neutron stars (BNSs) emit signals in both the gravitational-wave (GW) and electromagnetic (EM) spectra. Famously, the 2017 multi-messenger observation of GW170817 led to scientific discoveries across cosmology, nuclear physics, and gravity. Central to these results were the sky localization and distance obtained from GW data, which, in the case of GW170817, helped to identify the associated EM transient, AT 2017gfo, 11 hours after the GW signal. Fast analysis of GW data is critical for directing time-sensitive EM observations; however, due to challenges arising from the length and complexity of signals, it is often necessary to make approximations that sacrifice accuracy. Here, we develop a machine learning approach that performs complete BNS inference in just one second without making any such approximations. This is enabled by a new method for explicit integration of physical domain knowledge into neural networks. Our approach enhances multi-messenger observations by providing (i) accurate localization even before the merger; (ii) improved localization precision by $\sim30\%$ compared to approximate low-latency methods; and (iii) detailed information on luminosity distance, inclination, and masses, which can be used to prioritize expensive telescope time. Additionally, the flexibility and reduced cost of our method open new opportunities for equation-of-state and waveform systematics studies. Finally, we demonstrate that our method scales to extremely long signals, up to an hour in length, thus serving as a blueprint for data analysis for next-generation ground- and space-based detectors.
△ Less
Submitted 12 July, 2024;
originally announced July 2024.
-
Evidence for eccentricity in the population of binary black holes observed by LIGO-Virgo-KAGRA
Authors:
Nihar Gupte,
Antoni Ramos-Buades,
Alessandra Buonanno,
Jonathan Gair,
M. Coleman Miller,
Maximilian Dax,
Stephen R. Green,
Michael Pürrer,
Jonas Wildberger,
Jakob Macke,
Bernhard Schölkopf
Abstract:
Binary black holes (BBHs) in eccentric orbits produce distinct modulations the emitted gravitational waves (GWs). The measurement of orbital eccentricity can provide robust evidence for dynamical binary formation channels. We analyze 57 GW events from the first, second and third observing runs of the LIGO-Virgo-KAGRA (LVK) Collaboration using a multipolar aligned-spin inspiral-merger-ringdown wave…
▽ More
Binary black holes (BBHs) in eccentric orbits produce distinct modulations the emitted gravitational waves (GWs). The measurement of orbital eccentricity can provide robust evidence for dynamical binary formation channels. We analyze 57 GW events from the first, second and third observing runs of the LIGO-Virgo-KAGRA (LVK) Collaboration using a multipolar aligned-spin inspiral-merger-ringdown waveform model with two eccentric parameters: eccentricity and relativistic anomaly. This is made computationally feasible with the machine-learning code DINGO which accelerates inference by 2-3 orders of magnitude compared to traditional inference. First, we find eccentric aligned-spin versus quasi-circular aligned-spin $\log_{10}$ Bayes factors of 1.84 to 4.75 (depending on the glitch mitigation) for GW200129, 3.0 for GW190701 and 1.77 for GW200208_22. We measure $e_{\text{gw}, 10Hz}$ to be $0.27_{-0.12}^{+0.10}$ to $0.17_{-0.13}^{+0.14}$ for GW200129, $0.35_{-0.11}^{+0.32}$ for GW190701 and $0.35_{-0.21}^{+0.18}$ for GW200208_22. Second, we find $\log_{10}$ Bayes factors between the eccentric aligned-spin versus quasi-circular precessing-spin hypothesis between 1.43 and 4.92 for GW200129, 2.61 for GW190701 and 1.23 for GW200208_22. Third, our analysis does not show evidence for eccentricity in GW190521, which has an eccentric aligned-spin against quasi-circular aligned-spin $\log_{10}$ Bayes factor of 0.04. Fourth, we estimate that if we neglect the spin-precession and use an astrophysical prior, the probability of one out of the 57 events being eccentric is greater than 99.5% or $(100 - 8.4 \times 10^{-4})$% (depending on the glitch mitigation). Fifth, we study the impact on parameter estimation when neglecting either eccentricity or higher modes in eccentric models. These results underscore the importance of including eccentric parameters in the characterization of BBHs for GW detectors.
△ Less
Submitted 22 April, 2024;
originally announced April 2024.
-
Ultralight vector dark matter search using data from the KAGRA O3GK run
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
the KAGRA Collaboration,
A. G. Abac,
R. Abbott,
H. Abe,
I. Abouelfettouh,
F. Acernese,
K. Ackley,
C. Adamcewicz,
S. Adhicary,
N. Adhikari,
R. X. Adhikari,
V. K. Adkins,
V. B. Adya,
C. Affeldt,
D. Agarwal,
M. Agathos,
O. D. Aguiar,
I. Aguilar,
L. Aiello,
A. Ain,
P. Ajith,
T. Akutsu,
S. Albanesi
, et al. (1778 additional authors not shown)
Abstract:
Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we prese…
▽ More
Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we present the result of a search for $U(1)_{B-L}$ gauge boson DM using the KAGRA data from auxiliary length channels during the first joint observation run together with GEO600. By applying our search pipeline, which takes into account the stochastic nature of ultralight DM, upper bounds on the coupling strength between the $U(1)_{B-L}$ gauge boson and ordinary matter are obtained for a range of DM masses. While our constraints are less stringent than those derived from previous experiments, this study demonstrates the applicability of our method to the lower-mass vector DM search, which is made difficult in this measurement by the short observation time compared to the auto-correlation time scale of DM.
△ Less
Submitted 5 March, 2024;
originally announced March 2024.
-
Search for Eccentric Black Hole Coalescences during the Third Observing Run of LIGO and Virgo
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
the KAGRA Collaboration,
A. G. Abac,
R. Abbott,
H. Abe,
F. Acernese,
K. Ackley,
C. Adamcewicz,
S. Adhicary,
N. Adhikari,
R. X. Adhikari,
V. K. Adkins,
V. B. Adya,
C. Affeldt,
D. Agarwal,
M. Agathos,
O. D. Aguiar,
I. Aguilar,
L. Aiello,
A. Ain,
P. Ajith,
T. Akutsu,
S. Albanesi,
R. A. Alfaidi
, et al. (1750 additional authors not shown)
Abstract:
Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effect…
▽ More
Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that were already identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total mass $M>70$ $M_\odot$) binaries covering eccentricities up to 0.3 at 15 Hz orbital frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place an upper limit for the merger rate density of high-mass binaries with eccentricities $0 < e \leq 0.3$ at $0.33$ Gpc$^{-3}$ yr$^{-1}$ at 90\% confidence level.
△ Less
Submitted 7 August, 2023;
originally announced August 2023.
-
EinsteinPy: A Community Python Package for General Relativity
Authors:
Shreyas Bapat,
Ritwik Saha,
Bhavya Bhatt,
Shilpi Jain,
Akshita Jain,
Sofía Ortín Vela,
Priyanshu Khandelwal,
Jyotirmaya Shivottam,
Jialin Ma,
Gim Seng Ng,
Pratyush Kerhalkar,
Hrishikesh Sudam Sarode,
Rishi Sharma,
Manvi Gupta,
Divya Gupta,
Tushar Tyagi,
Tanmay Rustagi,
Varun Singh,
Saurabh Bansal,
Naman Tayal,
Abhijeet Manhas,
Raphael Reyna,
Gaurav Kumar,
Govind Dixit,
Ratin Kumar
, et al. (24 additional authors not shown)
Abstract:
This paper presents EinsteinPy (version 0.3), a community-developed Python package for gravitational and relativistic astrophysics. Python is a free, easy to use a high-level programming language which has seen a huge expansion in the number of its users and developers in recent years. Specifically, a lot of recent studies show that the use of Python in Astrophysics and general physics has increas…
▽ More
This paper presents EinsteinPy (version 0.3), a community-developed Python package for gravitational and relativistic astrophysics. Python is a free, easy to use a high-level programming language which has seen a huge expansion in the number of its users and developers in recent years. Specifically, a lot of recent studies show that the use of Python in Astrophysics and general physics has increased exponentially. We aim to provide a very high level of abstraction, an easy to use interface and pleasing user experience. EinsteinPy is developed keeping in mind the state of a theoretical gravitational physicist with little or no background in computer programming and trying to work in the field of numerical relativity or trying to use simulations in their research. Currently, EinsteinPy supports simulation of time-like and null geodesics and calculates trajectories in different background geometries some of which are Schwarzschild, Kerr, and KerrNewmann along with coordinate inter-conversion pipeline. It has a partially developed pipeline for plotting and visualization with dependencies on libraries like Plotly, matplotlib, etc. One of the unique features of EinsteinPy is a sufficiently developed symbolic tensor manipulation utilities which are a great tool in itself for teaching yourself tensor algebra which for many beginner students can be overwhelmingly tricky. EinsteinPy also provides few utility functions for hypersurface embedding of Schwarzschild spacetime which further will be extended to model gravitational lensing simulation.
△ Less
Submitted 22 May, 2020;
originally announced May 2020.
-
Optimal Gravitational-wave Follow-up Tiling Strategies Using a Genetic Algorithm
Authors:
Nihar Gupte,
Imre Bartos
Abstract:
The identification of electromagnetic emission from gravitational-wave sources typically requires multiple follow-up observations due to the limited fields-of-view of follow-up observatories compared to the poorly localized direction of gravitational waves. Gravitational-wave localization regions are typically covered with multiple telescope pointings using a "honeycomb" structure, which is optima…
▽ More
The identification of electromagnetic emission from gravitational-wave sources typically requires multiple follow-up observations due to the limited fields-of-view of follow-up observatories compared to the poorly localized direction of gravitational waves. Gravitational-wave localization regions are typically covered with multiple telescope pointings using a "honeycomb" structure, which is optimal only on an infinite, flat surface. Here we present a machine-learning algorithm which uses genetic algorithms along with Broyden-Fletcher-Goldfarb-Shanno (BFGS) optimization to find an optimal configuration of tiles to cover the gravitational-wave sky localization area on a spherical surface.
△ Less
Submitted 27 May, 2020; v1 submitted 10 March, 2020;
originally announced March 2020.
-
Gravitational-wave follow-up with CTA after the detection of GRBs in the TeV energy domain
Authors:
Imre Bartos,
K. Rainer Corley,
Nihar Gupte,
Neil Ash,
Zsuzsa. Márka,
Szabolcs Márka
Abstract:
The recent discovery of TeV emission from gamma-ray bursts (GRBs) by the MAGIC and H.E.S.S. Cherenkov telescopes confirmed that emission from these transients can extend to very high energies. The TeV energy domain reaches the most sensitive band of the Cherenkov Telescope Array (CTA). This newly anticipated, improved sensitivity will enhance the prospects of gravitational-wave follow-up observati…
▽ More
The recent discovery of TeV emission from gamma-ray bursts (GRBs) by the MAGIC and H.E.S.S. Cherenkov telescopes confirmed that emission from these transients can extend to very high energies. The TeV energy domain reaches the most sensitive band of the Cherenkov Telescope Array (CTA). This newly anticipated, improved sensitivity will enhance the prospects of gravitational-wave follow-up observations by CTA to probe particle acceleration and high-energy emission from binary black hole and neutron star mergers, and stellar core-collapse events. Here we discuss the implications of TeV emission on the most promising strategies of choice for the gravitational-wave follow-up effort for CTA and Cherenkov telescopes more broadly. We find that TeV emission (i) may allow more than an hour of delay between the gravitational-wave event and the start of CTA observations; (ii) enables the use of CTA's small size telescopes that have the largest fields of view. We characterize the number of pointings needed to find a counterpart. (iii) We compute the annual follow-up time requirements and find that prioritization will be needed. (iv) Even a few telescopes could detect sufficiently nearby counterparts, raising the possibility of adding a handful of small-size or medium-size telescopes to the network at diverse geographic locations taking into account the positions of CTA and the LIGO-Virgo-KAGRA network. (v) The continued operation of VERITAS/H.E.S.S./MAGIC would be a useful compliment to CTA's follow-up capabilities by increasing the sky area that can be rapidly covered, especially for directions above and 'below' the United States in which the present network of gravitational-wave detectors is more sensitive.
△ Less
Submitted 8 October, 2019; v1 submitted 26 August, 2019;
originally announced August 2019.
-
How would a nearby kilonova look on camera?
Authors:
Nihar Gupte,
Imre Bartos
Abstract:
Kilonovae are cosmic optical flashes produced in the aftermath of the merger of two neutron stars. While the typical radiant flux of a kilonova can be as high as $10^{34}$\,W, they typically occur at cosmological distances, requiring meter-class or larger telescopes for their observation. Here we explore how a kilonova would look like from Earth if it occurred in the Solar System's backyard, 1000…
▽ More
Kilonovae are cosmic optical flashes produced in the aftermath of the merger of two neutron stars. While the typical radiant flux of a kilonova can be as high as $10^{34}$\,W, they typically occur at cosmological distances, requiring meter-class or larger telescopes for their observation. Here we explore how a kilonova would look like from Earth if it occurred in the Solar System's backyard, 1000 light years from Earth. This is a small distance on cosmological scales, with only one neutron-star merger expected to occur this close every 100,000,000 years. While humans will likely only see kilonovae at much greater distances, showing how such a nearby event would look on camera can help visualize these events, and demonstrate their unique spectral evolution.
△ Less
Submitted 4 May, 2019;
originally announced May 2019.
-
Observational consequences of structured jets from neutron star mergers in the local Universe
Authors:
Nihar Gupte,
Imre Bartos
Abstract:
The recent multi-messenger discovery of binary neutron star (BNS) merger GW170817 showed that $γ$-ray emission in short GRBs is wider than the central energetic narrow cone, and weakly expands out to tens of degrees. Here we explore some of the observational consequences of this structured emission, taking the reconstructed angular emission profile of gamma-ray burst GRB170817A to be typical. We c…
▽ More
The recent multi-messenger discovery of binary neutron star (BNS) merger GW170817 showed that $γ$-ray emission in short GRBs is wider than the central energetic narrow cone, and weakly expands out to tens of degrees. Here we explore some of the observational consequences of this structured emission, taking the reconstructed angular emission profile of gamma-ray burst GRB170817A to be typical. We calculate the expected fraction of gravitational-wave observations from BNS mergers that will have an observed GRB counterpart to be $\sim30\%$, implying that joint gravitational-wave and GRB observations will be common. Further, we find that $\sim10\%$ of observed short GRBs occur within 200 Mpc. Finally, we estimate a BNS merger rate of $\sim500$ Gpc$^{-3}$yr$^{-1}$.
△ Less
Submitted 19 August, 2018;
originally announced August 2018.