Search

Quasinormal modes (QNMs) are usually characterized by their time dependence; oscillations at specific frequencies predicted by black hole (BH) perturbation theory. QNMs are routinely identified in the ringdown of numerical relativity waveforms, are widely used in waveform modeling, and underpin key tests of general relativity and of the nature of compact objects; a program sometimes called BH spectroscopy. Perturbation theory also predicts a specific spatial shape for each QNM perturbation. For the Kerr metric, these are the ($s=-2$) spheroidal harmonics. Spatial information can be extracted from numerical relativity by fitting a feature with known time dependence to all of the spherical harmonic modes, allowing the shape of the feature to be reconstructed; a program initiated here and that we call BH cartography. Accurate spatial reconstruction requires fitting to many spherical harmonics and is demonstrated using highly accurate Cauchy-characteristic numerical relativity waveforms. The loudest QNMs are mapped, and their reconstructed shapes are found to match the spheroidal harmonic predictions. The cartographic procedure is also applied to the quadratic QNMs - nonlinear features in the ringdown - and their reconstructed shapes are compared with predictions from second-order perturbation theory. BH cartography allows us to determine the viewing angles that maximize the amplitude of the quadratic QNMs, an important guide for future searches, and is expected to lead to an improved understanding of nonlinearities in BH ringdown., Comment: 12 pages plus appendices, 14 figures

Upcoming space-based gravitational-wave detectors will be sensitive to millions and resolve tens of thousands of stellar-mass binary systems at mHz frequencies. The vast majority of these will be double white dwarfs in our Galaxy. The greatest part will remain unresolved, forming an incoherent stochastic foreground signal. Using state-of-the-art Galactic models for the formation and evolution of binary white dwarfs and accurate LISA simulated signals, we introduce a test for foreground Gaussianity and stationarity. We explain the former with a new analytical modulation induced by the LISA constellation motion and the intrinsic anisotropy of the source distribution. By demodulating the foreground signal, we reveal a deviation from Gaussianity in the 2-10 mHz frequency band. Our finding is crucial to design faithful data models, i.e. unbiased likelihoods for both individual sources and astrophysical foregrounds parameter estimation, and ultimately for an accurate interpretation of the LISA data., Comment: 17 pages, 8 figures

Searching for gravitational waves from stellar-mass binary black holes with LISA remains a challenging open problem. Conventional template-bank approaches to the search are impossible due to the prohibitive number of templates that would be required. This paper continues the development of a hierarchical semi-coherent stochastic search, extending it to a full end-to-end pipeline that is then applied to multiple mock LISA data streams which include simulated noise. Particle swarm optimization is used as a stochastic search algorithm, tracking multiple maxima of a semi-coherent search statistic defined over source parameter space. The pipeline is accelerated by the use of graphical processing units (GPUs). No prior information from observations by ground-based detectors is used; this is necessary in order to provide advance warning of the merger. We find that the pipeline is able to detect sources with signal-to-noise ratios as low as $\rho\sim 17$, we demonstrate that these searches can directly seed coarse parameter estimation in cases where the search trigger is loud. An example of how the false-alarm probability can be estimated for this type of GW search is also included., Comment: 15 pages + appendices, 7 figures, submitted to PRD

The nature of progenitors of Type Ia supernovae has long been debated, primarily due to the elusiveness of the progenitor systems to traditional electromagnetic observation methods. We argue that gravitational wave observations with the upcoming Laser Interferometer Space Antenna (LISA) offer the most promising way to test one of the leading progenitor scenarios - the double-degenerate scenario, which involves a binary system of two white dwarf stars. In this study we review published results, supplementing them with additional calculations for the context of Type Ia supernovae. We discuss the fact that LISA will be able to provide a complete sample of double white dwarf Type Ia supernova progenitors with orbital periods shorter than 16-11 minutes (gravitational wave frequencies above 2-3 millihertz). Such a sample will enable a statistical validation of the double-degenerate scenario by simply counting whether LISA detects enough double white dwarf binaries to account for the measured Type Ia merger rate in Milky Way-like galaxies. Additionally, we illustrate how LISA's capability to measure the chirp mass will set lower bounds on the primary mass, revealing whether detected double white dwarf binaries will eventually end up as a Type Ia supernova. We estimate that the expected LISA constraints on the Type Ia merger rate for the Milky Way will be 4-9%. We also discuss the potential gravitational wave signal from a Type Ia supernova assuming a double-detonation mechanism and explore how multi-messenger observations could significantly advance our understanding of these transient phenomena., Comment: accepted version

Stellar-mass binary black holes may have circumbinary disks if formed through common-envelope evolution or within gaseous environments. Disks can drive binaries into wider and more eccentric orbits, while gravitational waves harden and circularise them. We combine cutting-edge evolution prescriptions for disk-driven binaries with well-known equations for gravitational-wave-driven evolution, and study the evolution of stellar-mass binary black holes. We find that binaries are driven by their disk to an equilibrium eccentricity, $0.2\lesssim e_\mathrm{eq} \lesssim0.5$, that dominates their evolution. Once they transition to the GW-dominated regime their eccentricity decreases rapidly; we find that stellar-mass binary black holes with long-lived disks will likely be observed in LISA with detectable eccentricities $\sim 10^{-2}$ at $0.01$ Hz, with the precise value closely correlating with the binary's initial mass ratio. This may lead stellar-mass binary black holes with CBDs observed in LISA to be confused with dynamically-formed binary black holes., Comment: 7 pages, 3 figures. Accepted for publication in MNRAS Letters

