LIGO Scientific Collaboration Research Articles

Constraints on compact binary merger evolution from spin-orbit misalignment in gravitational-wave observations

ABSTRACT The identification of the first confirmed neutron star–black hole (NS-BH) binary mergers by the LIGO, Virgo, and KAGRA collaboration provides the opportunity to investigate the properties of the early sample of confirmed and candidate events. Here, we focus primarily on the tilt angle of the BH’s spin relative to the orbital angular momentum vector of the binary, and the implications for the physical processes that determine this tilt. The posterior tilt distributions of GW200115 and the candidate events GW190426_152155 and GW190917_114630 peak at significantly anti-aligned orientations (though display wide distributions). Producing these tilts through isolated binary evolution would require stronger natal kicks than are typically considered (and preferentially polar kicks would be ruled out), and/or an additional source of tilt such as stable mass transfer. The early sample of NS-BH events are less massive than expected for classical formation channels, and may provide evidence for efficient mass transfer that results in the merger of more massive NS-BH binaries before their evolution to the compact phase is complete. We predict that future gravitational-wave detections of NS-BH events will continue to display total binary masses of ≈7 M⊙ and mass ratios of q ∼ 3 if this interpretation is correct. Conversely, the high mass of the candidate GW191219_163120 suggests a dynamical capture origin. Large tilts in a significant fraction of merging NS-BH systems would weaken the prospects for electromagnetic detection. However, EM observations, including non-detections, can significantly tighten the constraints on spin and mass ratio.

Primordial black hole production during first-order phase transitions

Primordial black holes (PBHs) produced in the early Universe have attracted wide interest for their ability to constitute dark matter and explain the compact binary coalescence. We propose a new mechanism of PBH production during first-order phase transitions (PTs) and find that PBHs are naturally produced during PTs model-independently. Because of the randomness of the quantum tunneling, there always exists some probability that the vacuum decay is postponed in a whole Hubble volume. Since the vacuum energy density remains constant while radiation is quickly redshifted in the expanding Universe, the postponed vacuum decay then results in overdense regions, which finally collapse into PBHs as indicated by numerical simulations. Utilizing this result one can obtain mutual predictions and constraints between PBHs and GWs from PTs. The predicted mass function of PBHs is nearly monochromatic. We investigate two typical cases and find that 1) PBHs from a PT constitute all dark matter and GWs peak at $1$Hz, 2) PBHs from a PT can explain the coalescence events observed by LIGO-Virgo collaboration, and meanwhile GWs can explain the common-spectrum process detected by NANOGrav collaboration.

Dark matter component in hadronic models with short-range correlations

A relativistic mean-field hadronic model with a dark matter (DM) particle coupled to nucleons including short-range correlations (SRC) is applied to study neutron stars (NS). The lightest neutralino is chosen as the dark particle candidate, which interacts with nucleons by the exchange of Higgs bosons. A detailed thermodynamical analysis shows that the contribution of the DM fermions to the energy density of the matter composed by these particles and nucleons is completely dominated by the DM kinetic terms. The model reproduces satisfactorily the constraints on the mass-radius diagram obtained from the analysis of the combined data from the NICER mission, LIGO collaboration, and mass measurements from radio observations. We show that the SRC balance the reduction of the neutron star mass due to the DM component, and because of that the model is able to present more massive NS. We also present a study of the effect, in the NS mass-radius profiles, of the uncertainties in some bulk parameters related to the hadronic sector. We find that it is possible to generate parametrizations, with DM content, compatible with the recent astrophysical constraints and with the uncertainty in the symmetry energy slope obtained from the results reported by the updated lead radius experiment (PREX-2).

On the Gravitational Waves Detected by the LIGO-Virgo Collaboration

On the Gravitational Waves Detected by the LIGO-Virgo Collaboration

On the Angular Momentum Transport Efficiency within the Star Constrained from Gravitational-wave Observations

