Gravitational Wave Events Research Articles

Identifying glitches near gravitational-wave signals from compact binary coalescences using the Q-transform

We present a computational method to identify glitches in gravitational-wave data that occur nearby gravitational-wave signals from compact binary coalescences. The Q-transform, an established tool in LIGO-Virgo-KAGRA data analysis, computes the probability of any excess in the data surrounding a signal against the assumption of a Gaussian noise background, flagging any significant glitches. Subsequently, we perform validation tests on this computational method to ensure self-consistency in colored Gaussian noise, as well as data that contains a gravitational-wave event after subtracting the signal using the best-fit template. Finally, a comparison of our glitch identification results from real events in LIGO-Virgo’s third observing run against the list of events which required glitch mitigation shows that this tool will be useful in providing precise information about data quality to improve astrophysical analyses of these events.

Cocoon breakout and escape from the ejecta of neutron star mergers

ABSTRACT The cocoon is an inevitable product of a jet propagating through ambient matter, and takes a fair fraction of the jet energy. In short gamma-ray bursts (sGRBs), the ambient matter is the ejecta from the merger of neutron stars, expanding with a high velocity ∼0.2c, in contrast to the static stellar envelope in collapsars. Using 2D relativistic hydrodynamic simulations with the ejecta density profile as ρ ∝ r−2, we find that the expansion makes a big difference; only 0.5–5 per cent of the cocoon mass escapes from (faster than) the ejecta, with an opening angle 20°–30°, while it is $\sim 100{{\ \rm per\ cent}}$ and spherical in collapsars. We also analytically obtain the shares of mass and energies for the escaped and trapped cocoons. Because the mass of the escaped cocoon is small and the trapped cocoon is concealed by the ejecta and the escaped cocoon, we suggest that it is unlikely that cooling emission from the sGRB-jet heated cocoon was observed as a counterpart to the gravitational wave event GW170817.

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.

A Standard Siren Measurement of the Hubble Constant Using Gravitational-wave Events from the First Three LIGO/Virgo Observing Runs and the DESI Legacy Survey

We present a new constraint on the Hubble constant H 0 using a sample of well-localized gravitational-wave (GW) events detected during the first three LIGO/Virgo observing runs as dark standard sirens. In the case of dark standard sirens, a unique host galaxy is not identified, and the redshift information comes from the distribution of potential host galaxies. From the third LIGO/Virgo observing run detections, we add the asymmetric-mass binary black hole GW190412 and the high-confidence GW candidates S191204r, S200129m, and S200311bg to the sample of dark standard sirens analyzed in Palmese et al. Our sample contains the top 20% (based on localization) GW events and candidates to date with significant coverage by the Dark Energy Spectroscopic Instrument Legacy Survey. We combine the H 0 posterior for eight dark siren events, finding (68% highest density interval) for a prior in H 0 uniform between [20, 140] km s−1 Mpc−1. This result shows that a combination of eight well-localized dark sirens combined with an appropriate galaxy catalog is able to provide an H 0 constraint that is competitive (∼20% versus 18% precision) with a single bright standard siren analysis (i.e., assuming the electromagnetic counterpart) using GW170817. When combining the posterior with that from GW170817, we obtain . This result is broadly consistent with recent H 0 estimates from both the cosmic microwave background and supernovae.

A physics-based neural network reconstruction of the dense matter equation of state from neutron star observables

We introduce a novel technique that utilizes a physics-driven deep learning method to reconstruct the dense matter equation of state from neutron star observables, particularly the masses and radii. The proposed framework involves two neural networks: one to optimize the EoS using Automatic Differentiation in the unsupervised learning scheme; and a pre-trained network to solve the Tolman–Oppenheimer–Volkoff (TOV) equations. The gradient-based optimization process incorporates a Bayesian picture into the proposed framework. The reconstructed EoS is proven to be consistent with the results from conventional methods. Furthermore, the resulting tidal deformation is in agreement with the limits obtained from the gravitational wave event, GW170817.

Searching for neutral particles at the highest energies at the Pierre Auger Observatory

