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 LISA can provide strong, complementary bounds on bosonic self-interactions, while the constraining power of current detectors is marginal.

Within General Relativity, the unique stationary solution of an isolated black hole is the Kerr spacetime, which has a peculiar multipolar structure depending only on its mass and spin. We develop a general method to extract the multipole moments of arbitrary stationary spacetimes and apply it to a large family of horizonless microstate geometries. The latter can break the axial and equatorial symmetry of the Kerr metric and have a much richer multipolar structure, which provides a portal to constrain fuzzball models phenomenologically. We find numerical evidence that all multipole moments are always larger (in absolute value) than those of a Kerr black hole with the same mass and spin. Current measurements of the quadrupole moment of black-hole candidates could place only mild constraints on fuzzballs, while future gravitational-wave detections of extreme mass-ratio inspirals with the space mission LISA will improve these bounds by orders of magnitude.

We develop a new method to measure neutron star parameters and derive constraints on the equation of state of dense matter by fitting the frequencies of simultaneous Quasi Periodic Oscillation modes observed in the X-ray flux of accreting neutron stars in low mass X-ray binaries. To this aim we calculate the fundamental frequencies of geodesic motion around rotating neutron stars based on an accurate general-relativistic approximation for their external spacetime. Once the fundamental frequencies are related to the observed frequencies through a QPO model, they can be fit to the data to obtain estimates of the three parameters describing the spacetime, namely the neutron star mass, angular momentum and quadrupole moment. From these parameters we derive information on the neutron star structure and equation of state. We present a proof of principle of our method applied to pairs of kHz QPO frequencies observed from three systems (4U1608-52, 4U0614+09 and 4U1728-34). We identify the kHz QPOs with the azimuthal and the periastron precession frequencies of matter orbiting the neutron star, and via our Bayesian inference technique we derive constraints on the neutrons stars' masses and radii. This method is applicable to other geodesic-frequency-based QPO models.

A generic feature of nearly out-of-equilibrium dissipative systems is that they resonate through a set of quasinormal modes. Black holes - the absorbing objects par excellence - are no exception. When formed in a merger, black holes vibrate in a process called "ringdown", which leaves the gravitational-wave footprint of the event horizon. In some models of quantum gravity which attempt to solve the information-loss paradox and the singularities of General Relativity, black holes are replaced by regular, horizonless objects with a tiny effective reflectivity. Motivated by these scenarios, here we develop a generic framework to the study of the ringdown of a compact object with various shades of darkness. By extending the black-hole membrane paradigm, we map the interior of any compact object in terms of the bulk and shear viscosities of a fictitious fluid located at the surface, with the black-hole limit being a single point in a three-dimensional parameter space. We unveil some remarkable features of the ringdown and some universal properties of the light ring in this framework. We also identify the region of the parameter space which can be probed by current and future gravitational-wave detectors. A general feature is the appearance of mode doublets which are degenerate only in the black-hole limit. We argue that the merger event GW150914 already imposes a strong lower bound on the compactness of the merger remnant of approximately 99% of the black-hole compactness. This places model-independent constraints on black-hole alternatives such as diffuse "fuzzballs" and nonlocal stars.

We study the impact of gas accretion on the orbital evolution of black-hole binaries initially at large separation in the band of the planned Laser Interferometer Space Antenna (LISA). We focus on two sources: (i)~stellar-origin black-hole binaries~(SOBHBs) that can migrate from the LISA band to the band of ground-based gravitational-wave observatories within weeks/months; and (ii) intermediate-mass black-hole binaries~(IMBHBs) in the LISA band only. Because of the large number of observable gravitational-wave cycles, the phase evolution of these systems needs to be modeled to great accuracy to avoid biasing the estimation of the source parameters. Accretion affects the gravitational-wave phase at negative ($-4$) post-Newtonian order, and is therefore dominant for binaries at large separations. If accretion takes place at the Eddington or at super-Eddington rate, it will leave a detectable imprint on the dynamics of SOBHBs. In optimistic astrophysical scenarios, a multiwavelength strategy with LISA and a ground-based interferometer can detect about $10$ (a few) SOBHB events for which the accretion rate can be measured at $50\%$ ($10\%$) level. In all cases the sky position can be identified within much less than $0.4\,{\rm deg}^2$ uncertainty. Likewise, accretion at $\gtrsim 10\%$ ($\gtrsim 100\%$) of the Eddington rate can be measured in IMBHBs up to redshift $z\approx 0.1$ ($z\approx 0.5$), and the position of these sources can be identified within less than $0.01\,{\rm deg}^2$ uncertainty. Altogether, a detection of SOBHBs or IMBHBs would allow for targeted searches of electromagnetic counterparts to black-hole mergers in gas-rich environments with future X-ray detectors (such as Athena) and radio observatories (such as SKA).

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.

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.

Ultralight bosonic fields are compelling dark-matter candidates and arise in a variety of beyond-Standard-Model scenarios. These fields can tap energy and angular momentum from spinning black holes through superradiant instabilities, during which a macroscopic bosonic condensate develops around the black hole. Striking features of this phenomenon include gaps in the spin-mass distribution of astrophysical black holes and a continuous gravitational-wave~(GW) signal emitted by the condensate. So far these processes have been studied in great detail for scalar fields and, more recently, for vector fields. Here we take an important step forward in the black-hole superradiance program by computing, analytically, the instability time scale, the direct GW emission, and the stochastic background, in the case of massive tensor (i.e., spin-$2$) fields. Our analysis is valid for any black hole spin and for small boson masses. The instability of massive spin-$2$ fields shares some properties with the scalar and vector cases, but its phenomenology is much richer, for example there exist multiple modes with comparable instability time scales, and the dominant GW signal is hexadecapolar rather than quadrupolar. Electromagnetic and GW observations of spinning black holes in the mass range $M\in(1,10^{10})M_\odot$ can constrain the mass of a putative spin-$2$ field in the range $10^{-22} \lesssim m_b\,{\rm c^2/eV} \lesssim 10^{-10}$. For $10^{-17}\lesssim m_b\,{\rm c^2/eV}\lesssim 10^{-15}$, the space mission LISA could detect the continuous GW signal for sources at redshift $z=20$, or even larger.

We study novel solitonic solutions to Einstein-Klein-Gordon theory in the presence of a periodic scalar potential arising in models of axion-like particles. The potential depends on two parameters: the mass of the scalar field $m_a$ and the decay constant $f_a$; the standard case of the QCD axion is recovered when $m_a\propto1/f_a$. When $f_a\to\infty$ the solutions reduce to the standard case of "mini" boson stars supported by a massive free scalar field. As the energy scale $f_a$ of the scalar self-interactions decreases we unveil several novel features of the solution: new stability branches emerge at high density, giving rise to very compact, radially stable, boson stars. Some of the most compact configurations acquire a photon sphere. When $f_a$ is at the GUT scale, a boson star made of QCD axions can have a mass up to ten solar masses and would be more compact than a neutron star. Gravitational-wave searches for these exotic compact objects might provide indirect evidence for ultralight axion-like particles in a region not excluded by the black-hole superradiant instability.

The LIGO and Virgo Interferometers have so far provided 11 gravitational-wave (GW) observations of black-hole binaries. Similar detections are bound to become very frequent in the near future. With the current and upcoming wealth of data, it is possible to confront specific formation models with observations. We investigate here whether current data are compatible with the hypothesis that LIGO/Virgo black holes are of primordial origin. We compute in detail the mass and spin distributions of primordial black holes (PBHs), their merger rates, the stochastic background of unresolved coalescences, and confront them with current data from the first two observational runs, also including the recently discovered GW190412. We compute the best-fit values for the parameters of the PBH mass distribution at formation that are compatible with current GW data. In all cases, the maximum fraction of PBHs in dark matter is constrained by these observations to be $f_{\text{PBH}}\approx {\rm few}\times 10^{-3}$. We discuss the predictions of the PBH scenario that can be directly tested as new data become available. In the most likely formation scenarios where PBHs are born with negligible spin, the fact that at least one of the components of GW190412 is moderately spinning is incompatible with a primordial origin for this event, unless accretion or hierarchical mergers are significant. In the absence of accretion, current non-GW constraints already exclude that LIGO/Virgo events are all of primordial origin, whereas in the presence of accretion the GW bounds on the PBH abundance are the most stringent ones in the relevant mass range. A strong phase of accretion during the cosmic history would favour mass ratios close to unity, and a redshift-dependent correlation between high masses, high spins and nearly-equal mass binaries, with the secondary component spinning faster than the primary.

The black hole uniqueness and the no-hair theorems imply that the quasinormal spectrum of any astrophysical black hole is determined solely by its mass and spin. The countably infinite number of quasinormal modes of a Kerr black hole are thus related to each other and any deviations from these relations provide a strong hint for physics beyond the general theory of relativity. To test the no-hair theorem using ringdown signals, it is necessary to detect at least two quasinormal modes. In particular, one can detect the fundamental mode along with a subdominant overtone or with another angular mode, depending on the mass ratio and the spins of the progenitor binary. Also in the light of the recent discovery of GW190412, studying how the mass ratio affects the prospect of black hole spectroscopy using overtones or angular modes is pertinent, and this is the major focus of our study. First, we provide ready-to-use fits for the amplitudes and phases of both the angular modes and overtones as a function of mass ratio $q\in[0,10]$. Using these fits we estimate the minimum signal-to-noise ratio for detectability, resolvability, and measurability of subdominant modes/tones. We find that performing black-hole spectroscopy with angular modes is preferable when the binary mass ratio is larger than $q\approx 1.2$ (provided that the source is not located at a particularly disfavoured inclination angle). For nonspinning, equal-mass binary black holes, the overtones seem to be the only viable option to perform a spectroscopy test of the no-hair theorem. However this would require a large ringdown signal-to-noise ratio ($\approx 100$ for a $5\%$ accuracy test with two overtones) and the inclusion of more than one overtone to reduce modelling errors, making black-hole spectroscopy with overtones impractical in the near future.

We consider the constraints on the fraction of dark matter in the universe in the form of primordial black holes taking into account the crucial role of accretion which may change both their mass and mass function. We show that accretion may drastically weaken the constraints at the present epoch for primordial black holes with masses larger than a few solar masses.

Primordial black holes in the mass range of ground-based gravitational-wave detectors can comprise a significant fraction of the dark matter. Mass and spin measurements from coalescences can be used to distinguish between an astrophysical or a primordial origin of the binary black holes. In standard scenarios the spin of primordial black holes is very small at formation. However, the mass and spin can evolve through the cosmic history due to accretion. We show that the mass and spin of primordial black holes are correlated in a redshift-dependent fashion, in particular primordial black holes with masses below ${\cal O}(30)M_\odot$ are likely non-spinning at any redshift, whereas heavier black holes can be nearly extremal up to redshift $z\sim10$. The dependence of the mass and spin distributions on the redshift can be probed with future detectors such as the Einstein Telescope. The mass and spin evolution affect the gravitational waveform parameters, in particular the distribution of the final mass and spin of the merger remnant, and that of the effective spin of the binary. We argue that, compared to the astrophysical-formation scenario, a primordial origin of black hole binaries might better explain the spin distribution of merger events detected by LIGO-Virgo, in which the effective spin parameter of the binary is compatible to zero except possibly for few high-mass events. Upcoming results from LIGO-Virgo third observation run might reinforce or weaken these predictions.

