Search

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

Gravitational lensing offers unique opportunities to learn about the astrophysical origin of distant sources, the abundance of intervening objects acting as lenses, and gravity and cosmology in general. However, all this information can only be retrieved as long as one can disentangle each effect from the finite number of observables. In the geometric optics regime, typical of electromagnetic radiation, when the wavelength of the lensed signal is small compared to the size of the lens, there are invariance transformations that change the mass of the lens and the source-lens configuration but leave the observables unchanged. Neglecting this ``mass-sheet degeneracy'' can lead to biased lens parameters or unrealistic low uncertainties, which could then transfer to an incorrect cosmography study. This might be different for gravitational waves as their long wavelengths can be comparable to the lens size and lensing enters into the wave-optics limit. We explore the existence of invariance transformations in the wave-optics regime of gravitational-wave lensing, extending previous work and examining the implications for astrophysical and cosmological studies. We study these invariance transformations using three different methods of increasing level of complexity: template mismatch, Fisher Matrix, and Bayesian parameter estimation. We find that, for a sufficiently loud signal, the degeneracy is partially broken and the lens and cosmological parameters, e.g. $H_0$, can be retrieved independently and unbiased. In current ground-based detectors, though, considering also population studies, a strong constraint on these parameters seems quite remote and the prevailing degeneracy implies a larger uncertainty in the lens model reconstruction. However, with better sensitivity of the third-generation ground-based detectors, a meaningful constraint on $H_0$ is possible to obtain., Comment: 23 pages, 24 figures. Comments are welcome

The LIGO-Virgo-KAGRA (LVK) Collaboration has made breakthrough discoveries in gravitational-wave astronomy, a new field that provides a different means of observing our Universe. Gravitational-wave discoveries are possible thanks to the work of thousands of people from across the globe working together. In this article, we discuss the range of engagement activities used to communicate LVK gravitational-wave discoveries and the stories of the people behind the science, using the activities surrounding the release of the third Gravitational-Wave Transient Catalog as a case study., Comment: 15 pages, 7 figures, published in JCOM: https://doi.org/10.22323/2.23070803

For the first detection of a novel astrophysical phenomenon, scientific standards are particularly high. Especially in a multi-messenger context, there are also opportunity costs to follow-up observations on any detection claims. So in searching for the still elusive lensed gravitational waves, care needs to be taken in controlling false positives. In particular, many methods for identifying strong lensing rely on some form of parameter similarity or waveform consistency, which under rapidly growing catalog sizes can expose them to false positives from coincident but unlensed events if proper care is not taken. And searches for waveform deformations in all lensing regimes are subject to degeneracies we need to mitigate between lensing, intrinsic parameters, insufficiently modelled effects such as orbital eccentricity, or even deviations from general relativity. Robust lensing studies also require understanding and mitigating glitches and non-stationarities in the detector data. This article reviews sources of possible false positives (and their flip side: false negatives) in gravitational-wave lensing searches and the main approaches the community is pursuing to mitigate them., Comment: 14 pages, 2 figures, submitted to Philosophical Transactions of the Royal Society special issue on Multi-Messenger Gravitational Lensing

'$ler$' is a statistics-based Python package specifically designed for computing detectable rates of both lensed and unlensed GW events, catering to the requirements of the LIGO-Virgo-KAGRA Scientific Collaboration and astrophysics research scholars. The core functionality of '$ler$' intricately hinges upon the interplay of various components which include sampling the properties of compact-binary sources, lens galaxies characteristics, solving lens equations to derive properties of resultant images, and computing detectable GW rates. This comprehensive functionality builds on the leveraging of array operations and linear algebra from the $numpy$ library, enhanced by interpolation methods from $scipy$ and Python's $multiprocessing$ capabilities. Efficiency is further boosted by the $numba$ library's Just-In-Time ($njit$) compilation, optimizing extensive numerical computations and employing the inverse transform sampling method to replace more cumbersome rejection sampling. The modular design of '$ler$' not only optimizes speed and functionality but also ensures adaptability and upgradability, supporting the integration of additional statistics as research evolves. Currently, '$ler$' is an important tool in generating simulated GW events, both lensed and unlensed, and provides astrophysically accurate distributions of event-related parameters for both detectable and non-detectable events. This functionality aids in event validation and enhances the forecasting of detection capabilities across various GW detectors to study such events. The architecture of the '$ler$' API facilitates seamless compatibility with other software packages, allowing researchers to integrate and utilize its functionalities based on specific scientific requirements., Comment: 5 pages, 1 Logo in each of the pages, this is for the JOSS publication

Gravitational lensing has been extensively observed for electromagnetic signals, but not yet for gravitational waves (GWs). Detecting lensed GWs will have many astrophysical and cosmological applications, and becomes more feasible as the sensitivity of the ground-based detectors improves. One of the missing ingredients to robustly identify lensed GWs is to ensure that the statistical tests used are robust under the choice of underlying waveform models. We present the first systematic study of possible waveform systematics in identifying candidates for strongly lensed GW event pairs, focusing on the posterior overlap method. To this end, we compare Bayes factors from all posteriors using different waveforms included in GWTC data releases from the first three observing runs (O1-O3). We find that waveform choice yields a wide spread of Bayes factors in some cases. However, it is likely that no event pairs from O1 to O3 were missed due to waveform choice. We also perform parameter estimation with additional waveforms for interesting cases, to understand the observed differences. We also briefly explore if computing the overlap from different runs for the same event can be a useful metric for waveform systematics or sampler issues, independent of the lensing scenario., Comment: 31 pages, 5 figures, matches version accepted by CQG

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