Searching for gravitational-wave signals is a challenging and computationally intensive endeavor undertaken by multiple independent analysis pipelines. While detection depends only on observed noisy data, it is sometimes inconsistently defined in terms of source parameters that in reality are unknown, e.g., by placing a threshold on the optimal signal-to-noise ratio (SNR). We present a method to calibrate unphysical thresholds to search results by performing Bayesian inference on real observations using a model that simultaneously parametrizes the intrinsic network optimal SNR distribution and the effect of search sensitivity on it. We find consistency with a fourth-order power law and detection thresholds of $10.5_{-2.4}^{+2.1}$, $11.2_{-1.4}^{+1.2}$, and $9.1_{-0.5}^{+0.5}$ (medians and $90\%$ credible intervals) for events with false-alarm rates less than $1\,\mathrm{yr}^{-1}$ in the first, second, and third LIGO-Virgo-KAGRA observing runs, respectively. Though event selection can only be self-consistently reproduced by physical searches, employing our inferred thresholds allows approximate observation-calibrated selection criteria to be applied when efficiency is required and injection campaigns are infeasible., Comment: 10 pages, 5 main figures (2 appendix figures)

Detecting stochastic gravitational wave backgrounds (SGWBs) with The Laser Interferometer Space Antenna (LISA) is among the mission science objectives. Disentangling SGWBs of astrophysical and cosmological origin is a challenging task, further complicated by the noise level uncertainties. In this study, we introduce a Bayesian methodology to infer upon SGWBs, taking inspiration from Gaussian stochastic processes. We investigate the suitability of the approach for signal of unknown spectral shape. We do by discretely exploring the model hyperparameters, a first step towards a more efficient transdimensional exploration. We apply the proposed method to a representative astrophysical scenario: the inference on the astrophysical foreground of Extreme Mass Ratio Inspirals, recently estimated in~\cite{Pozzoli2023}. We find the algorithm to be capable of recovering the injected signal even with large priors, while simultaneously providing estimate of the noise level., Comment: 16 pages, 10 figures, 0 tables

LISA will detect $\sim \! 10^4$ Galactic binaries, the majority being double white dwarfs. However, approximately $\sim \! 1 \textrm{--} 5 \%$ of these systems will contain neutron stars which, if they can be correctly identified, will provide new opportunities for studying binary evolution pathways involving mass reversal and supernovae as well as being promising targets for multi-messenger observations. Eccentricity, expected from neutron star natal kicks, will be a key identifying signature for binaries containing a neutron star. Eccentric binaries radiate at widely-spaced frequency harmonics that must first be identified as originating from a single source and then analysed coherently. A multi-harmonic heterodyning approach for this type of data analysis is used to perform Bayesian parameter estimation on a range of simulated eccentric LISA signals. This is used to: (i) investigate LISA's ability to measure orbital eccentricity and to quantify the minimum detectable eccentricity; (ii) demonstrate how eccentricity and periastron precession help to break the mass degeneracy allowing the individual component masses to be inferred, potentially confirming the presence of a neutron star; (iii) investigate the possibility of source misidentification when the individual harmonics of an eccentric binary masquerade as separate circular binaries; and (iv) investigate the possibility of source reclassification, where parameter estimation results of multiple circular analyses are combined in postprocessing to quickly infer the parameters of an eccentric source. The broader implications of this for the ongoing design of the LISA global fit are also discussed., Comment: 10 pages + appendices, 9 figures, published in MNRAS

Neutron star-white dwarf (NS+WD) binaries offer a unique opportunity for studying NS-specific phenomena with gravitational waves. In this paper, we employ the binary population synthesis technique to study the Galactic population of NS+WDs with the future Laser Interferometer Space Antenna (LISA). We anticipate approximately $\mathcal{O}(10^2)$ detectable NS+WDs by LISA, encompassing both circular and eccentric binaries formed via different pathways. Despite the challenge of distinguishing NS+WDs from more prevalent double white dwarfs in the LISA data (especially at frequencies below 2 mHz), we show that their eccentricity and chirp mass distributions may provide avenues to explore the NS natal kicks and common envelope evolution. Additionally, we investigate the spatial distribution of detectable NS+WDs relative to the Galactic plane and discuss prospects for identifying electromagnetic counterparts at radio wavelengths. Our results emphasise LISA's capability to detect and characterise NS+WDs and to offer insights into the properties of the underlying population. Our conclusions carry significant implications for shaping LISA data analysis strategies and future data interpretation., Comment: Accepted for publication by MNRAS