In viable models of minicharged dark matter, astrophysical black holes might be charged under a hidden $U(1)$ symmetry and are formally described by the same Kerr-Newman solution of Einstein-Maxwell theory. These objects are unique probes of minicharged dark matter and dark photons. We show that the recent gravitational-wave detection of a binary black-hole coalescence by aLIGO provides various observational bounds on the black hole's charge, regardless of its nature. The pre-merger inspiral phase can be used to constrain the dipolar emission of (ordinary and dark) photons, whereas the detection of the quasinormal modes set an upper limit on the final black hole's charge. By using a toy model of a point charge plunging into a Reissner-Nordstrom black hole, we also show that in dynamical processes the (hidden) electromagnetic quasinormal modes of the final object are excited to considerable amplitude in the gravitational-wave spectrum only when the black hole is nearly extremal. The coalescence produces a burst of low-frequency dark photons which might provide a possible electromagnetic counterpart to black-hole mergers in these scenarios.

Superradiance is a radiation enhancement process that involves dissipative systems. With a 60 year-old history, superradiance has played a prominent role in optics, quantum mechanics and especially in relativity and astrophysics. In General Relativity, black-hole superradiance is permitted by the ergoregion, that allows for energy, charge and angular momentum extraction from the vacuum, even at the classical level. Stability of the spacetime is enforced by the event horizon, where negative energy-states are dumped. Black-hole superradiance is intimately connected to the black-hole area theorem, Penrose process, tidal forces, and even Hawking radiation, which can be interpreted as a quantum version of black-hole superradiance. Various mechanisms (as diverse as massive fields, magnetic fields, anti-de Sitter boundaries, nonlinear interactions, etc...) can confine the amplified radiation and give rise to strong instabilities. These "black-hole bombs" have applications in searches of dark matter and of physics beyond the Standard Model, are associated to the threshold of formation of new black hole solutions that evade the no-hair theorems, can be studied in the laboratory by devising analog models of gravity, and might even provide a holographic description of spontaneous symmetry breaking and superfluidity through the gauge-gravity duality. This work is meant to provide a unified picture of this multifaceted subject. We focus on the recent developments in the field, and work out a number of novel examples and applications, ranging from fundamental physics to astrophysics.

Einstein's theory of general relativity predicts that the only stationary configuration of an isolated black hole is the Kerr spacetime, which has a unique multipolar structure and a spherical shape when non-spinning. This is in striking contrast to the case of other self-gravitating objects, which instead can in principle have arbitrary deformations even in the static case. Here we develop a general perturbative framework to construct stationary stars with small axisymmetric deformations, and study explicitly compact stars with an intrinsic quadrupole moment. The latter can be sustained, for instance, by crust stresses or strong magnetic fields. While our framework is general, we focus on quadrupolar deformations of neutron stars induced by an anisotropic crust, which continuously connect to spherical neutron stars in the isotropic limit. Deformed neutron stars might provide a more accurate description for stellar remnants formed in supernovae and in binary mergers, and can be used to improve constraints on the neutron-star equation of state through gravitational-wave detections and through the observation of low-mass X-ray binaries. We argue that, if the (dimensionless) intrinsic quadrupole moment is of a few percent or higher, the effect of the deformation is stronger than that of tidal interactions in coalescing neutron-star binaries, and might also significantly affect the electromagnetic signal from accreting neutron stars.

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.

Validating the black-hole no-hair theorem with gravitational-wave observations of compact binary coalescences provides a compelling argument that the remnant object is indeed a black hole as described by the general theory of relativity. This requires performing a spectroscopic analysis of the post-merger signal and resolving the frequencies of either different angular modes or overtones (of the same angular mode). For a nearly-equal mass binary black-hole system, only the dominant angular mode ($l=m=2$) is sufficiently excited and the overtones are instrumental to perform this test. Here we investigate the robustness of modelling the post-merger signal of a binary black hole coalescence as a superposition of overtones. Further, we study the bias expected in the recovered frequencies as a function of the start time of a spectroscopic analysis and provide a computationally cheap procedure to choose it based on the interplay between the expected statistical error due to the detector noise and the systematic errors due to waveform modelling. Moreover, since the overtone frequencies are closely spaced, we find that resolving the overtones is particularly challenging and requires a loud ringdown signal. Rayleigh's resolvability criterion suggests that in an optimistic scenario a ringdown signal-to-noise ratio larger than $\sim 30$ (achievable possibly with LIGO at design sensitivity and routinely with future interferometers such as Einstein Telescope, Cosmic Explorer, and LISA) is necessary to resolve the overtone frequencies. We then conclude by discussing some conceptual issues associated with black-hole spectroscopy with overtones.

The defining feature of a classical black hole is being a perfect absorber. Any evidence showing otherwise would indicate a departure from the standard black-hole picture. Energy and angular momentum absorption by the horizon of a black hole is responsible for tidal heating in a binary. This effect is particularly important in the latest stages of an extreme mass ratio inspiral around a spinning supermassive object, one of the main targets of the future LISA mission. We study how this effect can be used to probe the nature of supermassive objects in a model independent way. We compute the orbital dephasing and the gravitational-wave signal emitted by a point particle in circular, equatorial motion around a spinning supermassive object to the leading order in the mass ratio. Absence of absorption by the central object can affect the gravitational-wave signal dramatically, especially at high spin. This effect will make it possible to put an unparalleled upper bound on the reflectivity of exotic compact objects, at the level of ${\cal O}(0.01)\%$. This stringent bound would exclude the possibility of observing echoes in the ringdown of a supermassive binary merger.

Gravitational-wave echoes in the post-merger signal of a binary coalescence are predicted in various scenarios, including near-horizon quantum structures, exotic states of matter in ultracompact stars, and certain deviations from general relativity. The amplitude and frequency of each echo is modulated by the photon-sphere barrier of the remnant, which acts as a spin- and frequency-dependent high-pass filter, decreasing the frequency content of each subsequent echo. Furthermore, a major fraction of the energy of the echo signal is contained in low-frequency resonances corresponding to the quasi-normal modes of the remnant. Motivated by these features, in this work we provide an analytical gravitational-wave template in the low-frequency approximation describing the post-merger ringdown and the echo signal of a spinning ultracompact object. Besides the standard ringdown parameters, the template is parametrized in terms of only two physical quantities: the reflectivity coefficient and the compactness of the remnant. We discuss novel effects related to the spin and to the complex reflectivity, such as a more involved modulation of subsequent echoes, the mixing of two polarizations, and the ergoregion instability in case of perfectly-reflecting spinning remnants. Finally, we compute the errors in the estimation of the template parameters with current and future gravitational-wave detectors using a Fisher matrix framework. Our analysis suggests that models with almost perfect reflectivity can be excluded/detected with current instruments, whereas probing values of the reflectivity smaller than $80\%$ at $3\sigma$ confidence level requires future detectors (Einstein Telescope, Cosmic Explorer, LISA). The template developed in this work can be easily implemented to perform a matched-filter based search for echoes and to constrain models of exotic compact objects.

Ultralight bosonic fields in the mass range $\sim (10^{-20}-10^{-11})\,{\rm eV}$ can trigger a superradiant instability that extracts energy and angular momentum from an astrophysical black hole with mass $M\sim(5,10^{10})M_\odot$, forming a nonspherical, rotating condensate around it. So far, most studies of the evolution and end-state of the instability have been limited to initial data containing only the fastest growing superradiant mode. By studying the evolution of multimode data in a quasi-adiabatic approximation, we show that the dynamics is much richer and depend strongly on the energy of the seed, on the relative amplitude between modes, and on the gravitational coupling. If the seed energy is a few percent of the black-hole mass, a black hole surrounded by a mixture of superradiant and nonsuperradiant modes with comparable amplitudes might not undergo a superradiant unstable phase, depending on the value of the boson mass. If the seed energy is smaller, as in the case of an instability triggered by quantum fluctuations, the effect of nonsuperradiant modes is negligible. We discuss the implications of these findings for current constraints on ultralight fields with electromagnetic and gravitational-wave observations.

Future gravitational-wave observations will enable unprecedented and unique science in extreme gravity and fundamental physics answering questions about the nature of dynamical spacetimes, the nature of dark matter and the nature of compact objects.

Black holes are unique among astrophysical sources: they are the simplest macroscopic objects in the Universe, and they are extraordinary in terms of their ability to convert energy into electromagnetic and gravitational radiation. Our capacity to probe their nature is limited by the sensitivity of our detectors. The LIGO/Virgo interferometers are the gravitational-wave equivalent of Galileo's telescope. The first few detections represent the beginning of a long journey of exploration. At the current pace of technological progress, it is reasonable to expect that the gravitational-wave detectors available in the 2035-2050s will be formidable tools to explore these fascinating objects in the cosmos, and space-based detectors with peak sensitivities in the mHz band represent one class of such tools. These detectors have a staggering discovery potential, and they will address fundamental open questions in physics and astronomy. Are astrophysical black holes adequately described by general relativity? Do we have empirical evidence for event horizons? Can black holes provide a glimpse into quantum gravity, or reveal a classical breakdown of Einstein's gravity? How and when did black holes form, and how do they grow? Are there new long-range interactions or fields in our universe, potentially related to dark matter and dark energy or a more fundamental description of gravitation? Precision tests of black hole spacetimes with mHz-band gravitational-wave detectors will probe general relativity and fundamental physics in previously inaccessible regimes, and allow us to address some of these fundamental issues in our current understanding of nature.

We derive novel constraints on millicharged dark matter and ultralight axion-like particles using pulsar timing and fast radio burst observations. Millicharged dark matter affects the dispersion measure of the time of arrival of radio pulses in a way analogous to free electrons. Light pseudo-scalar dark matter, on the other hand, causes the polarization angle of radio signals to oscillate. We show that current and future data can set strong constraints in both cases. For dark matter particles of charge $\epsilon e$, these constraints are ${\epsilon}/{m_{\rm milli}} \lesssim 10^{-8}{\rm eV}^{-1}$, for masses $m_{\rm milli}\gtrsim 10^{-6}\,$eV. For axion-like particles, the analysis of signals from pulsars yields constraints in the axial coupling of the order of $g/m_a\lesssim 10^{-13} {\rm GeV}^{-1}/(10^{-22}{\rm eV})$. Both bounds scale as $(\rho/\rho_{\rm dm})^{1/2}$ if the energy density $\rho$ of the components is a fraction of the total dark matter energy density $\rho_{\rm dm}$. We do a detailed study of both effects using data from two samples of pulsars in the galaxy and in globular clusters, as well as data from FRB 121102 and PSR J0437$-$4715. We show that in both cases actual pulsar data constrain a new region of the parameter space for these models, and will improve with future pulsar-timing observations.