Machine learning can be a powerful tool to discover new signal types in astronomical data. We here apply it to search for long-duration transient gravitational waves triggered by pulsar glitches, which could yield physical insight into the mostly unknown depths of the pulsar. Current methods to search for such signals rely on matched filtering and a brute-force grid search over possible signal durations, which is sensitive but can become very computationally expensive. We develop a method to search for post-glitch signals on combining matched filtering with convolutional neural networks, which reaches similar sensitivities to the standard method at false-alarm probabilities relevant for practical searches, while being significantly faster. We specialize to the Vela glitch during the LIGO-Virgo O2 run, and set upper limits on the gravitational-wave strain amplitude from the data of the two LIGO detectors for both constant-amplitude and exponentially decaying signals., Comment: 19 pages, 9 figures. Comments welcome

Improved observational constraints on the orbital parameters of the low-mass X-ray binary Scorpius~X-1 were recently published in Killestein et al (2023). In the process, errors were corrected in previous orbital ephemerides, which have been used in searches for continuous gravitational waves from Sco~X-1 using data from the Advanced LIGO detectors. We present the results of a re-analysis of LIGO detector data from the third observing run of Advanced LIGO and Advanced Virgo using a model-based cross-correlation search. The corrected region of parameter space, which was not covered by previous searches, was about 1/3 as large as the region searched in the original O3 analysis, reducing the required computing time. We have confirmed that no detectable signal is present over a range of gravitational-wave frequencies from $25\textrm{Hz}$ to $1600\textrm{Hz}$, analogous to the null result of Abbott et al (2022). Our search sensitivity is comparable to that of Abbott et al (2022), who set upper limits corresponding, between $100\textrm{Hz}$ and $200\textrm{Hz}$, to an amplitude $h_0$ of about $10^{-25}$ when marginalized isotropically over the unknown inclination angle of the neutron star's rotation axis, or less than $4\times 10^{-26}$ assuming the optimal orientation., Comment: 8 pages, 3 figures, 2 tables. Typeset with AASTeX 6.3.1. Accepted for publication in The Astrophysical Journal. arXiv admin note: text overlap with arXiv:2209.02863

Pulsars are rotating neutron stars that emit periodic electromagnetic radiation. While pulsars generally slow down as they lose energy, some also experience glitches: spontaneous increases of their rotational frequency. According to several models, these glitches can also lead to the emission of long-duration transient gravitational waves (GWs). We present detection prospects for such signals by comparing indirect energy upper limits on GW strain for known glitches with the sensitivity of current and future ground-based GW detectors. We first consider the optimistic case of generic constraints based on the glitch size and find that realistic matched-filter searches in the fourth LIGO-Virgo-KAGRA observing run(O4) could make a detection, or set constraints below these indirect upper limits, for equivalents of 36 out of 726 previously observed glitches, and 74 in the O5 run. With the third-generation Einstein Telescope or Cosmic Explorer, 35-40% of glitches would be accessible. When specialising to a scenario where transient mountains produce the post-glitch GW emission, following the Yim & Jones model, the indirect upper limits are stricter. Out of the smaller set of 119 glitches with measured healing parameter, as needed for predictions under that model, only 6 glitches would have been within reach for O4 and 14 for O5, with a similar percentage as before with third generation detectors. We also discuss how this model matches the observed glitch population., Comment: 16 pages, 9 figures, matches version accepted by MNRAS

Neutron stars are known to show accelerated spin-up of their rotational frequency called a glitch. Highly magnetized rotating neutron stars (pulsars) are frequently observed by radio telescopes (and in other frequencies), where the glitch is observed as irregular arrival times of pulses which are otherwise very regular. A glitch in an isolated neutron star can excite the fundamental (f)-mode oscillations which can lead to gravitational wave generation. Electromagnetic observations of pulsars (and hence pulsar glitches) require the pulsar to be oriented so that the jet is pointed toward the detector, but this is not a requirement for gravitational wave emission which is more isotropic and not jetlike. Hence, gravitational wave observations have the potential to uncover nearby neutron stars where the jet is not pointed towards the Earth. In this work, we study the prospects of finding glitching neutron stars using a generic all-sky search for short-duration gravitational wave transients. The analysis covers the high-frequency range from $1-4$ kHz of LIGO-Virgo detectors for signals up to a few seconds. We set upper limits for the third observing run of the LIGO-Virgo detectors and present the prospects for upcoming observing runs of LIGO, Virgo, KAGRA, and LIGO India. We find the detectable glitch size will be around $10^{-5}$ Hz for the fifth observing run for pulsars with spin frequencies and distances comparable to the Vela pulsar. We also present the prospects of localizing the direction in the sky of these sources with gravitational waves alone, which can facilitate electromagnetic follow-up. We find that for the five detector configuration, the localization capability for a glitch size of $10^{-5}$ Hz is around $132\,\mathrm{deg}^{2}$ at $1\sigma$ confidence for $50\%$ of events with distance and spin frequency as that of Vela.

