Search

Dynamical interactions in dense star clusters could significantly influence the properties of black holes, leaving imprints on their gravitational-wave signatures. While previous studies have mostly focused on repeated black hole mergers for spin and mass growth, this work examines the impact of physical collisions and close encounters between black holes and (non-compact) stars. Using Monte Carlo N-body models of dense star clusters, we find that a large fraction of black holes retained upon formation undergo collisions with stars. Within our explored cluster models, the proportion of binary black hole mergers affected by stellar collisions ranges from $10\%$ to $60\%$. If all stellar-mass black holes are initially non-spinning, we find that up to $40\%$ of merging binary black holes may have components with dimensionless spin parameter $\chi\gtrsim 0.2$ because of prior stellar collisions, while typically about $10\%$ have spins near $\chi = 0.7$ from prior black hole mergers. We demonstrate that young star clusters are especially important environments as they can produce collisions of black holes with very massive stars, allowing significant spin up of the black holes through accretion. Our predictions for black hole spin distributions from these stellar collisions highlight their sensitivity to accretion efficiency, underscoring the need for detailed hydrodynamic calculations to better understand the accretion physics following these interactions., Comment: 16 pages, 8 figures, 2 tables. Submitted to ApJ

On the 29th of May 2023, the LIGO-Virgo-KAGRA Collaboration observed a compact binary coalescence event consistent with a neutron star-black hole merger, though the heavier object of mass 2.5-4.5 $M_\odot$ would fall into the purported lower mass gap. An alternative explanation for apparent observations of events in this mass range has been suggested as strongly gravitationally lensed binary neutron stars. In this scenario, magnification would lead to the source appearing closer and heavier than it really is. Here, we investigate the chances and possible consequences for the GW230529 event to be gravitationally lensed. We find this would require high magnifications and we obtain low rates for observing such an event, with a relative fraction of lensed versus unlensed observed events of $2 \times 10^{-3}$ at most. When comparing the lensed and unlensed hypotheses accounting for the latest rates and population model, we find a 1/58 chance of lensing, disfavoring this option. Moreover, when the magnification is assumed to be strong enough to bring the mass of the heavier binary component below the standard limits on neutron star masses, we find high probability for the lighter object to have a sub-solar mass, making the binary even more exotic than a mass-gap neutron star-black hole system. Even when the secondary is not sub-solar, its tidal deformability would likely be measurable, which is not the case for GW230529. Finally, we do not find evidence for extra lensing signatures such as the arrival of additional lensed images, type-II image dephasing, or microlensing. Therefore, we conclude it is unlikely for GW230529 to be a strongly gravitationally lensed binary neutron star signal., Comment: 15 pages, 11 figures

On 29 May 2023, the LIGO Livingston observatory detected the gravitational-wave signal GW230529_181500 from the merger of a neutron star with a lower mass-gap compact object. Its long inspiral signal provides a unique opportunity to test General Relativity (GR) in a parameter space previously unexplored by strong-field tests. In this work, we performed parameterized inspiral tests of GR with GW230529_181500. Specifically, we search for deviations in the frequency-domain GW phase by allowing for agnostic corrections to the post-Newtonian coefficients. We performed tests with the Flexible Theory Independent (FTI) and Test Infrastructure for General Relativity (TIGER) frameworks using several quasi-circular waveform models that capture different physical effects (higher modes, spins, tides). We find that the signal is consistent with GR for all deviation parameters. Assuming the primary object is a black hole, we obtain particularly tight constraints on the dipole radiation at $-1$PN order of $|\delta\hat{\varphi}_{-2}| \lesssim 8 \times 10^{-5}$, which is a factor $\sim17$ times more stringent than previous bounds from the neutron star--black hole merger GW200115_042309, as well as on the 0.5PN and 1PN deviation parameters. We discuss some challenges that arise when analyzing this signal, namely biases due to correlations with tidal effects and the degeneracy between the 0PN deviation parameter and the chirp mass. To illustrate the importance of GW230529_181500 for tests of GR, we mapped the agnostic $-1$PN results to a class of Einstein-scalar-Gauss-Bonnet (ESGB) theories of gravity. We also conducted an analysis probing the specific phase deviation expected in ESGB theory and obtain an upper bound on the Gauss-Bonnet coupling of $\ell_{\rm GB} \lesssim 0.51~\rm{M}_\odot$ ($\sqrt{\alpha_{\rm GB}} \lesssim 0.28$ km), which is better than any previously reported constraint., Comment: 19 pages, 9 figures