Supermassive binaries detectable by the future space gravitational-wave interferometer LISA might allow to distinguish black holes from ultracompact horizonless objects, even when the latter are motivated by quantum-gravity considerations. We show that a measurement of very small tidal Love numbers at the level of $10\%$ accuracy (as achievable with "golden binaries") may also allow to distinguish between different models of these exotic compact objects, even when taking into account an intrinsic uncertainty in the object radius putatively due to quantum mechanics. We argue that there is no conceptual obstacle in performing these measurements, the main challenge remains the detectability of small tidal effects and an accurate waveform modelling. Our analysis uses only coordinate-independent quantities related to the proper radial distance and the total mass of the object.

Very compact objects probe extreme gravitational fields and may be the key to understand outstanding puzzles in fundamental physics. These include the nature of dark matter, the fate of spacetime singularities, or the loss of unitarity in Hawking evaporation. The standard astrophysical description of collapsing objects tells us that massive, dark and compact objects are black holes. Any observation suggesting otherwise would be an indication of beyond-the-standard-model physics. Null results strengthen and quantify the Kerr black hole paradigm. The advent of gravitational-wave astronomy and precise measurements with very long baseline interferometry allow one to finally probe into such foundational issues. We overview the physics of exotic dark compact objects and their observational status, including the observational evidence for black holes with current and future experiments.

Anisotropic stresses are ubiquitous in nature, but their modeling in General Relativity is poorly understood and frame dependent. We introduce the first study on the dynamical properties of anisotropic self-gravitating fluids in a covariant framework. Our description is particularly useful in the context of tests of the black hole paradigm, wherein ultracompact objects are used as black hole mimickers but otherwise lack a proper theoretical framework. We show that: (i) anisotropic stars can be as compact and as massive as black holes, even for very small anisotropy parameters; (ii) the nonlinear dynamics of the 1+1 system is in good agreement with linearized calculations, and shows that configurations below the maximum mass are nonlinearly stable; (iii) strongly anisotropic stars have vanishing tidal Love numbers in the black-hole limit; (iv) their formation will usually be accompanied by gravitational-wave echoes at late times.

The tidal deformability of a self-gravitating object leaves an imprint on the gravitational-wave signal of an inspiral which is paramount to measure the internal structure of the binary components. We unveil here a surprisingly unnoticed effect: in the extreme-mass ratio limit the tidal Love number of the central object (i.e. the quadrupole moment induced by the tidal field of its companion) affects the gravitational waveform at the leading order in the mass ratio. This effect acts as a magnifying glass for the tidal deformability of supermassive objects but was so far neglected, probably because the tidal Love numbers of a black hole (the most natural candidate for a compact supermassive object) are identically zero. We argue that extreme-mass ratio inspirals detectable by the future LISA mission might place constraints on the tidal Love numbers of the central object which are roughly 8 orders of magnitude more stringent than current ones on neutron stars, potentially probing all models of black hole mimickers proposed so far.

The existence of black holes and of spacetime singularities is a fundamental issue in science. Despite this, observations supporting their existence are scarce, and their interpretation unclear. We overview how strong a case for black holes has been made in the last few decades, and how well observations adjust to this paradigm. Unsurprisingly, we conclude that observational evidence for black holes is impossible to come by. However, just like Popper's black swan, alternatives can be ruled out or confirmed to exist with a single observation. These observations are within reach. In the next few years and decades, we will enter the era of precision gravitational-wave physics with more sensitive detectors. Just as accelerators require larger and larger energies to probe smaller and smaller scales, more sensitive gravitational-wave detectors will be probing regions closer and closer to the horizon, potentially reaching Planck scales and beyond. What may be there, lurking?

Motivated by the lack of a general parametrization for exotic compact objects, we construct a class of perturbative solutions valid for small (but otherwise generic) multipolar deviations from a Schwarzschild metric in general relativity. We introduce two classes of exotic compact objects, with "soft" and "hard" hair, for which the curvature at the surface is respectively comparable to or much larger than that at the corresponding black-hole horizon. We extend the Hartle-Thorne formalism to relax the assumption of equatorial symmetry and to include deformations induced by multipole moments higher than the spin, thus constructing the most general, axisymmetric quasi-Schwarzschild solution to Einstein's vacuum equations. We explicitly construct several particular solutions of objects with soft hair, which might be useful for tests of quasi-black-hole metrics, and also to study deformed neutron stars. We show that the more compact a soft exotic object is, the less hairy it will be. All its multipole moments can approach their corresponding Kerr values only in two ways as their compactness increases: either logarithmically (or faster) if the moments are spin-induced, or linearly (or faster) otherwise. Our results suggest that it is challenging (but possibly feasible with next-generation gravitational-wave detectors) to distinguish Kerr black holes from a large class of ultracompact exotic objects on the basis of their different multipolar structure.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

Spinning horizonless compact objects may be unstable against an 'ergoregion instability'. We investigate this mechanism for electromagnetic perturbations of ultracompact Kerr-like objects with a reflecting surface, extending previous (numerical and analytical) work limited to the scalar case. We derive an analytical result for the frequency and the instability time scale of unstable modes which is valid at small frequencies. We argue that our analysis can be directly extended to gravitational perturbations of exotic compact objects in the black-hole limit. The instability for electromagnetic and gravitational perturbations is generically stronger than in the scalar case and it requires larger absorption to be quenched. We argue that exotic compact objects with spin $\chi\lesssim 0.7$ ($\chi\lesssim 0.9$) should have an absorption coefficient of at least $0.3\%$ ($6\%$) to remain linearly stable, and that an absorption coefficient of at least $\approx60\%$ would quench the instability for any spin. We also show that - in the static limit - the scalar, electromagnetic, and gravitatonal perturbations of the Kerr metric are related to one another through Darboux transformations.

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.

In this brief note, we clarify certain aspects related to the magnetic (i.e., odd parity or axial) tidal Love numbers of a star in general relativity. Magnetic tidal deformations of a compact star had been computed in 2009 independently by Damour and Nagar and by Binnington and Poisson. More recently, Landry and Poisson showed that the magnetic tidal Love numbers depend on the assumptions made on the fluid, in particular they are different (and of opposite sign) if the fluid is assumed to be in static equilibrium or if it is irrotational. We show that the zero-frequency limit of the Regge-Wheeler equation forces the fluid to be irrotational. For this reason, the results of Damour and Nagar are equivalent to those of Landry and Poisson for an irrotational fluid, and are expected to be the most appropriate to describe realistic configurations.

In this paper we describe the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to accretion flows in the strong field gravity regime around both stellar-mass and supermassive black-holes. eXTP has the unique capability of using advanced 'spectral-timing-polarimetry' techniques to analyze the rapid variations with three orthogonal diagnostics of the flow and its geometry, yielding unprecedented insight into the inner accreting regions, the effects of strong field gravity on the material within them and the powerful outflows which are driven by the accretion process.

We compute the spin-tidal couplings that affect the dynamics of two orbiting bodies at the leading order in the post-Newtonian (PN) framework and to linear order in the spin. These corrections belong to two classes: (i) terms arising from the coupling between the ordinary tidal terms and the point-particle terms, which depend on the standard tidal Love numbers of order $l$ and affect the gravitational-wave (GW) phase at $(2l+5/2)$PN order and (ii) terms depending on the rotational tidal Love numbers, recently introduced in previous work, that affect the GW phase at $(2l+1/2+\delta_{2l})$PN order. For circular orbits and spins orthogonal to the orbital plane, all leading-order spin-tidal terms enter the GW phase at $1.5$PN order relative to the standard, quadrupolar, tidal deformability term (and, thus, before the standard octupolar tidal deformability terms). We present the GW phase that includes all tidal terms up to $6.5$PN order and to linear order in the spin. We comment on a conceptual issue related to the inclusion of the rotational tidal Love numbers in a Lagrangian formulation and on the relevance of spin-tidal couplings for parameter estimation in coalescing neutron-star binaries and for tests of gravity.

We study the nonlinear dynamics of black holes that carry scalar hair and binaries composed of such black holes. The scalar hair is due to a linear or exponential coupling between the scalar and the Gauss--Bonnet invariant. We work perturbatively in the coupling constant of that interaction but nonperturbatively in the fields. We first consider the dynamical formation of hair for isolated black holes of arbitrary spin and determine the final state. This also allows us to compute for the first time the scalar quasinormal modes of rotating black holes in the presence of this coupling. We then study the evolution of nonspinning black-hole binaries with various mass ratios and produce the first scalar waveform for a coalescence. An estimate of the energy loss in scalar radiation and the effect this has on orbital dynamics and the phase of the GWs (entering at quadratic order in the coupling) show that GW detections can set the most stringent constraint to date on theories that exhibit a coupling between a scalar field and the Gauss--Bonnet invariant.

Gravitational-wave astronomy has the potential to explore one of the deepest and most puzzling aspects of Einstein's theory: the existence of black holes. A plethora of ultracompact, horizonless objects have been proposed to arise in models inspired by quantum gravity. These objects may solve Hawking's information-loss paradox and the singularity problem associated with black holes, while mimicking almost all of their classical properties. They are, however, generically unstable on relatively short timescales. Here, we show that this "ergoregion instability" leads to a strong stochastic background of gravitational waves, at a level detectable by current and future gravitational-wave detectors. The absence of such background in the first observation run of Advanced LIGO already imposes the most stringent limits to date on black-hole alternatives, showing that certain models of "quantum-dressed" stellar black holes can be at most a small percentage of the total population. The future LISA mission will allow for similar constraints on supermassive black-hole mimickers.

Gravitational-wave echoes in the post-merger ringdown phase are under intense scrutiny as probes of near-horizon quantum structures and as signatures of exotic states of matter in ultracompact stars. We present an analytical template that describes the ringdown and the echo signal for nonspinning objects in terms of two physical parameters: the reflectivity and the redshift at the surface of the object. We characterize the properties of the template and adopt it in a preliminary parameter estimation with current (aLIGO) and future (Cosmic Explorer, Einstein Telescope, LISA) gravitational-wave detectors. For fixed signal-to-noise ratio in the post-merger phase, the constraints on the model parameters depend only mildly on the details of the detector sensitivity curve, but depend strongly on the reflectivity. Our analysis suggests that it might be possible to detect or rule out Planckian corrections at the horizon scale for perfectly-reflecting ultracompact objects at $5\sigma$ confidence level with Advanced LIGO/Virgo. On the other hand, signal-to-noise ratios in the ringdown phase equal to $\approx 100$ (as achievable with future interferometers) might allow us to probe near-horizon quantum structures with reflectivity $\gtrsim30\%$ ($\gtrsim85\%$) at $2\sigma$ ($3\sigma$) level.

A tentative detection of gravitational-wave echoes in the post-merger signal of GW170817 has been recently claimed at $4.2\sigma$ significance level. It has been speculated that the signal might provide evidence for near-horizon quantum structures in the remnant exotic object. We point out that if the remnant object is an ultracompact neutron star, echoes are expected for objects with radius only slightly smaller than that of an ordinary neutron star. The reported echoes at ~72 Hz are compatible with a toy model of incompressible star with mass approximately $M\in(2,3) M_\odot$ and radius close to the Buchdahl limit, R~9GM/(4c^2). If confirmed, low-frequency gravitational-wave echoes would be in tension with all current neutron-star models and would have dramatic implications for nuclear physics and gravity.