This paper considers the problem of searching for quiet, long-duration and broadband gravitational wave signals, such as stellar-mass binary black hole binaries, in mock LISA data. We propose a method that combines a semi-coherent likelihood with the use of a particle swarm optimizer capable of efficiently exploring a large parameter space. The semi-coherent analysis is used to widen the peak of the likelihood distribution over parameter space, congealing secondary peaks and thereby assisting in localizing the posterior bulk. An iterative strategy is proposed, using particle swarm methods to initially explore a wide, loosely-coherent likelihood and then progressively constraining the signal to smaller regions in parameter space by increasing the level of coherence. The properties of the semi-coherent likelihood are first demonstrated using the well-studied binary neutron star signal GW170817. As a proof of concept, the method is then successfully applied to a simplified search for a stellar-mass binary black hole in zero-noise LISA data. Finally, we conclude by discussing what remains to be done to develop this into a fully-capable search and how the method might also be adapted to tackle the EMRI search problem in LISA., Comment: 15 pages + appendices, 9 figures. Accepted in Physical Review D on 8 September 2023

Pulsar timing arrays (PTAs) and the Laser Interferometer Space Antenna (LISA) will open complementary observational windows on massive black-hole binaries (MBHBs), i.e., with masses in the range $\sim 10^6 - 10^{10}\,$ M$_{\odot}$. While PTAs may detect a stochastic gravitational-wave background from a population of MBHBs, during operation LISA will detect individual merging MBHBs. To demonstrate the profound interplay between LISA and PTAs, we estimate the number of MBHB mergers that one can expect to observe with LISA by extrapolating direct observational constraints on the MBHB merger rate inferred from PTA data. For this, we postulate that the common signal observed by PTAs (and consistent with the increased evidence recently reported) is an astrophysical background sourced by a single MBHB population. We then constrain the LISA detection rate, $\mathcal{R}$, in the mass-redshift space by combining our Bayesian-inferred merger rate with LISA's sensitivity to spin-aligned, inspiral-merger-ringdown waveforms. Using an astrophysically-informed formation model, we predict a 95$\%$ upper limit on the detection rate of $\mathcal{R} < 134\,{\rm yr}^{-1}$ for binaries with total masses in the range $10^7 - 10^8\,$ M$_{\odot}$. For higher masses, i.e., $>10^8\,$ M$_{\odot}$, we find $\mathcal{R} < 2\,(1)\,\mathrm{yr}^{-1}$ using an astrophysically-informed (agnostic) formation model, rising to $11\,(6)\,\mathrm{yr}^{-1}$ if the LISA sensitivity bandwidth extends down to $10^{-5}$ Hz. Forecasts of LISA science potential with PTA background measurements should improve as PTAs continue their search.

In massive scalar-tensor gravity, core-collapse supernovae are strong sources of scalar-polarized gravitational waves. These can be detectable out to large distance. The dispersive nature of the propagation of waves in the massive scalar field mean the gravitational wave signals are long lived and many such signals can overlap to form a stochastic background. Using different models for the population of supernova events in the nearby universe, we compute predictions for the energy-density in the stochastic scalar-polarized gravitational wave background from core-collapse events in massive scalar-tensor gravity for theory parameters that facilitate strong scalarization. The resulting energy density is below the current constraints on a Gaussian stochastic gravitational wave background but large enough to be detectable with the current generation of detectors when they reach design sensitivity, indicating that it will soon be possible to place new constraints on the parameter space of massive scalar-tensor gravity., Comment: to match published version in Phys.Rev.D

The observation of massive black hole binary systems is one of the main science objectives of the Laser Interferometer Space Antenna (LISA). The instrument's design requirements have recently been revised: they set a requirement at $0.1\,\mathrm{mHz}$, with no additional explicit requirements at lower frequencies. This has implications for observations of the short-lived signals produced by the coalescence of massive and high-redshift binaries. Here we consider the most pessimistic scenario: the (unlikely) case in which LISA has no sensitivity below $0.1\,\mathrm{mHz}$. We show that the presence of higher multipoles (beyond the dominant $\ell = |m| = 2$ mode) in the gravitational radiation from these systems, which will be detectable with a total signal-to-noise ratio $\sim 10^3$, allows LISA to retain the capability to accurately measure the physical parameters, the redshift, and to constrain the sky location. To illustrate this point, we consider a few select binaries in a total (redshifted) mass range of $4 \times10^6 - 4 \times 10^7\,M_\odot$ whose ($\ell = |m| = 2$) gravitational-wave signals last between $\approx 12$ hours and $\approx 20$ days in band. We model the emitted gravitational radiation using the highly accurate (spin-aligned) waveform approximant IMRPhenomXHM and carry out a fully coherent Bayesian analysis on the LISA noise-orthogonal time-delay-interferometry channels., Comment: 13 pages, 6 figures, comments welcome!