GW230529 is the first compact binary coalescence detected by the LIGO-Virgo-KAGRA collaboration with at least one component mass confidently in the lower mass-gap, corresponding to the range 3-5$M_{\odot}$. If interpreted as a neutron star-black hole merger, this event has the most symmetric mass ratio detected so far and therefore has a relatively high probability of producing electromagnetic (EM) emission. However, no EM counterpart has been reported. At the merger time $t_0$, Swift-BAT and Fermi-GBM together covered 100$\%$ of the sky. Performing a targeted search in a time window $[t_0-20 \text{s},t_0+20 \text{s}]$, we report no detection by the Swift-BAT and the Fermi-GBM instruments. Combining the position-dependent $\gamma-$ray flux upper limits and the gravitational-wave posterior distribution of luminosity distance, sky localization and inclination angle of the binary, we derive constraints on the characteristic luminosity and structure of the jet possibly launched during the merger. Assuming a top-hat jet structure, we exclude at 90$\%$ credibility the presence of a jet which has at the same time an on-axis isotropic luminosity $\gtrsim 10^{48}$ erg s$^{-1}$, in the bolometric band 1 keV-10 MeV, and a jet opening angle $\gtrsim 15$ deg. Similar constraints are derived testing other assumptions about the jet structure profile. Excluding GRB 170817A, the luminosity upper limits derived here are below the luminosity of any GRB observed so far., Comment: 18 pages, 1 table, 11 figures

We describe an implementation of the relative binning technique to speed up parameter estimation of gravitational-wave signals. We first give a pedagogical overview of relative binning, discussing also the expressions for the likelihood marginalized over phase and distance. Then, we describe the details of the code in \texttt{Bilby}, an open-source software package commonly used for parameter estimation of gravitational-wave sources. Our code is able to reproduce the parameters of GW170817 in 14 hours on a single-core CPU, performs well on simulated signals, and passes the percentile-percentile (p-p) tests. We also illustrate that relative binning is an ideal technique to estimate the parameters of signals in next-generation gravitational wave detectors., Comment: 19 pages, 11 figures

The astrophysical formation channels of binary black hole systems predict correlations between their mass, spin, and redshift distributions, which can be probed with gravitational-wave observations. Population-level analysis of the latest LIGO-Virgo-KAGRA catalog of binary black hole mergers has identified evidence for such correlations assuming linear evolution of the mean and width of the effective spin distribution as a function of the binary mass ratio and merger redshift. However, the complex astrophysical processes at play in compact binary formation do not necessarily predict linear relationships between the distributions of these parameters. In this work, we relax the assumption of linearity and instead search for correlations using a more flexible cubic spline model. Our results suggest a nonlinear correlation between the width of the effective spin distribution and redshift. We also show that the LIGO-Virgo-Kagra collaborations may find convincing Bayesian evidence for nonlinear correlations by the end of the fourth observing run, O4. This highlights the valuable role of flexible models in population analyses of compact-object binaries in the era of growing catalogs., Comment: 16 pages, 11 figures

Gravitational-wave observations by the Laser Interferometer Gravitational-Wave Observatory (LIGO) and Virgo have provided us a new tool to explore the Universe on all scales from nuclear physics to the cosmos and have the massive potential to further impact fundamental physics, astrophysics, and cosmology for decades to come. In this paper we have studied the science capabilities of a network of LIGO detectors when they reach their best possible sensitivity, called A#, given the infrastructure in which they exist and a new generation of observatories that are factor of 10 to 100 times more sensitive (depending on the frequency), in particular a pair of L-shaped Cosmic Explorer observatories (one 40 km and one 20 km arm length) in the US and the triangular Einstein Telescope with 10 km arms in Europe. The presence of one or two A# observatories in a network containing two or one next generation observatories, respectively, will provide good localization capabilities for facilitating multimessenger astronomy and precision measurement of the Hubble parameter. Two Cosmic Explorer observatories are indispensable for achieving precise localization of binary neutron star events, facilitating detection of electromagnetic counterparts and transforming multimessenger astronomy. Their combined operation is even more important in the detection and localization of high-redshift sources, such as binary neutron stars, beyond the star-formation peak, and primordial black hole mergers, which may occur roughly 100 million years after the Big Bang. The addition of the Einstein Telescope to a network of two Cosmic Explorer observatories is critical for accomplishing all the identified science metrics. For most metrics the triple network of next generation terrestrial observatories are a factor 100 better than what can be accomplished by a network of three A# observatories., Comment: 48 pages, 20 figures, 14 tables

Compact-object binary mergers consisting of one neutron star and one black hole (NSBHs) have long been considered promising progenitors for gamma-ray bursts, whose central engine remains poorly understood. Using gravitational-wave constraints on the population-level NSBH mass and spin distributions we find that at most $20~\mathrm{Gpc}^{-3}\mathrm{yr}^{-1}$ of gamma-ray bursts in the local universe can have NSBH progenitors., Comment: Submitted to RNAAS

