Proca stars (arXiv:1508.05395) are everywhere regular, asymptotically flat self-gravitating solitons composed of complex, massive vector bosons, forming a macroscopic, star-sized, Bose-Einstein condensate. They have been suggested as a classical, dynamical model of exotic compact objects, black hole mimickers and dark matter candidates. In spherical symmetry, they are qualitatively similar to their scalar cousins, the standard (scalar) boson stars in all respects studied so far. Here we study the tidal deformability and the quadrupolar tidal Love numbers, both electric-type and magnetic-type, of spherically symmetric Proca stars. The equations for the perturbations are derived and the numerical values of the Love numbers for some concrete background solutions are computed. We observe that both Love numbers are qualitatively similar to the ones of the scalar boson stars; in particular, the electric-type (magnetic-type) quadrupolar Love numbers are positive (negative) for Proca stars. Quantitatively, for the same compactness, the electric and magnetic Love numbers in the Proca case are closer (in magnitude) than in the scalar case, with the electric Love numbers having a slightly larger magnitude than the magnetic ones.

The ringdown and shadow of the astrophysically significant Kerr Black Hole (BH) are both intimately connected to a special set of bound null orbits known as Light Rings (LRs). Does it hold that a generic equilibrium BH must possess such orbits? In this letter we prove the following theorem. A stationary, axi-symmetric, asymptotically flat black hole spacetime in 1+3 dimensions, with a non-extremal, topologically spherical, Killing horizon admits, at least, one standard LR outside the horizon for each rotation sense. The proof relies on a topological argument and assumes $C^2$-smoothness and circularity, but makes no use of the field equations. The argument is also adapted to recover a previous theorem establishing that a horizonless ultra-compact object must admit an even number of non-degenerate LRs, one of which is stable.

Hypothetical ultralight bosonic fields will spontaneously form macroscopic bosonic halos around Kerr black holes, via superradiance, transferring part of the mass and angular momentum of the black hole into the halo. Such process, however, is only efficient if resonant: when the Compton wavelength of the field approximately matches the gravitational scale of the black hole. For a complex-valued field, the process can form a stationary, bosonic field-black hole equilibrium state - a black hole with synchronised hair. For sufficiently massive black holes, such as the one at the centre of the M87 supergiant elliptic galaxy, the hairy black hole can be robust against its own superradiant instabilities, within a Hubble time. Studying the shadows of such scalar hairy black holes, we constrain the amount of hair which is compatible with the Event Horizon Telescope (EHT) observations of the M87 supermassive black hole, assuming the hair is a condensate of ultralight scalar particles of mass $\mu\sim 10^{-20}$ eV, as to be dynamically viable. We show the EHT observations set a weak constraint, in the sense that typical hairy black holes that could develop their hair dynamically, are compatible with the observations, when taking into account the EHT error bars and the black hole mass/distance uncertainty.

We consider a static, axially symmetric spacetime describing the superposition of a Schwarzschild black hole (BH) with a thin and heavy accretion disk. The BH-disk configuration is a solution of the Einstein field equations within the Weyl class. The disk is sourced by a distributional energy-momentum tensor and it is located at the equatorial plane. It can be interpreted as two streams of counter-rotating particles, yielding a total vanishing angular momentum. The phenomenology of the composed system depends on two parameters: the fraction of the total mass in the disk, $m$, and the location of the inner edge of the disk, $a$. We start by determining the sub-region of the space of parameters wherein the solution is physical, by requiring the velocity of the disk particles to be sub-luminal and real. Then, we study the null geodesic flow by performing backwards ray-tracing under two scenarios. In the first scenario the composed system is illuminated by the disk and in the second scenario the composed system is illuminated by a far-away celestial sphere. Both cases show that, as $m$ grows, the shadow becomes more prolate. Additionally, the first scenario makes clear that as $m$ grows, for fixed $a$, the geometrically thin disk appears optically enlarged, i.e., thicker, when observed from the equatorial plane. This is to due to light rays that are bent towards the disk, when backwards ray traced. In the second scenario, these light rays can cross the disk (which is assumed to be transparent) and may oscillate up to a few times before reaching the far away celestial sphere. Consequently, an almost equatorial observer sees different patches of the sky near the equatorial plane, as a chaotic "mirage". As $m\rightarrow 0$ one recovers the standard test, i.e., negligible mass, disk appearance.

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.

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.

