Network Of Gravitational Wave Detectors Research Articles

Calibrating the global network of gravitational wave observatories via laser power calibration at NIST and PTB

Abstract Current gravitational wave observatories rely on Photon Calibrators (Pcals) that use laser radiation pressure to generate displacement fiducials used to calibrate detector output signals. Reducing calibration uncertainty enables optimal extraction of astrophysical information such as source distance and sky position from detected signals. For the ongoing O4 observation run that started on May 24, 2023, the global gravitational wave detector network is employing a new calibration scheme with transfer standards calibrated at both the National Institute of Standards and Technology (NIST) and the Physikalisch-Technische Bundesanstalt (PTB). These transfer standards will circulate between the observatories and the metrology institutes to provide laser power calibration traceable to the International System of Units (SI) and enable assessment and reduction of relative calibration errors for the observatory network. The Laser Interferometer Gravitational-Wave Observatory (LIGO) project and the Virgo project are currently participating in the new calibration scheme. The Large-scale Cryogenic Gravitational-wave Telescope project (KAGRA) is expected to join in 2024, with the LIGO Aundha Observatory (LAO) in India joining later. Before implementing this new scheme, a NIST-PTB bilateral comparison was conducted. The results of this comparison, with significantly lower uncertainty than previous studies, are reported. We also describe the transfer of power sensor calibration, including detailed uncertainty estimates, from the transfer standards calibrated by NIST and PTB to the sensors operating continuously at the interferometer end stations. Finally, we discuss the ongoing calibration of Pcal-induced displacement fiducials for the O4 observing run. Achieved combined standard uncertainty levels as low as 0.3 % facilitate calibrating the interferometer output signals with sub-percent accuracy.

Multi-messenger prospects for black hole – neutron star mergers in the O4 and O5 runs

The existence of merging black hole-neutron star (BHNS) binaries has been ascertained through the observation of their gravitational wave (GW) signals. However, to date, no definitive electromagnetic (EM) emission has been confidently associated with these mergers. Such an association could help unravel crucial information on these systems, for example, their BH spin distribution, the equation of state (EoS) of the neutron star and the rate of heavy element production. We modeled the multi-messenger (MM) emission from BHNS mergers detectable during the fourth (O4) and fifth (O5) observing runs of the LIGO-Virgo-KAGRA (LVK) GW detector network in order to provide detailed predictions that can help enhance the effectiveness of observational efforts and extract the highest possible scientific information from such remarkable events. Our methodology is based on a population synthesis approach, which includes the modeling of the signal-to-noise ratio of the GW signal in the detectors, the GW-inferred sky localization of the source, the kilonova (KN) optical and near-infrared light curves, the relativistic jet gamma-ray burst (GRB) prompt emission peak photon flux, and the GRB afterglow light curves in the radio, optical, and X-ray bands. The resulting prospects for BHNS MM detections during O4 are not promising, with an LVK GW detection rate of 15.0−8.8+15.4 yr−1, but joint MM rates of ∼10−1 yr−1 for the KN and ∼10−2 yr−1 for the jet-related emission. In O5, we found an overall increase in expected detection rates by around an order of magnitude, owing to both the enhanced sensitivity of the GW detector network and the coming online of future EM facilities. Considering variations in the NS EoS and BH spin distribution, we find that the detection rates can increase further by up to a factor of several tens. Finally, we discuss direct searches for the GRB radio afterglow with large field-of-view instruments during O5 and beyond as a new possible follow-up strategy in the context of ever-dimming prospects for KN detection due to the recession of the GW horizon.

Spin Doctors: How to Diagnose a Hierarchical Merger Origin