Gravitational-wave astronomy has revolutionized humanity's view of the universe, a revolution driven by observations that no other field can make. This white paper describes an observatory that builds on decades of investment by the National Science Foundation and that will drive discovery for decades to come: Cosmic Explorer. Major discoveries in astronomy are driven by three related improvements: better sensitivity, higher precision, and opening new observational windows. Cosmic Explorer promises all three and will deliver an order-of-magnitude greater sensitivity than LIGO. Cosmic Explorer will push the gravitational-wave frontier to almost the edge of the observable universe using technologies that have been proven by LIGO during its development. With the unprecedented sensitivity that only a new facility can deliver, Cosmic Explorer will make discoveries that cannot yet be anticipated, especially since gravitational waves are both synergistic with electromagnetic observations and can reach into regions of the universe that electromagnetic observations cannot explore. With Cosmic Explorer, scientists can use the universe as a laboratory to test the laws of physics and study the nature of matter. Cosmic Explorer allows the United States to continue its leading role in gravitational-wave science and the international network of next-generation observatories. With its extraordinary discovery potential, Cosmic Explorer will deliver revolutionary observations across astronomy, physics, and cosmology including: Black Holes and Neutron Stars Throughout Cosmic Time, Multi-Messenger Astrophysics and Dynamics of Dense Matter, New Probes of Extreme Astrophysics, Fundamental Physics and Precision Cosmology, Dark Matter and the Early Universe.

Along their path from source to observer, gravitational waves may be gravitationally lensed by massive objects. This results in distortions of the observed signal which can be used to extract new information about fundamental physics, astrophysics, and cosmology. Searches for these distortions amongst the observed signals from the current detector network have already been carried out, though there have as yet been no confident detections. However, predictions of the observation rate of lensing suggest detection in the future is a realistic possibility. Therefore, preparations need to be made to thoroughly investigate the candidate lensed signals. In this work, we present some of the follow-up analyses and strategies that could be applied to assess the significance of such events and ascertain what information may be extracted about the lens-source system from such candidate signals by applying them to a number of O3 candidate events, even if these signals did not yield a high significance for any of the lensing hypotheses. For strongly-lensed candidates, we verify their significance using a background of simulated unlensed events and statistics computed from lensing catalogs. We also look for potential electromagnetic counterparts. In addition, we analyse in detail a candidate for a strongly-lensed sub-threshold counterpart that is identified by a new method. For microlensing candidates, we perform model selection using a number of lens models to investigate our ability to determine the mass density profile of the lens and constrain the lens parameters. We also look for millilensing signatures in one of the lensed candidates. Applying these additional analyses does not lead to any additional evidence for lensing in the candidates that have been examined. However, it does provide important insight into potential avenues to deal with high-significance candidates in future observations., Comment: 29 pages, 27 figures

The collection of gravitational waves (GWs) that are either too weak or too numerous to be individually resolved is commonly referred to as the gravitational-wave background (GWB). A confident detection and model-driven characterization of such a signal will provide invaluable information about the evolution of the Universe and the population of GW sources within it. We present a new, user-friendly Python--based package for gravitational-wave data analysis to search for an isotropic GWB in ground--based interferometer data. We employ cross-correlation spectra of GW detector pairs to construct an optimal estimator of the Gaussian and isotropic GWB, and Bayesian parameter estimation to constrain GWB models. The modularity and clarity of the code allow for both a shallow learning curve and flexibility in adjusting the analysis to one's own needs. We describe the individual modules which make up {\tt pygwb}, following the traditional steps of stochastic analyses carried out within the LIGO, Virgo, and KAGRA Collaboration. We then describe the built-in pipeline which combines the different modules and validate it with both mock data and real GW data from the O3 Advanced LIGO and Virgo observing run. We successfully recover all mock data injections and reproduce published results., Comment: 32 pages, 14 figures

While much has been learned about black holes by analyzing the latest LVK catalog, GWTC-3, a measurement of the astrophysical distribution of the black hole spin orientations remains elusive. This is usually probed by measuring the cosine of the tilt angle ($\cos\tau$) between each black hole spin and the orbital angular momentum, $\cos\tau=+1$ being perfect alignment. Abbott et al. has modeled the $\cos\tau$ distribution as a mixture of an isotropic component and a Gaussian component with mean fixed at +1 and width measured from the data. We want to verify if the data require the existence of such a peak at $\cos\tau=+1$. We use various alternative models for the astrophysical tilt distribution and measure their parameters using the LVK GWTC-3 catalog. We find that a) Augmenting the LVK model such that the mean $\mu$ of the Gaussian is not fixed at +1 returns results that strongly depend on priors. If we allow $\mu>+1$ then the resulting astrophysical $\cos\tau$ distribution peaks at +1 and looks linear, rather than Gaussian. If we constrain $-1\leq \mu\leq+1$ the Gaussian component peaks at $\mu=0.48^{+0.46}_{-0.99}$ (median and 90% symmetric credible interval). Two other 2-component mixture models yield $\cos\tau$ distributions that either have a broad peak centered at $0.19^{+0.22}_{-0.18}$ or a plateau that spans the range [-0.5, +1], without a clear peak at +1. b) All of the models we considered agree on the fact that there is no excess of black hole tilts at around -1. c) While yielding quite different posteriors, the models considered in this work have Bayesian evidences that are the same within error bars. We conclude that the current dataset is not sufficiently informative to draw any model-independent conclusions on the astrophysical distribution of spin tilts, except that there is no excess of spins with negatively aligned tilts., Comment: 6 pages body + 10 pages appendices. Following Ref suggestion, now includes differential rate plots per unit costilt. Data release on Zenodo linked from Data availability section. Note: the abstract in the metadata is shorter than the paper's due to Arxiv restriction to 1920 characters

