Gravitational-wave Transients Research Articles

An evaluation of the BALROG and RoboBA algorithms for determining the position of Fermi/GBM GRBs

ABSTRACT The Fermi/GBM instrument is a vital source of detections of gamma-ray bursts and has an increasingly important role to play in understanding gravitational-wave transients. In both cases, its impact is increased by accurate positions with reliable uncertainties. We evaluate the RoboBA and BALROG algorithms for determining the position of gamma-ray bursts detected by the Fermi/GBM instrument. We construct a sample of 54 bursts with detections both by Swift/BAT and by Fermi/GBM. We then compare the positions predicted by RoboBA and BALROG with the positions measured by BAT, which we can assume to be the true position. We find that RoboBA and BALROG are similarly precise for bright bursts whose uncertainties are dominated by systematic errors, but RoboBA performs better for faint bursts whose uncertainties are dominated by statistical noise. We further find that the uncertainties in the positions predicted by RoboBA are consistent with the distribution of position errors, whereas BALROG seems to be underestimating the uncertainties by a factor of about 2. Additionally, we consider the implications of these results for the follow-up of the optical afterglows of Fermi/GBM bursts. In particular, for the DDOTI wide-field imager we conclude that a single pointing is best. Our sample would allow a similar study to be carried out for other telescopes.

Search for Gravitational-wave Transients Associated with Magnetar Bursts in Advanced LIGO and Advanced Virgo Data from the Third Observing Run

Gravitational waves are expected to be produced from neutron star oscillations associated with magnetar giant flares and short bursts. We present the results of a search for short-duration (milliseconds to seconds) and long-duration (∼100 s) transient gravitational waves from 13 magnetar short bursts observed during Advanced LIGO, Advanced Virgo, and KAGRA’s third observation run. These 13 bursts come from two magnetars, SGR 1935+2154 and Swift J1818.0−1607. We also include three other electromagnetic burst events detected by Fermi-GBM which were identified as likely coming from one or more magnetars, but they have no association with a known magnetar. No magnetar giant flares were detected during the analysis period. We find no evidence of gravitational waves associated with any of these 16 bursts. We place upper limits on the rss of the integrated incident gravitational-wave strain that reach 3.6 × 10−23 /Hz at 100 Hz for the short-duration search and 1.1 × 10−22 /Hz at 450 Hz for the long-duration search. For a ringdown signal at 1590 Hz targeted by the short-duration search the limit is set to 2.3 × 10−22 /Hz . Using the estimated distance to each magnetar, we derive upper limits on the emitted gravitational-wave energy of 1.5 × 1044 erg (1.0 × 1044 erg) for SGR 1935+2154 and 9.4 × 1043 erg (1.3 × 1044 erg) for Swift J1818.0−1607, for the short-duration (long-duration) search. Assuming isotropic emission of electromagnetic radiation of the burst fluences, we constrain the ratio of gravitational-wave energy to electromagnetic energy for bursts from SGR 1935+2154 with the available fluence information. The lowest of these ratios is 4.5 × 103.

The high energy X-ray probe (HEX-P): sensitive broadband X-ray observations of transient phenomena in the 2030s

HEX-P is a probe-class mission concept that will combine high spatial resolution X-ray imaging (<10′′ FWHM) and broad spectral coverage (0.2–80 keV) with an effective area superior to NuSTAR above 10 keV to enable revolutionary new insights into a variety of astrophysical problems, especially those related to compact objects, accretion and outflows. HEX-P will launch at a time when the sky is being routinely scanned for transient gravitational wave, electromagnetic and neutrino phenomena that will require the capabilities of a sensitive, broadband X-ray telescope for follow up studies. These include the merger of compact objects such as neutron stars and black holes, stellar explosions, and the birth of new compact objects. A response time to target of opportunity observation requests of <24 hours and a field of regard of 3π steradians will allow HEX-P to probe the accretion and ejecta from these transient phenomena through the study of relativistic outflows and reprocessed emission, provide unique capabilities for understanding jet physics, and potentially revealing the nature of the central engine.

QoQ: a Q-transform based test for gravitational wave transient events