Pulsars are spinning neutron stars which emit an electromagnetic beam. We expect pulsars to slowly decrease their rotational frequency. However, sudden increases of the rotational frequency have been observed from different pulsars. These events are called "glitches" and they are followed by a relaxation phase with timescales from days to months. Gravitational wave (GW) emission may follow these peculiar events. We give an overview of the setup for an analysis of GW data from the Advanced LIGO and Virgo third observing run (O3) arXiv:2112.10990 searching for transient GW signals lasting hours to months after glitches in known pulsars. The search method consists of placing a template grid in frequency-spindown space with fixed grid spacings. Then, for each point we compute the transient F-statistic which is maximized over a set of transient parameters like the duration and start time of the potential signals. A threshold on the detection statistic is then set, and we search for peaks over the parameter space for each candidate., Comment: 5 pages, submitted to Proceedings of TAUP 2021

The general-relativistic phenomenon of spin-induced orbital precession has not yet been observed in strong-field gravity. Gravitational-wave observations of binary black holes (BBHs) are prime candidates, since we expect the astrophysical binary population to contain precessing binaries. Imprints of precession have been investigated in several signals, but no definitive identification of orbital precession has been reported in any one of the 84 BBH observations to date by the Advanced LIGO and Virgo detectors. Here we report the measurement of strong-field precession in the LIGO-Virgo-Kagra (LVK) gravitational-wave signal GW200129. The binary's orbit precesses at a rate ten orders of magnitude faster than previous weak-field measurements from binary pulsars. We also find that the primary black hole is likely highly spinning. According to current binary population estimates a GW200129-like signal is extremely unlikely, and therefore presents a direct challenge to many current binary formation models., Comment: 16 pages, 6 figures. Update to match published version

Continuous gravitational waves are long-lasting forms of gravitational radiation produced by persistent quadrupolar variations of matter. Standard expected sources for ground-based interferometric detectors are neutron stars presenting non-axisymmetries such as crustal deformations, r-modes or free precession. More exotic sources could include decaying ultralight boson clouds around spinning black holes. A rich suite of data-analysis methods spanning a wide bracket of thresholds between sensitivity and computational efficiency has been developed during the last decades to search for these signals. In this work, we review the current state of searches for continuous gravitational waves using ground-based interferometer data, focusing on searches for unknown sources. These searches typically consist of a main stage followed by several post-processing steps to rule out outliers produced by detector noise. So far, no continuous gravitational wave signal has been confidently detected, although tighter upper limits are placed as detectors and search methods are further developed., Comment: 31 pages, 2 figures, matches published version

Searches for gravitational-wave signals are often based on maximizing a detection statistic over a bank of waveform templates, covering a given parameter space with a variable level of correlation. Results are often evaluated using a noise-hypothesis test, where the background is characterized by the sampling distribution of the loudest template. In the context of continuous gravitational-wave searches, properly describing said distribution is an open problem: current approaches focus on a particular detection statistic and neglect template-bank correlations. We introduce a new approach using extreme value theory to describe the distribution of the loudest template's detection statistic in an arbitrary template bank. Our new proposal automatically generalizes to a wider class of detection statistics, including (but not limited to) line-robust statistics and transient continuous-wave signal hypotheses, and improves the estimation of the expected maximum detection statistic at a negligible computing cost. The performance of our proposal is demonstrated on simulated data as well as by applying it to different kinds of (transient) continuous-wave searches using O2 Advanced LIGO data. We release an accompanying Python software package, distromax, implementing our new developments., Comment: 24 pages, 23 figures, comments welcome. Package freely available in https://github.com/Rodrigo-Tenorio/distromax

We present the first application of a hierarchical Markov Chain Monte Carlo (MCMC) follow-up on continuous gravitational-wave candidates from real-data searches. The follow-up uses an MCMC sampler to draw parameter-space points from the posterior distribution, constructed using the matched-filter as a log-likelihood. As outliers are narrowed down, coherence time increases, imposing more restrictive phase-evolution templates. We introduce a novel Bayes factor to compare results from different stages: The signal hypothesis is derived from first principles, while the noise hypothesis uses extreme value theory to derive a background model. The effectiveness of our proposal is evaluated on fake Gaussian data and applied to a set of 30 outliers produced by different continuous wave searches on O2 Advanced LIGO data. The results of our analysis suggest all but five outliers are inconsistent with an astrophysical origin under the standard continuous wave signal model. We successfully ascribe four of the surviving outliers to instrumental artifacts and a strong hardware injection present in the data. The behavior of the fifth outlier suggests an instrumental origin as well, but we could not relate it to any known instrumental cause., Comment: 23 pages, 15 figures, matches PRD published version