We perform numerical evolutions of the fully non-linear Einstein-(complex, massive)Klein-Gordon and Einstein-(complex)Proca systems, to assess the formation and stability of spinning bosonic stars. In the scalar/vector case these are known as boson/Proca stars. Firstly, we consider the formation scenario. Starting with constraint-obeying initial data, describing a dilute, axisymmetric cloud of spinning scalar/Proca field, gravitational collapse towards a spinning star occurs, via gravitational cooling. In the scalar case the formation is transient, even for a non-perturbed initial cloud; a non-axisymmetric instability always develops ejecting all the angular momentum from the scalar star. In the Proca case, by contrast, no instability is observed and the evolutions are compatible with the formation of a spinning Proca star. Secondly, we address the stability of an existing star, a stationary solution of the field equations. In the scalar case, a non-axisymmetric perturbation develops collapsing the star to a spinning black hole. No such instability is found in the Proca case, where the star survives large amplitude perturbations; moreover, some excited Proca stars decay to, and remain as, fundamental states. Our analysis suggests bosonic stars have different stability properties in the scalar/vector case, which we tentatively relate to their toroidal/spheroidal morphology. A parallelism with instabilities of spinning fluid stars is briefly discussed.

We construct asymptotically flat, spinning, regular on and outside an event horizon, scalarised black holes (SBHs) in extended scalar-tensor-Gauss-Bonnet models. They reduce to Kerr BHs when the scalar field vanishes. For an illustrative choice of non-minimal coupling, we scan the domain of existence. For each value of spin, SBHs exist in an interval between two critical masses, with the lowest one vanishing in the static limit. Non-uniqueness with Kerr BHs of equal global charges is observed; the SBHs are entropically favoured. This suggests SBHs form dynamically from the spontaneous scalarisation of Kerr BHs, which are prone to a scalar-triggered tachyonic instability, below the largest critical mass. Phenomenologically, the introduction of BH spin damps the maximal observable difference between comparable scalarised and vacuum BHs. In the static limit, (perturbatively stable) SBHs can store over 20% of the spacetime energy outside the event horizon; in comparison with Schwarzschild BHs, their geodesic frequency at the ISCO can differ by a factor of 2.5 and deviations in the shadow areal radius may top 40%. As the BH spin grows, low mass SBHs are excluded, and the maximal relative differences decrease, becoming of order $\sim$ few % for dimensionless spin $j\gtrsim 0.5$. This reveals a spin selection effect: non-GR effects are only significant for low spin. We discuss if and how the recently measured shadow size of the M87 supermassive BH, constrains the length scale of the Gauss-Bonnet coupling.

Classical phenomenological aspects of acoustic perturbations on a draining bathtub geometry where a surface with reflectivity $\mathcal{R}$ is set at a small distance from the would-be acoustic horizon, which is excised, are addressed. Like most exotic compact objects featuring an ergoregion but not a horizon, this model is prone to instabilities when $|\mathcal{R}|^2\approx 1$. However, stability can be attained for sufficiently slow drains when $|\mathcal{R}|^2\lesssim70\%$. It is shown that the superradiant scattering of acoustic waves is more effective when their frequency approaches one of the system's quasi-normal mode frequencies.

Recently, two simple criteria were proposed to assess if vacua emerging from an effective scalar field theory are part of the string "landscape" or "swampland". The former are the vacua that emerge from string compactifications; the latter are not obtained by any such compactification and hence may not survive in a UV completed theory of gravity. So far, these criteria have been applied to inflationary and dark energy models. Here we consider them in the context of solitonic compact objects made up of scalar fields: boson stars. Analysing several models (static, rotating, with and without self-interactions), we find that, in this context, the criteria are not independent. Furthermore, we find the universal behaviour that in the region wherein the boson stars are expected to be perturbatively stable, the compact objects may be part of the landscape. By contrast, in the region where they may be faithful black hole mimickers, in the sense they possess a light ring, the criteria fail (are obeyed) for static (rotating) ultracompact boson stars, which should thus be part of the swampland (landscape). We also consider hairy black holes interpolating between these boson stars and the Kerr solution and establish the part of the domain of existence where the swampland criteria are violated. In interpreting these results one should bear in mind, however, that the swampland criteria are not quantitatively strict.

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.

We consider the status of black hole solutions with non-trivial scalar fields but no gauge fields, in four dimensional asymptotically flat space-times, reviewing both classical results and recent developments. We start by providing a simple illustration on the physical difference between black holes in electro-vacuum and scalar-vacuum. Next, we review no-scalar-hair theorems. In particular, we detail an influential theorem by Bekenstein and stress three key assumptions: 1) the type of scalar field equation; 2) the spacetime symmetry inheritance by the scalar field; 3) an energy condition. Then, we list regular (on and outside the horizon), asymptotically flat BH solutions with scalar hair, organizing them by the assumption which is violated in each case and distinguishing primary from secondary hair. We provide a table summary of the state of the art.

