Search

Extreme Mass Ratio Inspirals (EMRIs) are key sources for the future space-based gravitational wave detector LISA, and are considered promising probes of fundamental physics. Here, we present the first complete Bayesian analysis of EMRI signals in theories with an additional massless scalar, which could arise in an extension of General Relativity or of the Standard Model of Particle Physics. We develop a waveform model accurate at adiabatic order for equatorial eccentric orbits around spinning black holes. Using full Bayesian inference, we forecast LISA's ability to probe the presence of new fundamental fields with EMRI observations.

General relativity (GR) has proven to be a highly successful theory of gravity since its inception. The theory has thrivingly passed numerous experimental tests, predominantly in weak gravity, low relative speeds, and linear regimes, but also in the strong-field and very low-speed regimes with binary pulsars. Observable gravitational waves (GWs) originate from regions of spacetime where gravity is extremely strong, making them a unique tool for testing GR, in previously inaccessible regions of large curvature, relativistic speeds, and strong gravity. Since their first detection, GWs have been extensively used to test GR, but no deviations have been found so far. Given GR's tremendous success in explaining current astronomical observations and laboratory experiments, accepting any deviation from it requires a very high level of statistical confidence and consistency of the deviation across GW sources. In this paper, we compile a comprehensive list of potential causes that can lead to a false identification of a GR violation in standard tests of GR on data from current and future ground-based GW detectors. These causes include detector noise, signal overlaps, gaps in the data, detector calibration, source model inaccuracy, missing physics in the source and in the underlying environment model, source misidentification, and mismodeling of the astrophysical population. We also provide a rough estimate of when each of these causes will become important for tests of GR for different detector sensitivities. We argue that each of these causes should be thoroughly investigated, quantified, and ruled out before claiming a GR violation in GW observations., Comment: Review article; 1 figure; 1 table; comments welcome

Phase transitions can play an important role in the cosmological constant problem, allowing the underlying vacuum energy, and therefore the value of the cosmological constant, to change. Deep within the core of neutron stars, the local pressure may be sufficiently high to trigger the QCD phase transition, thus generating a shift in the value of the cosmological constant. The gravitational effects of such a transition should then be imprinted on the properties of the star. Working in the framework of General Relativity, we provide a new model of the stellar interior, allowing for a QCD and a vacuum energy phase transition. We determine the impact of a vacuum energy jump on mass-radius relations, tidal deformability-radius relations, I-Love-Q relations and on the combined tidal deformability measured in neutron star binaries., Comment: 21 pages, 15 figures

The Lunar Gravitational-wave Antenna (LGWA) is a proposed array of next-generation inertial sensors to monitor the response of the Moon to gravitational waves (GWs). Given the size of the Moon and the expected noise produced by the lunar seismic background, the LGWA would be able to observe GWs from about 1 mHz to 1 Hz. This would make the LGWA the missing link between space-borne detectors like LISA with peak sensitivities around a few millihertz and proposed future terrestrial detectors like Einstein Telescope or Cosmic Explorer. In this article, we provide a first comprehensive analysis of the LGWA science case including its multi-messenger aspects and lunar science with LGWA data. We also describe the scientific analyses of the Moon required to plan the LGWA mission.

Gravitational wave observations can test the validity of General Relativity (GR) in the strong field regime. Certain classes of scalar-tensor theories indeed predict that compact objects can exhibit significant deviations from their GR counterparts. Here we explore the quasinormal modes of axial perturbations in spherically symmetric black holes in three such classes: (i) dilatonic black holes with an additional scalar-Ricci coupling (EdRGB), (ii) spontaneously scalarized black holes (EsRGB) with a quadratic coupling to the Gauss-Bonnet invariant and the Ricci scalar, (iii) spontaneously scalarized black holes with a quadratic and a quartic coupling to the Gauss-Bonnet invariant., Comment: 12 pages, 7 figures

In the aftermath of a binary black hole merger event, the gravitational wave signal emitted by the remnant black hole is modeled as a superposition of damped sinusoids known as quasinormal modes. While the dominant quasinormal modes originating from linear black hole perturbation theory have been studied extensively in this post-merger "ringdown" phase, more accurate models of ringdown radiation include the nonlinear modes arising from higher-order perturbations of the remnant black hole spacetime. We explore the detectability of quadratic quasinormal modes with both ground- and space-based next-generation detectors. We quantify how predictions of the quadratic mode detectability depend on the quasinormal mode starting times. We then calculate the signal-to-noise ratio of quadratic modes for several detectors and binary black hole populations, focusing on the ($220\times220$) mode - i.e., on the quadratic term sourced by the square of the linear $(220)$ mode. For the events with the loudest quadratic mode signal-to-noise ratios, we additionally compute statistical errors on the mode parameters in order to further ascertain the distinguishability of the quadratic mode from the linear quasinormal modes. The astrophysical models used in this paper suggest that while the quadratic mode may be detectable in at most a few events with ground-based detectors, the prospects for detection with the Laser Interferometer Space Antenna (LISA) are more optimistic., Comment: 21 pages, 16 figures, 2 tables; revised to match published version; new plots in App. D

We develop a fully analytical waveform model for precessing binaries with arbitrary spin vectors using post-Newtonian~(PN) theory in the extreme mass-ratio limit and a hierarchical multi-scale analysis. The analytical model incorporates leading PN order spin precession dynamics from spin-orbit, spin-spin, and quadrupole-monopole couplings, and 2PN order dissipative dynamics truncated to first order in the mass ratio $q \ll 1$. Due to the pure analytic nature of the model, the framework developed herein can readily be extended to both higher PN and higher-$q$ order. Although the PN series is asymptotic to this limit, our results can be used to estimate how precession affects the measurability of certain binary parameters, and to inform and compare with other waveform approximants, such as effective-one-body models, hybrid waveforms, and self-force calculations., Comment: 29 pages, 10 figures, published version