In this paper we present an extensive analysis of the GW190521 gravitational wave event with the current (fourth) generation of phenomenological waveform models for binary black hole coalescences. GW190521 stands out from other events since only a few wave cycles are observable. This leads to a number of challenges, one being that such short signals are prone to not resolve approximate waveform degeneracies, which may result in multi-modal posterior distributions. The family of waveform models we use includes a new fast time-domain model IMRPhenomTPHM, which allows us extensive tests of different priors and robustness with respect to variations in the waveform model, including the content of spherical harmonic modes. We clarify some issues raised in a recent paper [Nitz&Capano], associated with possible support for a high-mass ratio source, but confirm their finding of a multi-modal posterior distribution, albeit with important differences in the statistical significance of the peaks. In particular, we find that the support for both masses being outside the PISN mass-gap, and the support for an intermediate mass ratio binary are drastically reduced with respect to what Nitz&Capano found. We also provide updated probabilities for associating GW190521 to the potential electromagnetic counterpart from ZTF., Comment: 20 pages, 13 figures, 7 tables. Data release at https://zenodo.org/record/4751655

We present IMRPhenomTPHM, a phenomenological model for the gravitational wave signals emitted by the coalescence of quasi-circular precessing binary black holes systems. The model is based on the "twisting up" approximation, which maps non-precessing signals to precessing ones in terms of a time dependent rotation described by three Euler angles, and which has been utilized in several frequency domain waveform models that have become standard tools in gravitational wave data analysis. Our model is however constructed in the time domain, which allows several improvements over the frequency domain models: we do not use the stationary phase approximation, we employ a simple approximation for the precessing Euler angles for the ringdown signal, and we implement a new method for computing the Euler angles through the evolution of the spin dynamics of the system, which is more accurate and also computationally efficient., Comment: 18 pages, 10 figures, 4 tables

We consider the ten confidently detected gravitational-wave signals in the GWTC-1 catalog which are consistent with mergers of binary black hole systems, and perform a thorough parameter estimation re-analysis. This is made possible by using computationally efficient waveform models of the current (fourth) generation of the IMRPhenom family of phenomenological waveform models, which consists of the IMRPhenomX frequency-domain modelsand the IMRPhenomT time-domain models. The analysis is performed with both precessing and non-precessing waveform models with and without subdominant spherical harmonic modes. Results for all events are validated with convergence tests, discussing in particular the events GW170729 and GW151226. For the latter and the other two lowest-mass events, we also compare results between two independent sampling codes, Bilby and LALInference. We find overall consistent results with the original GWTC-1 results, with all Jensen-Shannon divergences between the previous results using IMRPhenomPv2 and our default IMRPhenomXPHM posteriors below 0.045 bits, but we also discuss cases where including subdominant harmonics and/or precession influences the posteriors., Comment: 24 pages, 19 figures, 9 tables matches version accepted by MNRAS

Gravitational waves in the sensitivity band of ground-based detectors can be emitted by a number of astrophysical sources, including not only binary coalescences, but also individual spinning neutron stars. The most promising signals from such sources, although not yet detected, are long-lasting, quasi-monochromatic Continuous Waves (CWs). The PyFstat package provides tools to perform a range of CW data-analysis tasks. It revolves around the F-statistic, a matched-filter detection statistic for CW signals that has been one of the standard methods for LIGO-Virgo CW searches for two decades. PyFstat is built on top of established routines in LALSuite but through its more modern Python interface it enables a flexible approach to designing new search strategies. Hence, it serves a dual function of (i) making LALSuite CW functionality more easily accessible through a Python interface, thus facilitating the new user experience and, for developers, the exploratory implementation of novel methods; and (ii) providing a set of production-ready search classes for use cases not yet covered by LALSuite itself, most notably for MCMC-based followup of promising candidates from wide-parameter-space searches., Comment: 4 pages, updated to match published version. Software repository: https://github.com/PyFstat/PyFstat/

In this work we present an extension of the time domain phenomenological model IMRPhenomT for gravitational wave signals from binary black hole coalescences to include subdominant harmonics, specifically the $(l=2, m=\pm 1)$, $(l=3, m=\pm 3)$, $(l=4, m=\pm 4)$ and $(l=5, m=\pm 5)$ spherical harmonics. We also improve our model for the dominant $(l=2, m=\pm 2)$ mode and discuss mode mixing for the $(l=3, m=\pm 2)$ mode. The model is calibrated to numerical relativity solutions of the full Einstein equations up to mass ratio 18, and to numerical solutions of the Teukolsky equations for higher mass ratios. This work complements the latest generation of traditional frequency domain phenomenological models (IMRPhenomX), and provides new avenues to develop computationally efficient models for gravitational wave signals from generic compact binaries., Comment: 22 pages, 18 figures. https://dcc.ligo.org/P2000524

Searches for continuous gravitational waves from unknown sources attempt to detect long-lasting gravitational radiation by identifying Doppler-modulated signatures in the data. Semicoherent methods allow for wide parameter space surveys, identifying interesting regions to be followed up using more sensitive (and computationally expensive) tools. Thus, it is required to properly understand the parameter space structure under study, as failing to do so could significantly affect the effectiveness of said strategies. We introduce a new measure for distances in parameter space suited for semicoherent continuous wave searches. This novel approach, based on comparing time-frequency tracks, can be applied to any kind of quasi-monochromatic continuous wave signals and adapts itself to the underlying structure of the parameter space under study. In a first application to the post-processing stage of an all-sky search for continuous waves from neutron stars in binary systems, we demonstrate a search sensitivity improvement by solely replacing previous ad hoc distance measures in the candidate clustering procedure by the new proposal., Comment: 13 pages, 7 figures. Matches PRD published version