Testing the true nature of black holes - the no-hair hypothesis - will become increasingly more precise in the next few years as new observational data is collected in both the gravitational wave channel and the electromagnetic channel. In this paper we consider numerically generated spacetimes of Kerr black holes with synchronised scalar hair and build stationary models of magnetized thick disks (or tori) around them. Our approach assumes that the disks are not self-gravitating, they obey a polytropic equation of state, the distribution of their specific angular momentum is constant, and they are marginally stable, i.e. the disks completely fill their Roche lobe. Moreover, contrary to existing approaches in the literature, our models are thermodinamically relativist, as the specific enthalpy of the fluid can adopt values significantly larger than unity. We study the dependence of the morphology and properties of the accretion tori on the type of black hole considered, from purely Kerr black holes with varying degrees of spin parameter, namely from a Schwarzschild black hole to a nearly extremal Kerr case, to Kerr black holes with scalar hair with different ADM mass and horizon angular velocity. Comparisons between the disk properties for both types of black holes are presented. The sequences of magnetized, equilibrium disks models discussed in this study can be used as initial data for numerical relativity codes to investigate their dynamical (non-linear) stability and used in tandem with ray-tracing codes to obtain synthetic images of black holes (i.e. shadows) in astrophysically relevant situations where the light source is provided by an emitting accretion disk.

Extended scalar-tensor-Gauss-Bonnet (eSTGB) gravity has been recently argued to exhibit spontaneous scalarisation of vacuum black holes (BHs). A similar phenomenon can be expected in a larger class of models, which includes e.g. Einstein-Maxwell-scalar (EMS) models, where spontaneous scalarisation of electrovacuum BHs should occur. EMS models have no higher curvature corrections, a technical simplification over eSTGB models that allows us to investigate, fully non-linearly, BH scalarisation in two novel directions. Firstly, numerical simulations in spherical symmetry show, dynamically, that Reissner-Nordstr\"om (RN) BHs evolve into a perturbatively stable scalarised BH. Secondly, we compute the non-spherical sector of static scalarised BH solutions bifurcating from the RN BH trunk. Scalarised BHs form an infinite (countable) number of branches, and possess a large freedom in their multipole structure. Unlike the case of electrovacuum, the EMS model admits static, asymptotically flat, regular on and outside the horizon BHs without spherical symmetry and even without any spatial isometries, which are thermodynamically preferred over the electrovacuum state. We speculate on a possible dynamical role of these non-spherical scalarised BHs.

Kerr black holes with synchronised hair [arXiv:1403.2757, arXiv:1603.02687] are a counter example to the no hair conjecture, in General Relativity minimally coupled to simple matter fields (with mass $\mu$) obeying all energy conditions. Since these solutions have, like Kerr, an ergoregion it has been a lingering possibility that they are afflicted by the superradiant instability, the same process that leads to their dynamical formation from Kerr. A recent breakthrough [arXiv:1711.08464] confirmed this instability and computed the corresponding timescales for a sample of solutions. We discuss how these results and other observations support two conclusions: $1)$ starting from the Kerr limit, the increase of hair for fixed coupling $\mu M$ (where $M$ is the BH mass) increases the timescale of the instability; $2)$ there are hairy solutions for which this timescale, for astrophysical black hole masses, is larger than the age of the Universe. The latter conclusion introduces the limited, but physically relevant concept of effective stability. The former conclusion, allows us to identify an astrophysically viable domain of such effectively stable hairy black holes, occurring, conservatively, for $M\mu \lesssim 0.25$. These are hairy BHs that form dynamically, from the superradiant instability of Kerr, within an astrophysical timescale, but whose own superradiant instability occurs only in a cosmological timescale.

Highly spinning Kerr black holes with masses $M = 1 - 100\ M_{\odot}$ are subject to an efficient superradiant instability in the presence of bosons with masses $\mu \sim 10^{-10} - 10^{-12}\ {\rm eV}$. We observe that this matches the effective plasma-induced photon mass in diffuse galactic or intracluster environments ($\omega_{\rm pl} \sim 10^{-10} - 10^{-12}\ {\rm eV}$). This suggests that bare Kerr black holes within galactic or intracluster environments, possibly even including the ones produced in recently observed gravitational wave events, are unstable to formation of a photon cloud that may contain a significant fraction of the mass of the original black hole. At maximal efficiency, the instability timescale for a massive vector is milliseconds, potentially leading to a transient rate of energy extraction from a black hole in principle as large as $\sim 10^{55} \ {\rm erg \, s}^{-1}$. We discuss possible astrophysical effects this could give rise to, including a speculative connection to Fast Radio Bursts.

There is an exciting prospect of obtaining the shadow of astrophysical black holes (BHs) in the near future with the Event Horizon Telescope. As a matter of principle, this justifies asking how much one can learn about the BH horizon itself from such a measurement. Since the shadow is determined by a set of special photon orbits, rather than horizon properties, it is possible that different horizon geometries yield similar shadows. One may then ask how sensitive is the shadow to details of the horizon geometry? As a case study, we consider the double Schwarzschild BH and analyse the impact on the lensing and shadows of the conical singularity that holds the two BHs in equilibrium -- herein taken to be a strut along the symmetry axis in between the two BHs. Whereas the conical singularity induces a discontinuity of the scattering angle of photons, clearly visible in the lensing patterns along the direction of the strut's location, it produces no observable effect on the shadows, whose edges remain everywhere smooth. The latter feature is illustrated by examples including both equal and unequal mass BHs. This smoothness contrasts with the intrinsic geometry of the (spatial sections of the) horizon of these BHs, which is not smooth, and provides a sharp example on how BH shadows are insensitive to some horizon geometry details. This observation, moreover, suggests that for the study of their shadows, this static double BH system may be an informative proxy for a dynamical binary.