Double white dwarfs (DWDs) will be the most numerous gravitational-wave (GW) sources for the Laser Interferometer Space Antenna (LISA). Most of the Galactic DWDs will be unresolved and will superpose to form a confusion noise foreground, the dominant LISA noise source around $\sim 0.5\mathrm{-}3\,\mathrm{mHz}$. A small fraction of these sources will stand out from the background and be individually detectable. Uniquely among GW sources, a handful of these binaries will be known in advance from electromagnetic (EM) observations and will be guaranteed sources of detectable GWs in the LISA band; these are known as verification binaries (VBs). High-cadence photometric surveys are continuously discovering new VB systems, and their number will continue to grow ahead of the launch of LISA. We analyse, in a fully Bayesian framework, all the currently known VB candidates with the latest design requirements for the LISA mission and find that 25 of the considered sources can be detected within a $4\,\mathrm{yr}$ observation time. We explore what can be expected from GW observations, both alone and in combination with EM observations, and estimate the VB's time to detection in the early months of LISA operations. We also show how VBs can be analysed in the case where their GW signals compete with many other unknown binary signals (both resolved and unresolved) from a realistic Galactic population of DWDs., Comment: 13 pages, 8 figures, plus appendices (accepted in Monthly Notices of the Royal Astronomical Society on 26 April 2023)

We reanalyze the GW150914 post-merger data searching for quasinormal modes beyond the fundamental, quadrupolar mode. There is currently an ongoing disagreement in the literature about whether, and to what extent, the data contains evidence for a quasinormal mode overtone. We use a frequency-domain approach to ringdown data analysis that was recently proposed by the authors. Our analysis has several advantages compared to other analyses performed mainly in the time domain; in particular, the source sky position and the ringdown start time are marginalized over (as opposed to simply being fixed) as part of a Bayesian ringdown analysis. We find tentative evidence for an overtone in GW150914, but at a lower significance than reported elsewhere. Our preferred analysis, marginalizing over the uncertainty in the time of peak strain amplitude, gives a posterior on the overtone amplitude peaked away from zero at $\sim 1.8\sigma$., Comment: 12 pages, 9 figures, plus appendices (accepted in Physical Review D on 19 July 2022)

Understanding the formation and evolution of the stellar-mass binary black holes discovered by LIGO and Virgo is a challenge that spans many areas of astrophysics, from stellar evolution, dynamics and accretion disks, to possible exotic early universe processes. Over the final years of their lives, stellar-mass binaries radiate gravitational waves that are first observable by space-based detectors (such as LISA) and then ground-based instruments (such as LIGO, Virgo and the next generation observatories Cosmic Explorer and the Einstein Telescope). Using state-of-the-art waveform models and parameter-estimation pipelines for both ground- and space-based observations, we show that (the expected handful of) these multiband observations will allow at least percent-level measurements of all 17 parameters that describe the binary, the possible identification of a likely host galaxy, and the forewarning of the merger days in advance allowing telescopes at multiple wavelengths to search for any electromagnetic signature associated to it. Multiband sources will therefore be a gold mine for astrophysics, but we also show that they could be less useful as laboratories for fundamental tests of general relativity than has been previously suggested., Comment: 21 pages, 13 figures (including supplementary material)

We propose a novel, frequency-domain approach to the analysis of the gravitational-wave ringdown signal of binary black holes and the identification of quasinormal mode frequencies of the remnant. Our approach avoids the issues of spectral leakage that would normally be expected (associated with the abrupt start of the ringdown) by modeling the inspiral and merger parts of the signal using a flexible sum of sine-Gaussian wavelets truncated at the onset of the ringdown. Performing the analysis in the frequency domain allows us to use standard (and by now well-established) Bayesian inference pipelines for gravitational wave data as well as giving us the ability to readily search over the sky position and the ringdown start time, although we find that it is necessary to use an informative prior for the latter. We test our method by using it to analyze several simulated signals with varying signal-to-noise ratios injected into two- and three-detector networks. We find that our frequency-domain approach is generally able to place tighter constraints on the remnant black-hole mass and spin than a standard time-domain analysis., Comment: 12 pages, 7 figures, plus appendix (accepted in Physical Review D on 16 November 2021)

We describe a Bayesian formalism for analyzing individual gravitational-wave events in light of the rest of an observed population. This analysis reveals how the idea of a "population-informed prior" arises naturally from a suitable marginalization of an underlying hierarchical Bayesian model which consistently accounts for selection effects. Our formalism naturally leads to the presence of "leave-one-out" distributions which include subsets of events. This differs from other approximations, also known as empirical Bayes methods, which effectively double count one or more events. We design a double-reweighting post-processing strategy that uses only existing data products to reconstruct the resulting population-informed posterior distributions. Although the correction we highlight is an important conceptual point, we find it has a limited impact on the current catalog of gravitational-wave events. Our approach further allows us to study, for the first time in the gravitational-wave literature, correlations between the parameters of individual events and those of the population., Comment: 11 pages, 5 figures + supplementary material

We present a Bayesian parameter-estimation pipeline to measure the properties of inspiralling stellar-mass black hole binaries with LISA. Our strategy (i) is based on the coherent analysis of the three noise-orthogonal LISA data streams, (ii) employs accurate and computationally efficient post-Newtonian waveforms accounting for both spin-precession and orbital eccentricity, and (iii) relies on a nested sampling algorithm for the computation of model evidences and posterior probability density functions of the full 17 parameters describing a binary. We demonstrate the performance of this approach by analyzing the LISA Data Challenge (LDC-1) dataset, consisting of 66 quasi-circular, spin-aligned binaries with signal-to-noise ratios ranging from 3 to 14 and times to merger ranging from 3000 to 2 years. We recover 22 binaries with signal-to-noise ratio higher than 8. Their chirp masses are typically measured to better than $0.02 M_\odot$ at $90\%$ confidence, while the sky-location accuracy ranges from 1 to 100 square degrees. The mass ratio and the spin parameters can only be constrained for sources that merge during the mission lifetime. In addition, we report on the successful recovery of an eccentric, spin-precessing source at signal-to-noise ratio 15 for which we can measure an eccentricity of $3\times 10^{-3}$., Comment: 15 pages, 9 figures, 2 tables (accepted to Physical Review D on 19 July 2021)