Neutron star-black hole (NSBH) mergers detected in gravitational waves have the potential to shed light on supernova physics, the dense matter equation of state, and the astrophysical processes that power their potential electromagnetic counterparts. We use the population of four candidate NSBH events detected in gravitational waves so far with a false alarm rate $\leq 1~\mathrm{yr}^{-1}$ to constrain the mass and spin distributions and multimessenger prospects of these systems. We find that the black holes in NSBHs are both less massive and have smaller dimensionless spins than those in black hole binaries. We also find evidence for a mass gap between the most massive neutron stars and least massive black holes in NSBHs at 98.6% credibility. Using an approach driven by gravitational-wave data rather than binary simulations, we find that fewer than 14% of NSBH mergers detectable in gravitational waves will have an electromagnetic counterpart. While the inferred presence of a mass gap and fraction of sources with a counterpart depend on the event selection and prior knowledge of source classification, the conclusion that the black holes in NSBHs have lower masses and smaller spin parameters than those in black hole binaries is robust. Finally, we propose a method for the multimessenger analysis of NSBH mergers based on the nondetection of an electromagnetic counterpart and conclude that, even in the most optimistic case, the constraints on the neutron star equation of state that can be obtained with multimessenger NSBH detections are not competitive with those from gravitational-wave measurements of tides in binary neutron star mergers and radio and X-ray pulsar observations., Comment: 13 pages, 9 figures, Data release URL in paper; matches published version

Compact-object binaries including a white dwarf component are unique among gravitational-wave sources because their evolution is governed not just by general relativity and tides, but also by mass transfer. While the black hole and neutron star binaries observed with ground-based gravitational-wave detectors are driven to inspiral due to the emission of gravitational radiation--manifesting as a "chirp-like" gravitational-wave signal--the astrophysical processes at work in double white dwarf (DWD) systems can cause the inspiral to stall and even reverse into an outspiral. The dynamics of the DWD outspiral thus encode information about tides, which tell us about the behaviour of electron-degenerate matter. We carry out a population study to determine the effect of the strength of tides on the distributions of the DWD binary parameters that the Laser Interferometer Space Antenna (LISA) will be able to constrain. We find that the strength of tidal coupling parameterized via the tidal synchronization timescale at the onset of mass transfer affects the distribution of gravitational-wave frequencies and frequency derivatives for detectably mass-transferring DWD systems. Using a hierarchical Bayesian framework informed by binary population synthesis simulations, we demonstrate how this parameter can be inferred using LISA observations. By measuring the population properties of DWDs, LISA will be able to probe the behavior of electron-degenerate matter., Comment: version accepted in ApJ

The population-level distributions of the masses, spins, and redshifts of binary black holes (BBHs) observed using gravitational waves can shed light on how these systems form and evolve. Because of the complex astrophysical processes shaping the inferred BBH population, models allowing for correlations among these parameters will be necessary to fully characterize these sources. We hierarchically analyze the BBH population detected by LIGO and Virgo with a model allowing for correlations between the effective aligned spin and the primary mass and redshift. We find that the width of the effective spin distribution grows with redshift at 98.6% credibility. We determine this trend to be robust under the application of several alternative models and additionally verify that such a correlation is unlikely to be spuriously introduced using a simulated population. We discuss the possibility that this correlation could be due to a change in the natal black hole spin distribution with redshift., Comment: Matches published version

It is expected that the orbital planes of gravitational-wave (GW) sources are isotropically distributed. However, both physical and technical factors, such as alternate theories of gravity with birefringence, catalog contamination, and search algorithm limitations, could result in inferring a non-isotropic distribution. Showing that the inferred astrophysical distribution of the orbital orientations is indeed isotropic can thus be used to rule out some violations of general relativity, as a null test about the purity of the GW catalog sample, and as a check that selection effects are being properly accounted for. We augment the default mass/spins/redshift model used by the LIGO-Virgo-KAGRA Collaboration in their most recent analysis to also measure the astrophysical distribution of orbital orientations. We show that the 69 binary black holes in GWTC-3 are consistent with having random orbital orientations. The inferred distribution is highly symmetric around $\pi/2$, with skewness $\mathcal{S}_{\rm{post}}=0.01^{+0.17}_{-0.17}$. Meanwhile, the median of the inferred distribution has a Jensen-Shannon divergence of $1.4\times 10^{-4}$ bits when compared to the expected isotropic distribution., Comment: 5 pages, 3 figures, a lot of fun