Several quantum-gravity models of compact objects predict microscopic or even Planckian corrections at the horizon scale. We explore the possibility of measuring two model-independent, smoking-gun effects of these corrections in the gravitational waveform of a compact binary, namely the absence of tidal heating and the presence of tidal deformability. For events detectable by the future space-based interferometer LISA, we show that the effect of tidal heating dominates and allows one to constrain putative corrections down to the Planck scale. The measurement of the tidal Love numbers with LISA is more challenging but, in optimistic scenarios, it allows to constrain the compactness of a supermassive exotic compact object down to the Planck scale. Our analysis suggests that highly-spinning, supermassive binaries at 1-20 Gpc provide unparalleled tests of quantum-gravity effects at the horizon scale.

Ultralight bosons and axion-like particles appear naturally in different scenarios and could solve some long-standing puzzles. Their detection is challenging, and all direct methods hinge on unknown couplings to the Standard Model of particle physics. However, the universal coupling to gravity provides model-independent signatures for these fields. We explore here the superradiant instability of spinning black holes triggered in the presence of such fields. The instability taps angular momentum from and limits the maximum spin of astrophysical black holes. We compute, for the first time, the spectrum of the most unstable modes of a massive vector (Proca) field for generic black-hole spin and Proca mass. The observed stability of the inner disk of stellar-mass black holes can be used to derive \emph{direct} constraints on the mass of dark photons in the mass range $ 10^{-13}\,{\rm eV}\lesssim m_V \lesssim 3\times 10^{-12}\,{\rm eV}$. By including also higher azimuthal modes, similar constraints apply to axion-like particles in the mass range $6\times10^{-13}\,{\rm eV}\lesssim m_{\rm ALP} \lesssim 10^{-11}\, {\rm eV}$. Likewise, mass and spin distributions of supermassive BHs --~as measured through continuum fitting, K$\alpha$ iron line, or with the future space-based gravitational-wave detector LISA~-- imply indirect bounds in the mass range approximately $10^{-19}\,{\rm eV}\lesssim m_V, m_{\rm ALP} \lesssim 10^{-13}\, {\rm eV}$, for both axion-like particles and dark photons. Overall, superradiance allows to explore a region of approximately $8$ orders of magnitude in the mass of ultralight bosons.

Gravitational-wave astronomy can give us access to the structure of black holes, potentially probing microscopic or even Planckian corrections at the horizon scale, as those predicted by some quantum-gravity models of exotic compact objects. A generic feature of these models is the replacement of the horizon by a reflective surface. Objects with these properties are prone to the so-called ergoregion instability when they spin sufficiently fast. We investigate in detail a simple model consisting of scalar perturbations of a Kerr geometry with a reflective surface near the horizon. The instability depends on the spin, on the compactness, and on the reflectivity at the surface. The instability time scale increases only logarithmically in the black-hole limit and, for a perfectly reflecting object, this is not enough to prevent the instability from occurring on dynamical time scales. However, we find that an absorption rate at the surface as small as 0.4% (reflectivity coefficient as large as $|{\cal R}|^2=0.996$) is sufficient to quench the instability completely. Our results suggest that exotic compact objects are not necessarily ruled out by the ergoregion instability.

Ultralight bosons can induce superradiant instabilities in spinning black holes, tapping their rotational energy to trigger the growth of a bosonic condensate. Possible observational imprints of these boson clouds include (i) direct detection of the nearly monochromatic (resolvable or stochastic) gravitational waves emitted by the condensate, and (ii) statistically significant evidence for the formation of "holes" at large spins in the spin versus mass plane (sometimes also referred to as "Regge plane") of astrophysical black holes. In this work, we focus on the prospects of LISA and LIGO detecting or constraining scalars with mass in the range $m_s\in [10^{-19},\,10^{-15}]$ eV and $m_s\in [10^{-14},\,10^{-11}]$ eV, respectively. Using astrophysical models of black-hole populations calibrated to observations and black-hole perturbation theory calculations of the gravitational emission, we find that, in optimistic scenarios, LIGO could observe a stochastic background of gravitational radiation in the range $m_s\in [2\times 10^{-13}, 10^{-12}]$ eV, and up to $10^4$ resolvable events in a $4$-year search if $m_s\sim 3\times 10^{-13}\,{\rm eV}$. LISA could observe a stochastic background for boson masses in the range $m_s\in [5\times 10^{-19}, 5\times 10^{-16}]$, and up to $\sim 10^3$ resolvable events in a $4$-year search if $m_s\sim 10^{-17}\,{\rm eV}$. LISA could further measure spins for black-hole binaries with component masses in the range $[10^3, 10^7]~M_\odot$, which is not probed by traditional spin-measurement techniques. A statistical analysis of the spin distribution of these binaries could either rule out scalar fields in the mass range $\sim [4 \times 10^{-18}, 10^{-14}]$ eV, or measure $m_s$ with ten percent accuracy if light scalars in the mass range $\sim [10^{-17}, 10^{-13}]$ eV exist.

Ultralight scalar fields around spinning black holes can trigger superradiant instabilities, forming a long-lived bosonic condensate outside the horizon. We use numerical solutions of the perturbed field equations and astrophysical models of massive and stellar-mass black hole populations to compute, for the first time, the stochastic gravitational-wave background from these sources. In optimistic scenarios the background is observable by Advanced LIGO and LISA for field masses $m_s$ in the range $\sim [2\times 10^{-13}, 10^{-12}]\,{\rm eV}$ and $\sim 5\times[ 10^{-19}, 10^{-16}]\,{\rm eV}$, respectively, and it can affect the detectability of resolvable sources. Our estimates suggest that an analysis of the stochastic background limits from LIGO O1 might already be used to marginally exclude axions with mass $\sim 10^{-12.5}{\rm eV}$. Semicoherent searches with Advanced LIGO (LISA) should detect $\sim 15~(5)$ to $200~(40)$ resolvable sources for scalar field masses $3\times 10^{-13}$ ($10^{-17}$) eV. LISA measurements of massive BH spins could either rule out bosons in the range $\sim [10^{-18}, 2\times 10^{-13}]$ eV, or measure $m_s$ with ten percent accuracy in the range $\sim[10^{-17}, 10^{-13}]$ eV.

The existence of black holes and of spacetime singularities is a fundamental issue in science. Despite this, observations supporting their existence are scarce, and their interpretation unclear. We overview how strong a case for black holes has been made in the last few decades, and how well observations adjust to this paradigm. Unsurprisingly, we conclude that observational proof for black holes is impossible to come by. However, just like Popper's black swan, alternatives can be ruled out or confirmed to exist with a single observation. These observations are within reach. In the next few years and decades, we will enter the era of precision gravitational-wave physics with more sensitive detectors. Just as accelerators require larger and larger energies to probe smaller and smaller scales, more sensitive gravitational-wave detectors will be probing regions closer and closer to the horizon, potentially reaching Planck scales and beyond. What may be there, lurking?

Future very-large-area X-ray instruments (for which the effective area is larger than $>3$m$^2$) will be able to measure the frequencies of quasi-periodic oscillations~(QPOs) observed in the X-ray flux from accreting compact objects with sub-percent precision. If correctly modeled, QPOs can provide a novel way to test the strong-field regime of gravity. By using the relativistic precession model and a modified version of the epicyclic resonance model, we develop a method to test general relativity against a generic class of theories with quadratic curvature corrections. With the instrumentation being studied for future missions such as eXTP, LOFT, or STROBE-X, a measurement of at least two QPO triplets from a stellar mass black hole can set stringent constraints on the coupling parameters of quadratic gravity.

Light bosonic fields are ubiquitous in extensions of the Standard Model. Even when minimally coupled to gravity, these fields might evade the assumptions of the black-hole no-hair theorems and give rise to spinning black holes which can be drastically different from the Kerr metric. Furthermore, they allow for self-gravitating compact solitons, known as (scalar or Proca) boson stars. The quasi-periodic oscillations (QPOs) observed in the X-ray flux emitted by accreting compact objects carry information about the strong-field region, thus providing a powerful tool to constrain deviations from Kerr's geometry and to search for exotic compact objects. By using the relativistic precession model as a proxy to interpret the QPOs in terms of geodesic frequencies, we investigate how the QPO frequencies could be used to test the no-hair theorem and the existence of light bosonic fields near accreting compact objects. We show that a detection of two QPO triplets with current sensitivity can already constrain these models, and that the future eXTP mission or a LOFT-like mission can set very stringent constraints on black holes with bosonic hair and on (scalar or Proca) boson stars. The peculiar geodesic structure of compact scalar/Proca boson stars implies that these objects can easily be ruled out as alternative models for X-ray source GRO J1655-40.

In the presence of massive bosonic degrees of freedom, rotational superradiance can trigger an instability that spins down black holes. This leads to peculiar gravitational-wave signatures and distribution in the spin-mass plane, which in turn can impose stringent constraints on ultralight fields. Here, we demonstrate that there is an analogous spindown effect for conducting stars. We show that rotating stars amplify low frequency electromagnetic waves, and that this effect is largest when the time scale for conduction within the star is of the order of a light crossing time. This has interesting consequences for dark photons, as massive dark photons would cause stars to spin down due to superradiant instabilities. The time scale of the spindown depends on the mass of the dark photon, and on the rotation rate, compactness, and conductivity of the star. Existing measurements of the spindown rate of pulsars place direct constraints on models of dark sectors. Our analysis suggests that dark photons of mass $m_V \sim 10^{-12}$ eV are excluded by pulsar-timing observations. These constraints also exclude superradiant instabilities triggered by dark photons as an explanation for the spin limit of observed pulsars.

The tidal Love numbers (TLNs) encode the deformability of a self-gravitating object immersed in a tidal environment and depend significantly both on the object's internal structure and on the dynamics of the gravitational field. An intriguing result in classical general relativity is the vanishing of the TLNs of black holes. We extend this result in three ways, aiming at testing the nature of compact objects: (i) we compute the TLNs of exotic compact objects, including different families of boson stars, gravastars, wormholes, and other toy models for quantum corrections at the horizon scale. In the black-hole limit, we find a universal logarithmic dependence of the TLNs on the location of the surface; (ii) we compute the TLNs of black holes beyond vacuum general relativity, including Einstein-Maxwell, Brans-Dicke and Chern-Simons gravity; (iii) We assess the ability of present and future gravitational-wave detectors to measure the TLNs of these objects, including the first analysis of TLNs with LISA. Both LIGO, ET and LISA can impose interesting constraints on boson stars, while LISA is able to probe even extremely compact objects. We argue that the TLNs provide a smoking gun of new physics at the horizon scale, and that future gravitational-wave measurements of the TLNs in a binary inspiral provide a novel way to test black holes and general relativity in the strong-field regime.