Gravitational waves (GWs) at ultra-low frequencies (${\lesssim 100\,\mathrm{nHz}}$) are key to understanding the assembly and evolution of astrophysical black hole (BH) binaries with masses $\sim 10^{6}-10^{9}\,M_\odot$ at low redshifts. These GWs also offer a unique window into a wide variety of cosmological processes. Pulsar timing arrays (PTAs) are beginning to measure this stochastic signal at $\sim 1-100\,\mathrm{nHz}$ and the combination of data from several arrays is expected to confirm a detection in the next few years. The dominant physical processes generating gravitational radiation at $\mathrm{nHz}$ frequencies are still uncertain. PTA observations alone are currently unable to distinguish a binary BH astrophysical foreground from a cosmological background due to, say, a first order phase transition at a temperature $\sim 1-100\,\mathrm{MeV}$ in a weakly-interacting dark sector. This letter explores the extent to which incorporating integrated bounds on the ultra-low frequency GW spectrum from any combination of cosmic microwave background, big bang nucleosynethesis or astrometric observations can help to break this degeneracy., Comment: 4 pages and 4 figures, plus methods sections

Gravitational-wave observations of binary black holes allow new tests of general relativity to be performed on strong, dynamical gravitational fields. These tests require accurate waveform models of the gravitational-wave signal, otherwise waveform errors can erroneously suggest evidence for new physics. Existing waveforms are generally thought to be accurate enough for current observations, and each of the events observed to date appears to be individually consistent with general relativity. In the near future, with larger gravitational-wave catalogs, it will be possible to perform more stringent tests of gravity by analyzing large numbers of events together. However, there is a danger that waveform errors can accumulate among events: even if the waveform model is accurate enough for each individual event, it can still yield erroneous evidence for new physics when applied to a large catalog. This paper presents a simple linearised analysis, in the style of a Fisher matrix calculation, that reveals the conditions under which the apparent evidence for new physics due to waveform errors grows as the catalog size increases. We estimate that, in the worst-case scenario, evidence for a deviation from general relativity might appear in some tests using a catalog containing as few as 10-30 events above a signal-to-noise ratio of 20. This is close to the size of current catalogs and highlights the need for caution when performing these sorts of experiments., Comment: 14 pages, 4 figures

Modelling the end point of binary black hole mergers is a cornerstone of modern gravitational-wave astronomy. Extracting multiple quasinormal mode frequencies from the ringdown signal allows the remnant black hole to be studied in unprecedented detail. Previous studies on numerical relativity simulations of aligned-spin binaries have found that it is possible to start the ringdown analysis much earlier than previously thought if overtones (and possibly mirror modes) are included. This increases the signal-to-noise ratio in the ringdown making identification of subdominant modes easier. In this paper we study, for the first time, black hole binaries with misaligned spins and find a much greater variation in the performance of ringdown fits than in the aligned-spin case. The inclusion of mirror modes and higher harmonics, along with overtones, improves the reliability of ringdown fits with an early start time; however, there remain cases with poor performing fits. While using overtones in conjunction with an early ringdown start time is an enticing possibility, it is necessary to proceed with caution. We also consider for the first time the use of numerical relativity surrogate models in this type of quasinormal mode study and address important questions of accuracy in the underlying numerical waveforms used for the fit., Comment: 13 pages, plus appendices, 13 figures, v3 matches published version, see accompanying resources at https://doi.org/10.5281/zenodo.4538194

White dwarf stars are a well-established tool for studying Galactic stellar populations. Two white dwarfs in a tight binary system offer us an additional messenger - gravitational waves - for exploring the Milky Way and its immediate surroundings. Gravitational waves produced by double white dwarf (DWD) binaries can be detected by the future Laser Interferometer Space Antenna (LISA). Numerous and widespread DWDs have the potential to probe shapes, masses and formation histories of the stellar populations in the Galactic neighbourhood. In this work we outline a method for estimating the total stellar mass of Milky Way satellite galaxies based on the number of DWDs detected by LISA. To constrain the mass we perform a Bayesian inference using binary population synthesis models and considering the number of detected DWDs associated with the satellite and the measured distance to the satellite as the only inputs. Using a fiducial binary population synthesis model we find that for large satellites the stellar masses can be recovered to within 1) a factor two if the star formation history is known and 2) an order of magnitude when marginalising over different star formation history models. For smaller satellites we can place upper limits on their stellar mass. Gravitational wave observations can provide mass measurements for large satellites that are comparable, and in some cases more precise, than standard electromagnetic observations., Comment: published in MNRAS Letters