The observation of transient gravitational waves (GWs) is hindered by the presence of transient noise, colloquially referred to as glitches. These glitches can often be misidentified as GWs by searches for unmodeled transients using the excess-power type of methods and sometimes even excite template waveforms for compact binary coalescences while using matched filter techniques. They thus create a significant background in the searches. This background is more critical in getting identified promptly and efficiently within the context of real-time searches for GW transients. Such searches are the ones that have enabled multi-messenger astrophysics with the start of the Advanced LIGO and Advanced Virgo data taking in 2015 and they will continue to enable the field for further discoveries. With this work we propose and demonstrate the use of a signal-based test that quantifies the fidelity of the time-frequency decomposition of the putative signal based on first principles on how astrophysical transients are expected to be registered in the detectors and empirically measuring the instrumental noise. It is based on the Q-transform and a measure of the occupancy of the corresponding time-frequency pixels over select time-frequency volumes; we call it ‘QoQ’. Our method shows a 40% reduction in the number of retraction of public alerts that were issued by the LIGO-Virgo-KAGRA collaborations during the third observing run with negligible loss in sensitivity. Receiver Operator Characteristic measurements suggest the method can be used in online and offline searches for transients, reducing their background significantly.

Continuous gravitational waves from trapped magnetar ejecta and the connection to glitches and antiglitches

ABSTRACT Gravitational waves from isolated sources have eluded detection so far. The upper limit of long-lasting continuous gravitational wave emission can now probe physically motivated models with the most optimistic being strongly constrained. Naturally, one might want to relax the assumption of the gravitational wave being quasi-infinite in duration, leading to the idea of transient continuous gravitational waves. In this paper, we outline how to get transient continuous waves from magnetars (or strongly magnetized neutron stars) that exhibit glitches and/or antiglitches and apply the model to magnetar SGR 1935+2154. The toy model hypothesizes that at a glitch or antiglitch, mass is ejected from the magnetar but becomes trapped on its outward journey through the magnetosphere. Depending on the height of the trapped ejecta and the magnetic inclination angle, we are able to reproduce both glitches and antiglitches from simple angular momentum arguments. The trapped ejecta causes the magnetar to precess leading to gravitational wave emission at once and twice the magnetar’s spin frequency, for a duration equal to however long the ejecta is trapped for. We find that the gravitational waves are more detectable when the magnetar is: closer, rotating faster, or has larger glitches/antiglitches. The detectability also improves when the ejecta height and magnetic inclination angle have values near their critical values, though this requires more mass to be ejected to remain consistent with the observed glitch/antiglitch. We find it unlikely that gravitational waves will be detected from SGR 1935+2154 when using the trapped ejecta model.

Deep Multimessenger Search for Compact Binary Mergers in LIGO, Virgo, and Fermi/GBM Data from 2016–2017

GW170817–GRB 170817A provided the first observation of gravitational waves from a neutron star merger with associated transient counterparts across the entire electromagnetic spectrum. This discovery demonstrated the long-hypothesized association between short gamma-ray bursts and neutron star mergers. More joint detections are needed to explore the relation between the parameters inferred from the gravitational wave and the properties of the gamma-ray burst signal. We developed a joint multimessenger analysis of LIGO, Virgo, and Fermi/GBM data designed for detecting weak gravitational-wave transients associated with weak gamma-ray bursts. As such, it does not start from confident (GWTC-1) events only. Instead, we take the full list of existing compact binary coalescence triggers generated with the PyCBC pipeline from the second Gravitational-Wave Observing Run (O2), and reanalyze the entire set of public Fermi/GBM data covering this observing run to generate a corresponding set of gamma-ray burst candidate triggers. We then search for coincidences between the gravitational-wave and gamma-ray burst triggers without requiring a confident detection in any channel. The candidate coincidences are ranked according to a statistic combining each candidate’s strength in gravitational-wave and gamma-ray data, their time proximity, and the overlap of their sky localization. The ranking is then converted to a false alarm rate using time shifts between the gravitational-wave and gamma-ray burst triggers. We present the results using O2 triggers, which allowed us to check the validity of our method against GW170817–GRB 170817A. We also discuss the different configurations tested to maximize the significance of the joint detection.

Search for Gravitational Waves Associated with Fast Radio Bursts Detected by CHIME/FRB during the LIGO–Virgo Observing Run O3a