The multipole moments and the tidal Love numbers of neutron stars and quark stars satisfy certain relations which are almost insensitive to the star's internal structure. A natural question is whether the same relations hold for different compact objects and how they possibly approach the black-hole limit. Here we consider "gravastars", which are hypothetical compact objects sustained by their internal vacuum energy. Such solutions have been proposed as exotic alternatives to the black-hole paradigm because they can be as compact as black holes and exist in any mass range. By constructing slowly-rotating, thin-shell gravastars to quadratic order in the spin, we compute the moment of inertia $I$, the mass quadrupole moment $Q$, and the tidal Love number $\lambda$ in exact form. The $I$-$\lambda$-$Q$ relations of a gravastar are dramatically different from those of an ordinary compact star, but the black-hole limit is continuous, i.e. these quantities approach their Kerr counterparts when the compactness is maximum. Therefore, such relations can be used to discern a gravastar from an ordinary compact star, but not to break the degeneracy with the black-hole case. Based on these results, we conjecture that the full multipolar structure and the tidal deformability of a spinning, ultracompact gravastar are identical to those of a Kerr black hole. The approach to the black-hole limit is nonpolynomial, thus differing from the critical behavior recently found for strongly anisotropic neutron stars.

Gravitational waves emitted by distorted black holes---such as those arising from the coalescence of two neutron stars or black holes---carry not only information about the corresponding spacetime but also about the underlying theory of gravity. Although general relativity remains the simplest, most elegant and viable theory of gravitation, there are generic and robust arguments indicating that it is not the ultimate description of the gravitational universe. Here we focus on a particularly appealing extension of general relativity, which corrects Einstein's theory through the addition of terms which are second order in curvature: the topological Gauss-Bonnet invariant coupled to a dilaton. We study gravitational-wave emission from black holes in this theory, and {\bf(i)} find strong evidence that black holes are linearly (mode) stable against both axial and polar perturbations; {\bf(ii)} discuss how the quasinormal modes of black holes can be excited during collisions involving black holes, and finally {\bf(iii)} show that future ringdown detections with large signal-to-noise ratio would improve current constraints on the coupling parameter of the theory.

Gravitational waves from binary coalescences provide one of the cleanest signatures of the nature of compact objects. It has been recently argued that the post-merger ringdown waveform of exotic ultracompact objects is initially identical to that of a black-hole, and that putative corrections at the horizon scale will appear as secondary pulses after the main burst of radiation. Here we extend this analysis in three important directions: (i) we show that this result applies to a large class of exotic compact objects with a photon sphere for generic orbits in the test-particle limit; (ii) we investigate the late-time ringdown in more detail, showing that it is universally characterized by a modulated and distorted train of "echoes" of the modes of vibration associated with the photon sphere; (iii) we study for the first time equal-mass, head-on collisions of two ultracompact boson stars and compare their gravitational-wave signal to that produced by a pair of black-holes. If the initial objects are compact enough as to mimic a binary black-hole collision up to the merger, the final object exceeds the maximum mass for boson stars and collapses to a black-hole. This suggests that - in some configurations - the coalescence of compact boson stars might be almost indistinguishable from that of black-holes. On the other hand, generic configurations display peculiar signatures that can be searched for in gravitational-wave data as smoking guns of exotic compact objects.

We report on the Parallel Session BH4 "Gravitational fields with sources: From compact objects to black holes" of the 14th Marcel Grossmann Meeting held at Sapienza University of Rome in 2015.

It is commonly believed that the ringdown signal from a binary coalescence provides a conclusive proof for the formation of an event horizon after the merger. This expectation is based on the assumption that the ringdown waveform at intermediate times is dominated by the quasinormal modes of the final object. We point out that this assumption should be taken with great care, and that very compact objects with a light ring will display a similar ringdown stage, even when their quasinormal-mode spectrum is completely different from that of a black hole. In other words, universal ringdown waveforms indicate the presence of light rings, rather than of horizons. Only precision observations of the late-time ringdown signal, where the differences in the quasinormal-mode spectrum eventually show up, can be used to rule out exotic alternatives to black holes and to test quantum effects at the horizon scale.

The moment of inertia, the spin-induced quadrupole moment, and the tidal Love number of neutron-star and quark-star models are related through some relations which depend only mildly on the stellar equation of state. These "I-Love-Q" relations have important implications for astrophysics and gravitational-wave astronomy. An interesting problem is whether similar relations hold for other compact objects and how they approach the black-hole limit. To answer these questions, here we investigate the deformation properties of a large class of thin-shell gravastars, which are exotic compact objects that do not possess an event horizon nor a spacetime singularity. Working in a small-spin and small-tidal field expansion, we calculate the moment of inertia, the quadrupole moment, and the (quadrupolar electric) tidal Love number of gravastars with a polytropic thin shell. The I-Love-Q relations of a thin-shell gravastar are drastically different from those of an ordinary neutron star. The Love number and quadrupole moment are negative for less compact models and the I-Love-Q relations continuously approach the black-hole limit. We consider a variety of polytropic equations of state for the matter shell, and find no universality in the I-Love-Q relations. However, we cannot deny the possibility that, similarly to the neutron-star case, an approximate universality might emerge for a limited class of equations of state. Finally, we discuss how a measurement of the tidal deformability from the gravitational-wave detection of a compact-binary inspiral can be used to constrain exotic compact objects like gravastars.

During the motion of a binary pulsar around the Galactic center, the pulsar and its companion experience a wind of dark-matter particles that can affect the orbital motion through dynamical friction. We show that this effect produces a characteristic seasonal modulation of the orbit and causes a secular change of the orbital period whose magnitude can be well within the astonishing precision of various binary-pulsar observations. Our analysis is valid for binary systems with orbital period longer than a day. By comparing this effect with pulsar-timing measurements, it is possible to derive model-independent upper bounds on the dark-matter density at different distances $D$ from the Galactic center. For example, the precision timing of J1713+0747 imposes $\rho_{\rm DM}\lesssim 10^5\,{\rm GeV/cm}^3$ at $D\approx7\,{\rm kpc}$. The detection of a binary pulsar at $D\lesssim 10\,{\rm pc}$ could provide stringent constraints on dark-matter halo profiles and on growth models of the central black hole. The Square Kilometer Array can improve current bounds by 2 orders of magnitude, potentially constraining the local density of dark matter to unprecedented levels.

By extending our recent framework to describe the tidal deformations of a spinning compact object, we compute for the first time the tidal Love numbers of a spinning neutron star to linear order in the angular momentum. The spin of the object introduces couplings between electric and magnetic distortions and new classes of spin-induced ("rotational") tidal Love numbers emerge. We focus on stationary tidal fields, which induce axisymmetric perturbations. We present the perturbation equations for both electric-led and magnetic-led rotational Love numbers for generic multipoles and explicitly solve them for various tabulated equations of state and for a tidal field with an electric (even parity) and magnetic (odd parity) component with $\ell=2,3,4$. For a binary system close to the merger, various components of the tidal field become relevant. In this case we find that an octupolar magnetic tidal field can significantly modify the mass quadrupole moment of a neutron star. Preliminary estimates, assuming a spin parameter $\chi\approx0.05$, show modifications $\gtrsim10\%$ relative to the static case, at an orbital distance of five stellar radii. Furthermore, the rotational Love numbers as functions of the moment of inertia are much more sensitive to the equation of state than in the static case, where approximate universal relations at the percent level exist. For a neutron-star binary approaching the merger, we estimate that the approximate universality of the induced mass quadrupole moment deteriorates from $1\%$ in the static case to roughly $6\%$ when $\chi\approx0.05$. Our results suggest that spin-tidal couplings can introduce important corrections to the gravitational waveforms of spinning neutron-star binaries approaching the merger.

One century after its formulation, Einstein's general relativity has made remarkable predictions and turned out to be compatible with all experimental tests. Most of these tests probe the theory in the weak-field regime, and there are theoretical and experimental reasons to believe that general relativity should be modified when gravitational fields are strong and spacetime curvature is large. The best astrophysical laboratories to probe strong-field gravity are black holes and neutron stars, whether isolated or in binary systems. We review the motivations to consider extensions of general relativity. We present a (necessarily incomplete) catalog of modified theories of gravity for which strong-field predictions have been computed and contrasted to Einstein's theory, and we summarize our current understanding of the structure and dynamics of compact objects in these theories. We discuss current bounds on modified gravity from binary pulsar and cosmological observations, and we highlight the potential of future gravitational wave measurements to inform us on the behavior of gravity in the strong-field regime.

Gravitational theories with multiple scalar fields coupled to the metric and each other --- a natural extension of the well studied single-scalar-tensor theories --- are interesting phenomenological frameworks to describe deviations from general relativity in the strong-field regime. In these theories, the $N$-tuple of scalar fields takes values in a coordinate patch of an $N$-dimensional Riemannian target-space manifold whose properties are poorly constrained by weak-field observations. Here we introduce for simplicity a non-trivial model with two scalar fields and a maximally symmetric target-space manifold. Within this model we present a preliminary investigation of spontaneous scalarization for relativistic, perfect fluid stellar models in spherical symmetry. We find that the scalarization threshold is determined by the eigenvalues of a symmetric scalar-matter coupling matrix, and that the properties of strongly scalarized stellar configurations additionally depend on the target-space curvature radius. In preparation for numerical relativity simulations, we also write down the $3+1$ decomposition of the field equations for generic tensor-multi-scalar theories.

The deformability of a compact object induced by a perturbing tidal field is encoded in the tidal Love numbers, which depend sensibly on the object's internal structure. These numbers are known only for static, spherically-symmetric objects. As a first step to compute the tidal Love numbers of a spinning compact star, here we extend powerful perturbative techniques to compute the exterior geometry of a spinning object distorted by an axisymmetric tidal field to second order in the angular momentum. The spin of the object introduces couplings between electric and magnetic deformations and new classes of induced Love numbers emerge. For example, a spinning object immersed in a quadrupolar, electric tidal field can acquire some induced mass, spin, quadrupole, octupole and hexadecapole moments to second order in the spin. The deformations are encoded in a set of inhomogeneous differential equations which, remarkably, can be solved analytically in vacuum. We discuss certain subtleties in defining the multipole moments of the central object, which are due to the difficulty in separating the tidal field from the linear response of the object in the solution. By extending the standard procedure to identify the linear response in the static case, we prove analytically that the Love numbers of a Kerr black hole remain zero to second order in the spin. As a by-product, we provide the explicit form for a slowly-rotating, tidally-deformed Kerr black hole to quadratic order in the spin, and discuss its geodesic and geometrical properties.

Quasi-Periodic Oscillations (QPOs) observed in the X-ray flux emitted by accreting black holes, are associated to phenomena occurring near the horizon. Future very large area X-ray instruments will be able to measure QPO frequencies with very high precision, thus probing this strong-field region. By using the relativistic precession model, we show the way in which QPO frequencies could be used to test general relativity against those alternative theories of gravity which predict deviations from the classical theory in the strong-field regime. We consider one of the best motivated strong-curvature corrections to general relativity, namely the Einstein-Dilaton-Gauss-Bonnet theory, and show that a detection of QPOs with the expected sensitivity of the proposed ESA M-class mission LOFT would set the most stringent constraints on the parameter space of this theory.