Gravitational wave (GW) transients from binary neutron star (BNS) coalescences can, in principle, be subject to gravitational lensing thereby increasing the amplitude and signal-to-noise ratio. We estimate the rate of lensed BNS events resolvable by LIGO and Virgo and find that it is constrained by the current non-detection of a stochastic GW background. Following closely the formalism we developed (10.1103/PhysRevLett.125.141102) in the context of binary black hole lensing, we show that at current sensitivities the fraction of BNS coalescences with lensing magnifications $\mu> 1.02$ is less than $\sim 7\times 10^{-8}$ and therefore such events should not be expected in the near future. We also make predictions for projected future sensitivities., Comment: 4 pages, 2 figures, published in Physical Review D (22 October 2020)

We compute families of spherically symmetric neutron-star models in two-derivative scalar-tensor theories of gravity with a massive scalar field. The numerical approach we present allows us to compute the resulting spacetimes out to infinite radius using a relaxation algorithm on a compactified grid. We discuss the structure of the weakly and strongly scalarized branches of neutron-star models thus obtained and their dependence on the linear and quadratic coupling parameters $\alpha_0$, $\beta_0$ between the scalar and tensor sectors of the theory, as well as the scalar mass $\mu$. For highly negative values of $\beta_0$, we encounter configurations resembling a "gravitational atom", consisting of a highly compact baryon star surrounded by a scalar cloud. A stability analysis based on binding-energ calculations suggests that these configurations are unstable and we expect them to migrate to models with radially decreasing baryon density {\it and} scalar field strength., Comment: Corrected Eq.(6)

Gravitational waves (GWs) are subject to gravitational lensing in the same way as electromagnetic radiation. However, to date, no unequivocal observation of a lensed GW transient has been reported. Independently, GW observatories continue to search for the stochastic GW signal which is produced by many transient events at high redshift. We exploit a surprising connection between the lensing of individual transients and limits to the background radiation produced by the unresolved population of binary back hole mergers: we show that it constrains the fraction of individually resolvable lensed binary black holes to less than $\sim 4\times 10^{-5}$ at present sensitivity. We clarify the interpretation of existing, low redshift GW observations (obtained assuming no lensing) in terms of their apparent lensed redshifts and masses and explore constraints from GW observatories at future sensitivity. Based on our results, recent claims of observations of lensed events are statistically disfavoured., Comment: 9 pages, 6 figures, includes Supplemental material, published in Physical Review Letters

This paper provides an extended exploration of the inverse-chirp gravitational-wave signals from stellar collapse in massive scalar-tensor gravity reported in [Phys. Rev. Lett. {\bf 119}, 201103]. We systematically explore the parameter space that characterizes the progenitor stars, the equation of state and the scalar-tensor theory of the core collapse events. We identify a remarkably simple and straightforward classification scheme of the resulting collapse events. For any given set of parameters, the collapse leads to one of three end states, a weakly scalarized neutron star, a strongly scalarized neutron star or a black hole, possibly formed in multiple stages. The latter two end states can lead to strong gravitational-wave signals that may be detectable in present continuous-wave searches with ground-based detectors. We identify a very sharp boundary in the parameter space that separates events with strong gravitational-wave emission from those with negligible radiation., Comment: 26 pages, 19 figures, to match published version in Phys.Rev.D

The population of Milky Way satellite galaxies is of great interest for cosmology, fundamental physics, and astrophysics. They represent the faint end of the galaxy luminosity function, are the most dark-matter dominated objects in the local Universe, and contain the oldest and most metal-poor stellar populations. Recent surveys have revealed around 60 satellites, but this could represent less than half of the total. Characterization of these systems remains a challenge due to their low luminosity. We consider the gravitational wave observatory LISA as a potential tool for studying these satellites through observations of their short-period double white dwarf populations. LISA will observe the entire sky without selection effects due to dust extinction, complementing optical surveys, and could potentially discover massive satellites hidden behind the disk of the galaxy., Comment: 7 pages, 2 figures

A number of open problems hinder our present ability to extract scientific information from data that will be gathered by the near-future gravitational-wave mission LISA. Many of these relate to the modeling, detection and characterization of signals from binary inspirals with an extreme component-mass ratio of $\lesssim10^{-4}$. In this paper, we draw attention to the issue of systematic error in parameter estimation due to the use of fast but approximate waveform models; this is found to be relevant for extreme-mass-ratio inspirals even in the case of waveforms with $\gtrsim90\%$ overlap accuracy and moderate ($\gtrsim30$) signal-to-noise ratios. A scheme that uses Gaussian processes to interpolate and marginalize over waveform error is adapted and investigated as a possible precursor solution to this problem. Several new methodological results are obtained, and the viability of the technique is successfully demonstrated on a three-parameter example in the setting of the LISA Data Challenge., Comment: Published version

