Bosonic fields can give rise to self-gravitating structures. These are interesting hypothetical new "dark matter stars" and good descriptions of dark matter haloes if the fields are very light. We study the dynamical response of Newtonian boson stars (NBS) when excited by external matter (stars, planets or black holes) in their vicinities. Our setup can describe the interaction between a massive black hole and the surrounding environment, shortly after the massive body has undergone a "kick", due to the collapse of baryonic matter at the galactic center, or dark matter depletion as a reaction to an inspiralling binary. We perform the first self-consistent calculation of dynamical friction acting on moving bodies in these backgrounds. Binaries close to coalescence "stir" the NBS core, and backreaction affects gravitational waveforms at leading $-6PN$ order with respect to the dominant quadrupolar term; the coefficient is too small to allow detection by next-generation interferometers. We also show that the gravitational collapse to a supermassive black hole at the center of a NBS is accompanied by only a small change in the surrounding core. The NBS eventually gets accreted, but for astrophysical parameters this occurs only after several Hubble times.

The Galactic Center black hole Sagittarius A* is a variable NIR source that exhibits bright flux excursions called flares. The low-flux density turnover of the flux distribution is below the sensitivity of current single-aperture telescopes. We use the unprecedented resolution of the GRAVITY instrument at the VLTI. Our light curves are unconfused, overcoming the confusion limit of previous photometric studies. We analyze the light curves using standard statistical methods and obtain the flux distribution. We find that the flux distribution of SgrA* turns over at a median flux density of (1.1\pm0.3)mJy. We measure the percentiles of the flux distribution and use them to constrain the NIR K-band SED. Furthermore, we find that the flux distribution is intrinsically right-skewed to higher flux density in log space. Flux densities below 0.1mJy are hardly ever observed. In consequence, a single powerlaw or lognormal distribution does not suffice to describe the observed flux distribution in its entirety. However, if one takes into account a power law component at high flux densities, a lognormal distribution can describe the lower end of the observed flux distribution. We confirm the RMS-flux relation for Sgr~A* and find it to be linear for all flux densities in our observation. We conclude that Sgr~A* has two states: the bulk of the emission is generated in a lognormal process with a well-defined median flux density and this quiescent emission is supplemented by sporadic flares that create the observed power law extension of the flux distribution.

The star S2 orbiting the compact radio source Sgr A* is a precision probe of the gravitational field around the closest massive black hole (candidate). Over the last 2.7 decades we have monitored the star's radial velocity and motion on the sky, mainly with the SINFONI and NACO adaptive optics (AO) instruments on the ESO VLT, and since 2017, with the four-telescope interferometric beam combiner instrument GRAVITY. In this paper we report the first detection of the General Relativity (GR) Schwarzschild Precession (SP) in S2's orbit. Owing to its highly elliptical orbit (e = 0.88), S2's SP is mainly a kink between the pre-and post-pericentre directions of motion ~ +- 1 year around pericentre passage, relative to the corresponding Kepler orbit. The superb 2017-2019 astrometry of GRAVITY defines the pericentre passage and outgoing direction. The incoming direction is anchored by 118 NACO-AO measurements of S2's position in the infrared reference frame, with an additional 75 direct measurements of the S2-Sgr A* separation during bright states ('flares') of Sgr A*. Our 14-parameter model fits for the distance, central mass, the position and motion of the reference frame of the AO astrometry relative to the mass, the six parameters of the orbit, as well as a dimensionless parameter f_SP for the SP (f_SP = 0 for Newton and 1 for GR). From data up to the end of 2019 we robustly detect the SP of S2, del phi = 12' per orbital period. From posterior fitting and MCMC Bayesian analysis with different weighting schemes and bootstrapping we find f_SP = 1.10 +- 0.19. The S2 data are fully consistent with GR. Any extended mass inside S2's orbit cannot exceed ~ 0.1% of the central mass. Any compact third mass inside the central arcsecond must be less than about 1000 M_sun.

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.

The tidal deformability of compact objects by a companion has a detectable imprint in the gravitational waves emitted by a binary system. This effect is governed by the so-called tidal Love numbers. For a particular theory of gravity, these depend solely on the object internal structure and they vanish for black holes in general relativity. A measurement compatible with nonzero tidal Love numbers could thus provide evidence of new physics in the strong-field regime. However, in realistic astrophysical scenarios, compact objects are surrounded by a nonvacuum environment. In this work, we study the tidal deformability of configurations of black holes immersed in matter, focusing on two analytical models representing an anisotropic fluid and a thin-shell of matter around a black hole. We then apply our results to the astrophysically relevant case of a black hole surrounded by an accretion disk, in the parameter region of interest of the upcoming LISA mission. Our results indicate that there are challenges to overcome concerning tests of strong-field gravity using tidal Love numbers.

How well is the vacuum Kerr geometry a good description of the dark, compact objects in our universe? Precision measurements of accreting matter in the deep infrared and gravitational-wave measurements of coalescing objects are finally providing answers to this question. Here, we study the possibility of resonant excitation of the modes of the central object -- taken to be very compact but horizonless -- during an extreme-mass-ratio inspiral. We show that for very compact objects resonances are indeed excited. However, the impact of such excitation on the phase of the gravitational-wave signal is negligible, since resonances are crossed very quickly during inspiral.

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.

LISA can detect higher harmonics of the ringdown gravitational-wave signal from massive black-hole binary mergers with large signal-to-noise ratio. The most massive black-hole binaries are more likely to have electromagnetic counterparts, and the inspiral will contribute little to their signal-to-noise ratio. Here we address the following question: can we extract the binary parameters and localize the source using LISA observations of the ringdown only? Modulations of the amplitude and phase due to LISA's motion around the Sun can be used to disentangle the source location and orientation when we detect the long-lived inspiral signal, but they can not be used for ringdown-dominated signals, which are very short-lived. We show that (i) we can still measure the mass ratio and inclination of high-mass binaries by carefully combining multiple ringdown harmonics, and (ii) we can constrain the sky location and luminosity distance by relying on the relative amplitudes and phases of various harmonics, as measured in different LISA channels.

On May 29, 1919, at Ro\c{c}a Sundy, Principe island, Eddington confirms Einstein's general relativity theory for the first time by photographing stars behind the obscured Sun during a total eclipse. History was made. At Sobral, Eddington's astronomer colleagues photograph the same eclipse and also conclude that light from distant stars suffers a deflection when passing by the gravitational field of the Sun, in accordance with general relativity. With the confirmation of general relativity, a theory of gravitation at a fundamental level, physics became, once and for all, relativist and its future was outlined. The first world war had finished a few months before and the deep wounds between nations were yet to heal. Science wanted to be above it all showing that people could be united by a common goal. This year, the 100 years of this achievement was commemorated with the scientific conference "From Einstein and Eddington to LIGO: 100 years of gravitational light deflection" in Principe to celebrate this landmark event. We report here on this conference of celebration.

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

For the celebrations of the 100 years of the observations undertaken by Eddington at the island of Principe and collaborators at Sobral during a total solar eclipse in May 29, 1919, which have confirmed Einstein's theory of general relativity through the deflection of the incoming light from distant stars due to the spacetime curvature caused by the Sun, we highlight the main aspects of the theory, its tests and applications, focusing on some of its outstanding consequences. These are black holes, the object par excellence of general relativity, and gravitational waves, the gravitational probe for the distant Universe. We also point out some open issues.

Because of the very definition of black holes --- no light escapes them and falling objects get infinitely faint when approaching --- it is impossible to ever prove that they exist. However, electromagnetic and gravitational-wave observatories have now `seen' black holes. Datasets from these observations, released in 2019 and late 2018, give important hints about the environment, origin and growth of black holes.

We study the ability of gravitational-wave detectors to constrain the nature of the environment in which compact binaries evolve. We show that the strong dephasing induced by accretion and dynamical friction can constraint the density of the surrounding medium to orders of magnitude below that of accretion disks. Planned detectors, such as LISA or DECIGO, will be able to probe densities typical of those of dark matter.

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

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 consider the interaction between a plane wave and a (counter-moving) black hole. We show that energy is transferred from the black hole to the wave, giving rise to a negative absorption cross-section. Moving black holes absorb radiation and deposit energy in external radiation. Due to this effect, a black hole hole of mass M moving at relativistic speeds in a cold medium will appear to be surrounded by a bright "ring" of diameter 3\sqrt{3}GM/c^2 and thickness GM/c^2.

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.

Black holes are the simplest macroscopic objects, and provide unique tests of General Relativity. They have been compared to the Hydrogen atom in quantum mechanics. Here, we establish a few facts about the simplest systems bound by gravity: black hole binaries. We provide strong evidence for the existence of `global' photosurfaces surrounding the binary, and of binary quasinormal modes leading to exponential decay of massless fields when the binary spacetime is slightly perturbed. These two properties go hand in hand, as they did for isolated black holes. The binary quasinormal modes have high quality factor and may be prone to resonant excitations. Finally, we show that energy extraction from binaries is generic and we find evidence of a new mechanism -- akin to the Fermi acceleration process -- whereby the binary transfers energy to its surroundings in a cascading process. The mechanism is conjectured to work when the individual components spin, or are made of compact stars.

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 evolution of planets, stars and even galaxies is driven, to a large extent, by dynamical friction of gravitational origin. There is now a good understanding of the friction produced by extended media, either collisionless of fluid-like. However, the physics of accretion or protoplanetary disks, for instance, is described by slab-like geometries instead, compact in one spatial direction. Here, we find, for the first time, the gravitational wake due to a massive perturber moving through a slab-like medium, describing e.g. accretion disks with sharp transitions. We show that dynamical friction in such environments can be substantially reduced relatively to spatially extended profiles. Finally, we provide simple and accurate expressions for the gravitational drag force felt by the perturber, in both the subsonic and supersonic regime.

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?