The Laser Interferometer Space Antenna (LISA) is the first scientific endeavour to detect and study gravitational waves from space. LISA will survey the sky for Gravitational Waves in the 0.1 mHz to 1 Hz frequency band which will enable the study of a vast number of objects ranging from Galactic binaries and stellar mass black holes in the Milky Way, to distant massive black-hole mergers and the expansion of the Universe. This definition study report, or Red Book, presents a summary of the very large body of work that has been undertaken on the LISA mission over the LISA definition phase., Comment: 155 pages, with executive summary and table of contents

We develop a numerical approach to compute polar parity perturbations within fully relativistic models of black hole systems embedded in generic, spherically symmetric, anisotropic fluids. We apply this framework to study gravitational wave generation and propagation from extreme mass-ratio inspirals in the presence of several astrophysically relevant dark matter models, namely the Hernquist, Navarro-Frenk-White, and Einasto profiles. We also study dark matter spike profiles obtained from a fully relativistic calculation of the adiabatic growth of a BH within the Hernquist profile, and provide a closed-form analytic fit of these profiles. Our analysis completes prior numerical work in the axial sector, yielding a fully numerical pipeline to study black hole environmental effects. We study the dependence of the fluxes on the DM halo mass and compactness. We find that, unlike the axial case, polar fluxes are not adequately described by simple gravitational-redshift effects, thus offering an exciting avenue for the study of black hole environments with gravitational waves., Comment: 11 pages, 5 figures, 1 table

Black hole spectroscopy is a clean and powerful tool to test gravity in the strong-field regime and to probe the nature of compact objects. Next-generation ground-based detectors, such as the Einstein Telescope and Cosmic Explorer, will observe thousands of binary black hole mergers with large signal-to-noise ratios, allowing for accurate measurements of the remnant black hole quasinormal mode frequencies and damping times. In previous work we developed an observable-based parametrization of the quasinormal mode spectrum of spinning black holes beyond general relativity (ParSpec). In this paper we use this parametrization to ask: can next-generation detectors detect or constrain deviations from the Kerr spectrum by stacking multiple observations of binary mergers from astrophysically motivated populations? We focus on two families of tests: (i) agnostic (null) tests, and (ii) theory-based tests, which make use of quasinormal frequency calculations in specific modified theories of gravity. We consider in particular two quadratic gravity theories (Einstein-scalar-Gauss-Bonnet and dynamical Chern-Simons gravity) and various effective field theory-based extensions of general relativity. We find that robust inference of hypothetical corrections to general relativity requires pushing the slow-rotation expansion to high orders. Even when high-order expansions are available, ringdown observations alone may not be sufficient to measure deviations from the Kerr spectrum for theories with dimensionful coupling constants. This is because the constraints are dominated by "light" black hole remnants, and only few of them have sufficiently high signal-to-noise ratio in the ringdown. Black hole spectroscopy with next-generation detectors may be able to set tight constraints on theories with dimensionless coupling, as long as we assume prior knowledge of the mass and spin of the remnant black hole., Comment: references updated

These notes cover part of the lectures presented by Andrea Maselli for the 59th Winter School of Theoretical Physics and third COST Action CA18108 Training School 'Gravity -- Classical, Quantum and Phenomenology'. The school took place at Palac Wojan\'ow, Poland, from February 12th to 21st, 2023. The lectures focused on some key aspects of black hole physics, and in particular on the dynamics of particles and on the scattering of waves in the Schwarzschild spacetime. The goal of the course was to introduce the students to the concept of black hole quasi normal modes, to discuss their properties, their connection with the geodesic motion of massless particles, and to provide numerical approaches to compute their actual values., Comment: Lectures at Third Training School of COST Action "Quantum gravity phenomenology in the multi-messenger approach", to be published in PoS. 39 pages, 22 figures

LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmological distances; from the beginnings of inspirals that will venture into the ground-based detectors' view to the death spiral of compact objects into massive black holes, and many sources in between. Central to realising LISA's discovery potential are waveform models, the theoretical and phenomenological predictions of the pattern of gravitational waves that these sources emit. This white paper is presented on behalf of the Waveform Working Group for the LISA Consortium. It provides a review of the current state of waveform models for LISA sources, and describes the significant challenges that must yet be overcome., Comment: 239 pages, 11 figures, white paper from the LISA Consortium Waveform Working Group, invited for submission to Living Reviews in Relativity, updated with comments from community

Using the self-force approach, we present the premier first-post-adiabatic accuracy formalism for modelling compact binaries in theories with a massless scalar field non-minimally coupled to gravity. We limit the binary secondary to being a non-spinning compact body with no scalar dipole (we will address the spinning and scalar dipole cases in an upcoming paper). By producing an ansatz for the scalar charged point particle action, we derive first- and second-order perturbative field equations and equations of motion for the secondary compact object. Under our assumptions, implementing this formalism will produce sufficiently accurate waveform templates for precision measurements of the scalar charge of the secondary with LISA data on extreme-mass-ratio inspirals. Our formalism is consistent with almost general scalar-tensor theories of gravity. Implementing our formalism builds on self-force models in General Relativity; we show the incorporation into the two-timescale formalism is straightforward. Excitingly, implementation poses no significantly more challenging barriers than computing first-post adiabatic waveforms in General Relativity., Comment: 15 pages, 0 figures