The upcoming detection of gravitational waves by terrestrial interferometers will usher in the era of gravitational-wave astronomy. This will be particularly true when space-based detectors will come of age and measure the mass and spin of massive black holes with exquisite precision and up to very high redshifts, thus allowing for better understanding of the symbiotic evolution of black holes with galaxies, and for high-precision tests of General Relativity in strong-field, highly dynamical regimes. Such ambitious goals require that astrophysical environmental pollution of gravitational-wave signals be constrained to negligible levels, so that neither detection nor estimation of the source parameters are significantly affected. Here, we consider the main sources for space-based detectors -- the inspiral, merger and ringdown of massive black-hole binaries and extreme mass-ratio inspirals -- and account for various effects on their gravitational waveforms, including electromagnetic fields, cosmological evolution, accretion disks, dark matter, "firewalls" and possible deviations from General Relativity. We discover that the black-hole quasinormal modes are sharply different in the presence of matter, but the ringdown signal observed by interferometers is typically unaffected. The effect of accretion disks and dark matter depends critically on their geometry and density profile, but is negligible for most sources, except for few special extreme mass-ratio inspirals. Electromagnetic fields and cosmological effects are always negligible. We finally explore the implications of our findings for proposed tests of General Relativity with gravitational waves, and conclude that environmental effects will not prevent the development of precision gravitational-wave astronomy.

Superradiant instabilities of spinning black holes can be used to impose strong constraints on ultralight bosons, thus turning black holes into effective particle detectors. However, very little is known about the development of the instability and whether its nonlinear time evolution accords to the linear intuition. For the first time, we attack this problem by studying the impact of gravitational-wave emission and gas accretion on the evolution of the instability. Our quasi-adiabatic, fully-relativistic analysis shows that: (i) gravitational-wave emission does not have a significant effect on the evolution of the black hole, (ii) accretion plays an important role and (iii) although the mass of the scalar cloud developed through superradiance can be a sizeable fraction of the black-hole mass, its energy-density is very low and backreaction is negligible. Thus, massive black holes are well described by the Kerr geometry even if they develop bosonic clouds through superradiance. Using Monte Carlo methods and very conservative assumptions, we provide strong support to the validity of the linearized analysis and to the bounds of previous studies.

We construct models of slowly rotating, perfect-fluid neutron stars by extending the classical Hartle-Thorne formalism to generic scalar-tensor theories of gravity. Working at second order in the dimensionless angular momentum, we compute the mass M, radius R, scalar charge q, moment of inertia I and spin-induced quadrupole moment Q, as well as the tidal and rotational Love numbers. Our formalism applies to generic scalar-tensor theories, but we focus in particular on theories that allow for spontaneous scalarization. It was recently discovered that the moment of inertia, quadrupole moment and Love numbers are connected by approximately universal (i.e., equation-of-state independent) "I-Love-Q" relations. We find that similar relations hold also for spontaneously scalarized stars. More interestingly, the I-Love-Q relations in scalar-tensor theories coincide with the general relativistic ones within less than a few percent, even for spontaneously scalarized stars with the largest couplings allowed by current binary-pulsar constraints. This implies that astrophysical measurements of these parameters cannot be used to discriminate between general relativity and scalar-tensor theories, even if spontaneous scalarization occurs in nature. Because of the well known equivalence between f(R) theories and scalar-tensor theories, the theoretical framework developed in this paper can be used to construct rotating compact stellar models in f(R) gravity. Our slow-rotation expansion can also be used as a benchmark for numerical calculations of rapidly spinning neutron stars in generic scalar-tensor theories.

Ultracompact objects are self-gravitating systems with a light ring. It was recently suggested that fluctuations in the background of these objects are extremely long-lived and might turn unstable at the nonlinear level, if the object is not endowed with a horizon. If correct, this result has important consequences: objects with a light ring are black holes. In other words, the nonlinear instability of ultracompact stars would provide a strong argument in favor of the "black hole hypothesis," once electromagnetic or gravitational-wave observations confirm the existence of light rings. Here we explore in some depth the mode structure of ultracompact stars, in particular constant-density stars and gravastars. We show that the existence of very long-lived modes -- localized near a second, stable null geodesic -- is a generic feature of gravitational perturbations of such configurations. Already at the linear level, such modes become unstable if the object rotates sufficiently fast to develop an ergoregion. Finally, we conjecture that the long-lived modes become unstable under fragmentation via a Dyson-Chandrasekhar-Fermi mechanism at the nonlinear level. Depending on the structure of the star, it is also possible that nonlinearities lead to the formation of small black holes close to the stable light ring. Our results suggest that the mere observation of a light ring is a strong evidence for the existence of black holes.

Magnetic fields surrounding spinning black holes can confine radiation and trigger superradiant instabilities. To investigate this effect, we perform the first fully-consistent linear analysis of the Ernst spacetime, an exact solution of the Einstein--Maxwell equations describing a black hole immersed in a uniform magnetic field $B$. In the limit in which the black-hole mass vanishes, the background reduces to the marginally stable Melvin spacetime. The presence of an event horizon introduces a small dissipative term, resulting in a set of long-lived -- or unstable -- modes. We provide a simple interpretation of the mode spectrum in terms of a small perfect absorber immersed in a confining box of size $\sim1/B$ and show that rotation triggers a superradiant instability. By studying scalar perturbations of a magnetized Kerr--Newman black hole, we are able to confirm and quantify the details of this instability. The instability time scale can be orders of magnitude shorter than that associated to massive bosonic fields. The instability extracts angular momentum from the event horizon, competing against accretion. This implies that strong magnetic fields set an upper bound on the black-hole spin. Conversely, observations of highly-spinning massive black holes impose an intrinsic limit to the strength of the surrounding magnetic field. We discuss the astrophysical implications of our results and the limitations of the Ernst spacetime to describe realistic astrophysical configurations.

In a close encounter with a neutron star, a primordial black hole can get gravitationally captured by depositing a considerable amount of energy into nonradial stellar modes of very high angular number $l$. If the neutron-star equation of state is sufficiently stiff, we show that the total energy loss in the point-particle approximation is formally divergent. Various mechanisms -including viscosity, finite-size effects and the elasticity of the crust- can damp high-$l$ modes and regularize the total energy loss. Within a short time, the black hole is trapped inside the star and disrupts it by rapid accretion. Estimating these effects, we predict that the existence of old neutron stars in regions where the dark-matter density rho_{DM}>10^2 sigma/(km/s) GeV/cm^3 (where sigma is the dark-matter velocity dispersion) limits the abundance of primordial black holes in the mass range 10^{17} g < m_{PBH} < 10^{24} g, which was previously unconstrained. In combination with existing limits, our results suggest that primordial black holes cannot be the dominant dark matter constituent.

[abridged abstract] No, within a broad class of scenarios. With the advent of gravitational-wave (GW) astronomy, environmental effects on the GW signal will eventually have to be quantified. Here we present a wide survey of the corrections due to these effects in two situations of great interest for GW astronomy: the black hole (BH) ringdown emission and the inspiral of two compact objects. We take into account various effects such as: electric charges, magnetic fields, cosmological evolution, possible deviations from General Relativity, firewalls, and various forms of matter such as accretion disks and dark matter halos. Our analysis predicts the existence of resonances dictated by the external mass distribution, which dominate the very late-time behavior of merger/ringdown waveforms. The mode structure can drastically differ from the vacuum case, yet the BH response to external perturbations is unchanged at the time scales relevant for detectors. This is because although the vacuum Schwarzschild resonances are no longer quasinormal modes of the system, they still dominate the response at intermediate times. Our results strongly suggest that both parametrized and ringdown searches should use at least two-mode templates. Our analysis of compact binaries shows that environmental effects are typically negligible for most eLISA sources, with the exception of very few special extreme mass ratio inspirals. We show in particular that accretion and hydrodynamic drag generically dominate over self-force effects for geometrically thin disks, whereas they can be safely neglected for geometrically thick disk environments, which are the most relevant for eLISA. Finally, we discuss how our ignorance of the matter surrounding compact objects implies intrinsic limits on the ability to constrain strong-field deviations from General Relativity.

In electrovacuum stationary, asymptotically flat black holes in scalar-tensor theories of gravity are described by the Kerr-Newman family of solutions, just as in general relativity. We show that there exist two mechanisms which can render Kerr black holes unstable when matter is present in the vicinity of the black hole, as this induces an effective mass for the scalar. The first mechanism is a tachyonic instability that appears when the effective mass squared is negative, triggering the development of scalar hair --- a black hole version of "spontaneous scalarization". The second instability is associated with superradiance and is present when the effective mass squared is positive and when the black hole spin exceeds a certain threshold. The second mechanism is also responsible for a resonant effect in the superradiant scattering of scalar waves, with amplification factors as large as 10^5 or more.

Modifications of General Relativity usually include extra dynamical degrees of freedom, which to date remain undetected. Here we explore the possibility of modifying Einstein's theory by adding solely nondynamical fields. With the minimal requirement that the theory satisfies the weak equivalence principle and admits a covariant Lagrangian formulation, we show that the field equations generically have to include higher-order derivatives of the matter fields. This has profound consequences for the viability of these theories. We develop a parametrization based on a derivative expansion and show that - to next to leading order - all theories are described by just two parameters. Our approach can be used to put stringent, theory-independent constraints on such theories, as we demonstrates using the Newtonian limit as an example.

We uncover two mechanisms that can render Kerr black holes unstable in scalar-tensor gravity, both associated to the presence of matter in the vicinity of the black hole and the fact that this introduces an effective mass for the scalar. Our results highlight the importance of understanding the structure of spacetime in realistic, astrophysical black holes in scalar-tensor theories.

We discuss recent developments in gravitational fields with sources, regular black holes, quasiblack holes, and analogue black holes, related to the talks presented at the corresponding Parallel Session AT3 of the 13th Marcel Grossmann Meeting.

No-hair theorems exclude the existence of nontrivial scalar and massive vector hair outside four-dimensional, static, asymptotically flat black-hole spacetimes. We show, by explicitly building nonlinear solutions, that black holes can support massive graviton hair in theories of massive gravity. These hairy solutions are, most likely, the generic end state of the recently discovered monopole instability of Schwarzschild black holes in massive graviton theories.

We consider the imprint of superradiant instabilities of nonevaporating primordial black holes (PBHs) on the spectrum of the cosmic microwave background (CMB). In the radiation dominated era, PBHs are surrounded by a roughly homogeneous cosmic plasma which endows photons with an effective mass through the plasma frequency. In this setting, spinning PBHs are unstable to a spontaneous spindown through the well-known "black-hole bomb" mechanism. At linear level, the photon density is trapped by the effective photon mass and grows exponentially in time due to superradiance. As the plasma density declines due to cosmic expansion, the associated energy around PBHs is released and dissipated in the CMB. We evaluate the resulting spectral distortions of the CMB in the redshift range 10^3 < z < 2x10^6. Using the existing COBE/FIRAS bounds on CMB spectral distortions, we derive upper limits on the fraction of dark matter that can be associated with spinning PBHs in the mass range 10^{-8}*Msun < M < 0.2*Msin. For maximally-spinning PBHs, our limits are much tighter than those derived from microlensing or other methods. Future data from the proposed PIXIE mission could improve our limits by several orders of magnitude.