The Pierre Auger Observatory, being the largest air-shower experiment in the world, offers an unprecedented exposure to neutral particles at the highest energies. Since the beginning of data collection more than 18 years ago, several searches for ultra-high-energy (UHE,E> 1017eV) photons and neutrinos have been performed. The upper limits on the diffuse flux of UHE photons and neutrinos derived from Auger data are among the most stringent in the world, severely constraining models for the origin of UHE cosmic rays. In addition, the Pierre Auger Observatory contributes to current efforts in multimessenger astronomy through follow-up searches for UHE photons and neutrinos in association with transient events, such as gravitational wave events. The various activities concerning searches for UHE photons and neutrinos in the data from the Pierre Auger Observatory are presented and the current results are summarized. In addition, future perspectives will be discussed.

Kilonovae and Optical Afterglows from Binary Neutron Star Mergers. II. Optimal Search Strategy for Serendipitous Observations and Target-of-opportunity Observations of Gravitational Wave Triggers

In the second work of this series, we explore the optimal search strategy for serendipitous and gravitational-wave-triggered target-of-opportunity (ToO) observations of kilonovae and optical short-duration gamma-ray burst (sGRB) afterglows from binary neutron star (BNS) mergers, assuming that cosmological kilonovae are AT2017gfo-like (but with viewing-angle dependence) and that the properties of afterglows are consistent with those of cosmological sGRB afterglows. A one-day cadence serendipitous search strategy with an exposure time of ∼30 s can always achieve an optimal search strategy of kilonovae and afterglows for various survey projects. We show that the optimal detection rates of the kilonovae (afterglows) are ∼0.3/0.6/1/20 yr−1 (∼50/60/100/800 yr−1) for Zwicky the Transient Facility (ZTF)/Multi-channel Photometric Survey Telescope (Mephisto)/Wide Field Survey Telescope (WFST)/Large Synoptic Survey Telescope (LSST), respectively. A better search strategy for SiTian than the current design is to increase the exposure time. In principle, a fully built SiTian can detect ∼7(2000) yr−1 kilonovae (afterglows). Population properties of electromagnetic (EM) signals detected by serendipitous observations are studied in detail. For ToO observations, we predict that one can detect ∼11 yr−1 BNS gravitational wave (GW) events during the fourth observing run (O4) by considering an exact duty cycle of the third observing run. The median GW sky localization area is expected to be ∼10 deg2 for detectable BNS GW events. For O4, we predict that ZTF/Mephisto/WFST/LSST can detect ∼5/4/3/3 kilonovae (∼1/1/1/1 afterglows) per year, respectively. The GW detection rates, GW population properties, GW sky localizations, and optimistic ToO detection rates of detectable EM counterparts for BNS GW events at the Advanced Plus, LIGO Voyager, and ET&CE eras are detailedly simulated in this paper.

DeepGlow: An efficient neural network emulator of physical afterglow models for gamma-ray bursts and gravitational-wave events

Abstract Gamma-ray bursts (GRBs) and double neutron star merger gravitational-wave events are followed by afterglows that shine from X-rays to radio, and these broadband transients are generally interpreted using analytical models. Such models are relatively fast to execute, and thus easily allow estimates of the energy and geometry parameters of the blast wave, through many trial-and-error model calculations. One problem, however, is that such analytical models do not capture the underlying physical processes as well as more realistic relativistic numerical hydrodynamic (RHD) simulations do. Ideally, those simulations are used for parameter estimation instead, but their computational cost makes this intractable. To this end, we present DeepGlow, a highly efficient neural network architecture trained to emulate a computationally costly RHD-based model of GRB afterglows, to within a few percent accuracy. As a first scientific application, we compare both the emulator and a different analytical model calibrated to RHD simulations, to estimate the parameters of a broadband GRB afterglow. We find consistent results between these two models, and also give further evidence for a stellar wind progenitor environment around this GRB source. DeepGlow fuses simulations that are otherwise too complex to execute over all parameters, to real broadband data of current and future GRB afterglows.

The Late Afterglow of GW170817/GRB 170817A: A Large Viewing Angle and the Shift of the Hubble Constant to a Value More Consistent with the Local Measurements

The multimessenger data of neutron star merger events are promising for constraining the Hubble constant. So far, GW170817 is still the unique gravitational wave event with multiwavelength electromagnetic counterparts. In particular, its radio and X-ray emission have been measured in the past ∼3–5 yr. In this work, we fit the long-lasting X-ray, optical, and radio afterglow light curves of GW170817/GRB 170817A, including the forward shock radiation from both the decelerating relativistic gamma-ray burst outflow and the subrelativistic kilonova outflow (though whether the second component contributes significantly is still uncertain), and find out a relatively large viewing angle (∼0.5 rad). Such a viewing angle has been taken as a prior in the gravitational wave data analysis, and the degeneracy between the viewing angle and the luminosity distance is broken. Finally, we have a Hubble constant , which is more consistent with that obtained by other local measurements. If rather similar values are inferred from multimessenger data of future neutron star merger events, it will provide critical support to the existence of the Hubble tension.