Gravitational-wave observations provide the unique opportunity of studying black hole formation channels and histories—but only if we can identify their origin. One such formation mechanism is the dynamical synthesis of black hole binaries in dense stellar systems. Given the expected isotropic distribution of component spins of binary black holes in gas-free dynamical environments, the presence of antialigned or in-plane spins with respect to the orbital angular momentum is considered a tell-tale sign of a merger’s dynamical origin. Even in the scenario where birth spins of black holes are low, hierarchical mergers attain large component spins due to the orbital angular momentum of the prior merger. However, measuring such spin configurations is difficult. Here, we quantify the efficacy of the spin parameters encoding aligned-spin (χ eff) and in-plane spin (χ p ) at classifying such hierarchical systems. Using Monte Carlo cluster simulations to generate a realistic distribution of hierarchical merger parameters from globular clusters, we can infer mergers’ χ eff and χ p . The cluster populations are simulated using Advanced LIGO-Virgo sensitivity during the detector network’s third observing period and projections for design sensitivity. Using a “likelihood-ratio”-based statistic, we find that ∼2% of the recovered population by the current gravitational-wave detector network has a statistically significant χ p measurement, whereas no χ eff measurement was capable of confidently determining a system to be antialigned with the orbital angular momentum at current detector sensitivities. These results indicate that measuring spin-precession through χ p is a more detectable signature of hierarchical mergers and dynamical formation than antialigned spins.

Joint gravitational wave-short GRB detection of binary neutron star mergers with existing and future facilities

ABSTRACT We explore the joint detection prospects of short gamma-ray bursts (sGRBs) and their gravitational wave (GW) counterparts by the current and upcoming high-energy GRB and GW facilities from binary neutron star (BNS) mergers. We consider two GW detector networks: (1) a four-detector network comprising LIGO Hanford, Livingston, Virgo, and Kagra (IGWN4) and (2) a future five-detector network including the same four detectors and LIGO India (IGWN5). For the sGRB detection, we consider existing satellites Fermi and Swift and the proposed all-sky satellite Daksha. Most of the events for the joint detection will be off-axis, hence, we consider a broad range of sGRB jet models predicting the off-axis emission. Also, to test the effect of the assumed sGRB luminosity function, we consider two different functions for one of the emission models. We find that for the different jet models, the joint sGRB and GW detection rates for Fermi and Swift with IGWN4 (IGWN5) lie within 0.07–0.62 yr−1 0.8–4.0 yr−1) and 0.02–0.14 yr−1 (0.15–1.0 yr−1), respectively, when the BNS merger rate is taken to be 320 Gpc−3 yr−1. With Daksha, the rates increase to 0.2–1.3 yr−1 (1.3–8.3 yr−1), which is 2–9 times higher than the existing satellites. We show that such a mission with higher sensitivity will be ideal for detecting a higher number of fainter events observed off-axis or at a larger distance. Thus, Daksha will boost the joint detections of sGRB and GW, especially for the off-axis events. Finally, we find that our detection rates with optimal SNRs are conservative, and noise in GW detectors can increase the rates further.

GW190425: Pan-STARRS and ATLAS coverage of the skymap and limits on optical emission associated with FRB 20190425A

ABSTRACT GW190425 is the second of two binary neutron star (BNS) merger events to be significantly detected by the Laser Interferometer Gravitational Wave (GW) Observatory (LIGO), Virgo and the Kamioka Gravitational Wave (KAGRA) detector network. With a detection only in LIGO Livingston, the skymap containing the source was large and no plausible electromagnetic counterpart was found in real-time searching in 2019. Here, we summarize Asteroid Terrestrial-Impact Last Alert System (ATLAS) and Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) wide-field optical coverage of the skymap beginning within 1 and 3 h, respectively, of the GW190425 merger time. More recently, a potential coincidence between GW190425 and a fast radio burst FRB 20190425A has been suggested, given their spatial and temporal coincidences. The smaller sky localization area of FRB 20190425A and its dispersion measure led to the identification of a likely host galaxy, UGC 10667 at a distance of 141 ± 10 Mpc. Our optical imaging covered the galaxy 6.0 h after GW190425 was detected and 3.5 h after the FRB 20190425A. No optical emission was detected and further imaging at +1.2 and +13.2 d also revealed no emission. If the FRB 20190425A and GW190425 association were real, we highlight our limits on kilonova emission from a BNS merger in UGC 10667. The model for producing FRB 20190425A from a BNS merger involves a supramassive magnetized neutron star spinning down by dipole emission on the time-scale of hours. We show that magnetar-enhanced kilonova emission is ruled out by optical upper limits. The lack of detected optical emission from a kilonova in UGC 10667 disfavours, but does not disprove, the FRB–GW link for this source.