The multipole moments of black holes in general relativity obey certain consistency relations known as the no-hair theorems. The details of this multipolar structure are imprinted into the gravitational waves emitted by binary black holes, particularly if the binary is precessing. If black holes do not obey the vacuum field equations of general relativity, then the no-hair theorems may be broken, and the observed gravitational waves will be modified, thus providing an important test of the no-hair theorems. Recently, analytic solutions to the precession dynamics and inspiral waveforms were computed within the context of binaries possessing non-axisymmetric mass quadrupole moments, which are parametrized by a modulus $q_{m}$ and phase $a_{m}$ with $m = 1,2$ the azimuthal spherical harmonic number. Here, we use a Fisher analysis to study plausible constraints one may obtain on generic, non-axisymmetry quadrupole configurations using current and future ground-based detectors. For non-precessing binaries, we generically find that no meaningful constraints can be placed with current detectors on the non-axisymmetry parameters $(q_{m}, a_{m})$ due to the presence of strong degeneracies with other waveform parameters, while with next generation detectors, only weak constraints are possible. For precessing configurations, the exact value of the uncertainty is strongly dependent on the sky location, system orientation relative to the line of sight, and initial inclination angle of the orbital angular momentum. After averaging over these parameters, we find that with GWTC-3-like events, one should be able to plausibly constraint non-axisymmetric mass quadrupole deviations to $\Delta q_{m} \sim 10^{-2}$ for LIGO at design sensitivity, and $\Delta q_{m} \sim 10^{-4}$ for the same sources with Einstein Telescope and Cosmic Explorer., Comment: 15 pages, 2 figures, published version

Next-generation ground-based gravitational-wave observatories such as the Einstein Telescope and Cosmic Explorer will detect $O(10^{5}-10^{6})$ signals from compact binary coalescences every year, the exact number depending on uncertainties in the binary merger rate. Several overlapping signals will be present in band at any given time, generating a confusion noise background. We study how this confusion noise affects constraints on possible deviations from general relativity induced by modified gravity and environmental effects. Confusion noise impacts only the signals that last longer in band. Even for a "golden" GW170817-like signal, the constraints broaden by a factor in the range $[10\%,40\%]$ $([70\%,110\%])$ for the fiducial (highest) value of the local binary neutron star merger rate. Our ability to test general relativity or constrain environmental effects will be limited by systematic errors, and not by confusion noise., Comment: 11 pages, 5 figures

We study extreme-mass-ratio systems in theories admitting the Schwarzschild solution and propagating a massive graviton. We show that, in addition to small corrections to the quadrupolar and higher-order modes, a dipolar mode is excited in these theories and we quantify its excitation. While LIGO-Virgo-KAGRA observations are not expected to impose meaningful constraints in the dipolar sector, future observations by the Einstein Telescope or by LISA, together with bounds from dispersion relations, can rule out theories of massive gravity admitting vacuum General Relativistic backgrounds. For the bound to be circumvented, one needs to move away from Ricci-flat solutions, and enter a territory where constraints based on wave propagation and dispersion relations are not reliable., Comment: 5 pages plus Supplemental Material

The Einstein Telescope (ET), the European project for a third-generation gravitational-wave detector, has a reference configuration based on a triangular shape consisting of three nested detectors with 10 km arms, where in each arm there is a `xylophone' configuration made of an interferometer tuned toward high frequencies, and an interferometer tuned toward low frequencies and working at cryogenic temperature. Here, we examine the scientific perspectives under possible variations of this reference design. We perform a detailed evaluation of the science case for a single triangular geometry observatory, and we compare it with the results obtained for a network of two L-shaped detectors (either parallel or misaligned) located in Europe, considering different choices of arm-length for both the triangle and the 2L geometries. We also study how the science output changes in the absence of the low-frequency instrument, both for the triangle and the 2L configurations. We examine a broad class of simple `metrics' that quantify the science output, related to compact binary coalescences, multi-messenger astronomy and stochastic backgrounds, and we then examine the impact of different detector designs on a more specific set of scientific objectives., Comment: 197 pages, 73 figures. v2: corrections in the part on the sensitivity to stochastic backgrounds. Accepted in JCAP

We develop a numerical approach to find asymptotically flat black hole solutions coupled to anisotropic fluids, described by generic density profiles. Our model allows for a variety of applications in realistic astrophysical scenarios, and is potentially able to describe the geometry of galaxies hosting supermassive black holes, dark matter environments and accretion phenomena. We apply our framework to a black hole surrounded by different families of dark matter profiles, namely the Hernquist, the Navarro-Frenk White and the Einasto models. We study the geodesic motion of light and of massive particles in such spacetimes. Moreover we compute gravitational axial perturbations induced by a small secondary on the numerical background, and determine the changes in the emitted gravitational wave fluxes compared to the vacuum case. Our analysis confirms and extend previous studies showing that modifications of orbital frequencies and axial fluxes can be described in terms of gravitational-redshift, regardless of the halo model., Comment: 8 pages, 5 figures

Compact boson star binaries are hypothetical sources for ground-based and space gravitational-wave detectors. Their signal would be a messenger for novel fundamental fields and could shed light on the dark matter. In this work, we further develop our analysis in Phys. Rev. D 102, 083002 (2020), aimed at constraining the properties of these objects with future observations. We use a coherent waveform template for the inspiral stage of boson star binaries with large quartic self interactions, including tidal deformability and the nonlinear dependence of the quadrupole moments on the spin in terms of the fundamental couplings of the scalar field theory. Performing a Bayesian analysis, we investigate the ability of a third-generation gravitational-wave detector such as the Einstein Telescope to distinguish these exotic sources from black holes and infer constraints on the fundamental couplings of the model.