Abstract The LIGO Scientific Collaboration and Virgo Collaboration (LIGO/Virgo) have recently reported in GWTC-2.1 eight additional candidate events with a probability of astrophysical origin greater than 0.5 in the LIGO/Virgo’s deeper search on O3a running. In GWTC-2.1, the majority of the effective inspiral spins (χ eff) show magnitudes consistent with zero, while two (GW190403-051519 and GW190805-211137) of the eight new events have χ eff > 0 (at 90% credibility). We note that GW190403-051519 was reported with χ eff = 0.70 − 0.27 + 0.15 and mass ratio q = 0.25 − 0.11 + 0.54 . Assuming a uniform prior probability between 0 and 1 for each black hole’s dimensionless spin magnitude, GW190403-051519 was reported with the dimensionless spin of the more massive black hole, χ 1 = 0.92 − 0.22 + 0.07 . This is the fastest black hole ever measured in all current gravitational-wave events. If the immediate progenitor of GW190403-051519 is a close binary system composed of a black hole and a helium star, which can be the natural outcome of the classical isolated binary evolution through the common envelope phase, this extremely high spin challenges, at least in that case, the existence of an efficient angular momentum transport mechanism between the stellar core and the radiative envelope of massive stars, as for instance predicted by the Tayler–Spruit dynamo or its revised version by Fuller et al.

Research contribution and impact of Indian institutes of science education and research (iisers) to physical sciences

The present study explores the research contribution of five selected Indian Institutes of Science Education and Research (IISERs) in the discipline of physical sciences over the last 15 years period, 2006–2020. The selected five IISERs published a total of 3,354 research articles which comprise 47.34% share of internationally co-authored articles. The publications are interpreted in terms of chronological growth, collaboration trend, focus sub-areas, scholarly communication channels, leading collaborating countries, keywords and citation impact. Among the selected IISERs, IISER-P followed by IISER-K and IISER-M were the major producers of research papers. Since 2016, a sharp rising trend in terms of publication output and international collaborative efforts was observed. Furthermore, the collaboration trend also witnessedthat there were a strong international collaboration linkages among IISERS and developed countries like the USA, Germany and Italy were the primary collaborating countries.In addition, the international experimental research consortia like CMS, LIGO, VIRGO and BELLE Collaboration played a crucial role in global cooperative research efforts of IISERs.

Testing General Relativity with Gravitational Waves: An Overview

The detections of gravitational-wave (GW) signals from compact binary coalescence by ground-based detectors have opened up the era of GW astronomy. These observations provide opportunities to test Einstein’s general theory of relativity at the strong-field regime. Here we give a brief overview of the various GW-based tests of General Relativity (GR) performed by the LIGO-Virgo collaboration on the detected GW events to date. After providing details for the tests performed in four categories, we discuss the prospects for each test in the context of future GW detectors. The four categories of tests include the consistency tests, parametrized tests for GW generation and propagation, tests for the merger remnant properties, and GW polarization tests.

The search of higher multipole radiation in gravitational waves from compact binary coalescences by a minimally-modelled pipeline

The coherent WaveBurst (cWB) pipeline implements a minimally-modelled search to find a coherent response in the network of gravitational wave detectors of the LIGO-Virgo Col-laboration in the time-frequency domain. In this manuscript, we provide a timely introduction to an extension of the cWB analysis to detect spectral features beyond the main quadrupolar emission of gravitational waves during the inspiral phase of compact binary coalescences; more detailed discussion will be provided in a forthcoming paper [1]. The search is performed by defining specific regions in the time-frequency map to extract the energy of harmonics of main quadrupole mode in the inspiral phase. This method has already been used in the GW190814 discovery paper (Astrophys. J. Lett. 896 L44). Here we show the procedure to detect the (3, 3) multipole in GW190814 within the cWB framework.

Evidence for lensing of gravitational waves from LIGO-Virgo data

