LIGO-Virgo Network Research Articles

Measuring neutron star radius with second and third generation gravitational wave detector networks

Abstract The next generation of ground-based interferometric gravitational wave detectors will observe mergers of black holes and neutron stars throughout cosmic time. A large number of the binary neutron star merger events will be observed with extreme high fidelity, and will provide stringent constraints on the equation of state of nuclear matter. In this paper, we investigate the systematic improvement in the measurability of the equation of state with increase in detector sensitivity by combining constraints obtained on the radius of a $1.4 \, \mathrm{M}_{\odot}$ neutron star from a simulated source population. Since the measurability of the equation of state depends on its stiffness, we consider a range of realistic equations of state that span the current observational constraints. We show that a single 40km Cosmic Explorer detector can pin down the neutron star radius for a soft, medium and stiff equation of state to an accuracy of 10m within a decade, whereas the current generation of ground-based detectors like the Advanced LIGO-Virgo network would take $\mathcal{O}(10^5)$ years to do so for a soft equation of state. 

Impact of Higher Harmonics of Gravitational Radiation on the Population Inference of Binary Black Holes

Templates modeling just the dominant mode of gravitational radiation are generally sufficient for the unbiased parameter inference of near-equal-mass compact binary mergers. However, neglecting the subdominant modes can bias the inference if the binary is significantly asymmetric, very massive, or has misaligned spins. In this work, we explore if neglecting these subdominant modes in the parameter estimation of nonspinning binary black hole mergers can bias the inference of their population-level properties such as mass and merger redshift distributions. Assuming the design sensitivity of the advanced LIGO-Virgo detector network, we find that neglecting subdominant modes will not cause a significant bias in the population inference, although including them will provide more precise estimates. This is primarily because asymmetric binaries are expected to be rarer in our detected sample, due to their intrinsic rareness and the observational selection effects. The increased precision in the measurement of the maximum black hole mass can help in better constraining the upper mass gap in the mass spectrum.

Search for Gravitational-lensing Signatures in the Full Third Observing Run of the LIGO–Virgo Network

Gravitational lensing by massive objects along the line of sight to the source causes distortions to gravitational wave (GW) signals; such distortions may reveal information about fundamental physics, cosmology, and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO-Virgo network. We search for repeated signals from strong lensing by (1) performing targeted searches for subthreshold signals, (2) calculating the degree of overlap among the intrinsic parameters and sky location of pairs of signals, (3) comparing the similarities of the spectrograms among pairs of signals, and (4) performing dual-signal Bayesian analysis that takes into account selection effects and astrophysical knowledge. We also search for distortions to the gravitational waveform caused by (1) frequency-independent phase shifts in strongly lensed images, and (2) frequency-dependent modulation of the amplitude and phase due to point masses. None of these searches yields significant evidence for lensing. Finally, we use the nondetection of GW lensing to constrain the lensing rate based on the latest merger-rate estimates and the fraction of dark matter composed of compact objects.

Prospects for a Precise Equation of State Measurement from Advanced LIGO and Cosmic Explorer

Gravitational-wave observations of neutron star mergers can probe the nuclear equation of state by measuring the imprint of the neutron star’s tidal deformability on the signal. We investigate the ability of future gravitational-wave observations to produce a precise measurement of the equation of state from binary neutron star inspirals. Because measurability of the tidal effect depends on the equation of state, we explore several equations of state that span current observational constraints. We generate a population of binary neutron stars as seen by a simulated Advanced LIGO–Virgo network, as well as by a planned Cosmic Explorer observatory. We perform Bayesian inference to measure the parameters of each signal, and we combine measurements across each population to determine R 1.4, the radius of a 1.4 M ⊙ neutron star. We find that, with 321 signals, the LIGO–Virgo network is able to measure R 1.4 to better than 2% precision for all equations of state we consider; however, we also find that achieving this precision could take decades of observation, depending on the equation of state and the merger rate. On the other hand, we find that with one year of observation, Cosmic Explorer will measure R 1.4 to better than 0.6% precision. In both cases, we find that systematic biases, such as from an incorrect mass prior, can significantly impact measurement accuracy, and efforts will be required to mitigate these effects.

Prospects of probing dark matter condensates with gravitational waves