Recent nonlinear completions of Fierz-Pauli theory for a massive spin-2 field include nonlinear massive gravity and bimetric theories. The spectrum of black-hole solutions in these theories is rich, and comprises the same vacuum solutions of Einstein's gravity enlarged to include a cosmological constant. It was recently shown that Schwarzschild (de Sitter) black holes in these theories are generically unstable against spherical perturbations. Here we show that a notable exception is partially massless gravity, where the mass of the graviton is fixed in terms of the cosmological constant by \mu^2=2\Lambda/3 and a new gauge invariance emerges. We find that general-relativity black holes are stable in this limit. Remarkably, the spectrum of massive gravitational perturbations is isospectral.

Neutron stars are extremely relativistic objects which abound in our universe and yet are poorly understood, due to the high uncertainty on how matter behaves in the extreme conditions which prevail in the stellar core. It has recently been pointed out that the moment of inertia I, the Love number lambda and the spin-induced quadrupole moment Q of an isolated neutron star, are related through functions which are practically independent of the equation of state. These surprising universal I-lambda-Q relations pave the way for a better understanding of neutron stars, most notably via gravitational-wave emission. Gravitational-wave observations will probe highly-dynamical binaries and it is important to understand whether the universality of the I-lambda-Q relations survives strong-field and finite-size effects. We apply a Post-Newtonian-Affine approach to model tidal deformations in compact binaries and show that the I-lambda relation depends on the inspiral frequency, but is insensitive to the equation of state. We provide a fit for the universal relation, which is valid up to a gravitational wave frequency of ~900 Hz and accurate to within a few percent. Our results strengthen the universality of I-lambda-Q relations, and are relevant for gravitational-wave observations with advanced ground-based interferometers. We also discuss the possibility of using the Love-compactness relation to measure the neutron-star radius with an uncertainty of about 10% or smaller from gravitational-wave observations.

Massive bosonic fields of arbitrary spin are predicted by general extensions of the Standard Model. It has been recently shown that there exists a family of bimetric theories of gravity - including massive gravity - which are free of Boulware-Deser ghosts at the nonlinear level. This opens up the possibility to describe consistently the dynamics of massive spin-2 particles in a gravitational field. Within this context, we develop the study of massive spin-2 fluctuations - including massive gravitons - around Schwarzschild and slowly-rotating Kerr black holes. Our work has two important outcomes. First, we show that the Schwarzschild geometry is linearly unstable for small tensor masses, against a spherically symmetric mode. Second, we provide solid evidence that the Kerr geometry is also generically unstable, both against the spherical mode and against long-lived superradiant modes. In the absence of nonlinear effects, the observation of spinning black holes bounds the graviton mass to be smaller than 5x10^{-23} eV.

The nature and properties of dark matter (DM) are both outstanding issues in physics. Besides clustering in halos, the universal character of gravity implies that self-gravitating compact DM configurations might be spread throughout the universe. The astrophysical signature of these objects may be used to probe fundamental particle physics, or even to provide an alternative description of compact objects in active galactic nuclei. Here we discuss the most promising dissection tool of these configurations: the inspiral of a compact stellar-size object and consequent gravitational-wave emission. The inward motion of this "test probe" encodes unique information about the nature of the central, supermassive DM configuration. When the probe travels through some compact DM profile we show that, within a Newtonian approximation, the quasi-adiabatic evolution of the inspiral is mainly driven by DM accretion into the small compact object and by dynamical friction, rather than by gravitational-wave radiation-reaction. These effects circularize the orbits and leave a peculiar imprint on the gravitational waves emitted at late time. When accretion dominates, the frequency and the amplitude of the gravitational-wave signal produced during the latest stages of the inspiral are nearly constant. In the exterior region we study a relativistic model in which the inspiral is driven by the emission of gravitational and scalar waves. Resonances in the energy flux appear whenever the orbital frequency matches the mass of the DM particle and they correspond to the excitation of the central object's quasinormal frequencies. Unexpectedly, these resonances can lead to large dephasing with respect to standard inspiral templates, to such an extent as to prevent detection with matched filtering techniques. We discuss some observational consequences of these effects for gravitational-wave detection.

Tidal effects have long ago locked the Moon in synchronous rotation with the Earth and progressively increase the Earth-Moon distance. This "tidal acceleration" hinges on dissipation. Binaries containing black holes may also be tidally accelerated, dissipation being caused by the event horizon - a flexible, viscous one-way membrane. In fact, this process is known for many years under a different guise: superradiance. In General Relativity, tidal acceleration is obscured by gravitational-wave emission. However, when coupling to light scalar degrees of freedom is allowed, an induced dipole moment produces a "polarization acceleration", which might be orders of magnitude stronger than tidal quadrupolar effects. Consequences for optical and gravitational-wave observations are intriguing and it is not impossible that imprints of such mechanism have already been observed.

Eddington-inspired Born-Infeld gravity was recently proposed as an alternative to general relativity that offers a resolution of spacetime singularities. The theory differs from Einstein's gravity only inside matter due to nondynamical degrees of freedom, and it is compatible with all current observations. We show that the theory is reminiscent of Palatini f(R) gravity and that it shares the same pathologies, such as curvature singularities at the surface of polytropic stars and unacceptable Newtonian limit. This casts serious doubts on its viability.

We discuss a general method to study linear perturbations of slowly rotating black holes which is valid for any perturbation field, and particularly advantageous when the field equations are not separable. As an illustration of the method we investigate massive vector (Proca) perturbations in the Kerr metric, which do not appear to be separable in the standard Teukolsky formalism. Working in a perturbative scheme, we discuss two important effects induced by rotation: a Zeeman-like shift of nonaxisymmetric quasinormal modes and bound states with different azimuthal number m, and the coupling between axial and polar modes with different multipolar index l. We explicitly compute the perturbation equations up to second order in rotation, but in principle the method can be extended to any order. Working at first order in rotation we show that polar and axial Proca modes can be computed by solving two decoupled sets of equations, and we derive a single master equation describing axial perturbations of spin s=0 and s=+-1. By extending the calculation to second order we can study the superradiant regime of Proca perturbations in a self-consistent way. For the first time we show that Proca fields around Kerr black holes exhibit a superradiant instability, which is significantly stronger than for massive scalar fields. Because of this instability, astrophysical observations of spinning black holes provide the tightest upper limit on the mass of the photon: mv<4x10^-20 eV under our most conservative assumptions. Spin measurements for the largest black holes could reduce this bound to mv<10^-22 eV or lower.

Generic extensions of the standard model predict the existence of ultralight bosonic degrees of freedom. Several ongoing experiments are aimed at detecting these particles or constraining their mass range. Here we show that massive vector fields around rotating black holes can give rise to a strong superradiant instability which extracts angular momentum from the hole. The observation of supermassive spinning black holes imposes limits on this mechanism. We show that current supermassive black hole spin estimates provide the tightest upper limits on the mass of the photon (mv<4x10^{-20} eV according to our most conservative estimate), and that spin measurements for the largest known supermassive black holes could further lower this bound to mv<10^{-22} eV. Our analysis relies on a novel framework to study perturbations of rotating Kerr black holes in the slow-rotation regime, that we developed up to second order in rotation, and that can be extended to other spacetime metrics and other theories.

A stellar-mass compact object spiraling into a supermassive black hole, an extreme-mass-ratio inspiral (EMRI), is one of the targets of future gravitational-wave detectors and it offers a unique opportunity to test General Relativity (GR) in the strong-field. We study whether generic scalar-tensor (ST) theories can be further constrained with EMRIs. We show that in the EMRI limit, all such theories universally reduce to massive or massless Brans-Dicke theory and that black holes do not emit dipolar radiation to all orders in post-Newtonian (PN) theory. For massless theories, we calculate the scalar energy flux in the Teukolsky formalism to all orders in PN theory and fit it to a high-order PN expansion. We derive the PN ST corrections to the Fourier transform of the gravitational wave response and map it to the parameterized post-Einsteinian framework. We use the effective-one-body framework adapted to EMRIs to calculate the ST modifications to the gravitational waveform. We find that such corrections are smaller than those induced in the early inspiral of comparable-mass binaries, leading to projected bounds on the coupling that are worse than current Solar System ones. Brans-Dicke theory modifies the weak-field, with deviations in the energy flux that are largest at small velocities. For massive theories, superradiance can lead to resonances in the scalar energy flux that can lead to floating orbits outside the innermost stable circular orbit and that last until the supermassive black hole loses enough mass and spin-angular momentum. If such floating orbits occur in the frequency band of LISA, they would lead to a large dephasing (~1e6 rads), preventing detection with GR templates. A detection that is consistent with GR would then rule out floating resonances at frequencies lower than the lowest observed frequency, allowing for the strongest constraints yet on massive ST theories.

In four-dimensional spacetime, moons around black holes generate low-amplitude tides, and the energy extracted from the hole's rotation is always smaller than the gravitational radiation lost to infinity. Thus, moons orbiting a black hole inspiral and eventually merge. However, it has been conjectured that in higher-dimensional spacetimes orbiting bodies generate much stronger tides, which backreact by tidally accelerating the body outwards. This effect, analogous to the tidal acceleration experienced by the Earth-Moon system, would determine the evolution of the binary. Here, we put this conjecture to the test, by studying matter coupled to a massless scalar field in orbit around a singly-spinning rotating black hole in higher dimensions. We show that in dimensions larger than five the energy extracted from the black hole through superradiance is larger than the energy carried out to infinity. Our numerical results are in excellent agreement with analytic approximations and lend strong support to the conjecture that tidal acceleration is the rule, rather than the exception, in higher dimensions. Superradiance dominates the energy budget and moons "outspiral"; for some particular orbital frequency, the energy extracted at the horizon equals the energy emitted to infinity and "floating orbits" generically occur. We give an interpretation of this phenomenon in terms of the membrane paradigm and of tidal acceleration due to energy dissipation across the horizon.

Viable corrections to the matter sector of Poisson's equation may result in qualitatively different astrophysical phenomenology, for example the gravitational collapse and the properties of compact objects can change drastically. We discuss a class of modified non-relativistic theories and focus on a relativistic completion, Eddington-inspired Born-Infeld gravity. This recently proposed theory is equivalent to General Relativity in vacuum, but its non-trivial coupling to matter prevents singularities in early cosmology and in the non-relativistic collapse of non-interacting particles. We extend our previous analysis, discussing further developments. We present a full numerical study of spherically symmetric non-relativistic gravitational collapse of dust. For any positive coupling, the final state of the collapse is a regular pressureless star rather than a singularity. We also argue that there is no Chandrasekhar limit for the mass of non-relativistic white dwarf in this theory. Finally, we extend our previous results in the fully relativistic theory by constructing static and slowly rotating compact stars governed by nuclear-physics inspired equations of state. In the relativistic theory, there exists an upper bound on the mass of compact objects, suggesting that black holes can still be formed in the relativistic collapse.