Testing time evolution of mass distribution of black hole mergers

The detection of gravitational-wave events revealed that there are numerous populations of black hole (BH) binaries that can merge within the age of the Universe. Although several formation channels of such binaries are known, considerable theoretical uncertainties associated with each channel defeat the robust prediction of how much each channel contributes to the total merger rate density. Given that the time evolution of the merger rate density in some channels is (exactly or nearly) independent of BH masses, clarifying this feature from observational data will shed some light on the nature of BH binaries. On the basis of this motivation, we formulate a methodology to perform a statistical test of whether the mass distribution of BH mergers evolves over time by hypothesis testing. Our statistical test requires neither a priori specification of the mass distribution, which is largely uncertain, nor that of the time dependence of merger rate. We then apply it to mock data for some concrete shapes of the merger rate density and show that the proposed method rejects/(does not reject) the null hypothesis correctly for a large sample size. After this verification, the method is applied to a catalog of the gravitational-wave events obtained during the LIGO-Virgo’s third observing run. We find that the selection bias degrades the effectiveness of our method for the O3 catalog owing to the reduction in the number of and the maximum redshifts of the merger events that we can explore. Within the range where the method can be applied, there is no indication of the time evolution of the mass distribution of merger rate density. This limitation will be eased in future observations containing more events, and our hypothesis testing will help determine whether the merger rate density evolves over time independently of BH masses.

From inflation to black hole mergers and back again: Gravitational-wave data-driven constraints on inflationary scenarios with a first-principle model of primordial black holes across the QCD epoch

Recent population studies have searched for a subpopulation of primordial black holes (PBHs) in the gravitational-wave (GW) events so far detected by LIGO/Virgo/KAGRA (LVK), in most cases adopting a phenomenological PBH mass distribution. When deriving such population from first principles in the standard scenario, however, the equation of state of the Universe at the time of PBH formation may strongly affect the PBH abundance and mass distribution, which ultimately depend on the power spectrum of cosmological perturbations. Here we improve on previous population studies on several aspects: (i) we adopt state-of-the-art PBH formation models describing the collapse of cosmological perturbations across the QCD epoch; (ii) we perform the first Bayesian multi-population inference on GW data including PBHs and directly using power spectrum parameters instead of phenomenological distributions; (iii) we critically confront the PBH scenario with LVK phenomenological models describing the GWTC-3 catalog both in the neutron-star and in the BH mass ranges, also considering PBHs as subpopulation of the total events. Our results confirm that LVK observations prevent the majority of the dark matter to be in the form of stellar mass PBHs. We find that the best fit PBH model can comprise a small fraction of the total events, in particular it can naturally explain events in the mass gaps. If the lower mass-gap event GW190814 is interpreted as a PBH binary, we predict that LVK should detect up to a few subsolar mergers and one to $\approx 30$ lower mass gap events during the upcoming O4 and O5 runs. Finally, mapping back the best-fit power spectrum into an ultra slow-roll inflationary scenario, we show that the latter predicts detectable PBH mergers in the LVK band, a stochastic GW background detectable by current and future instruments, and may include the entirety of dark matter in asteroid-mass PBHs.

Universal Nuclear Equation of State Introducing the Hypothetical X17 Boson

Within the scope of the Symmetry journal special issue on: “The Nuclear Physics of Neutron Stars”, we complemented the nuclear equation of state (EoS) with a hypothetical 17 MeV boson and observed that only instances with an admixture of 30%–40% satisfy all of the constraints. The successful EoS resulted in a radius of around 13 km for a neutron star with mass MNS≈1.4M⊙ and in a maximum mass of around MNS≈2.5M⊙. The value of the radius is in agreement with the recent measurement by NICER. The maximum mass is also in agreement with the mass of the remnant of the gravitational wave event GW190814. Thus, it appears that these EoSs satisfy all of the existing experimental constraints and can be considered as universal nuclear equations of state.

Wave effect of gravitational waves intersected with a microlens field: A new algorithm and supplementary study