The final black hole left behind after a binary black hole merger can attain a recoil velocity, or a "kick", reaching values up to 5000 km/s. This phenomenon has important implications for gravitational wave astronomy, black hole formation scenarios, testing general relativity, and galaxy evolution. We consider the gravitational wave signal from the binary black hole merger GW200129_065458 (henceforth referred to as GW200129), which has been shown to exhibit strong evidence of orbital precession. Using numerical relativity surrogate models, we constrain the kick velocity of GW200129 to $v_f \sim 1542^{+747}_{-1098}$ km/s or $v_f \gtrsim 698$ km/s (one-sided limit), at 90\% credibility. This marks the first identification of a large kick velocity for an individual gravitational wave event. Given the kick velocity of GW200129, we estimate that there is a less than $0.48\%$ ($7.7\%$) probability that the remnant black hole after the merger would be retained by globular (nuclear star) clusters. Finally, we show that kick effects are not expected to cause biases in ringdown tests of general relativity for this event, although this may change in the future with improved detectors., Comment: 5 pages. Visualizations available at https://vijayvarma392.github.io/GW200129. Matches published version

The binary neutron star (BNS) mass distribution measured with gravitational-wave observations has the potential to reveal information about the dense matter equation of state, supernova physics, the expansion rate of the universe, and tests of General Relativity. As most current gravitational-wave analyses measuring the BNS mass distribution do not simultaneously fit the spin distribution, the implied population-level spin distribution is the same as the spin prior applied when analyzing individual sources. In this work, we demonstrate that introducing a mismatch between the implied and true BNS spin distributions can lead to biases in the inferred mass distribution. This is due to the strong correlations between the measurements of the mass ratio and spin components aligned with the orbital angular momentum for individual sources. We find that applying a low-spin prior which excludes the true spin magnitudes of some sources in the population leads to significantly overestimating the maximum neutron star mass and underestimating the minimum neutron star mass at the population level with as few as six BNS detections. The safest choice of spin prior that does not lead to biases in the inferred mass distribution is one which allows for high spin magnitudes and tilts misaligned with the orbital angular momentum., Comment: version accepted by journal

The Wide-Field Infrared Transient Explorer (WINTER) is a new 1 $\text{deg}^2$ seeing-limited time-domain survey instrument designed for dedicated near-infrared follow-up of kilonovae from binary neutron star (BNS) and neutron star-black hole mergers. WINTER will observe in the near-infrared Y, J, and short-H bands (0.9-1.7 microns, to $\text{J}_{AB}=21$ magnitudes) on a dedicated 1-meter telescope at Palomar Observatory. To date, most prompt kilonova follow-up has been in optical wavelengths; however, near-infrared emission fades more slowly and depends less on geometry and viewing angle than optical emission. We present an end-to-end simulation of a follow-up campaign during the fourth observing run (O4) of the LIGO, Virgo, and KAGRA interferometers, including simulating 625 BNS mergers, their detection in gravitational waves, low-latency and full parameter estimation skymaps, and a suite of kilonova lightcurves from two different model grids. We predict up to five new kilonovae independently discovered by WINTER during O4, given a realistic BNS merger rate. Using a larger grid of kilonova parameters, we find that kilonova emission is $\approx$2 times longer-lived and red kilonovae are detected $\approx$1.5 times further in the infrared than in the optical. For 90% localization areas smaller than 150 (450) $\rm{deg}^{2}$, WINTER will be sensitive to more than 10% of the kilonova model grid out to 350 (200) Mpc. We develop a generalized toolkit to create an optimal BNS follow-up strategy with any electromagnetic telescope and present WINTER's observing strategy with this framework. This toolkit, all simulated gravitational-wave events, and skymaps are made available for use by the community., Comment: 19 pages, 11 figures, published in ApJ

This Horizon Study describes a next-generation ground-based gravitational-wave observatory: Cosmic Explorer. With ten times the sensitivity of Advanced LIGO, Cosmic Explorer will push gravitational-wave astronomy towards the edge of the observable universe ($z \sim 100$). The goals of this Horizon Study are to describe and evaluate design concepts for Cosmic Explorer; to plan for the United States' leadership in gravitational-wave astronomy; and to envisage the role of Cosmic Explorer in the international effort to build a "Third-Generation" (3G) observatory network that will make discoveries transformative across astronomy, physics, and cosmology., Comment: 173 pages