We search for gravitational-wave (GW) transients associated with fast radio bursts (FRBs) detected by the Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst Project, during the first part of the third observing run of Advanced LIGO and Advanced Virgo (2019 April 1 15:00 UTC–2019 October 1 15:00 UTC). Triggers from 22 FRBs were analyzed with a search that targets both binary neutron star (BNS) and neutron star–black hole (NSBH) mergers. A targeted search for generic GW transients was conducted on 40 FRBs. We find no significant evidence for a GW association in either search. Given the large uncertainties in the distances of our FRB sample, we are unable to exclude the possibility of a GW association. Assessing the volumetric event rates of both FRB and binary mergers, an association is limited to 15% of the FRB population for BNS mergers or 1% for NSBH mergers. We report 90% confidence lower bounds on the distance to each FRB for a range of GW progenitor models and set upper limits on the energy emitted through GWs for a range of emission scenarios. We find values of order 1051–1057 erg for models with central GW frequencies in the range 70–3560 Hz. At the sensitivity of this search, we find these limits to be above the predicted GW emissions for the models considered. We also find no significant coincident detection of GWs with the repeater, FRB 20200120E, which is the closest known extragalactic FRB.

Convolutional neural networks for the classification of glitches in gravitational-wave data streams

We investigate the use of convolutional neural networks (including the modern ConvNeXt network family) to classify transient noise signals (i.e. glitches) and gravitational waves (GWs) in data from the Advanced LIGO detectors. First, we use models with a supervised learning approach, both trained from scratch using the Gravity Spy dataset and employing transfer learning by fine-tuning pre-trained models in this dataset. Second, we also explore a self-supervised approach, pre-training models with automatically generated pseudo-labels. Our findings are very close to existing results for the same dataset, reaching values for the F1 score of 97.18% (94.15%) for the best supervised (self-supervised) model. We further test the models using actual GW signals from LIGO-Virgo’s O3 run. Although trained using data from previous runs (O1 and O2), the models show good performance, in particular when using transfer learning. We find that transfer learning improves the scores without the need for any training on real signals apart from the less than 50 chirp examples from hardware injections present in the Gravity Spy dataset. This motivates the use of transfer learning not only for glitch classification but also for signal classification.

Semianalytic sensitivity estimates for catalogs of gravitational-wave transients

Semianalytic sensitivity estimates for catalogs of gravitational-wave transients

An autoencoder neural network integrated into gravitational-wave burst searches to improve the rejection of noise transients

The gravitational-wave (GW) detector data are affected by short-lived instrumental or terrestrial transients, called ‘glitches’, which can simulate GW signals. Mitigation of glitches is particularly difficult for algorithms which target generic sources of short-duration GW transients (GWT), and do not rely on GW waveform models to distinguish astrophysical signals from noise, such as coherent WaveBurst (cWB). This work is part of the long-term effort to mitigate transient noises in cWB, which led to the introduction of specific estimators, and a machine-learning based signal-noise classification algorithm. Here, we propose an autoencoder neural network, integrated into cWB, that learns transient noises morphologies from GW time-series. We test its performance on the glitch family known as ‘blip’. The resulting sensitivity to generic GWT and binary black hole mergers significantly improves when tested on LIGO detectors data from the last observation period (O3b). At false alarm rate of one event per 50 years the sensitivity volume increases up to 30% for signal morphologies similar to blip glitches. In perspective, this tool can adapt to classify different transient noise classes that may affect future observing runs, enhancing GWT searches.

Accelerating tests of general relativity with gravitational-wave signals using hybrid sampling

The Advanced LIGO/Virgo interferometers have observed $\ensuremath{\sim}100$ gravitational-wave transients enabling new questions to be answered about relativity, astrophysics, and cosmology. However, many of our current procedures for computing these constraints will not scale well with the increased size of future transient catalogs. We introduce a novel hybrid sampling method in order to more efficiently perform parameterized tests of general relativity with gravitational-wave signals. Applying our method to the binary black hole merger GW150914 and simulated signals we find that our method is approximately an order of magnitude more efficient than the current method with conservative settings for our hybrid analysis. While we have focused on the specific problem of measuring potential deviations from relativity, our method is of much wider applicability to any problem that can be decomposed into a simple and more complex model(s).

Search for gravitational-wave bursts in the third Advanced LIGO-Virgo run with coherent WaveBurst enhanced by machine learning