Physics in curved spacetime describes a multitude of phenomena, ranging from astrophysics to high energy physics. The last few years have witnessed further progress on several fronts, including the accurate numerical evolution of the gravitational field equations, which now allows highly nonlinear phenomena to be tamed. Numerical relativity simulations, originally developed to understand strong field astrophysical processes, could prove extremely useful to understand high-energy physics processes like trans-Planckian scattering and gauge-gravity dualities. We present a concise and comprehensive overview of the state-of-the-art and important open problems in the field(s), along with guidelines for the next years. This writeup is a summary of the "NR/HEP Workshop" held in Madeira, Portugal from August 31st to September 3rd 2011.

We study the coupling of massive scalar fields to matter in orbit around rotating black holes. It is generally expected that orbiting bodies will lose energy in gravitational waves, slowly inspiralling into the black hole. Instead, we show that the coupling of the field to matter leads to a surprising effect: because of superradiance, matter can hover into "floating orbits" for which the net gravitational energy loss at infinity is entirely provided by the black hole's rotational energy. Orbiting bodies remain floating until they extract sufficient angular momentum from the black hole, or until perturbations or nonlinear effects disrupt the orbit. For slowly rotating and nonrotating black holes floating orbits are unlikely to exist, but resonances at orbital frequencies corresponding to quasibound states of the scalar field can speed up the inspiral, so that the orbiting body "sinks". These effects could be a smoking gun of deviations from general relativity.

We present, in closed analytic form, a general stationary, slowly rotating black hole, which is solution to a large class of alternative theories of gravity in four dimensions. In these theories, the Einstein-Hilbert action is supplemented by all possible quadratic, algebraic curvature invariants coupled to a scalar field. The solution is found as a deformation of the Schwarzschild metric in General Relativity. We explicitly derive the changes to the orbital frequency at the innermost stable circular orbit and at the light ring in closed form. These results could be useful when comparing General Relativity against alternative theories by (say) measurements of X-ray emission in accretion disks, or by stellar motion around supermassive black holes. When gravitational-wave astronomy comes into force, strong constraints on the coupling parameters can in principle be made.

We propose a new approach to test possible corrections to Newtonian gravity using solar physics. The high accuracy of current solar models and new precise observations allow one to constrain corrections to standard gravity at unprecedented levels. Our case study is Eddington-inspired gravity, an attractive modified theory of gravity which results in non-singular cosmology and collapse. The theory is equivalent to standard gravity in vacuum, but it sensibly differs from it within matter, for instance it affects the evolution and the equilibrium structure of the Sun, giving different core temperature profiles, deviations in the observed acoustic modes and in solar neutrino fluxes. Comparing the predictions from a modified solar model with observations, we constrain the coupling parameter of the theory, |kappa_g| < 3x10^5 m^5 s^2 / kg. Our results show that the Sun can be used to efficiently constraint alternative theories of gravity.

We develop a theoretical framework to study slowly rotating compact stars in a rather general class of alternative theories of gravity, with the ultimate goal of investigating constraints on alternative theories from electromagnetic and gravitational-wave observations of compact stars. Our Lagrangian includes as special cases scalar-tensor theories (and indirectly f(R) theories) as well as models with a scalar field coupled to quadratic curvature invariants. As a first application of the formalism, we discuss (for the first time in the literature) compact stars in Einstein-Dilaton-Gauss-Bonnet gravity. We show that compact objects with central densities typical of neutron stars cannot exist for certain values of the coupling constants of the theory. In fact, the existence and stability of compact stars sets more stringent constraints on the theory than the existence of black hole solutions. This work is a first step in a program to systematically rule out (possibly using Bayesian model selection) theories that are incompatible with astrophysical observations of compact stars.

A new, Eddington inspired theory of gravity was recently proposed by Banados and Ferreira. It is equivalent to General Relativity in vacuum, but differs from it inside matter. This viable, one parameter theory was shown to avoid cosmological singularities and turns out to lead to many other exciting new features that we report here. First, for a positive coupling parameter, the field equations have a dramatic impact on the collapse of dust, and do not lead to singularities. We further find that the theory supports stable, compact pressureless stars made of perfect fluid, which provide interesting models of self-gravitating dark matter. Finally, we show that the mere existence of relativistic stars imposes a strong, near optimal constraint on the coupling parameter, which can even be improved by observations of the moment of inertia of the double pulsar.

Dynamical Chern-Simons gravity is an interesting extension of General Relativity, which finds its way in many different contexts, including string theory, cosmological settings and loop quantum gravity. In this theory, the gravitational field is coupled to a scalar field by a parity-violating term, which gives rise to characteristic signatures. Here we investigate how Chern-Simons gravity would affect the quasi-circular inspiralling of a small, stellar-mass object into a large non-rotating supermassive black hole, and the accompanying emission of gravitational and scalar waves. We find the relevant equations describing the perturbation induced by the small object, and we solve them through the use of Green's function techniques. Our results show that for a wide range of coupling parameters, the Chern-Simons coupling gives rise to an increase in total energy flux, which translates into a fewer number of gravitational-wave cycles over a certain bandwidth. For space-based gravitational-wave detectors such as LISA, this effect can be used to constrain the coupling parameter effectively.

Quantum fields in compact stars can be amplified due to a semiclassical instability. This generic feature of scalar fields coupled to curvature may affect the birth and the equilibrium structure of relativistic stars. We point out that the semiclassical instability has a classical counterpart, which occurs exactly in the same region of the parameter space. For negative values of the coupling parameter the instability is equivalent to the well-known "spontaneous scalarization" effect: the plausible end-state of the instability is a static, asymptotically flat equilibrium configuration with nonzero expectation value for the quantum fields, which is compatible with experiments in the weak-field regime and energetically favored over stellar solutions in general relativity. For positive values of the coupling parameter the new configurations are energetically disfavored, and the end-point of the instability remains an open and interesting issue. The vacuum instability may provide a natural mechanism to produce spontaneous scalarization, leading to new experimental opportunities to probe the nature of vacuum energy via astrophysical observations of compact stars.

Superspinars are ultracompact objects whose mass M and angular momentum J violate the Kerr bound (cJ/GM^2>1). Recent studies analyzed the observable consequences of gravitational lensing and accretion around superspinars in astrophysical scenarios. In this paper we investigate the dynamical stability of superspinars to gravitational perturbations, considering either purely reflecting or perfectly absorbing boundary conditions at the "surface" of the superspinar. We find that these objects are unstable independently of the boundary conditions, and that the instability is strongest for relatively small values of the spin. Also, we give a physical interpretation of the various instabilities that we find. Our results (together with the well-known fact that accretion tends to spin superspinars down) imply that superspinars are very unlikely astrophysical alternatives to black holes.

Dynamical Chern-Simons gravity is an extension of General Relativity in which the gravitational field is coupled to a scalar field through a parity-violating Chern-Simons term. In this framework, we study perturbations of spherically symmetric black hole spacetimes, assuming that the background scalar field vanishes. Our results suggest that these spacetimes are stable, and small perturbations die away as a ringdown. However, in contrast to standard General Relativity, the gravitational waveforms are also driven by the scalar field. Thus, the gravitational oscillation modes of black holes carry imprints of the coupling to the scalar field. This is a smoking gun for Chern-Simons theory and could be tested with gravitational-wave detectors, such as LIGO or LISA. For negative values of the coupling constant, ghosts are known to arise, and we explicitly verify their appearance numerically. Our results are validated using both time evolution and frequency domain methods.

Gravitational waves from compact objects provide information about their structure, probing deep into strong-gravity regions. Here we illustrate how the presence or absence of an event horizon can produce qualitative differences in the gravitational waves emitted by ultra-compact objects. In order to set up a straw-man ultra-compact object with no event horizon, but which is otherwise almost identical to a black hole, we consider a nonrotating thin-shell model inspired by Mazur and Mottola's gravastar, which has a Schwarzschild exterior, a de Sitter interior and an infinitely thin shell with finite tension separating the two regions. As viewed from the external space-time, the shell can be located arbitrarily close to the Schwarzschild radius, so a gravastar might seem indistinguishable from a black hole when tests are only performed on its external metric. We study the linearized dynamics of the system, and in particular the junction conditions connecting internal and external gravitational perturbations. As a first application of the formalism we compute polar and axial oscillation modes of a thin-shell gravastar. We show that the quasinormal mode spectrum is completely different from that of a black hole, even in the limit when the surface redshift becomes infinite. Polar QNMs depend on the equation of state of matter on the shell and can be used to distinguish between different gravastar models. Our calculations suggest that low-compactness gravastars could be unstable when the sound speed on the shell vs/c>0.92.

It is generally accepted that Einstein's theory will get some as yet unknown corrections, possibly large in the strong field regime. An ideal place to look for these modifications is around the vicinities of compact objects such as black holes. Our case study here are Dilatonic Black Holes, which arise in the framework of Gauss-Bonnet couplings and one-loop corrected four-dimensional effective theory of heterotic superstrings at low energies. These are interesting objects as a prototype for alternative, yet well-behaved gravity theories: they evade the "no-hair" theorem of General Relativity but were proved to be stable against radial perturbations. We investigate the viability of these black holes as astrophysical objects and try to provide some means to distinguish them from black holes in General Relativity. We start by extending previous works and establishing the stability of these black holes against axial perturbations. We then look for solutions of the field equations describing slowly rotating black holes and study geodesic motion around this geometry. Depending on the values of mass, dilaton charge and angular momentum of the solution, one can have measurable differences in the ISCO location and orbital frequency, relatively to black holes in General Relativity. Such differences may be useful in future experiments, to discriminate between alternative theories of gravity.

Viable alternatives to astrophysical black holes include hyper-compact objects without horizon, such as gravastars, boson stars, wormholes and superspinars. The authors have recently shown that typical rapidly-spinning gravastars and boson stars develop a strong instability. That analysis is extended in this paper to a wide class of horizonless objects with approximate Kerr-like geometry. A detailed investigation of wormholes and superspinars is presented, using plausible models and mirror boundary conditions at the surface. Like gravastars and boson stars, these objects are unstable with very short instability timescales. This result strengthens previous conclusions that observed hyper-compact astrophysical objects with large rotation are likely to be black holes.

Most of the properties of black holes can be mimicked by horizonless compact objects such as gravastars and boson stars. We show that these ultra-compact objects develop a strong ergoregion instability when rapidly spinning. Instability timescales can be of the order of 0.1 seconds to 1 week for objects with mass M=1-10^6 solar masses and angular momentum J> 0.4 M^2. This provides a strong indication that ultra-compact objects with large rotation are black holes. Explosive events due to ergoregion instability have a well-defined gravitational-wave signature. These events could be detected by next-generation gravitational-wave detectors such as Advanced LIGO or LISA.