The increase in gravitational wave (GW) events has allowed receiving strong lensing image pairs of GWs. However, the wave effect (diffraction and interference) due to the microlens field contaminates the parameter estimation of the image pair, which may lead to a misjudgment of strong lensing signals. To quantify the influence of the microlens field, researchers need a large sample of statistical research. Nevertheless, due to the oscillation characteristic, the Fresnel-Kirchhoff diffraction integral's computational time hinders this aspect's study. Although many algorithms are available, most cannot be well applied to the case where the microlens field is embedded in galaxy/galaxy clusters. This work proposes a faster and more accurate algorithm for studying the wave optics effect of microlenses embedded in different types of strong lensing images. Additionally, we provide a quantitative estimation criterion for the lens plane boundary for the Fresnel-Kirchhoff diffraction integral. This algorithm can significantly facilitate the study of wave optics, particularly in the case of microlens fields embedded in galaxy/galaxy clusters.

Pulsar revival in neutron star mergers: multimessenger prospects for the discovery of pre-merger coherent radio emission

ABSTRACT We investigate pre-merger coherent radio emission from neutron star mergers arising due to the magnetospheric interaction between compact objects. We consider two plausible radiation mechanisms, and show that if one neutron star has a surface magnetic field Bs ≥ 1012G, coherent millisecond radio bursts with characteristic temporal morphology and inclination angle dependence are observable to Gpc distances with next-generation radio facilities. We explore multi-messenger and multi-wavelength methods of identification of a neutron star merger origin of radio bursts, such as in fast radio burst surveys, triggered observations of gamma-ray bursts and gravitational wave events, and optical/radio follow-up of fast radio bursts in search of kilonova and radio afterglow emission. We present our findings for current and future observing facilities, and make recommendations for verifying or constraining the model.

Deep learning model based on a bidirectional gated recurrent unit for the detection of gravitational wave signals

The application of deep learning to the detection of gravitational wave (GW) signals has attracted much attention in recent years, and nearly all of the existing deep learning-based GW detection algorithms are based on the convolutional neural network (CNN), which can be used to identify the GW merger moment and produce responding triggers. However, current CNN-based models are not best optimized tools in modeling the causal logic of sequential data in the time dimension to incorporate the information of the inspiral, merger and ringdown phases of the GW waveform. To address this issue we propose a deep learning model designed based on the bidirectional gated recurrent unit (BiGRU), a variant of the recurrent neural network (RNN) that is sensitive to the sequential logic of GW data. We train the neural network and implement the GW detection task by using the data obtained during the first, second and third observation runs (O1/O2/O3) of the LIGO-VIGRO scientific collaboration. We find that our BiGRU model achieves a false alarm probability (FAP) of $4.88\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}4}$ (corresponding false alarm rate is about 1 per 18.2 hours) and maintains a detection ratio of about 89.6%. We also construct a model ensemble containing two independently trained BiGRU models and find it does not trigger any false alarm events on our testing dataset. Meanwhile, our BiGRU model not only recovers all the reported GW events of O1-O2 runs, including the unique (as of today) confident (with observed electromagnetic counterparts) GW event originating from a binary neutron star (BNS) merger (GW170817), but also achieves a detection ratio of about 80% of the O3 run. By calculating and comparing the receiver operator characteristic (ROC) curves we demonstrate the superiority of the BiGRU model over traditional CNN-based models, and show that the BiGRU model, being optimized to deal with GW waveform that sequentially contains the inspiral, merger and ringdown phases, is a potentially powerful tool to study the GW data in the future.

Observational constraints on maximum mass limit and physical properties of anisotropic strange star models by gravitational decoupling in Einstein–Gauss–Bonnet gravity

ABSTRACT In this work, we are guided by the gravitational wave events GW 170817 and GW 190814 together with observations of neutron stars PSR J1614-2230, PSR J1903+6620, and LMC X-4 to model compact objects within the framework of Einstein–Gauss–Bonnet (EGB) gravity. In addition, we employ the extended gravitational decoupling (EGD) method to explore the impact of anisotropy by varying the decoupling parameter. We model strange quark stars in which the interior stellar fluid obeys the MIT Bag equation of state which represents a degenerated Fermi gas comprising of up, down, and strange quarks. In order to close the system of field equations describing the seed solution, we employ the Buchdahl ansatz for one of the metric functions. The θ sector is solved under the bifurcation: $\epsilon =\theta ^0_0$ and $P_r=\theta ^1_1$ leading to two new families of solutions. In order to test the physical viability of the models, we vary the EGB parameter (α) or the decoupling constant (β) to achieve the observed masses and radii of compact objects. Our models are able to account for low-mass stars for a range of β values while α is fixed. The present models mimic the secondary component of the GW 190814 with a mass range of 2.5–2.67 M⊙ and radii typically of the order of 11.76$^{+0.14}_{-0.19}$ km for large values of the EGB parameter and the decoupling constant. The energy exchange between fluids inside the stellar object is sensitive to model parameters which lead to stable configurations.