A passing gravitational wave causes a deflection in the apparent astrometric positions of distant stars. The effect of the speed of the gravitational wave on this astrometric shift is discussed. A stochastic background of gravitational waves would result in a pattern of astrometric deflections which are correlated on large angular scales. These correlations are quantified and investigated for backgrounds of gravitational waves with sub- and super-luminal group velocities. The statistical properties of the correlations are depicted in two equivalent and related ways: as correlation curves and as angular power spectra. Sub-(super-)luminal gravitational wave backgrounds have the effect of enhancing (suppressing) the power in low-order angular modes. Analytical representations of the redshift-redshift and redshift-astrometry correlations are also derived. The potential for using this effect for constraining the speed of gravity is discussed., Comment: 14 pages, 16 figures, accepted for publication in Physical Review D

Gravitational-wave (GW) observations provide a wealth of information on the nature and properties of black holes. Among these, tidal Love numbers or the multipole moments of the inspiralling and final objects are key to a number of constraints. Here, we consider these observations in the context of higher-dimensional scenarios, with flat large extra dimensions. We show that -- as might be anticipated, but not always appreciated in the literature -- physically motivated set-ups are unconstrained by gravitational-wave data. Dynamical processes that do not excite the Kaluza-Klein (KK) modes lead to a signal identical to that in four-dimensional general relativity in vacuum . In addition, any possible excitation of the KK modes is highly suppressed relative to the dominant quadrupolar term; given existing constraints on the extra dimensions and the masses of the objects seen in gravitational-wave observations, KK modes appear at post-Newtonian order $\sim 10^{11}$. Finally, we re-compute the tidal Love numbers of spherical black holes in higher dimensions. We confirm that these are different from zero, but comparing with previous computations we find a different magnitude and sign., Comment: 8 pages, RevTeX 4. v2: Minor corrections, version to appear in Physical Review D

The label switching problem arises in the Bayesian analysis of models containing multiple indistinguishable parameters with arbitrary ordering. Any permutation of these parameters is equivalent, therefore models with many such parameters have extremely multi-modal posterior distributions. It is difficult to sample efficiently from such posteriors. This paper discusses a solution to this problem which involves carefully mapping the input parameter space to a high dimensional hypertriangle. It is demonstrated that this solution is efficient even for large numbers of parameters and can be easily applied alongside any stochastic sampling algorithm. This method is illustrated using two example problems from the field of gravitational wave astronomy., Comment: Accepted to PRD, 8 pages, 4 figures plus Supplement. Data behind figures at: 10.5281/zenodo.3351629

The progenitors of the high-mass black-hole mergers observed by LIGO and Virgo are potential LISA sources and promising candidates for multiband GW observations. In this letter, we consider the minimum signal-to-noise ratio these sources must have to be detected by LISA. Our revised threshold of $\rho_{\rm thr}\sim 15$ is higher than previous estimates, which significantly reduces the expected number of events. We also point out the importance of the detector performance at high-frequencies and the duration of the LISA mission, which both influence the event rate substantially., Comment: 6 pages, 3 figures. Published in MNRAS letters. DOI 10.1093/mnrasl/slz104

We study how the gravitational wave signal from stellar collapse in scalar-tensor gravity varies under the influence of scalar self-interaction. To this end, we extract the gravitational radiation from numerical simulations of stellar collapse for a range of potentials with higher-order terms in addition to the quadratic mass term. Our study includes collapse to neutron stars and black holes and we find the strong inverse-chirp signals obtained for the purely quadratic potential to be exceptionally robust under changes in the potential at higher orders; quartic and sextic terms in the potential lead to noticeable differences in the wave signal only if their contribution is amplified, implying a relative fine-tuning to within 5 or more orders of magnitude between the mass and self-interaction parameters., Comment: 15 pages, 5 figures; submitted to CQG

Gravitational wave astrophysics has only just begun, and as current detectors are upgraded and new detectors are built, many new, albeit faint, features in the signals will become accessible. One such feature is the presence of time-dependent Doppler shifts, generated by the acceleration of the center of mass of the gravitational-wave emitting system. We here develop a generic method that takes a frequency-domain, gravitational-wave model devoid of Doppler shifts and introduces modifications that incorporate them. Building upon a perturbative expansion that assumes the Doppler-shift velocity is small relative to the speed of light, the method consists of the inclusion of a single term in the Fourier phase and two terms in the Fourier amplitude. We validate the method through matches between waveforms with a Doppler shift in the time domain and waveforms constructed with our method for two toy problems: constant accelerations induced by a distant third body and Gaussian accelerations that resemble a kick profile. We find mismatches below $\sim\!10^{-6}$ for all of the astrophysically relevant cases considered, and improve further at smaller velocities. The work presented here will allow for the use of future detectors to extract new, faint features in the signal from the noise., Comment: 11 pages, 5 figures, submitted to Phys. Rev. D

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress., Comment: White Paper for the COST action "Gravitational Waves, Black Holes, and Fundamental Physics", 272 pages, 12 figures; v4: updated references and author list. Overall improvements and corrections. To appear in Classical and Quantum Gravity