This paper presents a search for generic short-duration gravitational-wave (GW) transients (or GW bursts) in the data from the third observing run of Advanced LIGO and Advanced Virgo. We use a coherent WaveBurst (cWB) pipeline enhanced with a decision-tree classification algorithm for more efficient separation of GW signals from noise transients. The machine-learning (ML) algorithm is trained on a representative set of noise events and a set of simulated stochastic signals that are not correlated with any known signal model. This training procedure preserves the model-independent nature of the search. We demonstrate that the ML-enhanced cWB pipeline can detect GW signals at a larger distance than previous model-independent searches. The sensitivity improvements are achieved across the broad spectrum of simulated signals, with the goal of testing the robustness of this model-agnostic search. At a false-alarm rate of one event per century, the detectable signal amplitudes are reduced up to almost an order of magnitude, most notably for the single-cycle signal morphologies. This ML-enhanced pipeline also improves the detection efficiency of compact binary mergers in a wide range of masses, from stellar mass to intermediate-mass black holes, both with circular and elliptical orbits. After excluding previously detected compact binaries, no new gravitational-wave signals are observed for the twofold Hanford-Livingston and the threefold Hanford-Livingston-Virgo detector networks. With the improved sensitivity of the all-sky search, we obtain the most stringent constraints on the isotropic emission of gravitational-wave energy from short-duration burst sources.

Stripped-envelope stars in different metallicity environments

Stripped-envelope stars can be observed as Wolf-Rayet (WR) stars or as less luminous hydrogen-poor stars with low mass-loss rates and transparent winds. Both types are potential progenitors of Type I core-collapse supernovae (SNe). We used grids of core-collapse models obtained from single helium stars at different metallicities to study the effects of metallicity on the transients and remnants these stars produce. We characterised the surface and core properties of our core-collapse models and investigated their ‘explodability’ using three criteria. In the cases where explosions are predicted, we estimated the ejecta mass, explosion energy, nickel mass, and neutron star (NS) mass. Otherwise, we predicted the mass of the resulting black hole (BH). We constructed a simplified population model and find that the properties of SNe and compact objects depend strongly on metallicity. The ejecta masses and explosion energies for Type Ic SNe are best reproduced by models with Z = 0.04 that exhibit strong winds during core helium burning. This implies that either their mass-loss rates are underestimated or that Type Ic SN progenitors experience mass loss through other mechanisms before exploding. The distributions of ejecta masses, explosion energies, and nickel mass for Type Ib SNe are not well reproduced by progenitor models with WR mass loss, but are better reproduced if we assume no mass loss in progenitors with luminosities below the minimum WR star luminosity. We find that Type Ic SNe become more common as metallicity increases, and that the vast majority of progenitors of Type Ib SNe must be transparent-wind stripped-envelope stars. We find that several models with pre-collapse CO masses of up to ∼30 M⊙ may form ∼3 M⊙ BHs in fallback SNe. This may have important consequences for our understanding of SNe, binary BH and NS systems, X-ray binary systems, and gravitational wave transients.

The Influence of Disk Composition on the Evolution of Stars in the Disks of Active Galactic Nuclei

Disks of gas accreting onto supermassive black holes, powering active galactic nuclei (AGN), can capture stars from nuclear star clusters or form stars in situ via gravitational instability. The density and thermal conditions of these disks can result in rapid accretion onto embedded stars, dramatically altering their evolution in comparison to stars in the interstellar medium. Theoretical models predict that, when subjected to sufficiently rapid accretion, fresh gas replenishes hydrogen in the cores of these stars as quickly as it is burned into helium, reaching a quasi-steady state. Such massive, long-lived (“immortal”) stars may be capable of dramatically enriching AGN disks with helium, and would increase the helium abundance in AGN broad-line regions relative to that in the corresponding narrow-line regions and hosts. We investigate how the helium abundance of AGN disks alters the evolution of stars embedded therein. We find, in agreement with analytical arguments, that stars at a given mass are more luminous at higher helium mass fractions, and so undergo more radiation-driven mass loss. We further find that embedded stars tend to be less massive in disks with higher helium mass fractions, and that immortal stars are less common in such disks. Thus, disk composition can alter the rates of electromagnetic and gravitational wave transients as well as further chemical enrichment by embedded stars.