The Lambda-Cold Dark Matter model explains cosmological observations most accurately till date. However, it is still plagued with various shortcomings at galactic scales. Models of dark matter such as superfluid dark matter, Bose-Einstein Condensate(BEC) dark matter and fuzzy dark matter have been proposed to overcome some of these drawbacks. In this work, we probe these models using the current constraint on the gravitational wave (GW) propagation speed coming from the binary neutron star GW170817 detection by LIGO-Virgo detector network and use it to study the allowed parameter space for these three models for Advanced LIGO+Virgo, LISA, IPTA and SKA detection frequencies. The speed of GW has been shown to depend upon the refractive index of the medium, which in turn, depends on the dark matter model parameters through the density profile of the galactic halo. We constrain the parameter space for these models using the bounds coming from GW speed measurement and the Milky Way radius bound. Our findings suggest that with Advanced LIGO-Virgo detector sensitivity, the three models considered here remain unconstrained. A meaningful constraint can only be obtained for detection frequencies ≤ 10-9 Hz, which falls in the detection range of radio telescopes such as IPTA and SKA. Considering this best possible case, we find that out of the three condensate models, the fuzzy dark matter model is the most feasible scenario to be falsified/validated in near future.

Upper limits on the polarized isotropic stochastic gravitational-wave background from advanced LIGO-Virgo's first three observing runs

Parity violation is expected to generate an asymmetry between the amplitude of left and right-handed gravitational-wave modes which leads to a circularly polarized stochastic gravitational-wave background (SGWB). Due to the three independent baselines in the LIGO-Virgo network, we focus on the amplitude difference in strain power characterized by Stokes' parameters and do maximum-likelihood estimation to constrain the polarization degree of SGWB. Our results indicate that there is no evidence for the circularly polarized SGWB in the data. Furthermore, by modeling the SGWB as a power-law spectrum, we place upper limit on the normalized energy density Ωgw(25 Hz) < 5.3 × 10-9 at 95% confidence level after marginalizing over the polarization degree and spectral index.

Parametrized tests of post-Newtonian theory using principal component analysis

Searching for departures from general relativity (GR) in more than one post-Newtonian (PN) phasing coefficients, called a \emph{multi-parameter test}, is known to be ineffective given the sensitivity of the present generation of gravitational-wave (GW) detectors. Strong degeneracies in the parameter space make the outcome of the test uninformative. We argue that Principal Component Analysis (PCA) can remedy this problem by constructing certain linear combinations of the original PN parameters that are better constrained by gravitational-wave observations. By analyzing binary black hole events detected during the first and second observing runs (O1 and O2) of LIGO/Virgo, we show that the two dominant principal components can capture the essence of a multi-parameter test. Combining five binary black hole mergers during O1/O2, we find that the dominant linear combination of the PN coefficients obtained from PCA is consistent with GR within the 0.38 standard deviation of the posterior distribution. Furthermore, using a set of simulated \emph{non-GR} signals in the three-detector LIGO-Virgo network with designed sensitivities, we find that the method is capable of excluding GR with high confidence as well as recovering the injected values of the non-GR parameters with good precision.

Constraints on Kerr-Newman black holes from merger-ringdown gravitational-wave observations

We construct a template to model the post-merger phase of a binary black hole coalescence in the presence of a remnant $U(1)$ charge. We include the quasi-normal modes typically dominant during a binary black hole coalescence, $(\ell,m,n) = \{(2,2,0), (2,2,1)\}$ and also present analytical fits for the quasinormal mode frequencies of a Kerr-Newman black hole in terms of its spin and charge, here also including the $(3,3,0)$ mode. Aside from astrophysical electric charge, our template can accommodate extensions of the Standard Model, such as a dark photon. Applying the model to LIGO-Virgo detections, we find that we are unable to distinguish between the charged and uncharged hypotheses from a purely post-merger analysis of the current events. However, restricting the mass and spin to values compatible with the analysis of the full signal, we obtain a 90th percentile bound $\bar{q} < 0.33$ on the black hole charge-to-mass ratio, for the most favorable case of GW150914. Under similar assumptions, by simulating a typical loud signal observed by the LIGO-Virgo network at its design sensitivity, we assess that this model can provide a robust measurement of the charge-to-mass ratio only for values $\bar{q} \gtrsim 0.5$; here we also assume that the mode amplitudes are similar to the uncharged case in creating our simulated signal. Lower values, down to $\bar{q} \sim 0.3$, could instead be detected when evaluating the consistency of the pre-merger and post-merger emission.