Binary black hole spin measurements from gravitational wave observations can reveal the binary's evolutionary history. In particular, the spin orientations of the component black holes within the orbital plane, $\phi_1$ and $\phi_2$, can be used to identify binaries caught in the so-called spin-orbit resonances. In a companion paper, we demonstrate that $\phi_1$ and $\phi_2$ are best measured near the merger of the two black holes. In this work, we use these spin measurements to provide the first constraints on the full six-dimensional spin distribution of merging binary black holes. In particular, we find that there is a preference for $\Delta \phi = \phi_1 - \phi_2 \sim \pm \pi$ in the population, which can be a signature of spin-orbit resonances. We also find a preference for $\phi_1 \sim -\pi/4$ with respect to the line of separation near merger, which has not been predicted for any astrophysical formation channel. However, the strength of these preferences depends on our prior choices, and we are unable to constrain the widths of the $\phi_1$ and $\Delta \phi$ distributions. Therefore, more observations are necessary to confirm the features we find. Finally, we derive constraints on the distribution of recoil kicks in the population, and use this to estimate the fraction of merger remnants retained by globular and nuclear star clusters. We make our spin and kick population constraints publicly available., Comment: 5 pages plus supplement. Matches PRL version. Results available publicly at https://github.com/vijayvarma392/spin_kick_pop_GWTC2

Binary black hole spins are among the key observables for gravitational wave astronomy. Among the spin parameters, their orientations within the orbital plane, $\phi_1$, $\phi_2$ and $\Delta \phi=\phi_1-\phi_2$, are critical for understanding the prevalence of the spin-orbit resonances and merger recoils in binary black holes. Unfortunately, these angles are particularly hard to measure using current detectors, LIGO and Virgo. Because the spin directions are not constant for precessing binaries, the traditional approach is to measure the spin components at some reference stage in the waveform evolution, typically the point at which the frequency of the detected signal reaches 20 Hz. However, we find that this is a poor choice for the orbital-plane spin angle measurements. Instead, we propose measuring the spins at a fixed dimensionless time or frequency near the merger. This leads to significantly improved measurements for $\phi_1$ and $\phi_2$ for several gravitational wave events. Furthermore, using numerical relativity injections, we demonstrate that $\Delta \phi$ will also be better measured near the merger for louder signals expected in the future. Finally, we show that numerical relativity surrogate models are key for reliably measuring the orbital-plane spin orientations, even at moderate signal-to-noise ratios like $\sim 30-45$., Comment: 15 pages plus appendix. Matches PRD version

Time series analysis is ubiquitous in many fields of science including gravitational-wave astronomy, where strain time series are analyzed to infer the nature of gravitational-wave sources, e.g., black holes and neutron stars. It is common in gravitational-wave transient studies to apply a tapered window function to reduce the effects of spectral artifacts from the sharp edges of data segments. We show that the conventional analysis of tapered data fails to take into account covariance between frequency bins, which arises for all finite time series -- no matter the choice of window function. We discuss the origin of this covariance and show that as the number of gravitational-wave detections grows, and as we gain access to more high signal-to-noise ratio events, this covariance will become a non-negligible source of systematic error. We derive a framework that models the correlation induced by the window function and demonstrate this solution using both data from the first LIGO--Virgo transient catalog and simulated Gaussian noise., Comment: 10 pages, 9 figures in main body, accepted in PRR

An accurate and precise measurement of the spins of individual merging black holes is required to understand their origin. While previous studies have indicated that most of the spin information comes from the inspiral part of the signal, the informative spin measurement of the heavy binary black hole system GW190521 suggests that the merger and ringdown can contribute significantly to the spin constraints for such massive systems. We perform a systematic study into the measurability of the spin parameters of individual heavy binary black hole mergers using a numerical relativity surrogate waveform model including the effects of both spin-induced precession and higher-order modes. We find that the spin measurements are driven by the merger and ringdown parts of the signal for GW190521-like systems, but the uncertainty in the measurement increases with the total mass of the system. We are able to place meaningful constraints on the spin parameters even for systems observed at moderate signal-to-noise ratios, but the measurability depends on the exact six-dimensional spin configuration of the system. Finally, we find that the azimuthal angle between the in-plane projections of the component spin vectors at a given reference frequency cannot be well-measured for most of our simulated configurations even for signals observed with high signal-to-noise ratios., Comment: Version accepted in PRD

Primordial gravitational waves are expected to create a stochastic background encoding information about the early Universe that may not be accessible by other means. However, the primordial background is obscured by an astrophysical foreground consisting of gravitational waves from compact binaries. We demonstrate a Bayesian method for estimating the primordial background in the presence of an astrophysical foreground. Since the background and foreground signal parameters are estimated simultaneously, there is no subtraction step, and therefore we avoid astrophysical contamination of the primordial measurement, sometimes referred to as "residuals". Additionally, since we include the non-Gaussianity of the astrophysical foreground in our model, this method represents the statistically optimal approach to the simultaneous detection of a multi-component stochastic background.