We study the imprint of light scalar fields on gravitational waves from extreme mass ratio inspirals -- binary systems with a very large mass asymmetry. We first show that, to leading order in the mass ratio, any effects of the scalar on the waveform are captured fully by two parameters: the mass of the scalar and the scalar charge of the secondary compact object. We then use this theory-agnostic framework to show that the future observations by LISA will be able to simultaneously measure both of these parameters with enough accuracy to detect ultra-light scalars., Comment: 11 pages, 7 figures

Motivated by events in which black holes can lose their environment due to tidal interactions in a binary system, we develop a waveform model in which the tidal deformability interpolates between a finite value (dressed black hole) at relatively low frequency and a zero value (naked black hole) at high frequency. We then apply this model to the example case of a black hole dressed with an ultralight scalar field and investigate the detectability of the tidal Love number with the Einstein Telescope. We show that the parameters of the tidal deformability model could be measured with high accuracy, providing a useful tool to understand dynamical environmental effects taking place during the inspiral of a binary system., Comment: 7 pages, 3 figures. v2: matching published version

We establish a generic, fully-relativistic formalism to study gravitational-wave emission by extreme-mass-ratio systems in spherically-symmetric, non-vacuum black-hole spacetimes. The potential applications to astrophysical setups range from black holes accreting baryonic matter to those within axionic clouds and dark matter environments, allowing to assess the impact of the galactic potential, of accretion, gravitational drag and halo feedback on the generation and propagation of gravitational-waves. We apply our methods to a black hole within a halo of matter. We find fluid modes imparted to the gravitational-wave signal (a clear evidence of the black hole fundamental mode instability) and the tantalizing possibility to infer galactic properties from gravitational-wave measurements by sensitive, low-frequency detectors., Comment: v2. Supplemental Material added, minor improvements. Accepted for publication in Physical Review Letters

We estimate the accuracy in the measurement of the tidal Love number of a supermassive compact object through the detection of an extreme mass ratio inspiral~(EMRI) by the future LISA mission. A nonzero Love number would be a smoking gun for departures from the classical black hole prediction of General Relativity. We find that an EMRI detection by LISA could set constraints on the tidal Love number of a spinning central object with dimensionless spin $\hat a=0.9$ ($\hat a=0.99$) which are approximately four (six) orders of magnitude more stringent than what achievable with current ground-based detectors for stellar-mass binaries. Our approach is based on the stationary phase approximation to obtain approximate but accurate semi-analytical EMRI waveforms in the frequency-domain, which greatly speeds up high-precision Fisher-information matrix computations. This approach can be easily extended to several other tests of gravity with EMRIs and to efficiently account for multiple deviations in the waveform at the same time., Comment: 8 pages + appendices and references; 2 tables and 1 figure

Astrophysical observations of neutron stars have been widely used to infer the properties of the nuclear matter equation of state. Beside being a source of information on average properties of dense matter, however, the data provided by electromagnetic and gravitational wave (GW) facilities are reaching the accuracy needed to constrain, for the first time, nuclear dynamics in dense matter. In this work we assess the sensitivity of current and future neutron star observations to directly infer the strength of repulsive three-nucleon forces, which are key to determine the stiffness of the equation of state. Using a Bayesian approach we focus on the constraints that can be derived on three-body interactions from binary neutron star mergers observed by second and third-generation of gravitational wave interferometers. We consider both single and multiple observations. For current detectors at design sensitivity the analysis suggests that only low mass systems, with large signal-to-noise ratios (SNR), allow to reliably constrain the three-body forces. However, our results show that a single observation with a third-generation interferometer, such as the Einstein Telescope or Cosmic Explorer, will constrain the strength of the repulsive three-body potential with exquisite accuracy, turning third-generation GW detectors into new laboratories to study the nucleon dynamics., Comment: Minor changes to improve the discussion of the results. Technical details added to extend the equation of state description

An important step in the planning of future gravitational-wave (GW) detectors and of the networks they will form is the estimation of their detection and parameter-estimation capabilities, which is the basis of science-case studies. Several future GW detectors have been proposed or are under development, which might also operate and observe in parallel. These detectors include terrestrial, lunar, and space-borne detectors. In this paper, we present GWFish, a new software to simulate GW detector networks and to calculate measurement uncertainties based on the Fisher-matrix approximation. GWFish models the impact of detector motion on PE and makes it possible to analyze multiband scenarios, i.e., observation of a GW signal by different detectors in different frequency bands. We showcase a few examples for the Einstein Telescope (ET) including the sky-localization of binary neutron stars, and ET's capability to measure the polarization of GWs., Comment: 15 pages, 9 figures

The Laser Interferometer Space Antenna (LISA) has the potential to reveal wonders about the fundamental theory of nature at play in the extreme gravity regime, where the gravitational interaction is both strong and dynamical. In this white paper, the Fundamental Physics Working Group of the LISA Consortium summarizes the current topics in fundamental physics where LISA observations of GWs can be expected to provide key input. We provide the briefest of reviews to then delineate avenues for future research directions and to discuss connections between this working group, other working groups and the consortium work package teams. These connections must be developed for LISA to live up to its science potential in these areas., Comment: Accepted in: Living Reviews in Relativity