Recently, the LIGO-Virgo Collaboration (LVC) concluded that there is no evidence for lensed gravitational waves (GW) in the first half of the O3 run, claiming "We find the observation of lensed events to be unlikely, with the fractional rate at $\mu>2$ being $3.3\times 10^{-4}$". While we agree that the chance of an individual GW event being lensed at $\mu>2$ is smaller than $10^{-3}$, the number of observed events depends on the product of this small probability times the rate of mergers at high redshift. Observational constraints from the stochastic GW background indicate that the rate of conventional mass BBH mergers (8 < M (M$_{\odot}$) < 15) in the redshift range 1<z< 2 could be as high as O($10^7$) events per year, more than sufficient to compensate for the intrinsically low probability of lensing. To reach the LVC trigger threshold these events require high magnification, but would still produce up to 10 to 30 LVC observable events per year. Thus, all the LVC observed ordinary stellar mass BBH mergers from this epoch must be strongly lensed. By adopting low-rates at high redshift, LVC assumes that lensed events can not be taking place, thus incorrectly assigning them a closer distance and higher masses by a factor of a few (typically 2 to 5). The LVC adopted priors on time delay are in tension with the distribution of observed time delays in lensed quasars. Pairs of events like GW190421-GW190910 and GW190424-GW190910, which are directly assigned a probability of zero by LVC, should be instead considered as prime candidates to be strongly lensed GW pairs, since their separation in time is consistent with observations of time delays in lensed quasars. Correcting for the LVC wrong Bayesian priors, maximum merger rate of conventional mass BBH in 1<z<2, and gravitational lensing time-delay model, reverses the LVC conclusions and supports the strong gravitational lensing hypothesis.

Constraining Scalar-tensor Theories Using Neutron Star–Black Hole Gravitational Wave Events

Abstract With the continuous upgrade of detectors, greater numbers of gravitational wave (GW) events have been captured by the LIGO Scientific Collaboration and Virgo Collaboration (LVC), which offer a new avenue to test general relativity and explore the nature of gravity. Although various model-independent tests have been performed by LVC in previous works, it is still interesting to ask what constraints can be placed on specific models by current GW observations. In this work, we focus on three models of scalar-tensor theories, the Brans–Dicke theory (BD), the theory with scalarization phenomena proposed by Damour and Esposito-Farèse (DEF), and screened modified gravity (SMG). Of the four possible neutron star–black hole events that have occurred so far, we use two of them to place constraints. The other two are excluded in this work because of possible unphysical deviations. We consider the inspiral range with the cutoff frequency at the innermost stable circular orbit and add a modification of dipole radiation into the waveform template. The scalar charges of neutron stars in the dipole term are derived by solving the Tolman–Oppenheimer–Volkoff equations for different equations of state. The constraints are obtained by performing the full Bayesian inference with the help of the open source software Bilby. The results show that the constraints given by GWs are comparable to those given by pulsar timing experiments for DEF theory, but are not competitive with the current solar system constraints for BD and SMG theories.

Constraints from LIGO O3 Data on Gravitational-wave Emission Due to R-modes in the Glitching Pulsar PSR J0537–6910

Abstract We present a search for continuous gravitational-wave emission due to r-modes in the pulsar PSR J0537–6910 using data from the LIGO–Virgo Collaboration observing run O3. PSR J0537–6910 is a young energetic X-ray pulsar and is the most frequent glitcher known. The inter-glitch braking index of the pulsar suggests that gravitational-wave emission due to r-mode oscillations may play an important role in the spin evolution of this pulsar. Theoretical models confirm this possibility and predict emission at a level that can be probed by ground-based detectors. In order to explore this scenario, we search for r-mode emission in the epochs between glitches by using a contemporaneous timing ephemeris obtained from NICER data. We do not detect any signals in the theoretically expected band of 86–97 Hz, and report upper limits on the amplitude of the gravitational waves. Our results improve on previous amplitude upper limits from r-modes in J0537-6910 by a factor of up to 3 and place stringent constraints on theoretical models for r-mode-driven spin-down in PSR J0537–6910, especially for higher frequencies at which our results reach below the spin-down limit defined by energy conservation.

The Implications of High Black Hole Spins for the Origin of Binary Black Hole Mergers