The nature of dark matter is one of the longest-standing puzzles in science. Axions or axion-like particles are a key possibility, and arise in mechanisms to solve the strong CP problem but also in low-energy limits of string theory. Extensive experimental and observational efforts are actively looking for `axionic' imprints. Independently on their nature, their abundance, and on their contribution to the dark matter problem, axions form dense clouds around spinning black holes, grown by superradiant mechanisms. It was recently suggested that once couplings to photons are considered, an exponential (quantum) stimulated emission of photons ensues at large enough axion number. Here we solve numerically the classical problem in different setups. We show that laser-like emission from clouds exists at the classical level, and we provide the first quantitative description of the problem.

Gravitational-wave astronomy has the potential to substantially advance our knowledge of the cosmos, from the most powerful astrophysical engines to the initial stages of our universe. Gravitational waves also carry information about the nature of black holes. Here we investigate the potential of gravitational-wave detectors to test a proposal by Bekenstein and Mukhanov that the area of black hole horizons is quantized in units of the Planck area. Our results indicate that this quantization could have a potentially observable effect on the classical gravitational wave signals received by detectors. In particular, we find distorted gravitational-wave "echoes" in the post-merger waveform describing the inspiral and merger of two black holes. These echoes have a specific frequency content that is characteristic of black hole horizon area quantization.

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.

In general relativity, Maxwell's equations are embedded in curved spacetime through the minimal prescription, but this could change if strong-gravity modifications are present. We show that with a nonminimal coupling between gravity and a massless vector field, nonperturbative effects can arise in compact stars. We find solutions describing stars with nontrivial vector field configurations, some of which are associated to an instability, while others are not. The vector field can be interpreted either as the electromagnetic field, or as a hidden vector field weakly coupled with the standard model.

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.

The coupling between scalar and vector fields has a long and interesting history. Axions are one key possibility to solve the strong CP problem and axion-like particles could be one solution to the dark matter puzzle. Given the nature of the coupling, and the universality of free fall, nontrivial important effects are expected in regions where gravity is strong. Here, we show that i. A background EM field induces an axionic instability in flat space, for large enough electric fields. Conversely, a homogeneous harmonic axion field induces an instability in the Maxwell sector. When carried over to curved spacetime, this phenomena translates into generic instabilities of charged black holes (BHs). ii. In the presence of charge, BH uniqueness results are lost. We find solutions which are small deformations of the Kerr-Newman geometry and hairy stationary solutions without angular momentum, which are `dragged' by the axion. Axion fields must exist around spinning BHs if these are immersed in external magnetic fields. The axion profile can be obtained perturbatively from the electro-vacuum solution derived by Wald. iii. Ultralight axions trigger superradiant instabilities of spinning BHs and form an axionic cloud in the exterior geometry. The superradiant growth can be interrupted or suppressed through axionic or scalar couplings to EM. These couplings lead to periodic bursts of light, which occur throughout the history of energy extraction from the BH. We provide numerical and simple analytical estimates for the rates of these processes. iv. Finally, we discuss how plasma effects can affect the evolution of superradiant instabilities.

Effective field theory methods suggest that some rather-general extensions of General Relativity include, or are mimicked by, certain higher-order curvature corrections, with coupling constants expected to be small but otherwise arbitrary. Thus, the tantalizing prospect to test the fundamental nature of gravity with gravitational-wave observations, in a systematic way, emerges naturally. Here, we build black hole solutions in such a framework and study their main properties. Once rotation is included, we find the first purely gravitational example of geometries without $\mathbb{Z}_2$-symmetry. Despite the higher-order operators of the theory, we show that linearized fluctuations of such geometries obey second-order differential equations. We find nonzero tidal Love numbers. We study and compute the quasinormal modes of such geometries. These results are of interest to gravitational-wave science but also potentially relevant for electromagnetic observations of the galactic center or $X$-ray binaries.

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.

Closed photon orbits around isolated black holes are related to important aspects of black hole physics, such as strong lensing, absorption cross section of null particles and the way that black holes relax through quasinormal ringing. When two black holes are present -- such as during the inspiral and merger events of interest for gravitational-wave detectors -- the concept of closed photon orbits still exists, but its properties are basically unknown. With these applications in mind, we study here the closed photon orbits of two different static black hole binaries. The first one is the Majumdar-Papapetrou geometry describing two extremal, charged black holes in equilibrium, while the second one is the double sink solution of fluid dynamics, which describes (in a curved-spacetime language) two "dumb" holes. For the latter solution, we also characterize its dynamical response to external perturbations, and study how it relates to the photon orbits. In addition, we compute the ergoregion of such spacetime and show that it does not coincide with the event horizon.

It is well known that a spinning body moving in a fluid suffers a force orthogonal to its velocity and rotation axis --- it is called the Magnus effect. Recent simulations of spinning black holes and (indirect) theoretical predictions, suggest that a somewhat analogous effect may occur for purely gravitational phenomena. The magnitude and precise direction of this "gravitational Magnus effect" is still the subject of debate. Starting from the rigorous equations of motion for spinning bodies in General Relativity (Mathisson-Papapetrou equations), we show that indeed such an effect takes place and is a fundamental part of the spin-curvature force. The effect arises whenever there is a current of mass/energy, non-parallel to a body's spin. We compute the effect explicitly for some astrophysical systems of interest: a galactic dark matter halo, a black hole accretion disk, and the Friedmann-Lema\^itre-Robertson-Walker (FLRW) spacetime. It is seen to lead to secular orbital precessions potentially observable by future astrometric experiments and gravitational-wave detectors. Finally, we consider also the reciprocal problem: the "force" exerted by the body on the surrounding matter, and show that (from this perspective) the effect is due to the body's gravitomagnetic field. We compute it rigorously, showing the matching with its reciprocal, and clarifying common misconceptions in the literature regarding the action-reaction law in post-Newtonian gravity.

We consider motion in spherically symmetric but time-dependent backgrounds. This problem is of interest, for example, in the context of ultralight dark matter, where galactic haloes produce a time-dependent and periodic gravitational potential. We study the properties of motion of stars in such spacetimes, for different field strengths and frequency, and including dissipative effects. We show that orbital resonances may occur and that spectroscopic emission lines from stars in these geometries exhibit characteristic, periodic modulation patterns. In addition, we work out a fully relativistic and weak-field description of a special class of time-periodic geometries, that of scalar oscillatons. When applied to the galactic center, our results indicate that the motion of S2-like stars may carry distinguishable observational imprints of ultra-light dark matter.

Superradiant scattering is a radiation enhancement process that takes place in many contexts, and which has recently found exciting applications in astro and particle physics. In the framework of curved spacetime physics, it has been associated with the classical Penrose process for particles. Superradiance is usually also associated with bosonic fields around geometries with ergoregions and horizons. These notions are in clear tension however: the Penrose process occurs for horizonless geometries, and particles are composed of fermions. Here, we resolve the tension in its different aspects, by showing that (i) superradiance occurs for self-interacting fermions on flat spacetime; (ii) superradiance occurs also for horizonless geometries, where it leads to an ergoregion instability. Ultracompact, horizonless geometries will usually respond with echoes of growing amplitude, until rotational (or electrostatic) energy is extracted from the object; (iii) the Fourier-domain analysis leads to absence of superradiance when horizons are not present. We elucidate why this analysis fails to give meaningful results; (iv) finally, we show that superradiant, ergoregion instabilities have a particle analog of similar growth timescales and which can power the formation of a structure outside a compact, rotating star.

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.

The relaxation of a distorted black hole to its final state provides important tests of general relativity within the reach of current and upcoming gravitational wave facilities. In black hole perturbation theory, this phase consists of a simple linear superposition of exponentially damped sinusoids (the quasinormal modes) and of a power-law tail. How many quasinormal modes are necessary to describe waveforms with a prescribed precision? What error do we incur by only including quasinormal modes, and not tails? What other systematic effects are present in current state-of-the-art numerical waveforms? These issues, which are basic to testing fundamental physics with distorted black holes, have hardly been addressed in the literature. We use numerical relativity waveforms and accurate evolutions within black hole perturbation theory to provide some answers. We show that (i) a determination of the fundamental $l=m=2$ quasinormal mode to within $1\%$ or better requires the inclusion of at least the first overtone, and preferably of the first two or three overtones; (ii) a determination of the black hole mass and spin with precision better than $1\%$ requires the inclusion of at least two quasinormal modes for any given angular harmonic mode $(\ell,\,m)$. We also improve on previous estimates and fits for the ringdown energy radiated in the various multipoles. These results are important to quantify theoretical (as opposed to instrumental) limits in parameter estimation accuracy and tests of general relativity allowed by ringdown measurements with high signal-to-noise ratio gravitational wave detectors.

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?

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

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.