Rapid localization and inference on compact binary coalescences with the Advanced LIGO-Virgo-KAGRA gravitational-wave detector network

Rapid localization and inference on compact binary coalescences with the Advanced LIGO-Virgo-KAGRA gravitational-wave detector network

Subtraction of the confusion foreground and parameter uncertainty of resolvable galactic binaries on the networks of space-based gravitational-wave detectors

Subtraction of the confusion foreground and parameter uncertainty of resolvable galactic binaries on the networks of space-based gravitational-wave detectors

The short gamma-ray burst population in a quasi-universal jet scenario

We present a model of the short gamma-ray burst (SGRB) population under a ‘quasi-universal jet’ scenario in which jets can differ somewhat in their on-axis peak prompt emission luminosity, Lc, but share a universal angular luminosity profile, ℓ(θv) = L(θv)/Lc, as a function of the viewing angle, θv. The model was fitted, through a Bayesian hierarchical approach inspired by gravitational wave (GW) population analyses, to three observed SGRB samples simultaneously: the Fermi/GBM sample of SGRBs with spectral information available in the catalogue (367 events); a flux-complete sample of 16 Swift/BAT SGRBs that are also detected by the GBM and have a measured redshift; and a sample of SGRBs with a binary neutron star (BNS) merger counterpart, which only includes GRB 170817A at present. Particular care was put into modelling selection effects. The resulting model, which reproduces the observations, favours a narrow jet ‘core’ with half-opening angle θc = 2.1−1.4+2.4 deg (uncertainties hereon refer to 90% credible intervals from our fiducial ‘full sample’ analysis) whose peak luminosity, as seen on-axis, is distributed as a power law, p(Lc) ∝ Lc−A with A = 3.2−0.4+0.7, above a minimum isotropic-equivalent luminosity, Lc⋆ = 5−2+11 × 1051 erg s−1. For viewing angles larger than θc, the luminosity profile scales as a single power law, l ∝ θv−αL with αL = 4.7−1.4+1.2, with no evidence of a break, despite the model allowing for it. While the model implies an intrinsic ‘Yonetoku’ correlation between L and the peak photon energy, Ep, of the spectral energy distribution, its slope is somewhat shallower, Ep ∝ L0.4 ± 0.2, than the apparent one, and the normalisation is offset towards larger Ep due to selection effects. The implied local rate density of SGRBs (regardless of the viewing angle) is between about one hundred up to several thousand events per cubic gigaparsec per year, in line with the BNS merger rate density inferred from GW observations. Based on the model, we predict 0.2 to 1.3 joint GW+SGRB detections per year by the advanced GW detector network and Fermi/GBM during the O4 observing run.

Waves in a forest: a random forest classifier to distinguish between gravitational waves and detector glitches

The LIGO-Virgo-KAGRA (LVK) network of gravitational-wave (GW) detectors have observed many tens of compact binary mergers to date. Transient, non-Gaussian noise excursions, known as ‘glitches’, can impact signal detection in various ways. They can imitate true signals as well as reduce the confidence of real signals. In this work, we introduce a novel statistical tool to distinguish astrophysical signals from glitches, using their inferred source parameter posterior distributions as a feature set. By modelling both simulated GW signals and real detector glitches with a gravitational waveform model, we obtain a diverse set of posteriors which are used to train a random forest classifier. We show that random forests can identify differences in the posterior distributions for signals and glitches, aggregating these differences to tell apart signals from common glitch types with high accuracy of over 93%. We conclude with a discussion on the regions of parameter space where the classifier is prone to making misclassifications, and the different ways of implementing this tool into LVK analysis pipelines.

GWDALI: A Fisher-matrix based software for gravitational wave parameter-estimation beyond Gaussian approximation