The Arches cluster revisited – IV. Observational constraints on the binary properties of very massive stars

ABSTRACT Serving as the progenitors of electromagnetic and gravitational-wave transients, massive stars have received renewed interest in recent years. However, many aspects of their birth and evolution remain opaque, particularly in the context of binary interactions. The centre of our Galaxy hosts a rich cohort of very massive stars, which appear to play a prominent role in the ecology of the region. In this paper, we investigate the binary properties of the Arches cluster, which is thought to host a large number of very massive stars. A combination of multi-epoch near-IR spectroscopy and photometry was utilized to identify binaries. 13 from 36 cluster members meet our criteria to be classed as RV variable. Combining the spectroscopic data with archival radio and X-ray observations – to detect colliding wind systems – provides a lower limit to the binary fraction of ∼43 per cent; increasing to ≳50 per cent for the O-type hypergiants and WNLha. Dynamical and evolutionary masses reveal the primaries to be uniformly massive ($\gtrsim 50\, \mathrm{M}_{\odot }$). Where available, orbital analysis reveals a number of short period, highly eccentric binaries, which appear to be pre-interaction systems. Such systems are X-ray luminous, with 80 per cent above an empirical bound of (Lx/Lbol) ∼ 10−7 and their orbital configurations suggest formation and evolution via a single-star channel; however, we cannot exclude a binary formation channel for a subset. Qualitative comparison to surveys of lower mass OB-type stars confirms that the trend to an extreme binary fraction (≥60 per cent) extends to the most massive stars currently forming in the local Universe.

Correlated 1–1000 Hz magnetic field fluctuations from lightning over Earth-scale distances and their impact on gravitational wave searches

We report Earth-scale distance magnetic correlations from lightning strokes in the frequency range 1-1000 Hz at several distances ranging from 1100 to 9000 km. Noise sources which are correlated on Earth-scale distances can affect future searches for gravitational-wave signals with ground-based gravitational-wave interferometric detectors. We consider the impact of correlations from magnetic field fluctuations on gravitational-wave searches due to Schumann resonances ($<$50 Hz) as well as higher frequencies ($>$100 Hz). We demonstrate that individual lightning strokes are a likely source for the observed correlations in the magnetic field fluctuations at gravitational-wave observatories and discuss some of their characteristics. Furthermore, we predict their impact on searches for an isotropic gravitational-wave background, as well as for searches looking for short-duration transient gravitational waves, both unmodeled signals (bursts) as well as modeled signals (compact binary coalescence). Whereas the recent third observing run by LIGO and Virgo was free of an impact from correlated magnetic field fluctuations, future runs could be affected. For example, at current magnetic coupling levels, neutron star inspirals in third generation detectors are likely to be contaminated by multiple correlated lightning glitches. We suggest that future detector design should consider reducing lightning coupling by, for example, reducing the lightning-induced beam tube currents that pass through sensitive magnetic coupling regions in current detectors. We also suggest that the diurnal and seasonal variation in lightning activity may be useful in discriminating between detector correlations that are produced by gravitational waves and those produced by lightning.

Convolutional neural network to distinguish glitches from minute-long gravitational wave transients

Gravitational wave bursts are transient signals distinct from compact binary mergers that arise from a wide variety of astrophysical phenomena. Because most of these phenomena are poorly modeled, the use of traditional search methods such as matched filtering is excluded. Bursts include short ($<$10 seconds) and long (from 10 to a few hundreds of seconds) duration signals for which the detection is constrained by environmental and instrumental transient noises called glitches. Glitches contaminate burst searches, reducing the amount of useful data and limiting the sensitivity of current algorithms. It is therefore of primordial importance to locate and distinguish them from potential burst signals. In this paper, we propose to train a convolutional neural network to detect glitches in the time-frequency space of the cross-correlated LIGO noise. We show that our network is retrieving more than 95$\%$ of the glitches while being trained only on a subset of the existing glitch classes highlighting the sensitivity of the network to completely new glitch classes.

MWA rapid follow-up of gravitational wave transients: Prospects for detecting prompt radio counterparts