We prove the following theorem: axisymmetric, stationary solutions of the Einstein field equations formed from classical gravitational collapse of matter obeying the null energy condition, that are everywhere smooth and ultracompact (i.e., they have a light ring) must have at least two light rings, and one of them is stable. It has been argued that stable light rings generally lead to nonlinear spacetime instabilities. Our result implies that smooth, physically and dynamically reasonable ultracompact objects are not viable as observational alternatives to black holes whenever these instabilities occur on astrophysically short time scales. The proof of the theorem has two parts: (i) We show that light rings always come in pairs, one being a saddle point and the other a local extremum of an effective potential. This result follows from a topological argument based on the Brouwer degree of a continuous map, with no assumptions on the spacetime dynamics, and hence it is applicable to any metric gravity theory where photons follow null geodesics. (ii) Assuming Einstein's equations, we show that the extremum is a local minimum of the potential (i.e., a stable light ring) if the energy-momentum tensor satisfies the null energy condition.

Spherically symmetric bosonic stars are one of the few examples of gravitating solitons that are known to form dynamically, via a classical process of (incomplete) gravitational collapse. As stationary solutions of the Einstein--Klein-Gordon or the Einstein--Proca theory, bosonic stars may also become sufficiently compact to develop light rings and hence mimic, in principle, gravitational-wave observational signatures of black holes (BHs). In this paper, we discuss how these horizonless ultra-compact objects (UCOs) are actually distinct from BHs, both phenomenologically and dynamically. In the electromagnetic channel, the light ring associated phenomenology reveals remarkable lensing patterns, quite distinct from a standard BH shadow, with an infinite number of Einstein rings accumulating in the vicinity of the light ring, both inside and outside the latter. The strong lensing region, moreover, can be considerably smaller than the shadow of a BH with a comparable mass. Dynamically, we investigate the fate of such UCOs under perturbations, via fully non-linear numerical simulations and observe that, in all cases, they decay into a Schwarzschild BH within a time scale of $\mathcal{O}(M)$, where $M$ is the mass of the bosonic star. Both these studies reinforce how difficult it is for horizonless UCOs to mimic BH phenomenology and dynamics, in all its aspects.

East and Pretorius (arXiv:1704.04791) have successfully evolved, using fully non-linear numerical simulations, the superradiant instability of the Kerr black hole (BH) triggered by a massive, complex vector field. Evolutions terminate in stationary states of a vector field condensate synchronised with a rotating BH horizon. We show these end points are fundamental states of Kerr BHs with synchronised Proca hair. Motivated by the "experimental data" from these simulations we suggest a universal (i.e. field-spin independent), analytic model for the subset of BHs with sychronised hair that possess a quasi-Kerr horizon, applicable in the weak hair regime. Comparing this model with fully non-linear numerical solutions of BHs with synchronised scalar or Proca hair, we show the model is accurate for hairy BHs that may emerge dynamically from superradiance, whose domain we identify.

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

We continue our study on the capabilities of present and future X-ray missions to test the nature of astrophysical black hole candidates via X-ray reflection spectroscopy and distinguish Kerr black holes from other solutions of 4-dimensional Einstein's gravity in the presence of a matter field. Here we investigate the case of Kerr black holes with Proca hair [1]. The analysis of a sample of these configurations suggests that even extremely hairy black holes can mimic the iron line profile of the standard Kerr black holes, and, at least for the configurations of our study, we find that current X-ray missions cannot distinguish these objects from Kerr black holes. This contrasts with our previous findings for the case of Kerr black holes with scalar (rather than Proca) hair [2], even though such comparison may be biased by the limited sample. Future X-ray missions can detect the presence of Proca hair, but a theoretical knowledge of the expected intensity profile (currently missing) can be crucial to obtain strong constraints.

X-ray reflection spectroscopy can be a powerful tool to test the nature of astrophysical black holes. Extending previous work on Kerr black holes with scalar hair [1] and on boson stars [2], here we study whether astrophysical black hole candidates may be horizonless, self-gravitating, vector Bose-Einstein condensates, known as Proca stars [3]. We find that observations with current X-ray missions can only provide weak constraints and rule out solely Proca stars with low compactness. There are two reasons. First, at the moment we do not know the geometry of the corona, and therefore the uncertainty in the emissivity profile limits the ability to constrain the background metric. Second, the photon number count is low even in the case of a bright black hole binary, and we cannot have a precise measurement of the spectrum.