We introduce GWDALI, a new Fisher-matrix, python based software that computes likelihood gradients to forecast parameter-estimation precision of arbitrary network of terrestrial gravitational wave detectors observing compact binary coalescences. The main new feature with respect to analogous software is to assess parameter uncertainties beyond Fisher-matrix approximation, using the derivative approximation for Likelihood (DALI). The software makes optional use of the LSC algorithm library LAL and the stochastic sampling algorithm Bilby, which can be used to perform Monte-Carlo sampling of exact or approximate likelihood functions. As an example we show comparison of estimated precision measurement of selected astrophysical parameters for both the actual likelihood, and for a variety of its derivative approximations, which turn out particularly useful when the Fisher matrix is not invertible.

Overview of KAGRA : Data transfer and management

Abstract KAGRA is a newly built gravitational wave observatory, a laser interferometer with a 3 km arm length, located in Kamioka, Gifu prefecture, Japan. In this article, we describe the KAGRA data management system, i.e., recording of data, transfer from the KAGRA experiment site to computing resources, as well as data distribution to tier sites, including international sites in Taiwan and Korea. The amount of KAGRA data exceeded 1.0 PiB and increased by about 1.5 TB per day during operation in 2020. Our system has succeeded in data management, and has achieved performance that can withstand observations after 2023, that is, a transfer rate of 20 MB s-1 or more and file storage of sufficient capacity for petabyte class. We also discuss the sharing of data between the global gravitational-wave detector network with other experiments, namely LIGO and Virgo. The latency, which consists of calculation of calibrated strain data and transfer time within the global network, is very important from the view of multi-messenger astronomy using gravitational waves. Real-time calbrated data delivered from the KAGRA detector site and other detectors to our computing system arrive with about 4–15 seconds of latency. These latencies are sufficiently short compared to the time taken for gravitational wave event search computations. We also established a high-latency exchange of offline calibrated data that was aggregated with a better accuracy compared with real-time data.

GWSim: Python package for creating mock GW samples for different astrophysical populations and cosmological models of binary black holes

Precision cosmology with gravitational wave (GW) sources requires a deeper understanding of the interplay between GW source population and cosmological parameters governing the dynamics of the Universe. With the swift increase in GW detections, it is necessary to develop a tool for exploring many aspects of cosmology and fundamental physics; this tools allows to simulate GW mock samples for several populations and cosmological models. We have developed a new code called GWSim, which allows us to make a large sample of GW mock events from a broad range of configurations, while varying the cosmology, the merger rate, and the GW source parameters (i.e. mass and spin distributions in particular) for a given network of GW detectors. A large sample of simulated mock GW events will be useful for improving our understanding of the statistical properties of the distribution of GW sources, as long as it is detectable for a given detector noise and an astrophysical and cosmological model. It will also be useful to compare simulated samples with the observed distribution of the GW sources from data and infer the underlying population of the GW source parameters and cosmology. We restricted the cosmology to spatially flat universes, including models with varying dark energy equation of state. The GWSim code provides each mock event with a position in the sky and a redshift; these values can be those of random host galaxies coming from an isotropic and homogeneous simulated Universe or a user-supplied galaxy catalog. We used realistic detector configurations of the LIGO and Virgo network of detectors to demonstrate the performance of this code for the latest observation runs and the upcoming observation run.

Virgo detector characterization and data quality: tools

Detector characterization and data quality studies—collectively referred to as DetChar activities in this article—are paramount to the scientific exploitation of the joint dataset collected by the LIGO-Virgo-KAGRA global network of ground-based gravitational-wave (GW) detectors. They take place during each phase of the operation of the instruments (upgrade, tuning and optimization, data taking), are required at all steps of the dataflow (from data acquisition to the final list of GW events) and operate at various latencies (from near real-time to vet the public alerts to offline analyses). This work requires a wide set of tools which have been developed over the years to fulfill the requirements of the various DetChar studies: data access and bookkeeping; global monitoring of the instruments and of the different steps of the data processing; studies of the global properties of the noise at the detector outputs; identification and follow-up of noise peculiar features (whether they be transient or continuously present in the data); quick processing of the public alerts. The present article reviews all the tools used by the Virgo DetChar group during the third LIGO-Virgo Observation Run (O3, from April 2019 to March 2020), mainly to analyze the Virgo data acquired at EGO. Concurrently, a companion article focuses on the results achieved by the DetChar group during the O3 run using these tools.