Following the selection of The Gravitational Universe by ESA, and the successful flight of LISA Pathfinder, the LISA Consortium now proposes a 4 year mission in response to ESA's call for missions for L3. The observatory will be based on three arms with six active laser links, between three identical spacecraft in a triangular formation separated by 2.5 million km. LISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using Gravitational Waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the infant Universe at TeV energy scales, has known sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales near the horizons of black holes, all the way to cosmological scales. The LISA mission will scan the entire sky as it follows behind the Earth in its orbit, obtaining both polarisations of the Gravitational Waves simultaneously, and will measure source parameters with astrophysically relevant sensitivity in a band from below $10^{-4}\,$Hz to above $10^{-1}\,$Hz.

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 estimate the potential of present and future interferometric gravitational-wave detectors to test the Kerr nature of black holes through "gravitational spectroscopy," i.e. the measurement of multiple quasinormal mode frequencies from the remnant of a black hole merger. Using population synthesis models of the formation and evolution of stellar-mass black hole binaries, we find that Voyager-class interferometers will be necessary to perform these tests. Gravitational spectroscopy in the local Universe may become routine with the Einstein Telescope, but a 40-km facility like Cosmic Explorer is necessary to go beyond $z\sim 3$. In contrast, eLISA-like detectors should carry out a few - or even hundreds - of these tests every year, depending on uncertainties in massive black hole formation models. Many space-based spectroscopical measurements will occur at high redshift, testing the strong gravity dynamics of Kerr black holes in domains where cosmological corrections to general relativity (if they occur in nature) must be significant.

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.

Black holes in General Relativity are very simple objects. This property, that goes under the name of "no-hair," has been refined in the last few decades and admits several versions. The simplicity of black holes makes them ideal testbeds of fundamental physics and of General Relativity itself. Here we discuss the no-hair property of black holes, how it can be measured in the electromagnetic or gravitational window, and what it can possibly tell us about our universe.

We review recent progress in the application of numerical relativity techniques to astrophysics and high-energy physics. We focus on some developments that took place within the "Numerical Relativity and High Energy Physics" network, a Marie Curie IRSES action that we coordinated, namely: spin evolution in black hole binaries, high-energy black hole collisions, compact object solutions in scalar-tensor gravity, superradiant instabilities and hairy black hole solutions in Einstein's gravity coupled to fundamental fields, and the possibility to gain insight into these phenomena using analog gravity models.

We provide a detailed analysis of how bosonic dark matter "condensates" interact with compact stars, extending significantly the results of a recent Letter. We focus on bosonic fields with mass $m_B$, such as axions, axion-like candidates and hidden photons. Self-gravitating bosonic fields generically form "breathing" configurations, where both the spacetime geometry and the field oscillate, and can interact and cluster at the center of stars. We construct stellar configurations formed by a perfect fluid and a bosonic condensate, and which may describe the late stages of dark-matter accretion onto stars, in dark matter-rich environments. These composite stars oscillate at a frequency which is a multiple of $f=2.5\times 10^{14}\,\left(m_{B}c^2/eV\right)\,{\rm Hz}$. Using perturbative analysis and Numerical Relativity techniques, we show that these stars are generically stable, and we provide criteria for instability. Our results also indicate that the growth of the dark matter core is halted close to the Chandrasekhar limit. We thus dispel a myth concerning dark matter accretion by stars: dark mater accretion does not necessarily lead to the destruction of the star, nor to collapse to a black hole. Finally, we argue that stars with long-lived bosonic cores may also develop in other theories with effective mass couplings, such as (massless) scalar-tensor theories.

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.

The ``cosmic censorship conjecture'' asserts that all singularities arising from gravitational collapse are hidden within black holes. We investigate this conjecture in a setup of interest for tests of General Relativity: black hole solutions which are parametrically small deviations away from the Kerr solution. These solutions have an upper bound on rotation, beyond which a naked singularity is visible to outside observers. We study whether these (generic) spacetimes can be spun-up past extremality with point particles or accretion disks. Our results show that cosmic censorship is preserved for generic parameterizations. We also present examples of special geometries which can be spun-up past extremality.

The demand to obtain answers to highly complex problems within strong-field gravity has been met with significant progress in the numerical solution of Einstein's equations - along with some spectacular results - in various setups. We review techniques for solving Einstein's equations in generic spacetimes, focusing on fully nonlinear evolutions but also on how to benchmark those results with perturbative approaches. The results address problems in high-energy physics, holography, mathematical physics, fundamental physics, astrophysics and cosmology.

Searches for dark matter imprints are one of the most active areas of current research. We focus here on light fields with mass $m_B$, such as axions and axion-like candidates. Using perturbative techniques and full-blown nonlinear Numerical Relativity methods, we show that (i) dark matter can pile up in the center of stars, leading to configurations and geometries oscillating with frequency which is a multiple of f=$2.5 10^{14}$ $m_B c^2$/eV Hz. These configurations are stable throughout most of the parameter space, and arise out of credible mechanisms for dark-matter capture. Stars with bosonic cores may also develop in other theories with effective mass couplings, such as (massless) scalar-tensor theories. We also show that (ii) collapse of the host star to a black hole is avoided by efficient gravitational cooling mechanisms.

We establish that massive complex Abelian vector fields (mass $\mu$) can form gravitating solitons, when minimally coupled to Einstein's gravity. Such Proca stars (PSs) have a stationary, everywhere regular and asymptotically flat geometry. The Proca field, however, possesses a harmonic time dependence (frequency $w$), realizing Wheeler's concept of geons for an Abelian spin 1 field. We obtain PSs with both a spherically symmetric (static) and an axially symmetric (stationary) line element. The latter form a countable number of families labelled by an integer $m\in \mathbb{Z}^+$. PSs, like (scalar) boson stars, carry a conserved Noether charge, and are akin to the latter in many ways. In particular, both types of stars exist for a limited range of frequencies and there is a maximal ADM mass, $M_{max}$, attained for an intermediate frequency. For spherically symmetric PSs (rotating PSs with $m=1,2,3$), $M_{max}\simeq 1.058 M_{Pl}^2/\mu$ ($M_{max}\simeq 1.568,\, 2.337, \, 3.247 \, M_{Pl}^2/\mu$), slightly larger values than those for (mini-)boson stars. We establish perturbative stability for a subset of solutions in the spherical case and anticipate a similar conclusion for fundamental modes in the rotating case. The discovery of PSs opens many avenues of research, reconsidering five decades of work on (scalar) boson stars, in particular as possible dark matter candidates.

It has long been known that dissipation is a crucial ingredient in the superradiant amplification of wavepackets off rotating objects. We show that, once appropriate dissipation mechanisms are included, stars are also prone to superradiance and superradiant instabilities. In particular, ultra-light dark matter with small interaction cross section with the star material or self-annihilation can trigger a superradiant instability. On long timescales, the instability strips the star of most of its angular momentum. Whether or not new stationary configurations surrounded by scalar condensates exist, remains to be seen.

Soon after the discovery of the Kerr metric, Penrose realized that superradiance can be exploited to extract energy from black holes. The original idea (involving the breakup of a single particle) yields only modest energy gains. A variant of the Penrose process consists of particle collisions in the ergoregion. The collisional Penrose process has been explored recently in the context of dark matter searches, with the conclusion that the ratio $\eta$ between the energy of post-collision particles detected at infinity and the energy of the colliding particles should be modest ($\eta \lesssim 1.5$). Schnittman has shown that these studies underestimated the maximum efficiency by about one order of magnitude (i.e., $\eta \lesssim 15$). In this work we show that particle collisions in the vicinity of rapidly rotating black holes can produce high-energy ejecta and result in high efficiencies under much more generic conditions. The astrophysical likelihood of these events deserves further scrutiny, but our study hints at the tantalizing possibility that the collisional Penrose process may power gamma rays and ultra-high-energy cosmic rays.

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.

Scalar forces "screened" by the Vainshtein mechanism may hold the key to understanding the cosmological expansion of our universe, while predicting new and exciting features in the interaction between massive bodies. Here we explore the dynamics of the Vainshtein screening mechanism, focusing on the decoupling limit of the DGP braneworld scenario and dRGT massive gravity. We show that there is a vast set of initial conditions whose evolution is well defined and which are driven to the static, screening solutions of these theories. Screening solutions are stable and behave coherently under small fluctuations: they oscillate and eventually settle to an equilibrium configuration, the timescale for the oscillations and damping being dictated by the Vainshtein radius of the screening solutions. At very late times, a power-law decay ensues, in agreement with known analytical results. However, we also conjecture that physically interesting processes such as gravitational collapse of compact stars may not possess a well-posed initial value problem. Finally, we construct solutions with nontrivial multipolar structure describing the screening field of deformed, asymmetric bodies and show that higher multipoles are screened more efficiently than the monopole component.

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.