We re-analyse the gravitational-wave event GW190412 with state-of-the-art phenomenological waveform models. This event, which has been associated with a black hole merger, is interesting due to the significant contribution from subdominant harmonics. We use both frequency-domain and time-domain waveform models. The PhenomX waveform models constitute the fourth generation of frequency-domain phenomenological waveforms for black hole binary coalescence; they have more recently been complemented by the time-domain PhenomT models, which open up new strategies to model precession and eccentricity, and to perform tests of general relativity with the phenomenological waveforms approach. Both PhenomX and PhenomT have been constructed with similar techniques and accuracy goals, and due to their computational efficiency this "generation X" model family allows the routine use of subdominant spherical harmonics in Bayesian inference. We show the good agreement between these and other state-of-the-art waveform models for GW190412, and discuss the improvements over the previous generation of phenomenological waveform models. We also discuss practical aspects of Bayesian inference such as run convergence, variations of sampling parameters, and computational cost., Comment: 25 pages, 18 figures. Data release: https://doi.org/10.5281/zenodo.4079188

Machine learning has emerged as a popular and powerful approach for solving problems in astrophysics. We review applications of machine learning techniques for the analysis of ground-based gravitational-wave detector data. Examples include techniques for improving the sensitivity of Advanced LIGO and Advanced Virgo gravitational-wave searches, methods for fast measurements of the astrophysical parameters of gravitational-wave sources, and algorithms for reduction and characterization of non-astrophysical detector noise. These applications demonstrate how machine learning techniques may be harnessed to enhance the science that is possible with current and future gravitational-wave detectors.