We address the astrophysical imaging of a family of deformed Kerr black holes (BHs). These are stationary, asymptotically flat black hole (BH) spacetimes, that are solutions of General Relativity minimally coupled to a massive, complex scalar field: Kerr BHs with scalar hair (KBHsSH). Such BHs bifurcate from the vacuum Kerr solution and can be regarded as a horizon within a rotating boson star. In a recent letter, it was shown that KBHsSH can exhibit very distinct shadows from the ones of their vacuum counterparts. The setup therein, however, considered the light source to be a celestial sphere sufficiently far away from the BH. Here, we analyse KBHsSH surrounded by an emitting torus of matter, simulating a more realistic astrophysical environment, and study the corresponding lensing of light as seen by a very far away observer, to appropriately model ground-based observations of Sgr A*. We find that the differences in imaging between KBHsSH and comparable vacuum Kerr BHs remain, albeit less dramatic than those observed for the corresponding shadows in the previous setup. In particular, we highlight two observables that might allow differentiating KBHsSH and Kerr BHs. The first is the angular size of the photon ring (in a Kerr spacetime) or lensing ring (in a KBHSH spacetime), the latter being significantly smaller for sufficiently non-Kerr-like spacetimes. The second is the existence of an edge in the intensity distribution (the photon ring in Kerr spacetime). This edge can disappear for very non-Kerr-like KBHsSH. It is plausible, therefore, that sufficiently precise Very Long Baseline Interferometric observations of BH candidates can constrain this model.

The present paper is a sequel to our previous work [Y. Ni et al., JCAP 1607, 049 (2016)] in which we studied the iron K$\alpha$ line expected in the reflection spectrum of Kerr black holes with scalar hair. These metrics are solutions of Einstein's gravity minimally coupled to a massive, complex scalar field. They form a continuous bridge between a subset of Kerr black holes and a family of rotating boson stars depending on one extra parameter, the dimensionless scalar hair parameter $q$, ranging from 0 (Kerr black holes) to 1 (boson stars). Here we study the limiting case $q=1$, corresponding to rotating boson stars. For comparison, spherical boson stars are also considered. We simulate observations with XIS/Suzaku. Using the fact that current observations are well fit by the Kerr solution and thus requiring that acceptable alternative compact objects must be compatible with a Kerr fit, we find that some boson star solutions are relatively easy to rule out as potential candidates to explain astrophysical black holes, while other solutions, which are neither too dilute nor too compact are more elusive and we argue that they cannot be distinguished from Kerr black holes by the analysis of the iron line with current X-ray facilities.

In a recent letter, arXiv:1509.00021, it was shown that the lensing of light around rotating boson stars and Kerr black holes with scalar hair can exhibit chaotic patterns. Since no separation of variables is known (or expected) for geodesic motion on these backgrounds, we examine the 2D effective potentials for photon trajectories, to obtain a deeper understanding of this phenomenon. We find that the emergence of stable light rings on the background spacetimes, allows the formation of "pockets" in one of the effective potentials, for open sets of impact parameters, leading to an effective trapping of some trajectories, dubbed quasi-bound orbits. We conclude that pocket formation induces chaotic scattering, although not all chaotic orbits are associated to pockets. These and other features are illustrated in a gallery of examples, obtained with a new ray-tracing code, PYHOLE, which includes tools for a simple, simultaneous visualization of the effective potential together with the spacetime trajectory, for any given point in a lensing image. An analysis of photon orbits allows us to further establish a positive correlation between photon orbits in chaotic regions and those with more than one turning point in the radial direction; we recall that the latter is not possible around Kerr black holes. Moreover, we observe that the existence of several light rings around a horizon (several fundamental orbits, including a stable one), is a central ingredient for the existence of multiple shadows of a single hairy black hole. We also exhibit the lensing and shadows by Kerr black holes with scalar hair, observed away from the equatorial plane, obtained with PYHOLE.

Recently, a family of hairy black holes in 4-dimensional Einstein gravity minimally coupled to a complex, massive scalar field was discovered~\cite{hbh}. Besides the mass $M$ and spin angular momentum $J$, these objects are characterized by a Noether charge $Q$, measuring the amount of scalar hair, which is not associated to a Gauss law and cannot be measured at spatial infinity. Introducing a dimensionless scalar hair parameter $q$, ranging from 0 to 1, we recover (a subset of) Kerr black holes for $q=0$ and a family of rotating boson stars for $q=1$. In the present paper, we explore the possibility of measuring $q$ for astrophysical black holes with current and future X-ray missions. We study the iron K$\alpha$ line expected in the reflection spectrum of such hairy black holes and we simulate observations with Suzaku and eXTP. As a proof of concept, we point out, by analyzing a sample of hairy black holes, that current observations can already constrain the scalar hair parameter $q$, because black holes with $q$ close to 1 would have iron lines definitively different from those we observe in the available data. We conclude that a detailed scanning of the full space of solutions, together with data from the future X-ray missions, like eXTP, will be able to put relevant constraints on the astrophysical realization of Kerr black holes with scalar hair.

