
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
(MathissonPapapetrou equations), we show that indeed such an effect takes
place and is a fundamental part of the spincurvature force. The effect arises
whenever there is a current of mass/energy, nonparallel 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 FLWR spacetime.
It is seen to lead to secular orbital precessions potentially observable by
future astrometric experiments and gravitationalwave 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 actionreaction law in postNewtonian gravity.

The ability to directly detect gravitational waves has enabled us to
empirically probe the nature of ultracompact relativistic objects. Several
alternatives to the black holes of classical general relativity have been
proposed which do not have a horizon, in which case a newly formed object (e.g.
as a result of binary merger) may emit echoes: bursts of gravitational
radiation with varying amplitude and duration, but arriving at regular time
intervals. Unlike in previous templatebased approaches, we present a
morphologyindependent search method to find echoes in the data from
gravitational wave detectors, based on a decomposition of the signal in terms
of generalized wavelets consisting of multiple sineGaussians. The ability of
the method to discriminate between echoes and instrumental noise is assessed by
inserting into the noise two different signals: a train of sineGaussians, and
an echoing signal from an extreme massratio inspiral of a particle into a
Schwarzschild vacuum spacetime, with reflective boundary conditions close to
the horizon. We find that both types of signals are detectable for plausible
signaltonoise ratios in existing detectors and their nearfuture upgrades.
Finally, we show how the algorithm can provide a characterization of the echoes
in terms of the time between successive bursts, and damping and widening from
one echo to the next.

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

The ability to detect and scrutinize gravitational waves from the merger and
coalescence of compact binaries opens up the possibility to perform tests of
fundamental physics. One such test concerns the dark, nature of compact
objects: are they really black holes? It was recently pointed out that the
absence of horizons  while keeping the external geometry very close to that
of General Relativity  would manifest itself in a series of echoes in
gravitational wave signals. The observation of echoes by LIGO/Virgo or upcoming
facilities would likely inform us on quantum gravity effects or unseen types of
matter. Detection of such signals is in principle feasible with relatively
simple tools, but would benefit enormously from accurate templates. Here we
analytically individualize each echo waveform and show that it can be written
as a Dyson series, for arbitrary effective potential and boundary conditions.
We further apply the formalism to explicitly determine the echoes of a simple
toy model: the Dirac delta potential. Our results allow to read off a few known
features of echoes and may find application in the modelling for data analysis.

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 selfinteracting 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 Fourierdomain 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 quantumgravity models of compact objects predict microscopic or even
Planckian corrections at the horizon scale. We explore the possibility of
measuring two modelindependent, smokinggun 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 spacebased 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 highlyspinning, supermassive binaries at 120 Gpc
provide unparalleled tests of quantumgravity effects at the horizon scale.

The fate of Cauchy horizons, such as those found inside charged black holes,
is intrinsically connected to the decay of small perturbations exterior to the
event horizon. As such, the validity of the strong cosmic censorship (SCC)
conjecture is tied to how effectively the exterior damps fluctuations. Here, we
study massless scalar fields in the exterior of ReissnerNordstromde Sitter
black holes. Their decay rates are governed by quasinormal modes of the black
hole. We identify three families of modes in these spacetimes: one directly
linked to the photon sphere, well described by standard WKBtype tools; another
family whose existence and timescale is closely related to the de Sitter
horizon. Finally, a third family which dominates for nearextremallycharged
black holes and which is also present in asymptotically flat spacetimes. The
last two families of modes seem to have gone unnoticed in the literature. We
give a detailed description of linear scalar perturbations of such black holes,
and conjecture that SCC is violated in the near extremal regime.

Ultralight bosons and axionlike particles appear naturally in different
scenarios and could solve some longstanding 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
modelindependent 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 blackhole
spin and Proca mass. The observed stability of the inner disk of stellarmass
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 axionlike 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 spacebased
gravitationalwave 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 axionlike particles and dark photons. Overall,
superradiance allows to explore a region of approximately $8$ orders of
magnitude in the mass of ultralight bosons.

We study the ladder operator on scalar fields, mapping a solution of the
KleinGordon equation onto another solution with a different mass, when the
operator is at most first order in derivatives. Imposing the commutation
relation between the d'Alembertian, we obtain the general condition for the
ladder operator, which contains a nontrivial case which was not discussed in
the previous work [V. Cardoso, T. Houri and M. Kimura, Phys.Rev.D 96, 024044
(2017), arXiv:1706.07339]. We also discuss the relation with supersymmetric
quantum mechanics.

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 gravitationalwave physics with
more sensitive detectors. Just as accelerators require larger and larger
energies to probe smaller and smaller scales, more sensitive gravitationalwave
detectors will be probing regions closer and closer to the horizon, potentially
reaching Planck scales and beyond. What may be there, lurking?

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 powerlaw 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 stateoftheart 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 signaltonoise 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 blackhole populations calibrated to observations
and blackhole 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 blackhole
binaries with component masses in the range $[10^3, 10^7]~M_\odot$, which is
not probed by traditional spinmeasurement 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 longlived bosonic condensate outside the horizon. We
use numerical solutions of the perturbed field equations and astrophysical
models of massive and stellarmass black hole populations to compute, for the
first time, the stochastic gravitationalwave 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.