We present IMRPhenomXPHM, a phenomenological frequency-domain model for the gravitational-wave signal emitted by quasi-circular precessing binary black holes, which incorporates multipoles beyond the dominant quadrupole in the precessing frame. The model is a precessing extension of IMRPhenomXHM (Garc\'ia-Quir\'os 2020), based on approximate maps between aligned-spin waveform modes in the co-precessing frame and precessing waveform modes in the inertial frame, which is commonly referred to as "twisting up" the non-precessing waveforms. IMRPhenomXPHM includes IMRPhenomXP as a special case, the restriction to the dominant quadrupole contribution in the co-precessing frame. We implement two alternative mappings, one based on a single-spin PN approximation, as used in IMRPhenomPv2 (Hannam 2013), and one based on the double-spin MSA approach (Chatziioannou 2017). We include a detailed discussion of conventions used in the description of precessing binaries and of all choices made in constructing the model. The computational cost of \phXPHM is further reduced by extending the interpolation technique of (C. Garc\'ia-Quir\'os 2020) to the Euler angles. The accuracy, speed, robustness and modularity of the IMRPhenomX family will make these models productive tools for gravitational wave astronomy in the current era of greatly increased number and diversity of detected events., Comment: Version accepted for publication in Physical Review D

Strong gravitational lensing can produce multiple images of the same gravitational-wave signal, each arriving at different times and with different magnification. Previous work has explored if lensed pairs exist among the known high-significance events and found no evidence of this. However, the possibility remains that weaker counterparts of these events are present in the data, unrecovered by previous searches. We conduct a targeted search specifically looking for sub-threshold lensed images of known binary black hole (BBH) observations. We recover candidates matching three of the additional events first reported by Venumadhav et al. (2019), but find no evidence for additional BBH events. We also find no evidence that any of the Venumadhav et al. observations are lensed counterparts. We demonstrate how this type of counterpart search can constrain hypotheses about the overall source and lens populations and we rule out at very high confidence the extreme hypothesis that all heavy BBH detections are in fact lensed systems at high redshift with intrinsic masses < 15Msun., Comment: 16 pages. Comments welcome. Data release: https://github.com/icg-gravwaves/lensed-o1-o2-data-release. Version accepted for publication in PRD. Supplement merged into main paper following PRD style. Various updates/improvements during review process. Typos corrected

Gravitational waves (GWs) can offer a novel window into the structure and dynamics of neutron stars. Here we present the first search for long-duration quasi-monochromatic GW transients triggered by pulsar glitches. We focus on two glitches observed in radio timing of the Vela pulsar (PSR J0835-4510) on 12 December 2016 and the Crab pulsar (PSR J0534+2200) on 27 March 2017, during the Advanced LIGO second observing run (O2). We assume the GW frequency lies within a narrow band around twice the spin frequency as known from radio observatons. Using the fully-coherent transient-enabled F-statistic method to search for transients of up to four months in length. We find no credible GW candidates for either target, and through simulated signal injections we set 90% upper limits on (constant) GW strain as a function of transient duration. For the larger Vela glitch, we come close to beating an indirect upper limit for when the total energy liberated in the glitch would be emitted as GWs, thus demonstrating that similar post-glitch searches at improved detector sensitivity can soon yield physical constraints on glitch models., Comment: 10 pages, 3 figures (incl. appendices). Updated to match version accepted by PRD and fixed a few references

This paper describes a new semi-coherent method to search for transient gravitational waves of intermediate duration (hours to days). In order to search for newborn isolated neutron stars with their possibly very rapid spin-down, we model the frequency evolution as a power law. The search uses short Fourier transforms from the output of ground-based gravitational wave detectors and applies a weighted Hough transform, also taking into account the signal's amplitude evolution. We present the technical details for implementing the algorithm, its statistical properties, and a sensitivity estimate. A first example application of this method was in the search for GW170817 post-merger signals, and we verify the estimated sensitivity with simulated signals for this case., Comment: 13 pages, 14 figures

Van Putten & Della Valle (2018) have reported a possible detection of gravitational-wave 'extended emission' from a neutron star remnant of GW170817. Starting from the time-frequency evolution and total emitted energy of their reported candidate, we show that such an emission is not compatible with the current understanding of neutron stars. We explore the additional required physical assumptions to make a full waveform model, for example, taking the optimistic emission from a spining-down neutron star with fixed quadrupolar deformation, and study whether even an ideal single-template matched-filter analysis could detect an ideal, fully phase-coherent signal. We find that even in the most optimistic case an increase in energy and extreme parameters would be required for a confident detection with LIGO sensitivity as of 2018-08-17. The argument also holds for other waveform models following a similar time-frequency track and overall energy budget. Single-template matched filtering on the LIGO data around GW170817, and on data with added simulated signals, verifies the expected sensitivity scaling and the overall statistical expectation., Comment: 9 pages, 6 figures, updated version as accepted by MNRAS

Huang et al. [arXiv:1804.03101] have analysed the population of 15 known galactic Double Neutron Stars (DNSs) regarding the total masses of these systems. They suggest the existence of two sub-populations, and report likelihood-based preference for a two-component Gaussian mixture model over a single Gaussian distribution. This note offers a cautionary perspective on model selection for this data set: Especially for such a small sample size, a pure likelihood ratio test can encourage overfitting. This can be avoided by penalising models with a higher number of free parameters. Re-examining the DNS total mass data set within the class of Gaussian mixture models, this can be achieved through several simple and well-established statistical tests, including information criteria (AICc, BIC), cross-validation, Bayesian evidence ratios and a penalised EM-test. While this re-analysis confirms the basic finding that a two-component mixture is consistent with the data, the model selection criteria consistently indicate that there is no robust preference for it over a single-component fit. Additional DNS discoveries will be needed to settle the question of sub-populations., Comment: 9 pages and 10 figures including appendices, updated version as accepted by MNRAS

Gravitational waves are minute ripples in spacetime, first predicted by Einstein's general theory of relativity in 1916. Gravitational waves from rapidly-rotating neutron stars, whose shape deviates from perfect axisymmetry, are a potential astrophysical source of gravitational waves, but which so far have not been detected. The search for this type of signals, also known as continuous waves, presents a significant data analysis challenge, as their weak signatures are expected to be buried deep within the instrumental noise of the LIGO and Virgo detectors. The OctApps library provides various functions, written in Octave, intended to aid research scientists who perform searches for continuous gravitational waves., Comment: 3 pages

The F-statistic is an established method to search for continuous gravitational waves from spinning neutron stars. Prix et al. (2011) introduced a variant for transient quasi-monochromatic signals. Possible astrophysical scenarios for such transients include glitching pulsars, newborn neutron stars and accreting systems. Here we present a new implementation of the transient F-statistic, using pyCUDA to leverage the power of modern graphics processing units (GPUs). The obtained speedup allows efficient searches over much wider parameter spaces, especially when using more realistic transient signal models including time-varying (e.g. exponentially decaying) amplitudes. Hence, it can enable comprehensive coverage of glitches in known nearby pulsars, improve the follow-up of outliers from continuous-wave searches, and might be an important ingredient for future blind all-sky searches for unknown neutron stars., Comment: 13 pages, 3 figures; v2: updated reference to 1710.02327 and its erratum

For a brief moment, a binary black hole (BBH) merger can be the most powerful astrophysical event in the visible universe. Here we present a model fit for this gravitational-wave peak luminosity of nonprecessing quasicircular BBH systems as a function of the masses and spins of the component black holes, based on numerical relativity (NR) simulations and the hierarchical fitting approach introduced by X. Jim\'enez-Forteza et al. [Phys. Rev. D 95, 064024 (2017), arXiv:1611.00332]. This fit improves over previous results in accuracy and parameter-space coverage and can be used to infer posterior distributions for the peak luminosity of future astrophysical signals like GW150914 and GW151226. The model is calibrated to the l<=6 modes of 378 nonprecessing NR simulations up to mass ratios of 18 and dimensionless spin magnitudes up to 0.995, and includes unequal-spin effects. We also constrain the fit to perturbative numerical results for large mass ratios. Studies of key contributions to the uncertainty in NR peak luminosities, such as (i) mode selection, (ii) finite resolution, (iii) finite extraction radius, and (iv) different methods for converting NR waveforms to luminosity, allow us to use NR simulations from four different codes as a homogeneous calibration set. This study of systematic fits to combined NR and large-mass-ratio data, including higher modes, also paves the way for improved inspiral-merger-ringdown waveform models., Comment: Updated to match PRD version. 22 pages, 22 figures, 6 tables; data and example implementations provided as supplementary material

Numerical relativity is an essential tool in studying the coalescence of binary black holes (BBHs). It is still computationally prohibitive to cover the BBH parameter space exhaustively, making phenomenological fitting formulas for BBH waveforms and final-state properties important for practical applications. We describe a general hierarchical bottom-up fitting methodology to design and calibrate fits to numerical relativity simulations for the three-dimensional parameter space of quasicircular nonprecessing merging BBHs, spanned by mass ratio and by the individual spin components orthogonal to the orbital plane. Particular attention is paid to incorporating the extreme-mass-ratio limit and to the subdominant unequal-spin effects. As an illustration of the method, we provide two applications, to the final spin and final mass (or equivalently: radiated energy) of the remnant black hole. Fitting to 427 numerical relativity simulations, we obtain results broadly consistent with previously published fits, but improving in overall accuracy and particularly in the approach to extremal limits and for unequal-spin configurations. We also discuss the importance of data quality studies when combining simulations from diverse sources, how detailed error budgets will be necessary for further improvements of these already highly accurate fits, and how this first detailed study of unequal-spin effects helps in choosing the most informative parameters for future numerical relativity runs., Comment: 23 pages, 31 figures, 15 tables; data and example implementations as supplementary material. Updated to match PRD version

We report results of an all-sky search for periodic gravitational waves with frequency between 50 and 510 Hz from isolated compact objects, i.e. neutron stars. A new hierarchical multi-stage approach is taken, supported by the computing power of the Einstein@Home project, allowing to probe more deeply than ever before. 16 million sub-threshold candidates from the initial search [LVC,arXiv:1606.09619] are followed up in three stages. None of those candidates is consistent with an isolated gravitational wave emitter, and 90% confidence level upper limits are placed on the amplitudes of continuous waves from the target population. Between 170.5 and 171 Hz we set the most constraining 90% confidence upper limit on the strain amplitude h0 at 4.3x10-25 , while at the high end of our frequency range we achieve an upper limit of 7.6x10-25. These are the most constraining all-sky upper limits to date and constrain the ellipticity of rotating compact objects emitting at 300 Hz at a distance D to less than 6x10-7 [d/100pc]., Comment: 14 pages, 11 figures, extensive tables in the appendixes on pages 10-13

We report the results of a directed search for continuous gravitational-wave emission in a broad frequency range (between 50 and 1000 Hz) from the central compact object of the supernova remnant Cassiopeia A (Cas A). The data comes from the sixth science run of LIGO and the search is performed on the volunteer distributed computing network Einstein@Home. We find no significant signal candidate, and set the most constraining upper limits to date on the gravitational-wave emission from Cas A, which beat the indirect age-based upper limit across the entire search range. At around 170 Hz (the most sensitive frequency range), we set 90% confidence upper limits on the gravitational wave amplitude $h_0$ of $\sim\!\!~2.9\times 10^{-25}$, roughly twice as constraining as the upper limits from previous searches on Cas A. The upper limits can also be expressed as constraints on the ellipticity of Cas A; with a few reasonable assumptions, we show that at gravitational-wave frequencies greater than 300~Hz, we can exclude an ellipticity of $\gtrsim\!\!~10^{-5}$., Comment: 29 pages, 7 figures, 3 tables

The vulnerability to single-detector instrumental artifacts in standard detection methods for long-duration quasimonochromatic gravitational waves from nonaxisymmetric rotating neutron stars [continuous waves (CWs)] was addressed in past work [D. Keitel et al., Phys. Rev. D 89, 064023 (2014).] by a Bayesian approach. An explicit model of persistent single-detector disturbances led to a generalized detection statistic with improved robustness against such artifacts. Since many strong outliers in semicoherent searches of LIGO data are caused by transient disturbances that last only a few hours, we extend the noise model to cover such limited-duration disturbances, and demonstrate increased robustness in realistic simulated data. Besides long-duration CWs, neutron stars could also emit transient signals which, for a limited time, also follow the CW signal model (tCWs). As a pragmatic alternative to specialized transient searches, we demonstrate how to make standard semicoherent CW searches more sensitive to transient signals. Considering tCWs in a single segment of a semicoherent search, Bayesian model selection yields a new detection statistic that does not add significant computational cost. On simulated data, we find that it increases sensitivity towards tCWs, even of varying durations, while not sacrificing sensitivity to classical CW signals, and still being robust to transient or persistent single-detector instrumental artifacts., Comment: 16 pages, 6 figures, REVTeX4.1

The multi-detector F-statistic is close to optimal for detecting continuous gravitational waves (CWs) in Gaussian noise. However, it is susceptible to false alarms from instrumental artefacts, for example quasi-monochromatic disturbances ('lines'), which resemble a CW signal more than Gaussian noise. In a recent paper [Keitel et al 2014, PRD 89 064023], a Bayesian model selection approach was used to derive line-robust detection statistics for CW signals, generalising both the F-statistic and the F-statistic consistency veto technique and yielding improved performance in line-affected data. Here we investigate a generalisation of the assumptions made in that paper: if a CW analysis uses data from two or more detectors with very different sensitivities, the line-robust statistics could be less effective. We investigate the boundaries within which they are still safe to use, in comparison with the F-statistic. Tests using synthetic draws show that the optimally-tuned version of the original line-robust statistic remains safe in most cases of practical interest. We also explore a simple idea on further improving the detection power and safety of these statistics, which we however find to be of limited practical use., Comment: 21 pages, 11 figures, updated to match published version

The standard multidetector F-statistic for continuous gravitational waves is susceptible to false alarms from instrumental artifacts, for example monochromatic sinusoidal disturbances ('lines'). This vulnerability to line artifacts arises because the F-statistic compares the signal hypothesis to a Gaussian-noise hypothesis, and hence is triggered by anything that resembles the signal hypothesis more than Gaussian noise. Various ad-hoc veto methods to deal with such line artifacts have been proposed and used in the past. Here we develop a Bayesian framework that includes an explicit alternative hypothesis to model disturbed data. We introduce a simple line model that defines lines as signal candidates appearing only in one detector. This allows us to explicitly compute the odds between the signal hypothesis and an extended noise hypothesis, resulting in a new detection statistic that is more robust to instrumental artifacts. We present and discuss results from Monte-Carlo tests on both simulated data and on detector data from the fifth LIGO science run. We find that the line-robust statistic retains the detection power of the standard F-statistic in Gaussian noise. In the presence of line artifacts it is more sensitive, even compared to the popular F-statistic consistency veto, over which it improves by as much as a factor of two in detectable signal strength., Comment: 19 pages, 22 images in 8 figures, REVTeX4.1; updated to match published version

We study for the first time a three-dimensional octahedron constellation for a space-based gravitational wave detector, which we call the Octahedral Gravitational Observatory (OGO). With six spacecraft the constellation is able to remove laser frequency noise and acceleration disturbances from the gravitational wave signal without needing LISA-like drag-free control, thereby simplifying the payloads and placing less stringent demands on the thrusters. We generalize LISA's time-delay interferometry to displacement-noise free interferometry (DFI) by deriving a set of generators for those combinations of the data streams that cancel laser and acceleration noise. However, the three-dimensional configuration makes orbit selection complicated. So far, only a halo orbit near the Lagrangian point L1 has been found to be stable enough, and this allows only short arms up to 1400 km. We derive the sensitivity curve of OGO with this arm length, resulting in a peak sensitivity of about $2\times10^{-23}\,\mathrm{Hz}^{-1/2}$ near 100 Hz. We compare this version of OGO to the present generation of ground-based detectors and to some future detectors. We also investigate the scientific potentials of such a detector, which include observing gravitational waves from compact binary coalescences, the stochastic background and pulsars as well as the possibility to test alternative theories of gravity. We find a mediocre performance level for this short-arm-length detector, between those of initial and advanced ground-based detectors. Thus, actually building a space-based detector of this specific configuration does not seem very efficient. However, when alternative orbits that allow for longer detector arms can be found, a detector with much improved science output could be constructed using the octahedron configuration and DFI solutions demonstrated in this paper. (abridged), Comment: 17 pages, 5 figures, revtex4-1, matches version accepted by PRD

Continuous gravitational waves (CW) are expected from spinning neutron stars with non-axisymmetric deformations. A network of interferometric detectors (LIGO, Virgo and GEO600) is looking for these signals. They are predicted to be very weak and retrievable only by integration over long observation times. One of the standard methods of CW data analysis is the multi-detector F-statistic. In a typical search, the F-statistic is computed over a range in frequency, spin-down and sky position, and the candidates with highest F values are kept for further analysis. However, this detection statistic is susceptible to a class of noise artifacts, strong monochromatic lines in a single detector. By assuming an extended noise model - standard Gaussian noise plus single-detector lines - we can use a Bayesian odds ratio to derive a generalized detection statistic, the line veto (LV-) statistic. In the absence of lines, it behaves similarly to the F-statistic, but it is more robust against line artifacts. In the past, ad-hoc post-processing vetoes have been implemented in searches to remove these artifacts. Here we provide a systematic framework to develop and benchmark this class of vetoes. We present our results from testing this LV-statistic on simulated data., Comment: 2 pages, 1 figure, to be published in Proceedings of Statistical Challenges in Modern Astronomy V, Springer 2012

Context: Two-point correlation functions are used throughout cosmology as a measure for the statistics of random fields. When used in Bayesian parameter estimation, their likelihood function is usually replaced by a Gaussian approximation. However, this has been shown to be insufficient. Aims: For the case of Gaussian random fields, we search for an exact probability distribution of correlation functions, which could improve the accuracy of future data analyses. Methods: We use a fully analytic approach, first expanding the random field in its Fourier modes, and then calculating the characteristic function. Finally, we derive the probability distribution function using integration by residues. We use a numerical implementation of the full analytic formula to discuss the behaviour of this function. Results: We derive the univariate and bivariate probability distribution function of the correlation functions of a Gaussian random field, and outline how higher joint distributions could be calculated. We give the results in the form of mode expansions, but in one special case we also find a closed-form expression. We calculate the moments of the distribution and, in the univariate case, we discuss the Edgeworth expansion approximation. We also comment on the difficulties in a fast and exact numerical implementation of our results, and on possible future applications., Comment: 13 pages, 5 figures, updated to match version published in A&A (slightly expanded Sects. 5.3 and 6)