Using grey sirens to resolve the Hubble–Lemaître tension

ABSTRACT The measurement of the Hubble–Lemaître constant (H0) from the cosmic microwave background and the Type IA supernovae are at odds with each other. One way to resolve this tension is to use an independent way to measure H0. This can be accomplished by using gravitational-wave (GW) observations. Previous works have shown that with the onset of the next generation of GW detector networks, it will be possible to constrain H0 to better than 2 per cent precision (which is enough to resolve the tension) with binary black hole systems that are extremely well localized in the sky, also called golden dark sirens. Bright sirens like binary neutron star systems can also help resolve the tension if both the GW and the following electromagnetic counterpart are detected. In this work, we show that neutron star-black hole (NSBH) mergers can act both as golden dark sirens as well as bright sirens, thus, assigning them the term grey sirens. We assess the potential of using NSBH mergers to measure H0 and find that the Voyager network might be able to resolve the tension in an observation span of 5 yr. The next generation networks, which include the Cosmic Explorer detectors and the Einstein Telescope will be able to measure H0 to sub-per cent level just by using NSBH mergers.

Development of advanced photon calibrator for Kamioka gravitational wave detector (KAGRA).

The Kamioka Gravitational wave detector (KAGRA) cryogenic gravitational-wave observatory has commenced joint observations with the worldwide gravitational wave detector network. Precise calibration of the detector response is essential for accurately estimating parameters of gravitational wave sources. A photon calibrator is a crucial calibration tool used in laser interferometer gravitational-wave observatory, Virgo, and KAGRA, and it was utilized in joint observation3 with GEO600 in Germany in April 2020. In this paper, KAGRA implemented three key enhancements: a high-power laser, a power stabilization system, and remote beam position control. KAGRA employs a 20W laser divided into two beams that are injected onto the mirror surface. By utilizing a high-power laser, the response of the detector at kHz frequencies can be calibrated. To independently control the power of each laser beam, an optical follower servo was installed for power stabilization. The optical path of the photon calibrator's beam positions was controlled using pico-motors, allowing for the characterization of the detector's rotation response. Additionally, a telephoto camera and quadrant photodetectors were installed to monitor beam positions, and beam position control was implemented to optimize the mirror response. In this paper, we discuss the statistical errors associated with the measurement of relative power noise. We also address systematic errors related to the power calibration model of the photon calibrator and the simulation of elastic deformation effects using finite element analysis. Ultimately, we have successfully reduced the total systematic error from the photon calibrator to 2.0%.

Measuring the ringdown scalar polarization of gravitational waves in Einstein-scalar-Gauss-Bonnet gravity

We model the scalar waves produced during the ringdown stage of binary black hole coalescence in Einstein scalar Gauss-Bonnet (EsGB) gravity, using numerical relativity simulations of the theory in the decoupling limit. Through a conformal coupling of the scalar field to the metric in the matter-field action, we show that the gravitational waves in this theory can have a scalar polarization. We model the scalar quasi-normal modes of the ringdown signal in EsGB gravity, and quantify the extent to which current and future gravitational wave detectors could observe the spectrum of scalar radiation emitted during the ringdown phase of binary black hole coalescence. We find that within the limits of the theory's coupling parameters set by current theoretical and observational constraints, the scalar ringdown signal from black hole remnants in the $10^1 - 10^3 \, M_{\odot}$ mass range is expected to be well below the detectability threshold with the current network of gravitational-wave detectors (LIGO-Virgo-KAGRA), but is potentially measurable with next-generation detectors such as the Einstein Telescope.

Analyses of overlapping gravitational wave signals using hierarchical subtraction and joint parameter estimation

ABSTRACT In the coming years, third-generation detectors such as Einstein Telescope and Cosmic Explorer will enter the network of ground-based gravitational-wave detectors. Their current design predicts a significantly improved sensitivity band with a lower minimum frequency than existing detectors. This, combined with the increased arm length, leads to two major effects: the detection of more signals and the detection of longer signals. Both will result in a large number of overlapping signals. It has been shown that such overlapping signals can lead to biases in the recovered parameters, which would adversely affect the science extracted from the observed binary merger signals. In this work, we analyse overlapping binary black hole coalescences with two methods to analyse multisignal observations: hierarchical subtraction and joint parameter estimation. We find that these methods enable a reliable parameter extraction in most cases and that joint parameter estimation is usually more precise but comes with higher computational costs.