A Fisher matrix for gravitational-wave population inference

ABSTRACT We derive a Fisher matrix for the parameters characterizing a population of gravitational-wave events. This provides a guide to the precision with which population parameters can be estimated with multiple observations, which becomes increasingly accurate as the number of events and the signal-to-noise ratio of the sampled events increase. The formalism takes into account individual event measurement uncertainties and selection effects, and can be applied to arbitrary population models. We illustrate the framework with two examples: an analytical calculation of the Fisher matrix for the mean and variance of a Gaussian model describing a population affected by selection effects, and an estimation of the precision with which the slope of a power-law distribution of supermassive black hole masses can be measured using extreme-mass-ratio inspiral observations. We compare the Fisher predictions to results from Monte Carlo analyses, finding very good agreement.

Strategy and Results of MASTER Network Follow-Up Observations of LIGO and Virgo Gravitational Wave Events within the Observational Sets O1, O2, and O3

Strategy and Results of MASTER Network Follow-Up Observations of LIGO and Virgo Gravitational Wave Events within the Observational Sets O1, O2, and O3

Repeated Mergers of Black Hole Binaries: Implications for GW190521

The gravitational wave event GW190521 involves the merger of two black holes of ∼85 M ⊙ and ∼66 M ⊙ forming an intermediate-mass black hole (IMBH) of mass ∼142 M ⊙. Both progenitors are challenging to explain within standard stellar evolution as they are within the upper black hole mass gap. We propose a dynamical formation pathway for this IMBH based on multiple mergers in the core of a globular cluster. We identify such scenarios from analysis of a set of 58 N-body simulations using NBODY6-gpu. In one of our simulations, we observe a stellar black hole undergoing a chain of seven binary mergers within 6 Gyr, attaining a final mass of 97.8 M ⊙. We discuss the dynamical interactions that lead to the final IMBH product, as well as the evolution of the black hole population in that simulation. We explore statistically the effects of gravitational recoil on the viability of such hierarchical mergers. From the analysis of all 58 simulations we observe additional smaller chains, tentatively inferring that an IMBH formation through hierarchical mergers is expected in the lifetime of a median-mass globular cluster with probability 0.01 ≲ p ≲ 0.1 without gravitational merger recoil. Using this order-of-magnitude estimate we show that our results are broadly consistent with the rate implied by GW190521, assuming that gravitational recoil ejection of progenitors has a low probability. We discuss implications for future gravitational wave detections, emphasizing the importance of studying such formation pathways for black holes within the upper mass gap as a means to constrain such modeling.

Contracting and Expanding Binary Black Holes in 3D Low-mass AGN Disks: The Importance of Separation

LIGO/Virgo has detected several binary black hole (BBH) merger events that may have originated in the accretion disks of active galactic nuclei (AGNs). These events require individual black hole masses that fall within the pair instability supernova mass gap, and therefore these black holes may have been grown from hierarchical mergers. AGN disks are a prime environment for hierarchical mergers and thus a potential location for the progenitors of BBH gravitational wave events. Understanding how a BBH embedded in an AGN disk interacts with the surrounding environment is thus crucial for determining whether this interaction can lead to its merger. However, there are few high-fidelity simulations of this process, and almost all are two-dimensional. We present the results from 3D, high-resolution, local shearing-box simulations of an embedded BBH interacting with an AGN disk. In these first simulations of their kind, we focus on determining the mass accretion rate and the orbital evolution rate at different BBH separations. We find that circular, equal-mass BBHs with separations greater than 10% of their Hill radius contract while accreting at a super-Eddington rate. At smaller separations, however, our 3D simulations find that BBHs expand their orbits. This result suggests that it may be difficult for an AGN disk to push a BBH to merger, but we discuss several mechanisms, including MHD turbulence and radiative and mechanical feedback, which could alleviate this difficulty.