Prospects for Measuring Off-axis Spins of Binary Black Holes with Plus-era Gravitational-wave Detectors

Abstract The mass and spin properties of binary black holes (BBHs) inferred from their gravitational-wave signatures reveal important clues about how these systems form. BBHs originating from isolated binary evolution are expected to have spins preferentially aligned with their orbital angular momentum, whereas there is no such preference in binaries formed via dynamical assembly. The fidelity with which near-future gravitational-wave detectors can measure off-axis spins will have implications for the study of BBH formation channels. In this work, we examine the degree to which the Advanced LIGO Plus (A+) and Advanced Virgo Plus (AdV+) interferometric detectors can measure both aligned and misaligned spins. We compare spin resolution between the LIGO-Virgo network operating at either A+/AdV+ (“Plus”) sensitivity or Advanced-era design (“Design”) sensitivity using simulated BBH gravitational-wave signals injected into synthetic detector noise. The signals are distributed over the mass-spin parameter space of likely BBH systems, accounting for the effects of precession and higher-order modes. We find that the Plus upgrades yield significant improvements in spin estimation for systems with unequal masses and moderate or large spins. Using simulated signals modeled after different types of hierarchical BBH mergers, we also conclude that the Plus detector network will yield substantially improved spin estimates for 1G+2G binaries compared to the Design network.

Multimessenger analysis strategy for core-collapse supernova search: gravitational waves and low-energy neutrinos

Core-collapse supernovae are fascinating astrophysical objects for multimessenger studies. Gravitational waves are expected to play an important role in the supernova explosion mechanism. Unfortunately, their modeling is challenging, due to the stochastic nature of the dynamics and the vast range of possible progenitors. Therefore, the gravitational wave detection from these objects is still elusive with already advanced detectors. Low-energy neutrinos will be emitted copiously during the core-collapse explosion and can help the gravitational wave counterpart search. In this work, we develop a multimessenger strategy to look for such astrophysical objects. We exploit a global network of both low-energy neutrino and gravitational wave detectors. First, we discuss how to improve the detection potential of the neutrino sub-network by exploiting the time profile of a neutrino burst from a core-collapse supernova. We show that in the proposed approach, neutrino detectors can gain at least 10% of detection efficiency at the distance where their efficiency drops. Then, we combine the information provided by gravitational wave and neutrino signals in a multimessenger analysis. In particular, by using the clusters of low-energy neutrinos observed by LVD and KamLAND detectors in combination with the gravitational wave triggers from LIGO-Virgo detector network, we obtain an increase of the probability to detect the gravitational wave signal from a core-collapse supernova at 60 kpc, from zero to ∼33% for some specific gravitational wave emission model.

Deep learning for core-collapse supernova detection

The detection of gravitational waves from core-collapse supernova (CCSN) explosions is a challenging task, yet to be achieved, in which it is key the connection between multiple messengers, including neutrinos and electromagnetic signals. In this work, we present a method for detecting these kind of signals based on machine learning techniques. We tested its robustness by injecting signals in the real noise data taken by the Advanced LIGO-Virgo network during the second observation run, O2. We trained a newly developed Mini-Inception Resnet neural network using time-frequency images corresponding to injections of simulated phenomenological signals, which mimic the waveforms obtained in 3D numerical simulations of CCSNe. With this algorithm we were able to identify signals from both our phenomenological template bank and from actual numerical 3D simulations of CCSNe. We computed the detection efficiency versus the source distance, obtaining that, for signal to noise ratio higher than 15, the detection efficiency is 70 % at a false alarm rate lower than 5%. We notice also that, in the case of O2 run, it would have been possible to detect signals emitted at 1 kpc of distance, whilst lowering down the efficiency to 60%, the event distance reaches values up to 14 kpc.

A self-consistent method to estimate the rate of compact binary coalescences with a Poisson mixture model