Abstract The LIGO–Virgo collaboration has reported 50 black hole–black hole (BH–BH) mergers and 8 candidates recovered from digging deeper into the detector noise. The majority of these mergers have low effective spins pointing toward low BH spins and efficient angular momentum (AM) transport in massive stars as proposed by several models (e.g., the Tayler–Spruit dynamo). However, out of these 58 mergers, 7 are consistent with having high effective-spin parameter (χ eff &gt; 0.3). Additionally, two events seem to have high effective spins sourced from the spin of the primary (more massive) BH. These particular observations could be used to discriminate between the isolated binary and dynamical formation channels. It might seem that high BH spins point to a dynamical origin if AM in stars is efficient and forms low-spinning BHs. In such a case dynamical formation is required to produce second and third generations of BH–BH mergers with typically high spinning BHs. Here we show, however, that isolated binary BH–BH formation naturally reproduces such highly spinning BHs. Our models start with efficient AM in massive stars that is needed to reproduce the majority of BH–BH mergers with low effective spins. Later, some of the binaries are subject to a tidal spin-up allowing the formation of a moderate fraction (∼10%) of BH–BH mergers with high effective spins (χ eff ≳ 0.4–0.5). In addition, isolated binary evolution can produce a small fraction of BH–BH mergers with almost maximally spinning primary BHs. Therefore, the formation scenario of these atypical BH–BH mergers remains to be found.

Detecting gravitational waves with disparate detector responses: Two new binary black hole mergers