Recently, various examples of asymptotically flat, rotating black holes (BHs) with synchronized hair have been explicitly constructed, including Kerr BHs with scalar or Proca hair, and Myers-Perry BHs with scalar hair and a mass gap, showing there is a general mechanism at work. All these solutions have been found numerically, integrating the fully non-linear field equations of motion from the event horizon outwards. Here, we address the spacetime geometry of these solutions inside the event horizon. Firstly, we provide arguments, within linear theory, that there is no regular inner horizon for these solutions. Then, we address this question fully non-linearly, using as a tractable model five dimensional, equal spinning, Myers-Perry hairy BHs. We find that, for non-extremal solutions: $(1)$ the inside spacetime geometry in the vicinity of the event horizon is smooth and the equations of motion can be integrated inwards; $(2)$ before an inner horizon is reached, the spacetime curvature grows (apparently) without bound. In all cases, our results suggest the absence of a smooth Cauchy horizon, beyond which the metric can be extended, for hairy BHs with synchronized hair.

For an observer, the Black Hole (BH) shadow is the BH's apparent image in the sky due to the gravitational lensing of nearby radiation, emitted by an external source. A recent class of solutions dubbed Kerr BHs with scalar hair possess smaller shadows than the corresponding Kerr BHs and, under some conditions, novel exotic shadow shapes can arise. Thus, these hairy BHs could potentially provide new shadow templates for future experiments such as the Event Horizon Telescope. In order to obtain the shadows, the backward ray-tracing algorithm is briefly introduced, followed by numerical examples of shadows of Kerr BHs with scalar hair contrasting with the Kerr analogues. Additionally, an analytical solution for the Kerr shadow is derived in closed form for a ZAMO observer at an arbitrary position.

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.

Bekenstein proved that in Einstein's gravity minimally coupled to one (or many) real, Abelian, Proca field, stationary black holes (BHs) cannot have Proca hair. Dropping Bekenstein's assumption that matter inherits spacetime symmetries, we show this model admits asymptotically flat, stationary, axi-symmetric, regular on and outside an event horizon BHs with Proca hair, for an even number of real (or an arbitrary number of complex) Proca fields. To establish it, we start by showing that a test, complex Proca field can form bound states, with real frequency, around Kerr BHs: stationary Proca clouds. These states exist at the threshold of superradiance. It was conjectured in arXiv:1403.2757, that the existence of such clouds at the linear level implies the existence of a new family of BH solutions at the non-linear level. We confirm this expectation and explicitly construct examples of such Kerr black holes with Proca hair (KBHsPH). For a single complex Proca field, these BHs form a countable number of families with three continuous parameters (ADM mass, ADM angular momentum and Noether charge). They branch off from the Kerr solutions that can support stationary Proca clouds and reduce to Proca stars when the horizon size vanishes. We present the domain of existence of one family of KBHsPH, as well as its phase space in terms of ADM quantities. Some physical properties of the solutions are discussed; in particular, and in contrast with Kerr BHs with scalar hair, some spacetime regions can be counter-rotating with respect to the horizon. We further establish a no-Proca-hair theorem for static, spherically symmetric BHs but allowing the complex Proca field to have a harmonic time dependence, which shows BHs with Proca hair in this model require rotation and have no static limit. KBHsPH are also disconnected from Kerr-Newman BHs with a real, massless vector field.

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.

Using backwards ray tracing, we study the shadows of Kerr black holes with scalar hair (KBHsSH). KBHsSH interpolate continuously between Kerr BHs and boson stars (BSs), so we start by investigating the lensing of light due to BSs. Moving from the weak to the strong gravity region, BSs - which by themselves have no shadows - are classified, according to the lensing produced, as: $(i)$ non-compact, which yield no multiple images; $(ii)$ compact, which produce an increasing number of Einstein rings and multiple images of the whole celestial sphere; $(iii)$ ultra-compact, which possess light rings, yielding an infinite number of images with (we conjecture) a self-similar structure. The shadows of KBHsSH, for Kerr-like horizons and non-compact BS-like hair, are analogous to, but distinguishable from, those of comparable Kerr BHs. But for non-Kerr-like horizons and ultra-compact BS-like hair, the shadows of KBHsSH are drastically different: novel shapes arise, sizes are considerably smaller and multiple shadows of a single BH become possible. Thus, KBHsSH provide quantitatively and qualitatively new templates for ongoing (and future) very large baseline interferometry (VLBI) observations of BH shadows, such as those of the Event Horizon Telescope.

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.

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.

Kerr black holes have their angular momentum, $J$, bounded by their mass, $M$: $Jc\leqslant GM^2$. There are, however, known black hole solutions violating this Kerr bound. We propose a very simple universal bound on the rotation, rather than on the angular momentum, of four-dimensional, stationary and axisymmetric, asymptotically flat black holes, given in terms of an appropriately defined horizon linear velocity, $v_H$. The $v_H$ bound is simply that $v_H$ cannot exceed the velocity of light. We verify the $v_H$ bound for known black hole solutions, including some that violate the Kerr bound, and conjecture that only extremal Kerr black holes saturate the $v_H$ bound.