The recently published GWTC-1 (Abbott B P et al (LIGO Scientific Collaboration and Virgo Collaboration) 2019 Phys. Rev. X 9 031040)—a journal article summarizing the search for gravitational waves (GWs) from coalescing compact binaries in data produced by the LIGO-Virgo network of ground-based detectors during their first and second observing runs—quoted estimates for the rates of binary neutron star, neutron star black hole binary, and binary black hole mergers, as well as assigned probabilities of astrophysical origin for various significant and marginal GW candidate events. In this paper, we delineate the formalism used to compute these rates and probabilities, which assumes that triggers above a low ranking statistic threshold, whether of terrestrial or astrophysical origin, occur as independent Poisson processes. In particular, we include an arbitrary number of astrophysical categories by redistributing, via mass-based template weighting, the foreground probabilities of candidate events, across source classes. We evaluate this formalism on synthetic GW data, and demonstrate that this method works well for the kind of GW signals observed during the first and second observing runs.

The updated DESGW processing pipeline for the third LIGO/VIRGO observing run

The DESGW group seeks to identify electromagnetic counterparts of gravitational wave events seen by the LIGO-VIRGO network, such as those expected from binary neutron star mergers or neutron star-black hole mergers. DESGW was active throughout the first two LIGO observing seasons, following up several binary black hole mergers and the first binary neutron star merger, GW170817. This work describes the modifications to the observing strategy generation and image processing pipeline between the second (ending in August 2017) and third (beginning in April 2019) LIGO observing seasons. The modifications include a more robust observing strategy generator, further parallelization of the image reduction software and difference imaging processing pipeline, data transfer streamlining, and a web page listing identified counterpart candidates that updates in real time. Taken together, the additional parallelization steps enable the identification of potential electromagnetic counterparts within fully calibrated search images in less than one hour, compared to the 3-5 hours it would typically take during the first two seasons. These performance improvements are critical to the entire EM follow-up community, as rapid identification (or rejection) of candidates enables detailed and rapid spectroscopic follow-up by multiple instruments, leading to more information about the environment immediately following such gravitational wave events.

Discrepancy in tidal deformability of GW170817 between the Advanced LIGO twin detectors

We find that the Hanford and Livingston detectors of Advanced LIGO derive a distinct posterior probability distribution of binary tidal deformability tilde{Lambda} of the first binary-neutron-star merger GW170817. By analyzing public data of GW170817 with a nested-sampling engine and the default TaylorF2 waveform provided by the LALInference package, the probability distribution of the binary tidal deformability derived by the LIGO-Virgo detector network turns out to be determined dominantly by the Hanford detector. Specifically, by imposing the flat prior on tidal deformability of individual stars, symmetric 90% credible intervals of tilde{Lambda} are estimated to be 527^{+619}_{-345} with the Hanford detector, 927^{+522}_{-619} with the Livingston detector, and 455^{+668}_{-281} with the LIGO-Virgo detector network. Furthermore, the distribution derived by the Livingston detector changes irregularly when we vary the maximum frequency of the data used in the analysis. This feature is not observed for the Hanford detector. While they are all consistent, the discrepancy and irregular behavior suggest that an in-depth study of noise properties might improve our understanding of GW170817 and future events.

Wider look at the gravitational-wave transients from GWTC-1 using an unmodeled reconstruction method

In this paper, we investigate the morphology of the events from the GWTC-1 catalog of compact binary coalescences as reconstructed by a method based on coherent excess power: we use an open-source version of the coherent WaveBurst (cWB) analysis pipeline, which does not make use of waveform models. The coherent response of the LIGO-Virgo network of detectors is estimated by using loose bounds on the duration and bandwidth of the signal. This pipeline version reproduces the same results that are reported for cWB in recent publications by the LIGO and Virgo collaborations. In particular, the sky localization and waveform reconstruction are in a good agreement with those produced by methods which exploit the detailed theoretical knowledge of the expected waveform for compact binary coalescences. However, in some cases cWB also detects features in excess in well-localized regions of the time-frequency plane. Here we focus on such deviations and present the methods devised to assess their significance. Out of the eleven events reported in the GWTC-1, in two cases -- GW151012 and GW151226 -- cWB detects an excess of coherent energy after the merger ($\Delta t \simeq 0.2$ s and $\simeq 0.1$ s, respectively) with p-values that call for further investigations ($0.004$ and $0.03$, respectively), though they are not sufficient to exclude noise fluctuations. We discuss the morphological properties and plausible interpretations of these features. We believe that the methodology described in the paper shall be useful in future searches for compact binary coalescences.

Constraining the Inclinations of Binary Mergers from Gravitational-wave Observations