Abstract We present and evaluate the prospects for detecting coherent radio counterparts to gravitational wave (GW) events using Murchison Widefield Array (MWA) triggered observations. The MWA rapid-response system, combined with its buffering mode ( $\sim$ 4 min negative latency), enables us to catch any radio signals produced from seconds prior to hours after a binary neutron star (BNS) merger. The large field of view of the MWA ( $\sim$ $1\,000\,\textrm{deg}^2$ at 120 MHz) and its location under the high sensitivity sky region of the LIGO-Virgo-KAGRA (LVK) detector network, forecast a high chance of being on-target for a GW event. We consider three observing configurations for the MWA to follow up GW BNS merger events, including a single dipole per tile, the full array, and four sub-arrays. We then perform a population synthesis of BNS systems to predict the radio detectable fraction of GW events using these configurations. We find that the configuration with four sub-arrays is the best compromise between sky coverage and sensitivity as it is capable of placing meaningful constraints on the radio emission from 12.6% of GW BNS detections. Based on the timescales of four BNS merger coherent radio emission models, we propose an observing strategy that involves triggering the buffering mode to target coherent signals emitted prior to, during or shortly following the merger, which is then followed by continued recording for up to three hours to target later time post-merger emission. We expect MWA to trigger on $\sim$ $5-22$ BNS merger events during the LVK O4 observing run, which could potentially result in two detections of predicted coherent emission.

Prospects for detecting transient quasi-monochromatic gravitational waves from glitching pulsars with current and future detectors

ABSTRACT Pulsars are rotating neutron stars that emit periodic electromagnetic radiation. While pulsars generally slow down as they lose energy, some also experience glitches: spontaneous increases of their rotational frequency. According to several models, these glitches can also lead to the emission of long-duration transient gravitational waves (GWs). We present detection prospects for such signals by comparing indirect energy upper limits on GW strain for known glitches with the sensitivity of current and future ground-based GW detectors. We first consider the optimistic case of generic constraints based on the glitch size and find that realistic matched-filter searches in the fourth LIGO–Virgo–KAGRA observing run (O4) could make a detection, or set constraints below these indirect upper limits, for equivalents of 36 out of 726 previously observed glitches, and 74 in the O5 run. With the third-generation Einstein Telescope or Cosmic Explorer, 35–40 per cent of glitches would be accessible. When specializing to a scenario where transient mountains produce the post-glitch GW emission, following the Yim &amp; Jones model, the indirect upper limits are stricter. Out of the smaller set of 119 glitches with measured healing parameter, as needed for predictions under that model, only 6 glitches would have been within reach for O4 and 14 for O5, with a similar percentage as before with third-generation detectors. We also discuss how this model matches the observed glitch population.

Prospects for detecting and localizing short-duration transient gravitational waves from glitching neutron stars without electromagnetic counterparts

Neutron stars are known to show accelerated spin-up of their rotational frequency called a glitch. Highly magnetized rotating neutron stars (pulsars) are frequently observed by radio telescopes (and in other frequencies), where the glitch is observed as irregular arrival times of pulses which are otherwise very regular. A glitch in an isolated neutron star can excite the fundamental (f)-mode oscillations which can lead to gravitational wave generation. Electromagnetic observations of pulsars (and hence pulsar glitches) require the pulsar to be oriented so that the jet is pointed toward the detector, but this is not a requirement for gravitational wave emission which is more isotropic and not jetlike. Hence, gravitational wave observations have the potential to uncover nearby neutron stars where the jet is not pointed towards the Earth. In this work, we study the prospects of finding glitching neutron stars using a generic all-sky search for short-duration gravitational wave transients. The analysis covers the high-frequency range from $1-4$ kHz of LIGO-Virgo detectors for signals up to a few seconds. We set upper limits for the third observing run of the LIGO-Virgo detectors and present the prospects for upcoming observing runs of LIGO, Virgo, KAGRA, and LIGO India. We find the detectable glitch size will be around $10^{-5}$ Hz for the fifth observing run for pulsars with spin frequencies and distances comparable to the Vela pulsar. We also present the prospects of localizing the direction in the sky of these sources with gravitational waves alone, which can facilitate electromagnetic follow-up. We find that for the five detector configuration, the localization capability for a glitch size of $10^{-5}$ Hz is around $132\,\mathrm{deg}^{2}$ at $1\sigma$ confidence for $50\%$ of events with distance and spin frequency as that of Vela.