Ultralight scalars have been predicted in a variety of scenarios, and
advocated as a possible component of dark matter. These fields can form compact
regular structures known as boson stars, orin the presence of
horizonsgive rise to nontrivial timedependent scalar hair and a stationary
geometry. Because these fields can be coherent over large spatial extents,
their interaction with "regular" matter can lead to very peculiar effects, most
notably resonances. Here we study the motion of stars in a background
describing black holes surrounded by nonaxially symmetric scalar field
profiles. By analyzing the system in a weakfield approach, we find that the
presence of a scalar field gives rise to secular effects akin to ones existing
in planetary and accretion disks. Particularly, the existence of resonances
between the orbiting stars and the scalar field may enable angular momentum
exchange between them, providing mechanisms similar to planetary migration.
Additionally, these mechanisms may allow \textit{floating orbits}, which are
stable radiating orbits. We also show, in the full relativistic case, that
these effects also appear when there is a direct coupling between the scalar
field and the stellar matter, which can arise due to the presence of a scalar
core in the star or in alternative theories of gravity.

The collision of black holes and the emission of gravitational radiation in
higher dimensional spacetimes are of interest in various research areas,
including the gaugegravity duality, the TeV gravity scenarios evoked for the
explanation of the hierarchy problem, and the largedimensionality limit of
general relativity. We present numerical simulations of headon collisions of
nonspinning, unequalmass black holes starting from rest in general relativity
with $4 \leq D\leq 10$ spacetime dimensions. We compare the energy and linear
momentum radiated in gravitational waves with perturbative predictions in the
extreme mass ratio limit, demonstrating the strength and limitations of black
hole perturbation theory in this context.

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 gravitationalwave physics with
more sensitive detectors. Just as accelerators require larger and larger
energies to probe smaller and smaller scales, more sensitive gravitationalwave
detectors will be probing regions closer and closer to the horizon, potentially
reaching Planck scales and beyond. What may be there, lurking?

The dynamical evolution of selfinteracting scalars is of paramount
importance in cosmological settings, and can teach us about the content of
Einstein's equations. In flat space, nonlinear scalar field theories can give
rise to localized, nonsingular, timedependent, longlived solutions called
{\it oscillons}. Here, we discuss the effects of gravity on the properties and
formation of these structures, described by a scalar field with a double well
potential. We show that oscillons continue to exist even when gravity is turned
on, and we conjecture that there exists a sequence of critical solutions with
infinite lifetime. Our results suggest that a new type of critical behavior
appears in this theory, characterized by modulations of the lifetime of the
oscillon around the scaling law and the modulations of the amplitude of the
critical solutions.

Ladder operators can be useful constructs, allowing for unique insight and
intuition. In fact, they have played a special role in the development of
quantum mechanics and field theory. Here, we introduce a novel type of ladder
operators, which map a scalar field onto another massive scalar field. We
construct such operators, in arbitrary dimensions, from closed conformal
Killing vector fields, eigenvectors of the Ricci tensor. As an example, we
explicitly construct these objects in antide Sitter spacetime (AdS) and show
that they exist for masses above the BreitenlohnerFreedman (BF) bound.
Starting from a regular seed solution of the massive KleinGordon equation
(KGE), mass ladder operators in AdS allow one to build a variety of regular
solutions with varying boundary condition at spatial infinity. We also discuss
mass ladder operator in the context of spherical harmonics, and the relation
between supersymmetric quantum mechanics and socalled Aretakis constants in an
extremal black hole.

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 Earthlike
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 microarcsecond level. This is about 1000
times better than ESA/Gaia's accuracy for the brightest objects and represents
a factor 1030 improvement for the faintest stars (depending on the exact
observational program). In the version submitted to ESA, we proposed an optical
(3501000nm) onaxis 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.

Gravitational waves can teach us not only about sources and the environment
where they were generated, but also about the gravitational interaction itself.
Here we study the features of gravitational radiation produced during the
scattering of a pointlike mass by a black hole. Our results are exact (to
numerical error) at any order in a velocity expansion, and are compared against
various approximations. At large impact parameter and relatively small
velocities our results agree to within percent level with various
postNewtonian and weakfield results. Further, we find good agreement with
scaling predictions in the weakfield/highenergy regime. Lastly, we achieve
striking agreement with zerofrequency estimates.

In the presence of massive bosonic degrees of freedom, rotational
superradiance can trigger an instability that spins down black holes. This
leads to peculiar gravitationalwave signatures and distribution in the
spinmass 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 pulsartiming 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 selfgravitating
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 blackhole
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 EinsteinMaxwell, BransDicke and ChernSimons gravity;
(iii) We assess the ability of present and future gravitationalwave 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
gravitationalwave measurements of the TLNs in a binary inspiral provide a
novel way to test black holes and general relativity in the strongfield
regime.

Analogue systems are a powerful instrument to investigate and understand in a
controlled setting many generalrelativistic effects. Here, we focus on
superradianttriggered instabilities and quasinormal modes. We consider a
compressible hydrodynamic vortex characterized by a polytropic equation of
state, the {\it polytropic hydrodynamic vortex}, a purely circulating system
with an ergoregion but no event horizon. We compute the quasinormal modes of
this system numerically with different methods, finding excellent agreement
between them. When the fluid velocity is larger than the speed of sound, an
ergoregion appears in the effective spacetime, triggering an "ergoregion
instability." We study the details of the instability for the polytropic
vortex, and in particular find analytic expressions for the marginallystable
configuration.

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

Ultralight degrees of freedom coupled to matter lead to resonances, which can
be excited when the Compton wavelength of the field equals a dynamical scale in
the problem. For binaries composed of a star orbiting a supermassive black
hole, these resonances lead to a smokinggun effect: a periastron distance
which {\it stalls}, even in the presence of gravitationalwave dissipation.
This effect, also called a {\it floating orbit}, occurs for generic equatorial
but eccentric orbits and we argue that finitesize effects are not enough to
suppress it.