The gravitational wave signal from a binary neutron star merger carries the imprint of the deformability properties of the coalescing bodies, and then of the equation of state of neutron stars. In current models of the waveforms emitted in these events, the contribution of tidal deformation is encoded in a set of parameters, the tidal Love numbers. More refined models include tidal-rotation couplings, described by an additional set of parameters, the rotational tidal Love numbers, which appear in the waveform at $6.5$ post-Newtonian order. For neutron stars with spins as large as $\sim0.1$, we show that neglecting tidal-rotation couplings may lead to a significant error in the parameter estimation by third-generation gravitational wave detectors. By performing a Fisher matrix analysis we assess the measurability of rotational tidal Love numbers, showing that their contribution in the waveform could be measured by third-generation detectors. Our results suggest that current models of tidal deformation in late inspiral should be improved in order to avoid waveform systematics and extract reliable information from gravitational wave signals observed by next generation detectors., Comment: 11 pages, 4 figures

The relativistic multipole moments provide a key ingredient to characterize the gravitational field around compact astrophysical objects. They play a crucial role in the description of the orbital evolution of coalescing binary systems and encode valuable information on the nature of the binary's components, which leaves a measurable imprint in their gravitational-wave emission. We present a new study on the multipolar structure of a class of arbitrarily spinning boson stars with quartic self-interactions in the large coupling limit, where these solutions are expected to be stable. Our results strengthen and extend previous numerical analyses, showing that even for the most compact configurations the multipolar structure deviates significantly from that of a Kerr black hole. We provide accurate data for the multipole moments as functions of the object's mass and spin, which can be directly used to construct inspiral waveform approximants and to perform parameter estimations and searches for boson star binaries., Comment: 11+3 pages, 11 figures, data for the quadrupole and octupole moments available at https://web.uniroma1.it/gmunu

We study binary systems in which a stellar mass compact object spirals into a massive black hole, known as extreme mass ratio inspirals, in scenarios with a new fundamental scalar field. Earlier work has shown that, in most interesting such scenarios and to leading order in the mass ratio, the massive black holes can be adequately approximated by the Kerr metric and the imprint of the scalar field on the waveform is fully controlled by the scalar charge of the stellar mass object. Here we use this drastic simplification in the inspiral modelling and consider eccentric equatorial orbits. We study how the scalar charge affects the orbital evolution for different eccentricities and different values of the black hole spin. We then determine how changes in the orbital evolution get imprinted on the waveform and assess LISA's capability to detect or constrain the scalar charge., Comment: Minor changes to match the version published on Phys. Rev. D

Self-gravitating bodies can have an arbitrarily complex shape, which implies a much richer multipolar structure than that of a black hole in General Relativity. With this motivation, we study the corrections to the dynamics of a binary system due to generic, nonaxisymmetric mass quadrupole moments to leading post-Newtonian (PN) order. Utilizing the method of osculating orbits and a multiple scale analysis, we find analytic solutions to the precession and orbital dynamics of a (generically eccentric) binary in terms of the dimensionless modulus parameters $\epsilon_{m}$, corresponding to axial $m=1$ and polar $m=2$ corrections from oblateness/prolateness. The solutions to the precession dynamics are exact for $0 \le \epsilon_{2} < 1$, and perturbative in $\epsilon_{1} \ll 1$. We further compute the leading order corrections to the gravitational wave amplitude and phase for a quasi-circular binary due to mass quadrupole effects. Making use of the stationary phase approximation and shifted uniform asymptotics (SUA), the corrections to the phase enter at relative 2PN order, while the amplitude modulations enter at -0.5PN order with a SUA amplitude correction at 3.25PN order, relative 2PN order to the leading order SUA correction. By investigating the dephasing due to generic quadrupole moments, we find that a phase difference $\gtrsim 0.1$~radians is achievable for $\epsilon_{m} \gtrsim 10^{-3}$, which suggests that constraints with current and future ground-based gravitational wave detectors are possible. Our results can be implemented in parameter estimation studies to constrain generic multipolar deformations of the Kerr geometry and of neutron stars., Comment: 21 pages, 4 figures, journal version

The next generation of ground-based gravitational-wave detectors will observe coalescences of black holes and neutron stars throughout the cosmos, thousands of them with exceptional fidelity. The Science Book is the result of a 3-year effort to study the science capabilities of networks of next generation detectors. Such networks would make it possible to address unsolved problems in numerous areas of physics and astronomy, from Cosmology to Beyond the Standard Model of particle physics, and how they could provide insights into workings of strongly gravitating systems, astrophysics of compact objects and the nature of dense matter. It is inevitable that observatories of such depth and finesse will make new discoveries inaccessible to other windows of observation. In addition to laying out the rich science potential of the next generation of detectors, this report provides specific science targets in five different areas in physics and astronomy and the sensitivity requirements to accomplish those science goals. This report is the second in a six part series of reports by the GWIC 3G Subcommittee: i) Expanding the Reach of Gravitational Wave Observatories to the Edge of the Universe, ii) The Next Generation Global Gravitational Wave Observatory: The Science Book (this report), iii) 3G R&D: R&D for the Next Generation of Ground-based Gravitational Wave Detectors, iv) Gravitational Wave Data Analysis: Computing Challenges in the 3G Era, v) Future Ground-based Gravitational-wave Observatories: Synergies with Other Scientific Communities, and vi) An Exploration of Possible Governance Models for the Future Global Gravitational-Wave Observatory Network., Comment: 69 pages, 18 figures