Gravitational waves from binary black holes have the potential to yield information on both of the intrinsic parameters that characterize the compact objects: their masses and spins. While the component masses are usually resolvable, the component spins have proven difficult to measure. This limitation stems in great part from our choice to inquire about the spins of the most and least massive objects in each binary, a question that becomes ill-defined when the masses are equal. In this paper we show that one can ask a different question of the data: what are the spins of the objects with the highest and lowest dimensionless spins in the binary? We show that this can significantly improve estimates of the individual spins, especially for binary systems with comparable masses. When applying this parameterization to the first 13 gravitational-wave events detected by the LIGO-Virgo collaboration (LVC), we find that the highest-spinning object is constrained to have nonzero spin for most sources and to have significant support at the Kerr limit for GW151226 and GW170729. A joint analysis of all the confident binary black hole detections by the LVC finds that, unlike with the traditional parametrization, the distribution of spin magnitude for the highest-spinning object has negligible support at zero spin. Regardless of the parameterization used, the configuration where all of the spins in the population are aligned with the orbital angular momentum is excluded from the 90% credible interval for the first ten events and from the 99% credible interval for all current confident detections., Comment: version accepted in PRL

Gravitational waves emitted by neutron star black hole mergers encode key properties of neutron stars - such as their size, maximum mass and spins - and black holes. However, the presence of matter and the high mass ratio makes generating long and accurate waveforms from these systems hard to do with numerical relativity, and not much is known about systematic uncertainties due to waveform modeling. We simulate gravitational waves from neutron star black hole mergers by hybridizing numerical relativity waveforms produced with the SpEC code with a recent numerical relativity surrogate NRHybSur3dq8Tidal. These signals are analyzed using a range of available waveform families, and statistical and systematic errors are reported. We find that at a network signal-to-noise ratio (SNR) of 30, statistical uncertainties are usually larger than systematic offsets, while at an SNR of 70 the two become comparable. The individual black hole and neutron star masses, as well as the mass ratios, are typically measured very precisely, though not always accurately at high SNR. At a SNR of 30 the neutron star tidal deformability can only be bound from above, while for louder sources it can be measured and constrained away from zero. All neutron stars in our simulations are non-spinning, but in no case we can constrain the neutron star spin to be smaller than $\sim0.4$ (90% credible interval). Waveform families whose late inspiral has been tuned specifically for neutron star black hole signals typically yield the most accurate characterization of the source parameters. Their measurements are in tension with those obtained using waveform families tuned against binary neutron stars, even for mass ratios that could be relevant for both binary neutron stars and neutron star black holes mergers., Comment: 39 pages, 24 figures, 15 tables

The initial separation of massive star binaries sets the timescale over which their compact remnants merge through the emission of gravitational waves. We show that the delay time distribution (DTD) of binary neutron stars or black holes can be inferred from the stochastic gravitational wave background (SGWB). If the DTD of a population is long, most of the mergers take place at low redshifts and the background would be rather quiet compared to a scenario in which the DTD is short leading to few individual detections at low redshift but a rather loud background. We show that different DTDs predict a factor of 5 difference in the magnitude of the gravitational wave background energy density ($\Omega_{\rm GW}$) and have the dominant effect on $\Omega_{\rm GW}$ over other factors such as the mass function of the primary BH mass, $p(m_1)$, the maximum considered BH mass ($M_{\rm max}$), and the effective spin of the black hole ($\chi_{\rm eff}$). A non-detection of such a background can rule out the short DTD scenario. We show that SGWB searches can rule out the short DTD scenario for the BBHs within about four years of observing time at advanced LIGO design sensitivty for a local merger rate of 30 $\rm Gpc^{-3} yr^{-1}$ assuming $p(m_1)\propto m_1^{-1}$, and $M_{\rm max}=50 M_{\odot}$., Comment: Accepted for publication at ApJ

In order to perform Bayesian parameter estimation to infer the source properties of gravitational waves from compact binary coalescences (CBCs), the noise characteristics of the detector must be understood. It is typically assumed that the detector noise is stationary and Gaussian, characterized by a power spectral density (PSD) that is measured with infinite precision. We present a new method to incorporate the uncertainty in the power spectral density estimation into the Bayesian inference of the binary source parameters and apply it to the first 11 CBC detections reported by the LIGO- Virgo Collaboration. We find that incorporating the PSD uncertainty only leads to variations in the positions and widths of the binary parameter posteriors on the order of a few percent. Our results are publicly available for download on git [1].