Abstract Much of the information we hope to extract from the gravitational-wave signatures of compact binaries is only obtainable when we can accurately constrain the inclination of the orbital plane relative to the line of sight. In this paper, we discuss in detail a degeneracy between the measurement of the binary distance and inclination that limits our ability to accurately measure the inclination using gravitational waves alone. This degeneracy is exacerbated by the expected distribution of events in the universe, which leads us to prefer face-on systems at a greater distance. We use a simplified model that only considers the binary distance and orientation and show that this gives comparable results to the full parameter estimates obtained from the binary neutron star merger GW170817. For the advanced LIGO-Virgo network, it is only binaries that are close to edge-on, i.e., with inclinations ι ≳ 75°, that will be distinguishable from face-on systems. Extended networks that have good sensitivity to both gravitational-wave polarizations will only be able to constrain the inclination of a face-on binary at a signal-to-noise ratio of 20 to ι ≲ 45°. Even for loud signals with signal-to-noise ratios of 100, face-on signals will only be constrained to have inclinations . In the absence of observable higher modes or orbital precession, this degeneracy will dominate the mass measurements of binary black hole mergers at cosmological distances.

Fast localization of coalescing binaries with a heterogeneous network of advanced gravitational wave detectors

We present the expected performance regarding fast sky localization of coalescing binaries with a network of three gravitational wave detectors having heterogeneous sensitivities, such as the LIGO-Virgo network. A hierarchical approach can be used in order to make an effective use of information from the least sensitive detector. In this approach, the presence of an event seen in coincidence in the two more sensitive detectors triggers a focused search in the data of the third, less sensitive, detector with a lower signal-to-noise ratio (SNR) threshold. We investigate the benefit for sky localization that can be expected in such an approach, using simulated data and signals. We find that as the sensitivity of Virgo approaches one third of the LIGO sensitivity, the accuracy and precision of the localization can be improved by about a factor 3 when Virgo data are searched with an SNR threshold around 3.5.

Digging the population of compact binary mergers out of the noise

Coalescing compact binaries emitting gravitational wave (GW) signals, as recently detected by the Advanced LIGO-Virgo network, constitute a population over the multi-dimensional space of component masses and spins, redshift, and other parameters. Characterizing this population is a major goal of GW observations and may be approached via parametric models. We demonstrate hierarchical inference for such models with a method that accounts for uncertainties in each binary merger's individual parameters, for mass-dependent selection effects, and also for the presence of a second population of candidate events caused by detector noise. Thus, the method is robust to potential biases from a contaminated sample and allows us to extract information from events that have a relatively small probability of astrophysical origin.

Discovering intermediate-mass black hole lenses through gravitational wave lensing

Intermediate-mass black holes are the missing link that connects stellar-mass to supermassive black holes and are key to understanding galaxy evolution. Gravitational waves, like photons, can be lensed, leading to discernable effects such as diffraction or repeated signals. We investigate the detectability of intermediate-mass black hole deflectors in the LIGO-Virgo detector network. In particular, we simulate gravitational waves with variable source distributions lensed by an astrophysical population of intermediate-mass black holes, and use standard LIGO tools to infer the properties of these lenses. We find detections of intermediate-mass black holes at $98\%$ confidence level over a wide range of binary and lens parameters. Therefore, we conclude that intermediate-mass black holes could be detected through lensing of gravitational waves in the LIGO-Virgo detector network.

The Enigmatic Spin Evolution of PSR J0537–6910: r-modes, Gravitational Waves, and the Case for Continued Timing

Abstract We discuss the unique spin evolution of the young X-ray pulsar PSR J0537–6910, a system in which the regular spin down is interrupted by glitches every few months. Drawing on the complete timing data from the Rossi X-ray Timing Explorer (from 1999 to 2011), we argue that a trend in the interglitch behavior points to an effective braking index close to n = 7, which is much larger than expected. This value is interesting because it would accord with the neutron star spinning down due to gravitational waves from an unstable r-mode. We discuss to what extent this, admittedly speculative, scenario may be consistent and if the associated gravitational-wave signal would be within reach of ground-based detectors. Our estimates suggest that one may, indeed, be able to use future observations to test the idea. Further precision timing would help to enhance the achievable sensitivity, and we advocate a joint observing campaign between the Neutron Star Interior Composition Explorer and the LIGO-Virgo network.