Kerr black holes with scalar hair are solutions of the Einstein-Klein-Gordon field equations describing regular (on and outside an event horizon), asymptotically flat black holes with scalar hair (arXiv:1403.2757). These black holes interpolate continuously between the Kerr solution and rotating boson stars in D=4 spacetime dimensions. Here we provide details on their construction, discussing properties of the ansatz, the field equations, the boundary conditions and the numerical strategy. Then, we present an overview of the parameter space of the solutions, and describe in detail the space-time structure of the black holes exterior geometry and of the scalar field for a sample of reference solutions. Phenomenological properties of potential astrophysical interest are also discussed, and the stability properties and possible generalizations are commented on. As supplementary material to this paper we make available numerical data files for the sample of reference solutions discussed, for public use.

n-DBI gravity explicitly breaks Lorentz invariance by the introduction of a unit time-like vector field, thereby giving rise to an extra (scalar) degree of freedom. We look for observational consequences of this mode in two setups. Firstly, we compute the parametrized post-Newtonian (PPN) expansion of the metric to first post-Newtonian order. Surprisingly, we find that the PPN parameters are exactly the same as in General Relativity (GR), and no preferred-frame effects are produced. In particular this means that n-DBI gravity is consistent with all GR solar system experimental tests. We discuss the origin of such degeneracy between n-DBI gravity and GR, and suggest it may also hold in higher post-Newtonian order. Secondly, we study gravitational scalar perturbations of a Friedmann-Robertson-Walker space-time with a cosmological constant $\Lambda \geq 0$. In the case of de Sitter space, we show that the scalar mode grows as the universe expands and, in contrast with a canonical scalar field coupled to GR, it does not freeze on super horizon scales.

Q-balls are regular extended `objects' that exist for some non-gravitating, self-interacting, scalar field theories with a global, continuous, internal symmetry, on Minkowski spacetime. Here, analogous objects are also shown to exist around rotating (Kerr) black holes, as non-linear bound states of a test scalar field. We dub such configurations Q-clouds. We focus on a complex massive scalar field with quartic plus hexic self-interactions. Without the self-interactions, linear clouds have been shown to exist, in synchronous rotation with the black hole horizon, along 1-dimensional subspaces - existence lines - of the Kerr 2-dimensional parameter space. They are zero modes of the superradiant instability. Non-linear Q-clouds, on the other hand, are also in synchronous rotation with the black hole horizon; but they exist on a 2-dimensional subspace, delimited by a minimal horizon angular velocity and by an appropriate existence line, wherein the non-linear terms become irrelevant and the Q-cloud reduces to a linear cloud. Thus, Q-clouds provide an example of scalar bound states around Kerr black holes which, generically, are not zero modes of the superradiant instability. We describe some physical properties of Q-clouds, whose backreaction leads to a new family of hairy black holes, continuously connected to the Kerr family.

Massive fields can exist in long-lived configurations around black holes. We examine how the gravitational wave signal of a perturbed black hole is affected by such `dirtiness' within linear theory. As a concrete example, we consider the gravitational radiation emitted by the infall of a massive scalar field into a Schwarzschild black hole. Whereas part of the scalar field is absorbed/scattered by the black hole and triggers gravitational wave emission, another part lingers in long-lived quasi-bound states. Solving numerically the Teukolsky master equation for gravitational perturbations coupled to the massive Klein-Gordon equation, we find a characteristic gravitational wave signal, composed by a quasi-normal ringing followed by a late time tail. In contrast to `clean' black holes, however, the late time tail contains small amplitude wiggles with the frequency of the dominating quasi-bound state. Additionally, an observer dependent beating pattern may also be seen. These features were already observed in fully non-linear studies; our analysis shows they are present at linear level, and, since it reduces to a 1+1 dimensional numerical problem, allows for cleaner numerical data. Moreover, we discuss the power law of the tail and that it only becomes universal sufficiently far away from the `dirty' black hole. The wiggly tails, by constrast, are a generic feature that may be used as a smoking gun for the presence of massive fields around black holes, either as a linear cloud or as fully non-linear hair.

We have recently reported (arXiv:1403.2757) the existence of Kerr black holes with scalar hair in General Relativity minimally coupled to a massive, complex scalar field. These solutions interpolate between boson stars and Kerr black holes. The latter have a well known topologically S^2 ergo-surface (ergo-sphere) whereas the former develop a S^1*S^1 ergo-surface (ergo-torus) in a region of parameter space. We show that hairy black holes always have an ergo-region, and that this region is delimited by either an ergo-sphere or an ergo-Saturn -- i.e. a S^2\oplus (S^1*S^1) ergo-surface. In the phase space of solutions, the ergo-torus can either appear disconnected from the ergo-sphere or pinch off from it. We provide a heuristic argument, based on a measure of the size of the ergo-region, that superradiant instabilities - which are likely to be present - are weaker for hairy black holes than for Kerr black holes with the same global charges. We observe that Saturn-like, and even more remarkable ergo-surfaces, should also arise for other rotating `hairy' black holes.