We introduce a new technique to search for gravitational wave events from compact binary mergers that produce a clear signal only in a single gravitational wave detector, and marginal signals in other detectors. Such a situation can arise when the detectors in a network have different sensitivities, or when sources have unfavorable sky locations or orientations. We start with a short list of loud single-detector triggers from regions of parameter space that are empirically unaffected by glitches (after applying signal-quality vetoes). For each of these triggers, we compute evidence for astrophysical origin from the rest of the detector network by coherently combining the likelihoods from all detectors and marginalizing over extrinsic geometric parameters. We report the discovery of two new binary black hole (BBH) mergers in the second observing run of Advanced LIGO and Virgo (O2), in addition to the ones that were reported in [B. P. Abbott et al. (LIGO Scientific and Virgo Collaborations), Phys. Rev. X 9, 031040 (2019) and [T. Venumadhav et al., Phys. Rev. D 101, 083030 (2020)]. We estimate that the two events have false alarm rates of one in 19 years (60 O2) and one in 11 years (36 O2). One of the events, GW170817A, has primary and secondary masses ${m}_{1}^{\mathrm{src}}=5{6}_{\ensuremath{-}10}^{+16}\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$ and ${m}_{2}^{\mathrm{src}}=4{0}_{\ensuremath{-}11}^{+10}\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$ in the source frame. The existence of GW170817A should be very informative about the theoretically predicted upper mass gap for stellar mass black holes. Its effective spin parameter is measured to be ${\ensuremath{\chi}}_{\mathrm{eff}}=0.5\ifmmode\pm\else\textpm\fi{}0.2$, which is consistent with the tendency of the heavier detected BBH systems to have large and positive effective spin parameters. The other event, GWC170402, will be discussed thoroughly in future work.

Search for Neutrinos in Coincidence with Gravitational Wave Events from the LIGO–Virgo O3a Observing Run with the Super-Kamiokande Detector

Abstract The Super-Kamiokande detector can be used to search for neutrinos in time coincidence with gravitational waves detected by the LIGO–Virgo Collaboration (LVC). Both low-energy (7–100 MeV) and high-energy (0.1–105 GeV) samples were analyzed in order to cover a very wide neutrino spectrum. Follow-ups of 36 (out of 39) gravitational waves reported in the GWTC-2 catalog were examined; no significant excess above the background was observed, with 10 (24) observed neutrinos compared with 4.8 (25.0) expected events in the high-energy (low-energy) samples. A statistical approach was used to compute the significance of potential coincidences. For each observation, p-values were estimated using neutrino direction and LVC sky map; the most significant event (GW190602_175927) is associated with a post-trial p-value of 7.8% (1.4σ). Additionally, flux limits were computed independently for each sample and by combining the samples. The energy emitted as neutrinos by the identified gravitational wave sources was constrained, both for given flavors and for all flavors assuming equipartition between the different flavors, independently for each trigger and by combining sources of the same nature.

Hidden Markov model tracking of continuous gravitational waves from a binary neutron star with wandering spin. III. Rotational phase tracking

A hidden Markov model (HMM) solved recursively by the Viterbi algorithm can be configured to search for persistent, quasimonochromatic gravitational radiation from an isolated or accreting neutron star, whose rotational frequency is unknown and wanders stochastically. Here an existing HMM analysis pipeline is generalized to track rotational phase and frequency simultaneously, by modeling the intra-step rotational evolution according to a phase-wrapped Ornstein-Uhlenbeck process, and by calculating the emission probability using a phase-sensitive version of the Bayesian matched filter known as the $\mathcal{B}$-statistic. The generalized algorithm tracks signals from isolated and binary sources with characteristic wave strain $h_0 \geq 1.3\times 10^{-26}$ in Gaussian noise with amplitude spectral density $4\times 10^{-24}\,{\rm Hz^{-1/2}}$, for a simulated observation composed of $N_T=37$ data segments, each $T_{\rm drift}=10\,{\rm days}$ long, the typical duration of a search for the low-mass X-ray binary (LMXB) Sco X$-$1 with the Laser Interferometer Gravitational Wave Observatory (LIGO). It is equally sensitive to isolated and binary sources and $\approx 1.5$ times more sensitive than the previous pipeline. Receiver operating characteristic curves and errors in the recovered parameters are presented for a range of practical $h_0$ and $N_T$ values. The generalized algorithm successfully detects every available synthetic signal in Stage I of the Sco X$-$1 Mock Data Challenge convened by the LIGO Scientific Collaboration, recovering the frequency and orbital semimajor axis with accuracies of better than $9.5\times 10^{-7}\,{\rm Hz}$ and $1.6\times 10^{-3}\,{\rm lt\,s}$ respectively. The Viterbi solver runs in $\approx 2\times 10^3$ CPU-hr for an isolated source and $\sim 10^5$ CPU-hr for a LMXB source in a typical, broadband ($0.5$-${\rm kHz}$) search.

Holographic bound on area of a compact binary merger remnant

Using concomitantly the generalized second law of black hole thermodynamics and the holographic Bekenstein entropy bound embellished by loop quantum gravity corrections to quantum black hole entropy, we show that the boundary cross sectional area of the postmerger remnant formed from the compact binary merger in gravitational wave detection experiments like GW150914 et seq., by the LIGO-VIRGO collaboration, is bounded from below. This lower bound is more general than the bound obtained from application of Hawking's classical area theorem for black holes, since it does not depend on whether the inspiralling compact binary pair or the postmerger remnant consists of black holes or other exotic compact objects. The derivation of the bound entails an estimate of the entropy of the gravitational waves emitted during the binary merger, which adapts to gravitational waves an extant formalism proposed originally for particle ensembles. The results for the minimal cross sectional area of the merger remnant due to binary compact mergers observed recently by the LIGO-VIRGO collaboration are discussed. While accurate measurement of the mass of the remnant for the binary neutron star merger GW170817 remains a challenge, we provide a proof of principle that for binary neutron star mergers our lower bound on the cross sectional area of the remnant provides an alternative approach to probe the validity of neutron star equations of state, independent of the tidal deformabilities of the components.

Filling the black hole mass gap: Avoiding pair instability in massive stars through addition of nonnuclear energy

In standard stellar evolution, stars with masses ranging from approximately $150\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$ to $240\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$ are expected to evolve to a pair-instability supernova with no black hole (BH) remnant. This evolutionary behavior leads to a predicted gap in the black hole mass function from approximately $50\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$ to $140\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$. Yet the LIGO and Virgo collaborations [Phys. Rev. Lett. 125, 101102 (2020)] recently discovered black holes with masses of $66\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$ and $85\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$ in the gravitational-wave event GW190521. We propose a new method to populate the BH mass gap. If an energy source is added throughout the star in addition to nuclear fusion, it is possible for the altered evolution to avoid the complete destruction of a pair-instability supernova, and instead a BH remnant is left behind. An example of an extra energy source is dark matter annihilation within the star, but our results hold more generally. We show this phenomenon by exploring the effect of adding an energy source independent of temperature and density to a $180\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$ star, using the mesa one-dimensional stellar evolution software. $\text{If}\ensuremath{\sim}50%$ of the star's energy is due to this new source, the star is capable of avoiding the pair instability entirely, reaching a presupernova mass of $119\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$ before evolving into a BH remnant in the mass gap.

Improved gravitational-wave constraints on higher-order curvature theories of gravity

Gravitational wave observations of compact binaries allow us to test general relativity (and modifications thereof) in the strong and highly dynamical field regime of gravity. Here, we confront two extensions to general relativity, dynamical Chern-Simons, and Einstein-dilaton-Gauss-Bonnet theories, against the gravitational wave sources from the GWTC-1 and GWTC-2 catalogs by the LIGO-Virgo Collaboration. By stacking the posterior of individual events, we strengthen the constraint on the square root of the coupling parameter in Einstein-dilaton-Gauss-Bonnet gravity to $\sqrt{{\ensuremath{\alpha}}_{\mathrm{EdGB}}}&lt;1.7\text{ }\text{ }\mathrm{km}$, but we are unable to place meaningful constraints on dynamical Chern-Simons gravity. Importantly, we also show that our bounds are robust to (i) the choice of general-relativity base waveform model, upon which we add modifications, (ii) unknown higher post-Newtonian order terms in the modifications to general relativity, (iii) the small-coupling approximation, and (iv) uncertainties on the nature of the constituent compact objects.

Anomaly detection in gravitational waves data using convolutional autoencoders

As of this moment, 50 gravitational wave (GW) detections have been announced, thanks to the observational efforts of the LIGO-Virgo collaboration, working with the Advanced LIGO and the Advanced Virgo interferometers. The detection of signals is complicated by the noise-dominated nature of the data. Conventional approaches in GW detection procedures require either precise knowledge of the GW waveform in the context of matched filtering searches or coincident analysis of data from multiple detectors. Furthermore, the analysis is prone to contamination by instrumental or environmental artifacts called glitches which either mimic astrophysical signals or reduce the overall quality of data. In this paper, we propose an alternative generic method of studying GW data based on detecting anomalies. The anomalies we study are transient signals, different from the slow non-stationary noise of the detector. The anomalies presented in the manuscript are mostly based on the GW emitted by the mergers of binary black hole systems. However, the presented study of anomalies is not limited only to GW alone, but also includes glitches occurring in the real LIGO/Virgo dataset available at the Gravitational Waves Open Science Center. To search for anomalies we employ deep learning algorithms, namely convolutional autoencoders, which are trained on both simulated and real detector data. We demonstrate the capabilities of our deep learning implementation in the reconstruction of injected signals. We study the influence of the GW strength, defined in terms of matched filter signal-to-noise ratio, on the detection of anomalies. Moreover, we present the application of our method for the localization in time of anomalies in the studied time-series data. We validate the results of anomaly searches on real data containing confirmed gravitational wave detections; we thus prove the generalization capabilities of our method, towards detecting GWs unknown to our deep learning models during training.

Exploring gravitational-wave detection and parameter inference using deep learning methods

We explore machine learning methods to detect gravitational waves (GW) from binary black hole (BBH) mergers using deep learning (DL) algorithms. The DL networks are trained with gravitational waveforms obtained from BBH mergers with component masses randomly sampled in the range from 5 to 100 solar masses and luminosity distances from 100 Mpc to, at least, 2000 Mpc. The GW signal waveforms are injected in public data from the O2 run of the Advanced LIGO and Advanced Virgo detectors, in time windows that do not coincide with those of known detected signals. We demonstrate that DL algorithms, trained with GW signal waveforms at distances of 2000 Mpc, still show high accuracy when detecting closer signals, within the ranges considered in our analysis. Moreover, by combining the results of the three-detector network in a unique RGB image, the single detector performance is improved by as much as 70%. Furthermore, we train a regression network to perform parameter inference on BBH spectrogram data and apply this network to the events from the GWTC-1 and GWTC-2 catalogs. Without significant optimization of our algorithms we obtain results that are mostly consistent with published results by the LIGO–Virgo Collaboration. In particular, our predictions for the chirp mass are compatible (up to 3σ) with the official values for 90% of events. From these results we conclude that the combination of computer vision techniques and deep-learning methods put forward in this work is a worthy addition to the GW astronomer’s toolbox.