[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.

Scalar fields pervade theoretical physics and are a fundamental ingredient to solve the dark matter problem, to realize the Peccei-Quinn mechanism in QCD or the string-axiverse scenario. They are also a useful proxy for more complex matter interactions, such as accretion disks or matter in extreme conditions. Here, we study the collision between scalar "clouds" and rotating black holes. For the first time we are able to compare analytic estimates and strong field, nonlinear numerical calculations for this problem. As the black hole pierces through the cloud it accretes according to the Bondi-Hoyle prediction, but is deflected through a purely kinematic gravitational "anti-Magnus" effect, which we predict to be present also during the interaction of black holes with accretion disks. After the interaction is over, we find large recoil velocities in the transverse direction. The end-state of the process belongs to the vacuum Kerr family if the scalar is massless, but can be a hairy black hole when the fundamental scalar is massive.

Fundamental fields are a natural outcome in cosmology and particle physics and might therefore serve as a proxy for more complex interactions. The equivalence principle implies that all forms of matter gravitate, and one therefore expects relevant, universal imprints of new physics in strong field gravity, such as that encountered close to black holes. Fundamental fields in the vicinities of supermassive black holes give rise to extremely long-lived, or even unstable, configurations which slowly extract angular momentum from the black hole or simply evolve non-linearly over long timescales, with important implications for particle physics and gravitational-wave physics. Here, we perform a fully non-linear study of scalar-field condensates around rotating black holes. We provide novel ways to specify initial data for the Einstein-Klein-Gordon system, with potential applications in a variety of scenarios. Our numerical results confirm the existence of long-lived bar-modes which act as lighthouses for gravitational wave emission: the scalar field condenses outside the black hole geometry and acts as a constant frequency gravitational-wave source for very long timescales. This effect could turn out to be a potential signature of beyond standard model physics and also a promising source of gravitational waves for future gravitational wave detectors.

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.

We summarize the parallel session B4: 'Analytic approximations, perturbation theory effective field theory methods and their applications' and the joint session B2/B4: 'Approximate solutions to Einstein equations: Methods and Applications', of the GR20 & Amaldi10 conference in Warsaw, July 2013. The contributed talks reported significant advances on various areas of research in gravity.

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.

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.

Black holes are the elementary particles of gravity, the final state of sufficiently massive stars and of energetic collisions. With a forty-year long history, black hole physics is a fully-blossomed field which promises to embrace several branches of theoretical physics. Here I review the main developments in highly dynamical black holes with an emphasis on high energy black hole collisions and probes of particle physics via superradiance. This write-up, rather than being a collection of well known results, is intended to highlight open issues and the most intriguing results.

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.

Scalar-tensor theories are a compelling alternative to general relativity and one of the most accepted extensions of Einstein's theory. Black holes in these theories have no hair, but could grow "wigs" supported by time-dependent boundary conditions or spatial gradients. Time-dependent or spatially varying fields lead in general to nontrivial black hole dynamics, with potentially interesting experimental consequences. We carry out a numerical investigation of the dynamics of single and binary black holes in the presence of scalar fields. In particular we study gravitational and scalar radiation from black-hole binaries in a constant scalar-field gradient, and we compare our numerical findings to analytical models. In the single black hole case we find that, after a short transient, the scalar field relaxes to static configurations, in agreement with perturbative calculations. Furthermore we predict analytically (and verify numerically) that accelerated black holes in a scalar-field gradient emit scalar radiation. For a quasicircular black-hole binary, our analytical and numerical calculations show that the dominant component of the scalar radiation is emitted at twice the binary's orbital frequency.

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.

The analytical understanding of quasinormal mode ringing requires an accurate knowledge of the Green's function describing the response of the black hole to external perturbations. We carry out a comprehensive study of quasinormal mode excitation for Kerr black holes. Relying on the formalism developed by Mano, Suzuki and Takasugi, we improve and extend previous calculations of the quasinormal mode residues in the complex frequency plane ("excitation factors" B_q). Using these results we compute the "excitation coefficients" C_q (essentially the mode amplitudes) in the special case where the source of the perturbations is a particle falling into the black hole along the symmetry axis. We compare this calculation with numerical integrations of the perturbation equations, and we show quantitatively how the addition of higher overtones improves the agreement with the numerical waveforms. Our results should find applications in models of the ringdown stage and in the construction of semianalytical template banks for gravitational-wave detectors, especially for binaries with large mass ratios and/or fast-spinning black holes.

Light bosonic degrees of freedom have become a serious candidate for dark matter. The evolution of these fields around curved spacetimes is poorly understood but is expected to display interesting effects. In particular, the interaction of light bosonic fields with supermassive black holes, key players in most galaxies, could provide colourful examples of superradiance and nonlinear bosenova-like collapse. In turn, the observation of spinning black holes is expected to impose stringent bounds on the mass of putative massive bosonic fields in our universe. Our purpose here is to present a comprehensive study of the evolution of linearized massive scalar and vector fields in the vicinities of rotating black holes. For a certain boson field mass range, the field can become trapped in a potential barrier outside the horizon and transition to a bound state. Because there are a number of such quasi-bound states, the generic outcome is an amplitude modulated sinusoidal, or beating, signal. We believe that the appearance of such beatings has gone unnoticed in the past, and in fact mistaken for exponential growth. The amplitude modulation of the signal depends strongly on the relative excitation of the overtones, which in turn is strongly tied to the bound-state geography. For the first time we explore massive vector fields in generic BH background which are hard, if not impossible, to separate in the Kerr background. Our results show that spinning BHs are generically strongly unstable against massive vector fields.

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.

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.

We perform fully non-linear numerical simulations of charged-black-hole collisions, described by the Einstein-Maxwell equations, and contrast the results against analytic expectations. We focus on head-on collisions of non-spinning black holes, starting from rest and with the same charge to mass ratio, Q/M. The addition of charge to black holes introduces a new interesting channel of radiation and dynamics, most of which seem to be captured by Newtonian dynamics and flat-space intuition. The waveforms can be qualitatively described in terms of three stages; (i) an infall phase prior to the formation of a common apparent horizon; (ii) a nonlinear merger phase which corresponds to a peak in gravitational and electromagnetic energy; (iii) the ringdown marked by an oscillatory pattern with exponentially decaying amplitude and characteristic frequencies that are in good agreement with perturbative predictions. We observe that the amount of gravitational-wave energy generated throughout the collision decreases by about three orders of magnitude as the charge-to-mass ratio Q/M is increased from 0 to 0.98. We interpret this decrease as a consequence of the smaller accelerations present for larger values of the charge. In contrast, the ratio of energy carried by electromagnetic to gravitational radiation increases, reaching about 22% for the maximum Q/M ratio explored, which is in good agreement with analytic predictions.

We review recent progress in numerical relativity simulations of black-hole (BH) spacetimes. Following a brief summary of the methods employed in the modeling, we summarize the key results in three major areas of BH physics: (i) BHs as sources of gravitational waves (GWs), (ii) astrophysical systems involving BHs, and (iii) BHs in high-energy physics. We conclude with a list of the most urgent tasks for numerical relativity in these three areas.

Nonlinear dynamics in cosmological backgrounds has the potential to teach us immensely about our universe, and also to serve as prototype for nonlinear processes in generic curved spacetimes. Here we report on dynamical evolutions of black holes in asymptotically de Sitter spacetimes. We focus on the head-on collision of equal mass binaries and for the first time compare analytical and perturbative methods with full blown nonlinear simulations. Our results include an accurate determination of the merger/scatter transition (consequence of an expanding background) for small mass binaries and a test of the Cosmic Censorship conjecture, for large mass binaries. We observe that, even starting from small separations, black holes in large mass binaries eventually lose causal contact, in agreement with the conjecture.

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.

A classical thought-experiment to destroy black holes was envisaged by Wald in 1974: it consists of throwing particles with large angular momentum into an extremal black hole, checking whether their capture can over-spin the black hole past the extremal limit and create a naked singularity. Wald showed that in the test-particle limit, particles that would be otherwise capable of producing naked singularities are simply scattered. Recently Jacobson and Sotiriou showed that if one considers instead a black hole that is almost, but not exactly extremal, then in the absence of backreaction effects particle capture could indeed over-spin the spacetime above the Kerr limit. Here we analyze back-reaction effects and show that for some of the trajectories giving rise to naked singularities, radiative effects can be neglected. However, for these orbits the conservative self-force is important, and seems to have the right sign to prevent the formation of naked singularities.

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.

Numerical relativity has seen incredible progress in the last years, and is being applied with success to a variety of physical phenomena, from gravitational-wave research and relativistic astrophysics to cosmology and high-energy physics. Here we probe the limits of current numerical setups, by studying collisions of unequal mass, non-rotating black holes of mass-ratios up to 1:100 and making contact with a classical calculation in General Relativity: the infall of a point-like particle into a massive black hole. Our results agree well with the predictions coming from linearized calculations of the infall of point-like particles into non-rotating black holes. In particular, in the limit that one hole is much smaller than the other, and the infall starts from an infinite initial separation, we recover the point-particle limit. Thus, numerical relativity is able to bridge the gap between fully non-linear dynamics and linearized approximations, which may have important applications. Finally, we also comment on the "spurious" radiation content in the initial data and the linearized predictions.

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.

Jacobson and Sotiriou showed that rotating black holes could be spun-up past the extremal limit by the capture of non-spinning test bodies, if one neglects radiative and self-force effects. This would represent a violation of the Cosmic Censorship Conjecture in four-dimensional, asymptotically flat spacetimes. We show that for some of the trajectories giving rise to naked singularities, radiative effects can be neglected. However, for these orbits the conservative self-force is important, and seems to have the right sign to prevent the formation of naked singularities.

Black objects in higher dimensional space-times have a remarkably richer structure than their four dimensional counterparts. They appear in a variety of configurations (e.g. black holes, black branes, black rings, black Saturns), and display complex stability phase diagrams. They might also play a key role in high energy physics: for energies above the fundamental Planck scale, gravity is the dominant interaction which, together with the hoop-conjecture, implies that the trans-Planckian scattering of point particles should be well described by black hole scattering. Higher dimensional scenarios with a fundamental Planck scale of the order of TeV predict, therefore, black hole production at the LHC, as well as in future colliders with yet higher energies. In this setting, accurate predictions for the production cross-section and energy loss (through gravitational radiation) in the formation of black holes in parton-parton collisions is crucial for accurate phenomenological modelling in Monte Carlo event generators. In this paper, we use the formalism and numerical code reported in arXiv:1001.2302 to study the head-on collision of two black holes. For this purpose we provide a detailed treatment of gravitational wave extraction in generic D-dimensional space-times, which uses the Kodama-Ishibashi formalism. For the first time, we present the results of numerical simulations of the head-on collision in five space-time dimensions, together with the relevant physical quantities. We show that the total radiated energy, when two black holes collide from rest at infinity, is approximately (0.089\pm 0.006)% of the centre of mass energy, slightly larger than the 0.055% obtained in the four dimensional case, and that the ringdown signal at late time is in very good agreement with perturbative calculations.

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.

The numerical evolution of Einstein's field equations in a generic background has the potential to answer a variety of important questions in physics: from applications to the gauge-gravity duality, to modelling black hole production in TeV gravity scenarios, analysis of the stability of exact solutions and tests of Cosmic Censorship. In order to investigate these questions, we extend numerical relativity to more general space-times than those investigated hitherto, by developing a framework to study the numerical evolution of D dimensional vacuum space-times with an SO(D-2) isometry group for D\ge 5, or SO(D-3) for D\ge 6. Performing a dimensional reduction on a (D-4)-sphere, the D dimensional vacuum Einstein equations are rewritten as a 3+1 dimensional system with source terms, and presented in the Baumgarte, Shapiro, Shibata and Nakamura (BSSN) formulation. This allows the use of existing 3+1 dimensional numerical codes with small adaptations. Brill-Lindquist initial data are constructed in D dimensions and a procedure to match them to our 3+1 dimensional evolution equations is given. We have implemented our framework by adapting the LEAN code and perform a variety of simulations of non-spinning black hole space-times. Specifically, we present a modified moving puncture gauge which facilitates long term stable simulations in D=5. We further demonstrate the internal consistency of the code by studying convergence and comparing numerical versus analytic results in the case of geodesic slicing for D=5,6.

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 has been suggested that rotating black holes could serve as particle colliders with arbitrarily high center-of-mass energy. Astrophysical limitations on the maximal spin, back-reaction effects and sensitivity to the initial conditions impose severe limits on the likelihood of such collisions.

Dynamical Chern-Simons (DCS) modified gravity is an attractive, yet relatively unexplored, candidate to an alternative theory of gravity. The DCS correction couples a dynamical scalar field to the gravitational field. In this framework, we analyze the perturbation formalism and stability properties of spherically symmetric black holes. Assuming that no background scalar field is present, gravitational perturbations with polar and axial parities decouple. We find no effect of the Chern-Simons coupling on the polar sector, while axial perturbations couple to the Chern-Simons scalar field. The axial sector can develop strong instabilities if the coupling parameter beta, associated to the dynamical coupling of the scalar field, is small enough; this yields a constraint on beta which is much stronger than the constraints previously known in the literature.

Quasinormal modes are eigenmodes of dissipative systems. Perturbations of classical gravitational backgrounds involving black holes or branes naturally lead to quasinormal modes. The analysis and classification of the quasinormal spectra requires solving non-Hermitian eigenvalue problems for the associated linear differential equations. Within the recently developed gauge-gravity duality, these modes serve as an important tool for determining the near-equilibrium properties of strongly coupled quantum field theories, in particular their transport coefficients, such as viscosity, conductivity and diffusion constants. In astrophysics, the detection of quasinormal modes in gravitational wave experiments would allow precise measurements of the mass and spin of black holes as well as new tests of general relativity. This review is meant as an introduction to the subject, with a focus on the recent developments in the field.

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.

Geodesic motion determines important features of spacetimes. Null unstable geodesics are closely related to the appearance of compact objects to external observers and have been associated with the characteristic modes of black holes. By computing the Lyapunov exponent, which is the inverse of the instability timescale associated with this geodesic motion, we show that, in the eikonal limit, quasinormal modes of black holes in any dimensions are determined by the parameters of the circular null geodesics. This result is independent of the field equations and only assumes a stationary, spherically symmetric and asymptotically flat line element, but it does not seem to be easily extendable to anti-de Sitter spacetimes. We further show that (i) in spacetime dimensions greater than four, equatorial circular timelike geodesics in a Myers-Perry black hole background are unstable, and (ii) the instability timescale of equatorial null geodesics in Myers-Perry spacetimes has a local minimum for spacetimes of dimension d > 5.

We study the transition from inspiral to plunge in general relativity by computing gravitational waveforms of non-spinning, equal-mass black-hole binaries. We consider three sequences of simulations, starting with a quasi-circular inspiral completing 1.5, 2.3 and 9.6 orbits, respectively, prior to coalescence of the holes. For each sequence, the binding energy of the system is kept constant and the orbital angular momentum is progressively reduced, producing orbits of increasing eccentricity and eventually a head-on collision. We analyze in detail the radiation of energy and angular momentum in gravitational waves, the contribution of different multipolar components and the final spin of the remnant. We find that the motion transitions from inspiral to plunge when the orbital angular momentum L=L_crit is about 0.8M^2. For L<L_crit the radiated energy drops very rapidly. Orbits with L of about L_crit produce our largest dimensionless Kerr parameter for the remnant, j=J/M^2=0.724. Generalizing a model recently proposed by Buonanno, Kidder and Lehner to eccentric binaries, we conjecture that (1) j=0.724 is the maximal Kerr parameter that can be obtained by any merger of non-spinning holes, and (2) no binary merger (even if the binary members are extremal Kerr black holes with spins aligned to the orbital angular momentum, and the inspiral is highly eccentric) can violate the cosmic censorship conjecture.

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.

We study the transformation of multipolar decompositions of gravitational radiation under rotations and boosts. Rotations to the remnant black hole's frame simplify the waveforms from the merger of generic spinning black hole binaries. Boosts may be important to get an accurate gravitational-wave phasing, especially for configurations leading to large recoil velocities of the remnant. As a test of our formalism we revisit the classic problem of point particles falling into a Schwarzschild black hole. Then we highlight by specific examples the importance of choosing the right frame in numerical simulations of unequal-mass, spinning binary black-hole mergers.

Some years ago, a new powerful technique, known as the Classical Effective Field Theory, was proposed to describe classical phenomena in gravitational systems. Here we show how this approach can be useful to investigate theoretically important issues, such as gravitational radiation in any spacetime dimension. In particular, we derive for the first time the Einstein-Infeld-Hoffman Lagrangian and we compute Einstein's quadrupole formula for any number of flat spacetime dimensions.

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.

Using recent results from numerical relativity simulations of non-spinning binary black hole mergers we revisit the problem of detecting ringdown waveforms and of estimating the source parameters, considering both LISA and Earth-based interferometers. We find that Advanced LIGO and EGO could detect intermediate-mass black holes of mass up to about 1000 solar masses out to a luminosity distance of a few Gpc. For typical multipolar energy distributions, we show that the single-mode ringdown templates presently used for ringdown searches in the LIGO data stream can produce a significant event loss (> 10% for all detectors in a large interval of black hole masses) and very large parameter estimation errors on the black hole's mass and spin. We estimate that more than 10^6 templates would be needed for a single-stage multi-mode search. Therefore, we recommend a "two stage" search to save on computational costs: single-mode templates can be used for detection, but multi-mode templates or Prony methods should be used to estimate parameters once a detection has been made. We update estimates of the critical signal-to-noise ratio required to test the hypothesis that two or more modes are present in the signal and to resolve their frequencies, showing that second-generation Earth-based detectors and LISA have the potential to perform no-hair tests.

We study the inspiral, merger and ringdown of unequal mass black hole binaries by analyzing a catalogue of numerical simulations for seven different values of the mass ratio (from q=M2/M1=1 to q=4). We compare numerical and Post-Newtonian results by projecting the waveforms onto spin-weighted spherical harmonics, characterized by angular indices (l,m). We find that the Post-Newtonian equations predict remarkably well the relation between the wave amplitude and the orbital frequency for each (l,m), and that the convergence of the Post-Newtonian series to the numerical results is non-monotonic. To leading order the total energy emitted in the merger phase scales like eta^2 and the spin of the final black hole scales like eta, where eta=q/(1+q)^2 is the symmetric mass ratio. We study the multipolar distribution of the radiation, finding that odd-l multipoles are suppressed in the equal mass limit. Higher multipoles carry a larger fraction of the total energy as q increases. We introduce and compare three different definitions for the ringdown starting time. Applying linear estimation methods (the so-called Prony methods) to the ringdown phase, we find resolution-dependent time variations in the fitted parameters of the final black hole. By cross-correlating information from different multipoles we show that ringdown fits can be used to obtain precise estimates of the mass and spin of the final black hole, which are in remarkable agreement with energy and angular momentum balance calculations.

The ringdown phase following a binary black hole merger is usually assumed to be well described by a linear superposition of complex exponentials (quasinormal modes). In the strong-field conditions typical of a binary black hole merger, non-linear effects may produce mode coupling. Artificial mode coupling can also be induced by the black hole's rotation, if the radiation field is expanded in terms of spin-weighted spherical (rather than spheroidal) harmonics. Observing deviations from linear black hole perturbation theory requires optimal fitting techniques to extract ringdown parameters from numerical waveforms, which are inevitably affected by errors. So far, non-linear least-squares fitting methods have been used as the standard workhorse to extract frequencies from ringdown waveforms. These methods are known not to be optimal for estimating parameters of complex exponentials. Furthermore, different fitting methods have different performance in the presence of noise. The main purpose of this paper is to introduce the gravitational wave community to modern variations of a linear parameter estimation technique first devised in 1795 by Prony: the Kumaresan-Tufts and matrix pencil methods. Using "test" damped sinusoidal signals in Gaussian white noise we illustrate the advantages of these methods, showing that they have variance and bias at least comparable to standard non-linear least-squares techniques. Then we compare the performance of different methods on unequal-mass binary black hole merger waveforms. The methods we discuss should be useful both theoretically (to monitor errors and search for non-linearities in numerical relativity simulations) and experimentally (for parameter estimation from ringdown signals after a gravitational wave detection).

Distorted black holes radiate gravitational waves. In the so-called ringdown phase radiation is emitted in a discrete set of complex quasinormal frequencies, whose values depend only on the black hole's mass and angular momentum. Ringdown radiation could be detectable with large signal-to-noise ratio by the Laser Interferometer Space Antenna LISA. If more than one mode is detected, tests of the black hole nature of the source become possible. The detectability of different modes depends on their relative excitation, which in turn depends on the cause of the perturbation (i.e. on the initial data). A ``universal'', initial data-independent measure of the relative mode excitation is encoded in the poles of the Green's function that propagates small perturbations of the geometry (``excitation factors''). We compute for the first time the excitation factors for general-spin perturbations of Kerr black holes. We find that for corotating modes with $l=m$ the excitation factors tend to zero in the extremal limit, and that the contribution of the overtones should be more significant when the black hole is fast rotating. We also present the first analytical calculation of the large-damping asymptotics of the excitation factors for static black holes, including the Schwarzschild and Reissner-Nordstrom metrics. This is an important step to determine the convergence properties of the quasinormal mode expansion.

We study equilibrium shapes, stability and possible bifurcation diagrams of fluids in higher dimensions, held together by either surface tension or self-gravity. We consider the equilibrium shape and stability problem of self-gravitating spheroids, establishing the formalism to generalize the MacLaurin sequence to higher dimensions. We show that such simple models, of interest on their own, also provide accurate descriptions of their general relativistic relatives with event horizons. The examples worked out here hint at some model-independent dynamics, and thus at some universality: smooth objects seem always to be well described by both ``replicas'' (either self-gravity or surface tension). As an example, we exhibit an instability afflicting self-gravitating (Newtonian) fluid cylinders. This instability is the exact analogue, within Newtonian gravity, of the Gregory-Laflamme instability in general relativity. Another example considered is a self-gravitating Newtonian torus made of a homogeneous incompressible fluid. We recover the features of the black ring in general relativity.

We show that if the visible universe is a membrane embedded in a higher-dimensional space, particles in uniform motion radiate gravitational waves because of spacetime lumpiness. This phenomenon is analogous to the electromagnetic diffraction radiation of a charge moving near to a metallic grating. In the gravitational case, the role of the metallic grating is played by the inhomogeneities of the extra-dimensional space, such as a hidden brane. We derive a general formula for gravitational diffraction radiation and apply it to a higher-dimensional scenario with flat compact extra dimensions. Gravitational diffraction radiation may carry away a significant portion of the particle's initial energy. This allows to set stringent limits on the scale of brane perturbations. Physical effects of gravitational diffraction radiation are briefly discussed.

We show that algebraically special modes lead to the instability of naked singularity spacetimes with negative mass. Four-dimensional negative-mass Schwarzschild and Schwarzschild-de Sitter spacetimes are unstable. Stability of the Schwarzschild-anti-de Sitter spacetime depends on boundary conditions. We briefly discuss the generalization of these results to charged and rotating singularities.

The evidence for supermassive Kerr black holes in galactic centers is strong and growing, but only the detection of gravitational waves will convincingly rule out other possibilities to explain the observations. The Kerr spacetime is completely specified by the first two multipole moments: mass and angular momentum. This is usually referred to as the ``no-hair theorem'', but it is really a ``two-hair'' theorem. If general relativity is the correct theory of gravity, the most plausible alternative to a supermassive Kerr black hole is a rotating boson star. Numerical calculations indicate that the spacetime of rotating boson stars is determined by the first three multipole moments (``three-hair theorem''). LISA could accurately measure the oscillation frequencies of these supermassive objects. We propose to use these measurements to ``count their hair'', unambiguously determining their nature and properties.

In models with extra dimensions, a black hole evaporates both in the bulk and on the visible brane, where standard model fields live. The exact emissivities of each particle species are needed to determine how the black hole decay proceeds. We compute and discuss the absorption cross-sections, the relative emissivities and the total power output of all known fields in the evaporation phase. Graviton emissivity is highly enhanced as the spacetime dimensionality increases. Therefore, a black hole loses a significant fraction of its mass in the bulk. This result has important consequences for the phenomenology of black holes in models with extra dimensions and black hole detection in particle colliders.

Newly formed black holes are expected to emit characteristic radiation in the form of quasi-normal modes, called ringdown waves, with discrete frequencies. LISA should be able to detect the ringdown waves emitted by oscillating supermassive black holes throughout the observable Universe. We develop a multi-mode formalism, applicable to any interferometric detectors, for detecting ringdown signals, for estimating black hole parameters from those signals, and for testing the no-hair theorem of general relativity. Focusing on LISA, we use current models of its sensitivity to compute the expected signal-to-noise ratio for ringdown events, the relative parameter estimation accuracy, and the resolvability of different modes. We also discuss the extent to which uncertainties on physical parameters, such as the black hole spin and the energy emitted in each mode, will affect our ability to do black hole spectroscopy.

We provide some considerations on the excitation of black hole quasinormal modes (QNMs) in different physical scenarios. Considering a simple model in which a stream of particles accretes onto a black hole, we show that resonant QNM excitation by hyperaccretion requires a significant amount of fine-tuning, and is quite unlikely to occur in nature. Then we summarize and discuss present estimates of black hole QNM excitation from gravitational collapse, distorted black holes and head-on black hole collisions. We emphasize the areas that, in our opinion, are in urgent need of further investigation from the point of view of gravitational wave source modeling.

The evolution of small perturbations around rotating black branes and strings, which are low energy solutions of string theory, are investigated. For simplicity, we concentrate on the Kerr solution times transverse flat extra dimensions, possibly compactified, but one can also treat other branes composed of any rotating black hole and extra transverse dimensions, as well as analogue black hole models and rotating bodies in fluid mechanics systems. It is shown that such a rotating black brane is unstable against any massless (scalar, vectorial, tensorial or other) field perturbation for a wide range of wavelengths and frequencies in the transverse dimensions. Since it holds for any massless field it can be considered, in this sense, a stronger instability than the one studied by Gregory and Laflamme. Accordingly, it has also a totally different physical origin. The perturbations can be stabilized if the extra dimensions are compactified to a length smaller than the minimum wavelength for which the instability settles in, resembling in this connection the Gregory-Laflamme case. Likewise, this instability will have no effect for astrophysical black holes. However, in the large extra dimensions scenario, where TeV scale black holes can be produced, this instability should be important. It seems plausible that the endpoint of this instability is a static, or very slowly rotating, black brane and some outgoing radiation at infinity.

Higher-dimensional scenarios allow for the formation of mini-black holes from TeV-scale particle collisions. The purpose of this paper is to review and compare different methods for the estimate of the total gravitational energy emitted in this process. To date, black hole formation has mainly been studied using an apparent horizon search technique. This approach yields only an upper bound on the gravitational energy emitted during black hole formation. Alternative calculations based on instantaneous collisions of point particles and black hole perturbation theory suggest that the emitted gravitational energy may be smaller. New and more refined methods may be necessary to accurately describe black hole formation in high-energy particle collisions.

We discuss some general aspects of acoustic black holes. We begin by describing the associated formalism with which acoustic black holes are established, then we show how to model arbitrary geometries by using a de Laval nozzle. It is argued that even though the Hawking temperature of these black holes is too low to be detected, acoustic black holes have interesting classical properties, some of which are outlined here, that should be explored.

We investigate the stability of higher dimensional rotating black holes against scalar perturbations. In particular, we make a thorough numerical and analytical analysis of six-dimensional black holes, not only in the low rotation regime but in the high rotation regime as well. Our results suggest that higher dimensional Kerr black holes are stable against scalar perturbations, even in the ultra-spinning regime.

Many properties of black holes can be studied using acoustic analogues in the laboratory through the propagation of sound waves. We investigate in detail sound wave propagation in a rotating acoustic (2+1)-dimensional black hole, which corresponds to the ``draining bathtub'' fluid flow. We compute the quasinormal mode frequencies of this system and discuss late-time power-law tails. Due to the presence of an ergoregion, waves in a rotating acoustic black hole can be superradiantly amplified. We compute superradiant reflection coefficients and instability timescales for the acoustic black hole bomb, the equivalent of the Press-Teukolsky black hole bomb. Finally we discuss quasinormal modes and late-time tails in a non-rotating canonical acoustic black hole, corresponding to an incompressible, spherically symmetric (3+1)-dimensional fluid flow.

I discuss the range of validity of Detweiler's formula for the resonant frequencies of rapidly rotating Kerr black holes. While his formula is correct for extremal black holes, it has also been commonly accepted that it describes very well the resonant frequencies of near extremal black holes, and that therefore there is a large number of modes clustering on the real axis as the black hole becomes extremal. I will show that this last statement is not only incorrect, but that it also does not follow from Detweiler's formula, provided it is handled with due care. It turns out that only the first n << -log{(r_+-r_-)/r_+} modes are well described by that formula, which translates, for any astrophysical black hole, into one or two modes only. All existing numerical data gives further support to this claim. I also discuss some implications of this result for recent investigations on the late-time dynamics of rapidly rotating black holes.

The study of the quasinormal modes (QNMs) of the 2+1 dimensional rotating draining bathtub acoustic black hole, the closest analogue found so far to the Kerr black hole, is performed. Both the real and imaginary parts of the quasinormal (QN) frequencies as a function of the rotation parameter B are found through a full non-linear numerical analysis. Since there is no change in sign in the imaginary part of the frequency as B is increased we conclude that the 2+1 dimensional rotating draining bathtub acoustic black hole is stable against small perturbations.

The exact computation of asymptotic quasinormal frequencies is a technical problem which involves the analytic continuation of a Schrodinger-like equation to the complex plane and then performing a method of monodromy matching at the several poles in the plane. While this method was successfully used in asymptotically flat spacetime, as applied to both the Schwarzschild and Reissner-Nordstrom solutions, its extension to non-asymptotically flat spacetimes has not been achieved yet. In this work it is shown how to extend the method to this case, with the explicit analysis of Schwarzschild de Sitter and large Schwarzschild Anti-de Sitter black holes, both in four dimensions. We obtain, for the first time, analytic expressions for the asymptotic quasinormal frequencies of these black hole spacetimes, and our results match previous numerical calculations with great accuracy. We also list some results concerning the general classification of asymptotic quasinormal frequencies in d-dimensional spacetimes.

A wave impinging on a Kerr black hole can be amplified as it scatters off the hole if certain conditions are satisfied giving rise to superradiant scattering. By placing a mirror around the black hole one can make the system unstable. This is the black hole bomb of Press and Teukolsky. We investigate in detail this process and compute the growing timescales and oscillation frequencies as a function of the mirror's location. It is found that in order for the system black hole plus mirror to become unstable there is a minimum distance at which the mirror must be located. We also give an explicit example showing that such a bomb can be built. In addition, our arguments enable us to justify why large Kerr-AdS black holes are stable and small Kerr-AdS black holes should be unstable.

Superradiance in black hole spacetimes can trigger instabilities. Here we show that, due to superradiance, small Kerr-anti-de Sitter black holes are unstable. Our demonstration uses a matching procedure, in a long wavelength approximation.

Black holes play a fundamental role in modern physics. They have characteristic oscillation modes, called quasinormal modes. Past studies have shown that these modes are important to our understanding of the dynamics of astrophysical black holes. Recent studies indicate that they are important as a link between gravitation and quantum mechanics. Thus, the investigation of these modes is a timeliness topic. Quasinormal modes dominate almost every process involving black holes, in particular gravitational wave emission during, for example, the collision between two black holes. It may be possible to create black holes at future accelerators, according to recent theories advocating the existence of extra dimensions in our universe. It is therefore important to study in depth the gravitational radiation emitted in high energy collision between two black holes in several dimensions, and also to make a theoretical study of gravitational waves in higher dimensions. In this thesis we shall make a thorough study of the quasinormal modes of black holes in several kinds of background spacetimes. We shall investigate the gravitational radiation given away when highly energetic particles collide with black holes, and also when two black holes collide with each other. Finally, we shall study the properties of gravitational waves in higher dimensions, for instance, we generalize Einstein's quadrupole formula.

We compute for the first time very highly damped quasinormal modes of the (rotating) Kerr black hole. Our numerical technique is based on a decoupling of the radial and angular equations, performed using a large-frequency expansion for the angular separation constant_{s}A_{l m}. This allows us to go much further in overtone number than ever before. We find that the real part of the quasinormal frequencies approaches a non-zero constant value which does not depend on the spin s of the perturbing field and on the angular index l: \omega_R=m\varpi(a). We numerically compute \varpi(a). Leading-order corrections to the asymptotic frequency are likely to be of order 1/\omega_I. The imaginary part grows without bound, the spacing between consecutive modes being a monotonic function of a.

We calculate the quasinormal modes (QNMs) for gravitational perturbations of the Schwarzschild black hole in the five dimensional (5D) spacetime with a continued fraction method. For all the types of perturbations (scalar-gravitational, vector-gravitational, and tensor-gravitational perturbations), the QNMs associated with l=2, l=3, and l=4 are calculated. Our numerical results are summarized as follows: (i) The three types of gravitational perturbations associated with the same angular quantum number l have a different set of the quasinormal (QN) frequencies; (ii) There is no purely imaginary frequency mode; (iii) The three types of gravitational perturbations have the same asymptotic behavior of the QNMs in the limit of the large imaginary frequencies. In Hawking temperature units these frequencies are given by log3 + i2pi(n+1/2) as n goes to infinity, where n is the mode number.

Motivated by recent suggestions that highly damped black hole quasinormal modes (QNM's) may provide a link between classical general relativity and quantum gravity, we present an extensive computation of highly damped QNM's of Kerr black holes. We do not limit our attention to gravitational modes, thus filling some gaps in the existing literature. The frequency of gravitational modes with l=m=2 tends to \omega_R=2 \Omega, \Omega being the angular velocity of the black hole horizon. If Hod's conjecture is valid, this asymptotic behaviour is related to reversible black hole transformations. Other highly damped modes with m>0 that we computed do not show a similar behaviour. The real part of modes with l=2 and m<0 seems to asymptotically approach a constant value \omega_R\simeq -m\varpi, \varpi\simeq 0.12 being (almost) independent of a. For any perturbing field, trajectories in the complex plane of QNM's with m=0 show a spiralling behaviour, similar to the one observed for Reissner-Nordstrom (RN) black holes. Finally, for any perturbing field, the asymptotic separation in the imaginary part of consecutive modes with m>0 is given by 2\pi T_H (T_H being the black hole temperature). We conjecture that for all values of l and m>0 there is an infinity of modes tending to the critical frequency for superradiance (\omega_R=m) in the extremal limit. Finally, we study in some detail modes branching off the so--called ``algebraically special frequency'' of Schwarzschild black holes. For the first time we find numerically that QNM multiplets emerge from the algebraically special Schwarzschild modes, confirming a recent speculation.

We make a thorough investigation of the asymptotic quasinormal modes of the four and five-dimensional Schwarzschild black hole for scalar, electromagnetic and gravitational perturbations. Our numerical results give full support to all the analytical predictions by Motl and Neitzke, for the leading term. We also compute the first order corrections analytically, by extending to higher dimensions, previous work of Musiri and Siopsis, and find excellent agreement with the numerical results. For generic spacetime dimension number D the first-order corrections go as $\frac{1}{n^{(D-3)/(D-2)}}$. This means that there is a more rapid convergence to the asymptotic value for the five dimensional case than for the four dimensional case, as we also show numerically.

We study the late-time tails appearing in the propagation of massless fields (scalar, electromagnetic and gravitational) in the vicinities of a D-dimensional Schwarzschild black hole. We find that at late times the fields always exhibit a power-law falloff, but the power-law is highly sensitive to the dimensionality of the spacetime. Accordingly, for odd D>3 we find that the field behaves as t^[-(2l+D-2)] at late times, where l is the angular index determining the angular dependence of the field. This behavior is entirely due to D being odd, it does not depend on the presence of a black hole in the spacetime. Indeed this tails is already present in the flat space Green's function. On the other hand, for even D>4 the field decays as t^[-(2l+3D-8)], and this time there is no contribution from the flat background. This power-law is entirely due to the presence of the black hole. The D=4 case is special and exhibits, as is well known, the t^[-(2l+3)] behavior. In the extra dimensional scenario for our Universe, our results are strictly correct if the extra dimensions are infinite, but also give a good description of the late time behaviour of any field if the large extra dimensions are large enough.

We investigate the electromagnetic radiation released during the high energy collision of a charged point particle with a four-dimensional Schwarzschild black hole. We show that the spectra is flat, and well described by a classical calculation. We also compare the total electromagnetic and gravitational energies emitted, and find that the former is supressed in relation to the latter for very high energies. These results could apply to the astrophysical world in the case charged stars and small charged black holes are out there colliding into large black holes, and to a very high energy collision experiment in a four-dimensional world. In this latter scenario the calculation is to be used for the moments just after the black hole formation, when the collision of charged debris with the newly formed black hole is certainly expected. Since the calculation is four-dimensional, it does not directly apply to Tev-scale gravity black holes, as these inhabit a world of six to eleven dimensions, although our results should qualitatively hold when extrapolated with some care to higher dimensions.

We present a thorough analysis for the quasinormal (QN) behavior, associated with the decay of scalar, electromagnetic and gravitational perturbations, of Schwarzschild-anti-de Sitter black holes. As it is known the anti-de Sitter (AdS) QN spectrum crucially depends on the relative size of the black hole to the AdS radius. There are three different types of behavior depending on whether the black hole is large, intermediate, or small. The results of previous works, concerning lower overtones for large black holes, are completed here by obtaining higher overtones for all the three black hole regimes. There are two major conclusions that one can draw from this work: First, asymptotically for high overtones, all the modes are evenly spaced, and this holds for all three types of regime, large, intermediate and small black holes, independently of l, where l is the quantum number characterizing the angular distribution; Second, the spacing between modes is apparently universal, in that it does not depend on the field, i.e., scalar, electromagnetic and gravitational QN modes all have the same spacing for high overtones. We are also able to prove why scalar and gravitational perturbations are isospectral, asymptotically for high overtones, by introducing appropriate superpartner potentials.

Self-phase modulation of spherical gravitational wavepackets propagating in a flat space-time in the presence of a tenuous distribution of matter is considered. Analogies with respect to similar effects in nonlinear optics are explored. Self phase modulation of waves emitted from a single source can eventually lead to an efficient energy dilution of the gravitational wave energy over an increasingly large spectral range. An explicit criterium for the occurrence of a significant spectral energy dilution is established.

Gravitational wave solutions to Einstein's equations and their generation are examined in D-dimensional flat spacetimes. First the plane wave solutions are analyzed; then the wave generation is studied with the solution for the metric tensor being obtained with the help of retarded D-dimensional Green's function. Due to the difficulties in handling the wave tails in odd dimensions we concentrate our study in even dimensions. We compute the metric quantities in the wave zone in terms of the energy momentum tensor at retarded time. Some special cases of interest are studied: first the slow motion approximation, where the $D$-dimensional quadrupole formula is deduced. Within the quadrupole approximation, we consider two cases of interest, a particle in circular orbit and a particle falling radially into a higher dimensional Schwarzschild black hole. Then we turn our attention to the gravitational radiation emitted during collisions lasting zero seconds, i.e., hard collisions. We compute the gravitational energy radiated during the collision of two point particles, in terms of a cutoff frequency. In the case in which at least one of the particles is a black hole, we argue this cutoff frequency should be close to the lowest gravitational quasinormal frequency. In this context, we compute the scalar quasinormal frequencies of higher dimensional Schwarzschild black holes. Finally, as an interesting new application of this formalism, we compute the gravitational energy release during the quantum process of black hole pair creation. These results might be important in light of the recent proposal that there may exist extra dimensions in the Universe, one consequence of which may be black hole creation at the Large Hadron Collider at CERN.

Recently Mendonca and Cardoso [Phys. Rev. D, vol. 66, 104009 (2002)] considered nonlinear gravitational wave packets propagating in flat space-time. They concluded that the evolution equation - to third order in amplitude - takes a similar form to what arises in nonlinear optics. Based on this equation, the authors found that nonlinear gravitational waves exhibit self-phase modulation and high harmonic generation leading to frequency up-shifting and spectral energy dilution of the gravitational wave energy. In this Brief Report we point out the fact - a possibility that seems to have been overlooked by Mendonca and Cardoso - that the nonlinear terms in the evolution equation cancels and, hence, that there is no amplitude evolution of the pulse. Finally we discuss scenarios where these nonlinearities may play a role.

In this paper, we consider the gravitational radiation generated by the collision of highly relativistic particles with rotating Kerr black holes. We use the Sasaki-Nakamura formalism to compute the waveform, energy spectra and total energy radiated during this process. We show that the gravitational spectrum for high-energy collisions has definite characteristic universal features, which are independent of the spin of the colliding objects. We also discuss possible connections between these results and the black hole-black hole collision at the speed of light process. With these results at hand, we predict that during the high speed collision of a non-rotating hole with a rotating one, about 35% of the total energy can get converted into gravitational waves. Thus, if one is able to produce black holes at the Large Hadron Collider, as much as 35% of the partons' energy should be emitted during the so called balding phase. This energy will be missing, since we don't have gravitational wave detectors able to measure such amplitudes. The collision at the speed of light between one rotating black hole and a non-rotating one or two rotating black holes turns out to be the most efficient gravitational wave generator in the Universe.

Nonlinear wave interaction of low amplitude gravitational waves in flat space-time is considered. Analogy with optics is established. It is shown that the flat metric space-time is equivalent to a centro-symmetric optical medium, with no second order susceptibility. The lowest order nonlinear effects are those due to the third order nonlinearity and include self-phase modulation and high harmonic generation. These processes lead to an efficient energy dilution of the gravitational wave energy over an increasingly large spectral range.

In this Letter we consider the radial infall along the symmetry axis of an ultra-relativistic point particle into a rotating Kerr black hole. We use the Sasaki-Nakamura formalism to compute the waveform, energy spectra and total energy radiated during this process. We discuss possible connections between these results and the black hole-black hole collision at the speed of light process.

We compute spectra, waveforms, angular distribution and total gravitational energy of the gravitiational radiation emitted during the radial infall of a massless particle into a Schwarzschild black hole. Our fully relativistic approach shows that (i) less than 50% of the total energy radiated to infinity is carried by quadrupole waves, (ii) the spectra is flat, and (iii) the zero frequency limit agrees extremely well with a prediction by Smarr. This process may be looked at as the limiting case of collisions between massive particles traveling at nearly the speed of light, by identifying the energy $E$ of the massless particle with $m_0 \gamma$, $m_0$ being the mass of the test particle and $\gamma$ the Lorentz boost parameter. We comment on the implications for the two black hole collision at nearly the speed of light process, where we obtain a 13.3% wave generation efficiency, and compare our results with previous results by D'Eath and Payne.

We present a complete relativistic analysis for the scalar radiation emitted by a particle in circular orbit around a Schwarzschild-anti-de Sitter black hole. If the black hole is large, then the radiation is concentrated in narrow angles- high multipolar distribution- i.e., the radiation is synchrotronic. However, small black holes exhibit a totally different behavior: in the small black hole regime, the radiation is concentrated in low multipoles. There is a transition mass at $M=0.427 R$, where $R$ is the AdS radius. This behavior is new, it is not present in asymptotically flat spacetimes.

We calculate the quasinormal modes and associated frequencies of the Banados, Zanelli and Teitelboim (BTZ) non-rotating black hole. This black hole lives in 2+1-dimensions in an asymptotically anti-de Sitter spacetime. We obtain exact results for the wavefunction and quasi normal frequencies of scalar, electromagnetic and Weyl (neutrino) perturbations.