Primordial black holes possibly formed in the early universe could provide a significant fraction of the dark matter and would be unique probes of inflation. A smoking gun for their discovery would be the detection of a subsolar mass compact object. We argue that extreme mass-ratio inspirals will be ideal to search for subsolar-mass black holes not only with LISA but also with third-generation ground-based detectors such as Cosmic Explorer and the Einstein Telescope. These sources can provide unparalleled measurements of the mass of the secondary object at subpercent level for primordial black holes as light as ${\cal O}(0.01)M_\odot$ up to luminosity distances around hundred megaparsec and few gigaparsec for LISA and Einstein Telescope, respectively, in a complementary frequency range. This would allow claiming, with very high statistical confidence, the detection of a subsolar-mass black hole, which would also provide a novel (and currently undetectable) family of sources for third-generation detectors., Comment: 7 pages, 2 figures. v2: matching version to appear in PRL

We introduce a family of solutions of Einstein's gravity minimally coupled to an anisotropic fluid, describing asymptotically flat black holes with "hair" and a regular horizon. These spacetimes can describe the geometry of galaxies harboring supermassive black holes, and are extensions of Einstein clusters to include horizons. They are useful to constrain the environment surrounding astrophysical black holes, using electromagnetic or gravitational-wave observations. We compute the main properties of the geometry, including the corrections to the ringdown stage induced by the external matter and fluxes by orbiting particles. The leading order effect to these corrections is a gravitational-redshift, but gravitational-wave propagation is affected by the galactic potential in a nontrivial way, and may be characterized with future observatories., Comment: 6 pages, 3 figures. v2: overall improvements in text and discussion. Version accepted for publication in Physical Review D Letters

The science objectives of the LISA mission have been defined under the implicit assumption of a 4 yr continuous data stream. Based on the performance of LISA Pathfinder, it is now expected that LISA will have a duty cycle of $\approx 0.75$, which would reduce the effective span of usable data to 3 yr. This paper reports the results of a study by the LISA Science Group, which was charged with assessing the additional science return of increasing the mission lifetime. We explore various observational scenarios to assess the impact of mission duration on the main science objectives of the mission. We find that the science investigations most affected by mission duration concern the search for seed black holes at cosmic dawn, as well as the study of stellar-origin black holes and of their formation channels via multi-band and multi-messenger observations. We conclude that an extension to 6 yr of mission operations is recommended., Comment: 50 pages, 19 figures, 5 tables. Matches version published in GERG

The Laser Interferometer Space Antenna, LISA, will detect gravitational wave signals from Extreme Mass Ratio Inspirals, where a stellar mass compact object orbits a supermassive black hole and eventually plunges into it. Here we report on LISA's capability to detect whether the smaller compact object in an Extreme Mass Ratio Inspiral is endowed with a scalar field, and to measure its scalar charge -- a dimensionless quantity that acts as a measure of how much scalar field the object carries. By direct comparison of signals, we show that LISA will be able to detect and measure the scalar charge with an accuracy of the order of percent, which is an unprecedented level of precision. This result is independent of the origin of the scalar field and of the structure and other properties of the small compact object, so it can be seen as a generic assessment of LISA's capabilities to detect new fundamental fields., Comment: Typo corrected to match the version in press

Extreme mass-ratio inspirals~(EMRIs) detectable by the Laser Inteferometric Space Antenna~(LISA) are unique probes of astrophysics and fundamental physics. Parameter estimation for these sources is challenging, especially because the waveforms are long, complicated, known only numerically, and slow to compute in the most relevant regime, where the dynamics is relativistic. We perform a time-consuming Fisher-matrix error analysis of the EMRI parameters using fully-relativistic numerical waveforms to leading order in an adiabatic expansion on a Kerr background, taking into account the motion of the LISA constellation, higher harmonics, and also including the leading correction from the spin of the secondary in the post-adiabatic approximation. We pay particular attention to the convergence of the numerical derivatives in the Fisher matrix and to the numerical stability of the covariance matrix, which for some systems requires computing the numerical waveforms with approximately $90$-digit precision. Our analysis confirms previous results (obtained with approximated but much more computationally efficient waveforms) for the measurement errors on the binary's parameters. We also show that the inclusion of higher harmonics improves the errors on the luminosity distance and on the orbital angular momentum angles by one order and two orders of magnitude, respectively, which might be useful to identify the environments where EMRIs live. We particularly focus on the measurability of the spin of the secondary, confirming that it cannot be measured with sufficient accuracy. However, due to correlations, its inclusion in the waveform model can deteriorate the accuracy on the measurements of other parameters by orders of magnitude, unless a physically-motivated prior on the secondary spin is imposed., Comment: 11 pages + appendices and references; 4 tables. Supplemental Mathematica notebooks available at https://web.uniroma1.it/gmunu/resources

Gravitational-wave measurements of the tidal deformability in neutron-star binary coalescences can be used to infer the still unknown equation of state (EoS) of dense matter above the nuclear saturation density. By employing a Bayesian-ranking test we quantify the ability of current and future gravitational-wave observations to discriminate among families of nuclear-physics based EoS which differ in particle content and ab-initio microscopic calculations. While the constraining power of GW170817 is limited, we show that even twenty coalescences detected by LIGO-Virgo at design sensitivity are not enough to discriminate between EoS with similar softness but distinct microphysics. However, just a single detection with a third-generation detector such as the Einstein Telescope or Cosmic Explorer will rule out several families of EoS with very strong statistical significance, and can discriminate among models which feature similar softness, hence constraining the properties of nuclear matter to unprecedented levels., Comment: 9 pages, 4 figures; v3: improved analysis, matches version to appear in PRL