We show that scalar hair can be added to rotating, vacuum black holes of general relativity. These hairy black holes (HBHs) clarify a lingering question concerning gravitational solitons: if a black hole can be added at the centre of a boson star, as it typically can for other solitons. We argue that it can, but only if it is spinning. The existence of such HBHs is related to the Kerr superradiant instability triggered by a massive scalar field. This connection leads to the following conjecture: a (hairless) black hole which is afflicted by the superradiant instability of a given field must allow hairy generalizations with that field.

We present a family of solutions of Einstein's gravity minimally coupled to a complex, massive scalar field, describing asymptotically flat, spinning black holes with scalar hair and a regular horizon. These hairy black holes (HBHs) are supported by rotation and have no static limit. Besides mass M and angular momentum J, they carry a conserved, continuous Noether charge Q measuring the scalar hair. HBHs branch off from the Kerr metric at the threshold of the superradiant instability and reduce to spinning boson stars in the limit of vanishing horizon area. They overlap with Kerr black holes for a set of (M,J) values. A single Killing vector field preserves the solutions, tangent to the null geodesic generators of the event horizon. HBHs can exhibit sharp physical differences when compared to the Kerr solution, such as J/M^2>1, quadrupole moment larger than J^2/M and larger orbital angular velocity at the innermost stable circular orbit. Families of HBHs connected to the Kerr geometry should exist in scalar (and other) models with more general self interactions.

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.

n-DBI gravity is a gravitational theory which yields near de Sitter inflation spontaneously at the cost of breaking Lorentz invariance by a preferred choice of foliation. We show that this breakdown endows n-DBI gravity with one extra physical gravitational degree of freedom: a scalar graviton. Its existence is established by Dirac's theory of constrained systems. Firstly, studying scalar perturbations around Minkowski space-time, we show that there exists one scalar degree of freedom and identify it in terms of the metric perturbations. Then, a general analysis is made in the canonical formalism, using ADM variables. It is useful to introduce an auxiliary scalar field, which allows recasting n-DBI gravity in an Einstein-Hilbert form but in a Jordan frame. Identifying the constraints and their classes we confirm the existence of an extra degree of freedom in the full theory, besides the two usual tensorial modes of the graviton. We then argue that, unlike the case of (the original proposal for) Horava-Lifschitz gravity, there is no evidence that the extra degree of freedom originates pathologies, such as vanishing lapse, instabilities and strong self-coupling at low energy scales.

We propose a 3 + 1 dimensional model of gravity which results in inflation at early times, followed by radiation- and matter-dominated epochs and a subsequent acceleration at late times. Both the inflation and late time acceleration are nearly de Sitter with a large hierarchy between the effective cosmological constants. There is no scalar field agent of inflation, and the transition from the inflation to the radiation-dominated period is smooth. This model is designed so that it yields, at the cost of giving up on Lorentz invariance in the gravitational sector, the Dirac-Born-Infeld type conformal scalar theory when the universe is conformally flat. It, however, resembles Einstein's gravity with the Gibbons-Hawking-York boundary term in weakly curved space-times.

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.

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.

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.

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.

The similarities between linearized gravity and electromagnetism are known since the early days of General Relativity. Using an exact approach based on tidal tensors, we show that such analogy holds only on very special conditions and depends crucially on the reference frame. This places restrictions on the validity of the "gravito-electromagnetic" equations commonly found in the literature.

A detailed non-linear analysis of the internal structure of spherical, charged black holes that are accreting scalar matter is performed in the framework of the Brans-Dicke theory of gravity. We choose the lowest value of the Brans-Dicke parameter that is compatible with observational constraints. First, the homogeneous approximation is used. It indicates that mass inflation occurs and that the variations of the Brans-Dicke scalar inside the black hole, which could in principle be large in the absence of mass inflation, become small when mass inflation does occur. Then, a full non-linear numerical study of the black hole interior perturbed by a self-gravitating massless uncharged scalar-field is performed. We use an algorithm with adaptive mesh refinement capabilities. In this way, the changes in the internal structure of the black hole caused by mass inflation are determined, as well as the induced variations of the Brans-Dicke scalar, confirming, qualitatively, the indications given by the homogeneous approximation.

A proposal is made for a cosmological D3/D7 model with a constant magnetic flux along the D7 world-volume. It describes an N=2 gauge model with Fayet-Iliopoulos terms and the potential of the hybrid P-term inflation. The motion of the D3-brane towards D7 in a phase with spontaneously broken supersymmetry provides a period of slow-roll inflation in the de Sitter valley, the role of the inflaton being played by the distance between D3 and D7-branes. After tachyon condensation a supersymmetric ground state is formed: a D3/D7 bound state corresponding to an Abelian non-linear (non-commutative) instanton. In this model the existence of a non-vanishing cosmological constant is associated with the resolution of the instanton singularity. We discuss a possible embedding of this model into a compactified M-theory setup.