Gravitational waves carry energy, angular momentum, and linear momentum. In generic binary black hole mergers, the loss of linear momentum imparts a recoil velocity, or a "kick", to the remnant black hole. We exploit recent advances in gravitational waveform and remnant black hole modeling to extract information about the kick from the gravitational wave signal. Kick measurements such as these are astrophysically valuable, enabling independent constraints on the rate of second-generation mergers. Further, we show that kicks must be factored into future ringdown tests of general relativity with third-generation gravitational wave detectors to avoid systematic biases. We find that, although little information can be gained about the kick for existing gravitational wave events, interesting measurements will soon become possible as detectors improve. We show that, once LIGO and Virgo reach their design sensitivities, we will reliably extract the kick velocity for generically precessing binaries--including the so-called superkicks, reaching up to 5000 km/s., Comment: Matches published version. 10 pages, including supplement

Gamma-ray burst (GRB) prompt emission is highly beamed, and understanding the jet geometry and beaming configuration can provide information on the poorly understood central engine and circum-burst environment. Prior to the advent of gravitational-wave astronomy, astronomers relied on observations of jet breaks in the multi-wavelength afterglow to determine the GRB opening angle, since the observer's viewing angle relative to the system cannot be determined from the electromagnetic data alone. Gravitational-wave observations, however, provide an independent measurement of the viewing angle. We describe a Bayesian method for determining the geometry of short GRBs using coincident electromagnetic and gravitational-wave observations. We demonstrate how an ensemble of multi-messenger detections can be used to measure the distributions of the jet energy, opening angle, Lorentz factor, and angular profile of short GRBs; we find that for a population of 100 such observations, we can constrain the mean of the opening angle distribution to within 10 degrees regardless of the angular emission profile. Conversely, the constraint on the energy distribution depends on the shape of the profile, which can be distinguished., Comment: Matches published version

The LIGO and Virgo Collaborations currently conduct searches for gravitational waves from compact binary coalescences in real-time. For promising candidate events, a sky map and distance estimation are released in low-latency, to facilitate their electromagnetic follow-up. Currently, no information is released about the masses of the compact objects. Recently, Margalit and Metzger (2019) have suggested that knowledge of the chirp mass of the detected binary neutron stars could be useful to prioritize the electromagnetic follow-up effort, and have urged the LIGO-Virgo collaboration to release chirp mass information in low-latency. One might worry that low-latency searches for compact binaries make simplifying assumptions that could introduce biases in the mass parameters: neutron stars are treated as point particles with dimensionless spins below $0.05$ and perfectly aligned with the orbital angular momentum. Furthermore, the template bank used to search for them has a finite resolution. In this paper we show that none of these limitations can introduce chirp mass biases larger than $\sim 10^{-3}~M_\odot$. Even the total mass is usually accurately estimated, with biases smaller than 6%. The mass ratio and effective inspiral spins, on the other hand, can suffer from more severe biases.

Recent gravitational-wave observations from the LIGO and Virgo observatories have brought a sense of great excitement to scientists and citizens the world over. Since September 2015,10 binary black hole coalescences and one binary neutron star coalescence have been observed. They have provided remarkable, revolutionary insight into the "gravitational Universe" and have greatly extended the field of multi-messenger astronomy. At present, Advanced LIGO can see binary black hole coalescences out to redshift 0.6 and binary neutron star coalescences to redshift 0.05. This probes only a very small fraction of the volume of the observable Universe. However, current technologies can be extended to construct "$3^\mathrm{rd}$ Generation" (3G) gravitational-wave observatories that would extend our reach to the very edge of the observable Universe. The event rates over such a large volume would be in the hundreds of thousands per year (i.e.tens per hour). Such 3G detectors would have a 10-fold improvement in strain sensitivity over the current generation of instruments, yielding signal-to-noise ratios of 1000 for events like those already seen. Several concepts are being studied for which engineering studies and reliable cost estimates will be developed in the next 5 years., Comment: For the 2020 Astro decadal

Bayesian parameter estimation is fast becoming the language of gravitational-wave astronomy. It is the method by which gravitational-wave data is used to infer the sources' astrophysical properties. We introduce a user-friendly Bayesian inference library for gravitational-wave astronomy, Bilby. This python code provides expert-level parameter estimation infrastructure with straightforward syntax and tools that facilitate use by beginners. It allows users to perform accurate and reliable gravitational-wave parameter estimation on both real, freely-available data from LIGO/Virgo, and simulated data. We provide a suite of examples for the analysis of compact binary mergers and other types of signal model including supernovae and the remnants of binary neutron star mergers. These examples illustrate how to change the signal model, how to implement new likelihood functions, and how to add new detectors. Bilby has additional functionality to do population studies using hierarchical Bayesian modelling. We provide an example in which we infer the shape of the black hole mass distribution from an ensemble of observations of binary black hole mergers., Comment: Submitted for publication