Probing ensemble properties of vortex-avalanche pulsar glitches with a stochastic gravitational-wave background search

A stochastic gravitational-wave background (SGWB) is expected to be produced by the superposition of individually undetectable, unresolved gravitational-wave (GW) signals from cosmological and astrophysical sources. Such a signal can be searched with dedicated techniques using the data acquired by a network of ground-based GW detectors. In this work, we consider the astrophysical SGWB resulting from pulsar glitches, which are sudden increases in the rotational pulsar frequency, within our Galaxy. More specifically, we assume glitches to be associated with quantized, superfluid, vortex avalanches in the pulsars, and we model the SGWB from the superposition of GW bursts emitted during the glitching phase. We perform a cross-correlation search for this SGWB-like signal employing the data from the first three observation runs of Advanced LIGO and Virgo. Not having found any evidence for a SGWB signal, we set upper limits on the dimensionless energy density parameter ${\mathrm{\ensuremath{\Omega}}}_{\mathrm{gw}}(f)$ for two different power-law SGWBs, corresponding to two different glitch regimes. We obtain ${\mathrm{\ensuremath{\Omega}}}_{\mathrm{gw}}(f)\ensuremath{\le}7.5\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}10}$ at 25 Hz for a spectral index $5/2$, and ${\mathrm{\ensuremath{\Omega}}}_{\mathrm{gw}}(f)\ensuremath{\le}5.7\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}17}$ at 25 Hz for a spectral index $17/2$. We then use these results to set constraints on the average glitch duration and the average radial motion of the vortices during the glitches for the population of the glitching Galactic pulsars, as a function of the Galactic glitch rate.

Measuring properties of primordial black hole mergers at cosmological distances: Effect of higher order modes in gravitational waves

Primordial black holes (PBHs) may form from the collapse of matter overdensities shortly after the Big Bang. One may identify their existence by observing gravitational wave (GW) emissions from merging PBH binaries at high redshifts $z\gtrsim 30$, where astrophysical binary black holes (BBHs) are unlikely to merge. The next-generation ground-based GW detectors, Cosmic Explorer and Einstein Telescope, will be able to observe BBHs with total masses of $\mathcal{O}(10-100)~M_{\odot}$ at such redshifts. This paper serves as a companion paper of arXiv:2108.07276, focusing on the effect of higher-order modes (HoMs) in the waveform modeling, which may be detectable for these high redshift BBHs, on the estimation of source parameters. We perform Bayesian parameter estimation to obtain the measurement uncertainties with and without HoM modeling in the waveform for sources with different total masses, mass ratios, orbital inclinations and redshifts observed by a network of next-generation GW detectors. We show that including HoMs in the waveform model reduces the uncertainties of redshifts and masses by up to a factor of two, depending on the exact source parameters. We then discuss the implications for identifying PBHs with the improved single-event measurements, and expand the investigation of the model dependence of the relative abundance between the BBH mergers originating from the first stars and the primordial BBH mergers as shown in arXiv:2108.07276.

Gwfish: A simulation software to evaluate parameter-estimation capabilities of gravitational-wave detector networks

An important step in the planning of future gravitational-wave (GW) detectors and of the networks they will form is the estimation of their detection and parameter-estimation capabilities, which is the basis of science-case studies. Several future GW detectors have been proposed or are under development, which might also operate and observe in parallel. These detectors include terrestrial, lunar, and space-borne detectors. In this paper, we present gwfish,11github.com/janosch314/GWFish. a new software to simulate GW detector networks and to calculate measurement uncertainties based on the Fisher-matrix approximation. gwfish models the impact of detector motion on PE and makes it possible to analyze multiband scenarios, i.e., observation of a GW signal by different detectors in different frequency bands. We showcase a few examples for the Einstein Telescope (ET) including the sky-localization of binary neutron stars, and ET’s capability to measure the polarization of GWs.