The Gaia mission offers a new opportunity to search for the low frequency gravitational wave background using astrometric measurements. In this paper, the astrometric effect of gravitational waves is reviewed, with a particular focus on the effect of non-Einsteinian gravitational wave polarizations. A stochastic gravitational wave background generates a correlated vector field of astrometric deflections on the sky. A convenient decomposition for the correlation matrix is introduced, enabling it to be calculated for all possible gravitational wave polarizations and compared to the redshift correlations from the pulsar-timing literature; in the case of a GR background of transverse traceless gravitational waves, this also allows us to identify an astrometric analog of the famous Hellings-Downs curve. Finally, the cross-correlation between the redshift and astrometric signal is also calculated; this may form the basis for future joint pulsar-timing and astrometry searches for arbitrarily polarized gravitational wave backgrounds., Comment: 14 pages, 12 figures

Extreme-mass-ratio-inspiral observations from future space-based gravitational-wave detectors such as LISA will enable strong-field tests of general relativity with unprecedented precision, but at prohibitive computational cost if existing statistical techniques are used. In one such test that is currently employed for LIGO black-hole binary mergers, generic deviations from relativity are represented by $N$ deformation parameters in a generalised waveform model; the Bayesian evidence for each of its $2^N$ combinatorial submodels is then combined into a posterior odds ratio for modified gravity over relativity in a null-hypothesis test. We adapt and apply this test to a generalised model for extreme-mass-ratio inspirals constructed on deformed black-hole spacetimes, and focus our investigation on how computational efficiency can be increased through an evidence-free method of model selection. This method is akin to the algorithm known as product-space Markov chain Monte Carlo, but uses nested sampling and improved error estimates from a rethreading technique. We perform benchmarking and robustness checks for the method, and find order-of-magnitude computational gains over regular nested sampling in the case of synthetic data generated from the null model., Comment: Published version

The quest for fundamental limitations on physical processes is old and venerable. Here, we investigate the maximum possible power, or luminosity, that any event can produce. We show, via full nonlinear simulations of Einstein's equations, that there exist initial conditions which give rise to arbitrarily large luminosities. However, the requirement that there is no past horizon in the spacetime seems to limit the luminosity to below the Planck value, ${{\cal L}_\textrm{P}\!=\!c^5/G}$. Numerical relativity simulations of critical collapse yield the largest luminosities observed to date, ${\approx \! 0.2 {\cal L}_\textrm{P}}$. We also present an analytic solution to the Einstein equations which seems to give an unboundedly large luminosity; this will guide future numerical efforts to investigate super-Planckian luminosities., Comment: 10 pages (including appendices), 4 figures, 1 table, submitted to Phys. Rev. D

This letter considers stellar core collapse in massive scalar--tensor theories of gravity. The presence of a mass term for the scalar field allows for dramatic increases in the radiated gravitational wave signal. There are several potential smoking gun signatures of a departure from general relativity associated with this process. These signatures could show up within existing LIGO--Virgo searches., Comment: 6 pages, 3 figures; to match version published in Phys.Rev.Lett

Gravitational waves (GWs) cause the apparent position of distant stars to oscillate with a characteristic pattern on the sky. Astrometric measurements (e.g. those made by Gaia) therefore provide a new way to search for GWs. The main difficulty facing such a search is the large size of the data set; Gaia observes more than one billion stars. In this letter the problem of searching for GWs from individually resolvable supermassive black hole binaries using astrometry is addressed for the first time; it is demonstrated how the data set can be compressed by a factor of more than $10^6$, with a loss of sensitivity of less than $1\%$. This technique is successfully used to recover artificially injected GWs from mock Gaia data. Repeated injections are used to calculate the sensitivity of Gaia as a function of frequency, and Gaia's directional sensitivity variation, or antenna pattern. Throughout the letter the complementarity of Gaia and pulsar timing searches for GWs is highlighted., Comment: 6 pages, 5 figures

The space based interferometer LISA will be capable of detecting the gravitational waves emitted by stellar mass black holes or neutron stars slowly inspiralling into the supermassive black holes found in the centre of most galaxies. The gravitational wave signal from such an extreme mass ratio inspiral (EMRI) event will provide a unique opportunity to test whether the spacetime metric around the central black hole is well described by the Kerr solution. In this paper a well studied model for EMRIs around Kerr black holes is extended to a family of parametrically deformed bumpy black holes which preserve the basic symmetries of the Kerr metric. The new EMRI model is then used to quantify the constraints that LISA observations of EMRIs may be able to place on the deviations, or bumps, on the Kerr metric., Comment: 20 pages, 6 figures

The extreme-mass-ratio inspirals (EMRIs) of stellar-mass compact objects into massive black holes are an important class of source for the future space-based gravitational-wave detector LISA. Detecting signals from EMRIs will require waveform models that are both accurate and computationally efficient. In this paper, we present the latest implementation of an augmented analytic kludge (AAK) model, publicly available at github.com/alvincjk/EMRI_Kludge_Suite as part of an EMRI waveform software suite. This version of the AAK model has improved accuracy compared to its predecessors, with two-month waveform overlaps against a more accurate fiducial model exceeding 0.97 for a generic range of sources; it also generates waveforms 5-15 times faster than the fiducial model. The AAK model is well suited for scoping out data analysis issues in the upcoming round of mock LISA data challenges. A simple analytic argument shows that it might even be viable for detecting EMRIs with LISA through a semi-coherent template bank method, while the use of the original analytic kludge in the same approach will result in around 90% fewer detections., Comment: Published version