Monitoring of vibrational eigenmodes of an elastic body excited by gravitational waves was one of the first concepts proposed for the detection of gravitational waves. At laboratory scale, these experiments became known as resonant-bar detectors first developed by Joseph Weber in the 1960s. Due to the dimensions of these bars, the targeted signal frequencies were in the kHz range. Weber also pointed out that monitoring of vibrations of Earth or Moon could reveal gravitational waves in the mHz band. His Lunar Surface Gravimeter experiment deployed on the Moon by the Apollo 17 crew had a technical failure rendering the data useless. In this article, we revisit the idea and propose a Lunar Gravitational-Wave Antenna (LGWA). We find that LGWA could become an important partner observatory for joint observations with the space-borne, laser-interferometric detector LISA, and at the same time contribute an independent science case due to LGWA's unique features. Technical challenges need to be overcome for the deployment of the experiment, and development of inertial vibration sensor technology lays out a future path for this exciting detector concept., Comment: 29 pages, 17 figures

We report the results of a study aimed at inferring direct information on the repulsive three-nucleon potential $V^R_{ijk}$\textemdash driving the stiffness of the nuclear matter equation of state at supranuclear densities\textemdash from astrophysical observations. Using a Bayesian approach, we exploit the measurements of masses, radii and tidal deformabalities performed by the NICER satellite and the LIGO/Virgo collaboration, as well as the mass of the heaviest observed pulsar, to constrain the strength of $V^R_{ijk}$. The baseline of our analysis is the widely employed nuclear Hamiltonian comprising the Argonne $v_{18}$ nucleon-nucleon potential andthe Urbana IX model of three-nucleon potential. The numerical results, largely determined by the bound on the maximum mass, suggest that existing and future facilities have the potential to provide valuable new insight into microscopic nuclear dynamics at supranuclear densities., Comment: Discussion of the results extended. Various changes to match the version in press on Phys. Rev. C

Gravitational-wave (GW) detections of binary neutron star coalescences play a crucial role to constrain the microscopic interaction of matter at ultrahigh density. Similarly, if boson stars exist in the universe their coalescence can be used to constrain the fundamental coupling constants of a scalar field theory. We develop the first coherent waveform model for the inspiral of boson stars with quartic interactions. The waveform includes coherently spin-induced quadrupolar and tidal-deformability contributions in terms of the masses and spins of the binary and of a single coupling constant of the theory. We show that future instruments such as the Einstein Telescope and the Laser Interferometer Space Antenna can provide strong complementary bounds on bosonic self-interactions, while the constraining power of current detectors is marginal., Comment: 13 pages, 8 figures; v3: fig.s 2-3 and tab.s 1-2 changed, results unchanged

Gravitational waves from the merger of two neutron stars cannot be easily distinguished from those produced by a comparable-mass mixed binary in which one of the companions is a black hole. Low-mass black holes are interesting because they could form in the aftermath of the coalescence of two neutron stars, from the collapse of massive stars, from matter overdensities in the primordial Universe, or as the outcome of the interaction between neutron stars and dark matter. Gravitational waves carry the imprint of the internal composition of neutron stars via the so-called tidal deformability parameter, which depends on the stellar equation of state and is equal to zero for black holes. We present a new data analysis strategy powered by Bayesian inference and machine learning to identify mixed binaries, hence low-mass black holes, using the distribution of the tidal deformability parameter inferred from gravitational-wave observations., Comment: 13 pages, 6 figures - v2: matches the published version in Phys. Rev. D 102, 023025

We study Extreme Mass Ratio Inspirals (EMRIs), during which a small body spirals into a supermassive black hole, in gravity theories with additional scalar fields. We first argue that no-hair theorems and the properties of known theories that manage to circumvent them introduce a drastic simplification to the problem: the effects of the scalar on supermassive black holes, if any, are mostly negligible for EMRIs in vast classes of theories. We then exploit this simplification to model the inspiral perturbatively and we demonstrate that the scalar charge of the small body leaves a significant imprint on gravitational wave emission. Although much higher precision is needed for waveform modelling, our results strongly suggest that this imprint is observable with LISA, rendering EMRIs promising probes of scalar fields., Comment: Minor changes to match the version in press on Phys. Rev. Lett

Extreme mass-ratio inspirals detectable by the future Laser Interferometer Space Antenna provide a unique way to test general relativity and fundamental physics. Motivated by this possibility, here we study in detail the EMRI dynamics in the presence of a spinning secondary, collecting and extending various results that appeared in previous work and also providing useful intermediate steps and new relations for the first time. We present the results of a frequency-domain code that computes gravitational-wave fluxes and the adiabatic orbital evolution for the case of circular, equatorial orbits with (anti)aligned spins. The spin of the secondary starts affecting the gravitational-wave phase at the same post-adiabatic order as the leading-order self-force terms and introduces a detectable dephasing, which can be used to measure it at $5-25\%$ level, depending on individual spins. In a companion paper we discuss the implication of this effect for tests of the Kerr bound., Comment: Improvements in response to referee reports. Typos corrected. Longer companion paper of arXiv:2003.08448. Data and code available at the website https://web.uniroma1.it/gmunu

An outstanding prediction of general relativity is the fact that the angular momentum $S$ of an isolated black hole with mass $\mu$ is limited by the Kerr bound, $S\leq G\mu^2/c$. Testing this cornerstone is challenging due to the difficulty in modelling spinning compact objects that violate this bound. We argue that precise, model-independent tests can be achieved by measuring gravitational waves from an extreme mass ratio inspiral around a supermassive object, one of the main targets of the future LISA mission. In the extreme mass ratio limit, the dynamics of the small compact object depends only on its multipole moments, which are free parameters. At variance with the comparable-mass case, accurate waveforms are valid also when the spin of the small object greatly exceeds the Kerr bound. By computing the orbital dephasing and the gravitational-wave signal emitted by a spinning point particle in circular, nonprecessing, equatorial motion around a Kerr black hole, we estimate that LISA will be able to measure the spin of the small compact object at the level of $10\%$. Together with mass measurements, this will allow for theory-agnostic, unprecedented constraints on string-theory inspired objects such as "superspinars", almost in their entire parameter space., Comment: Improvements in response to referee reports. Typos corrected. Data and code available at https://web.uniroma1.it/gmunu

Black-hole spectroscopy is arguably the most promising tool to test gravity in extreme regimes and to probe the ultimate nature of black holes with unparalleled precision. These tests are currently limited by the lack of a ringdown parametrization that is both robust and accurate. We develop an observable-based parametrization of the ringdown of spinning black holes beyond general relativity, which we dub ParSpec (Parametrized Ringdown Spin Expansion Coefficients). This approach is perturbative in the spin, but it can be made arbitrarily precise (at least in principle) through a high-order expansion. It requires O(10) ringdown detections, which should be routinely available with the planned space mission LISA and with third-generation ground-based detectors. We provide a preliminary analysis of the projected bounds on parametrized ringdown parameters with LISA and with the Einstein Telescope, and discuss extensions of our model that can be straightforwardly included in the future., Comment: 15 pages, 8 figures. v2: matches published version

Coalescences of binary white dwarfs represent a copious source of information for gravitational wave interferometers operating in the decihertz band. Moreover, according to the double degenerate scenario, they have been suggested to be possible progenitors of supernovae (SNe) Type Ia events. In this paper we discuss the detectability of gravitational waves emitted by the inspiral of double white dwarfs. We focus on the constraints that can be derived on the source's luminosity distance, and on other binary's parameters, such as the angular momentum orientation. We explore the possibility of coincident detections of gravitational and electromagnetic signals; the latter comes from the observation of the supernova counterpart. Confirmation of the double degenerate scenario would allow one to use distances inferred in the gravitational wave channel to consistently calibrate SNe as standard candles. We find that decihertz gravitational wave interferometers can measure the luminosity distance with relative accuracy better than $1\%$ for binaries at 100 Mpc. We show how multimessenger observations can put strong constraints on the Hubble constant, which are tighter than current bounds at low redshift, and how they can potentially shed new light on the differences with early-universe measurements., Comment: Results updated to match the version in press on A&A

The dynamics of coalescing compact binaries can be affected by the environment in which the systems evolve, leaving detectable signatures into the emitted gravitational signal. In this paper we investigate the ability of gravitational-wave detectors to constrain the nature of the environment in which compact binaries merge. We parametrize a variety of environmental effects by modifying the phase of the gravitational signal emitted by black hole and neutron star binaries. We infer the bounds on such effects by current and future generation of interferometers, studying their dependence on the binary's parameters. We show that the strong dephasing induced by accretion and dynamical friction can constraint the density of the surrounding medium to orders of magnitude below that of accretion disks. Planned detectors, such as LISA or DECIGO, will be able to probe densities typical of those of dark matter., Comment: Corrected typos in the figures compared to the version published on A&A

Linear perturbations of spherically symmetric spacetimes in general relativity are described by radial wave equations, with potentials that depend on the spin of the perturbing field. In previous work we studied the quasinormal mode spectrum of spacetimes for which the radial potentials are slightly modified from their general relativistic form, writing generic small modifications as a power-series expansion in the radial coordinate. We assumed that the perturbations in the quasinormal frequencies are linear in some perturbative parameter, and that there is no coupling between the perturbation equations. In general, matter fields and modifications to the gravitational field equations lead to coupled wave equations. Here we extend our previous analysis in two important ways: we study second-order corrections in the perturbative parameter, and we address the more complex (and realistic) case of coupled wave equations. We highlight the special nature of coupling-induced corrections when two of the wave equations have degenerate spectra, and we provide a ready-to-use recipe to compute quasinormal modes. We illustrate the power of our parametrization by applying it to various examples, including dynamical Chern-Simons gravity, Horndeski gravity and an effective field theory-inspired model., Comment: 14 pages, 5 figures, 3 tables. Tabulated coefficients for generic QNM calculations are available online (see the conclusions for a computational recipe). v2: fixed minor typos to match version in press in PRD. v3: added a few references, matches version published in PRD

We study the constraints on alternative theories of gravity that can be determined by multi-band observations of gravitational wave signals emitted from binary black hole coalescences. We focus on three types of General Relativity modifications induced by a generalised Brans-Dicke theory, and two classes of quadratic gravity, Einstein-dilaton-Gauss-Bonnet and dynamical Chern-Simons. Considering a network of space and ground-based detectors, supplied by a population of spinning binaries black hole, we show how the multi-band analysis improves the existing bounds on the theory's parameters by several orders of magnitude, for both pre- and post-Newtonian deviations. Our results also show the fundamental role played by an interferometer in the frequency range between LISA and advanced detectors, in constraining possible deviations from General Relativity., Comment: Results updated to match the version in press on Phys. Rev. D

Asymmetric Dark Stars, i.e., compact objects formed from the collapse of asymmetric dark matter could potentially produce a detectable photon flux if dark matter particles self-interact via dark photons that kinetically mix with ordinary photons. The morphology of the emitted spectrum is significantly different and therefore distinguishable from a typical black-body one. Given the above and the fact that asymmetric dark stars can have masses outside the range of neutron stars, the detection of such a spectrum can be considered as a smoking gun signature for the existence of these exotic stars., Comment: Minor changes to match the version published on IJMPD