Astrophysical imaging of Kerr black holes with scalar hair
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.
Over the last half a century, electromagnetic observations have gathered a solid body of evidence for the existence of astrophysical black holes (BHs) – see . Recently, the detection of the gravitational wave event GW150914 , has strengthened the case for astrophysical BHs, now using a completely independent channel. Observations with ever increasing precision, both in the gravitational wave and electromagnetic channels, are expected over the next decade, leading to a more thorough understanding of the true nature of astrophysical BH candidates and, in particular, discriminating the paradigmatic BHs of General Relativity (GR) from more exotic alternatives . It is thus quite timely to obtain the appropriate phenomenology for any theoretically sound alternative model.
Steady-state BHs in vacuum GR have been recognised to be surprisingly simple since the 1970s: they have only two degrees of freedom [5, 6, 7] and are described by the elegant Kerr metric .
Unexpectedly, in the last two years, new families of deformed Kerr BHs could be obtained within simple matter models. The first (and simplest) such family that was discovered, corresponds to Kerr BHs with scalar hair (KBHsSH)  (see also [15, 16]). These are solutions of Einstein’s gravity minimally coupled to a massive, complex, scalar field. This matter content obeys all energy conditions, and the solutions are regular on and outside an event horizon, thus providing a theoretically consistent model. Besides the vacuum limit, wherein they reduce to a subset of vacuum Kerr BHs, KBHsSH have also a solitonic limit, in which they reduce to well-known gravitating solitons known as boson stars [17, 18]. KBHsSH have been found numerically, but a formal proof of their existence was given in . Generalizations with self-interactions were constructed in [20, 21] and with massive Proca fields in .
KBHsSH inherit a well-known property from boson stars: the maximal ADM (Arnowitt, Deser, Misner) mass possible is of the order of the Compton wavelength of the scalar field. This demands, for the existence of either boson stars or KBHsSH with ADM masses of the order of the solar mass, or higher, that the scalar field must be ultra-light, with a mass of order of eV, or smaller – see the discussion in  for models including self-interactions. Such ultralight particles are not present in the standard model of particle physics, but have been predicted in beyond standard model scenarios . Thus, evidence for boson stars or KBHsSH could be faced as evidence for an ultralight particle beyond the standard model. Claiming such evidence, however, will only be possible if the phenomenology of these exotic compact objects is understood in detail and smoking guns that distinguish them from standard compact objects are identified.
Within the electromagnetic channel, a promising class of observations pertains to the lensing of light in the neighbourhood of some supermassive BH candidates, together with their shadows . The shadow of a BH is the region in the observer’s sky comprising the directions of photons that asymptotically approach the event horizon in a backward ray-tracing computation. This type of observable is being targeted by the Event Horizon Telescope (EHT) , a millimeter-wavelength Very Long Baseline Interferometry (VLBI) network, which has the potential to test the Kerr paradigm and constrain alternative models . In particular, the EHT will be able to image the shadow of the supermassive compact object at the center of the Galaxy, Sgr A*. The resolution of the instrument will reach as, which is smaller than the typical size of Sgr A*’s shadow, as.
In a recent letter , it was shown that the shadows of KBHsSH can be sharply distinct from the ones of the vacuum Kerr BHs, in some regions of the parameter space. This study, however, used a setup (first introduced in ) which had a primarily goal of producing visually striking images of the lensing, and the corresponding differences with respect to Kerr, rather than simulating a realistic astrophysical environment. Thus, a natural question is whether introducing a more realistic astrophysical environment masks the peculiar shadows obtained in . The main purpose of the current paper is to address this question.
We shall consider the same configurations of KBHsSH studied in detail in  and study their lensing, but with three main differences with respect to the former study. Firstly, the ray tracing will be done with a well tested code (Gyoto ), completely independent from that used in . Secondly, the BHs are surrounded by an accretion torus, using the same model studied recently around Kerr BHs  and boson stars . This is the light source, rather than a faraway celestial sphere, as in the setup considered in . Thirdly, we consider the observer to be at a realistic distance to model a ground-based observation of Sgr A, whereas the observations collected in  were from the viewpoint of a much closer observer. Our main conclusion is that, even with this more realistic astrophysical setup, observable differences remain between the hairy BHs and their vacuum counterparts, even though the astrophysical environment partly masks the strikingly different shapes observed in . More concretely, the total flux difference in the imaging of KBHsSH and comparable Kerr BHs is small, for the sample of cases analysed; but the photon ring size difference can have an observable signature. In more extreme cases the shadow, in the astrophysical setup, is essentially erased, in sharp contrast to that of the comparable Kerr BH.
This paper is organized as follows. In Section 2 we briefly review KBHsSH and, in particular, the sample of backgrounds to be addressed in this paper. In Section 3 we benchmark the use of Gyoto by comparing its results with the ones previously obtained for the shadows of KBHsSH. In Section 4 we introduce the astrophysical environment around KBHsSH and produce the corresponding imaging. We conclude in Section 5 with a summary of our results and an outlook.
2 The black hole backgrounds
The simplest KBHsSH are solutions of the Einstein-Klein-Gordon model, where the scalar field is free (no self-interactions), complex and massive (see [20, 21] for generalizations including self-interactions). This model is described by the action:
Throughout we use units with . This model has two constants, Newton’s constant and the scalar field mass . We take the two associated natural scales to be Planck’s mass , and the mass scale for boson stars/KBHsSH, which reads
The hairy Kerr BH solutions are found with the following metric ansatz for a stationary axisymmetric and circular spacetime:
being a constant (representing the radial coordinate of the event horizon), and , and are functions of the spheroidal coordinates . The vacuum Kerr metric can be written in this coordinate system. The explicit form of the coefficients can be found in [16, 22]. Note that the parameters (in ) and (in ) relate as . In the following we shall dub these as spheroidal prolate (SP) coordinates, Appendix A of .
KBHsSH form a countable number of families, labeled by the azimuthal harmonic index and the number of nodes in the scalar field. Here we shall focus on a particular member of this family, with the lowest azimuthal harmonic index (equal to one) and the lowest number of nodes (no nodes – see  for a discussion of the general case). The latter defines fundamental states; excited solutions are likely unstable, towards decay into fundamental solutions, as in the case of boson stars . The corresponding scalar field ansatz is
where is a constant. In the following we shall address three particular solutions of the model (1), with the ansatz (3)–(5), which are dubbed Configurations I-III in . The numerical data for these solutions is publicly available . Table 1 provides a brief description of the physical parameters of these configurations (from ).
These three KBHsSH configurations are compared to Kerr BH solutions with the same values of ADM mass and total angular momentum. These comparable Kerr solutions are denoted Kerr in . They are comparable observation-wise, because observations typically give access to the parameters and .
Configuration I is a “rather Kerr-like” KBHSH. Only of the mass and of angular momentum are stored in the scalar field. The horizon is also Kerr-like, in the sense that the Kerr bound is obeyed in terms of horizon quantities – a property which is not mandatory for other KBHsSH . Its Kerr-like shadow’s average radius is only a few percent smaller than that of its comparable Kerr counterpart. The latter is an almost extremal BH, with .
Configuration II is a “not-so-Kerr-like” KBHSH. In this case, the majority of the mass () and angular momentum () are stored in the scalar field. The horizon is non-Kerr-like, in the sense that it violates the Kerr bound in terms of horizon quantities. Its shadow is not only 25 smaller than that of the comparable Kerr BH , but it has also a peculiar shape – more square – than that observed for any vacuum Kerr BH. The comparable Kerr BH is less extremal than the corresponding one for configuration I, with .
Finally, configuration III, is a “very non-Kerr-like” KBHSH. Almost all mass () and angular momentum () are stored in the scalar field. The horizon is very non-Kerr-like, violating the Kerr bound in terms of horizon quantities by a factor of 6. The horizon shape is quite exotic – it is delimeted by a non-convex curve and there are multiple disconnected shadows. The comparable Kerr BH has .
In the following two sections we will readdress these configurations, firstly recomputing their shadows using Gyoto, and comparing the results with those previously obtained, and then performing their imaging in the astrophysical setup.
3 Ray-tracing setup and shadow computations
We use in this article the ray-tracing code Gyoto, which is open-source software [35, 28]. We employ this code to integrate numerically null geodesics in different KBHsSH numerical spacetimes and integrate the radiative transfer equation inside an accretion structure surrounding the BH. The geodesic integration is performed backward in time from a distant observer. In our setup, the observer is located at a radial coordinate corresponding to the distance between the Earth and the Galactic Center, i.e. kpc . The KBHsSH and Kerr solutions that we use are expressed in SP coordinates , as defined in . In order to relate the SP radial coordinate , expressed in units of , to a physical distance of kpc, we need to fix the mass of Sgr A*. Throughout this article we use , where is the ADM mass of the BH. The integration is performed using a Runge-Kutta-Fehlberg adaptive-step integrator at order 7/8 as implemented in the boost C++ library. A recent study  has demonstrated the ability of Gyoto to ray trace accurately even over such very large distances. The observer is located at some fixed inclination , where is the angle between the axis of rotation and the observer. This particular value of inclination is inherited from a previous study  where it was shown to be able to reproduce well the observed features of Sgr A* modeled as a Kerr BH. The observer is modeled by a screen, with every pixel corresponding to some direction of photon incidence. The total computed field of view is typically of as, unless otherwise stated. The various pixels of the screen are assigned with the value of the specific intensity transported by the photon corresponding to the pixel’s direction of incidence. The output of the ray-tracing calculation is thus a map of specific intensity over some small field of view, which, in the following, will be called an image.
We use the numerical metrics corresponding to the KBHsSH and comparable Kerr solutions described in the previous section, which are publicly available . We extended the Lorene library  to make it able to read these metrics and translate them to multi-domain spectral grids. Gyoto is then able to perform ray tracing using such numerical spacetimes. The Lorene class ScalarBH used for generating Gyoto-compatible metrics from the KBHsSH raw data is publicly available in the latest version of Lorene. We remark that in preparing the framework for this study, we have compared rotating boson star solutions numerically generated by the Kadath code  in  with those generated by the code used in , Fidisol/Cadsol , finding a very good agreement.
In this section, we want to compute the shadows of the three KBHsSH, configurations I, II and III, and their comparable Kerr BHs. Our aim is first to determine what the shadows look like in an “astrophysically neutral” setup ( no emission of electromagnetic radiation), which will be useful for discussing the astrophysically realistic images later on. We also want to compare our computations with the previous calculations developed in , to ensure consistency between the two completely independent ray-tracing codes. As a consequence, we do not consider any source of radiation and simply trace null geodesics until they approach the event horizon. The ray-traced images show only two intensity values: either when the backtraced photon came arbirary close to the event horizon, or otherwise.
Fig. 1 shows the shadows of the KBHSH and comparable Kerr BH of our configuration I. They look very similar to Fig. 5, top-left and top-middle panels of . Note, however, that therein an inclination of was used, whereas here we take . The inclination impact is very mild, for this configuration, and the change in is barely noticible in the shadow. Observe that in all solutions presented in this paper the BHs are rotating such that the left hand side of the image is moving towards the reader.
Fig. 2 shows the shadows of the KBHSH and comparable Kerr BH of our configuration II. They are very similar to the corresponding one in Fig. 5 of  (second line from the top). It is clear that the angular size of the shadow is smaller (by ) for the KBHSH as compared to the comparable Kerr setup.
Fig. 3 shows the shadows of the KBHSH and comparable Kerr BH of configuration III. The two shadows are extremely different, as already showed in . The KBHSH shadow is rather different from the lower-left panel of Fig. 5 of , but the dominating difference is due to the change in inclination since our shadow is not symmetric with respect to the horizontal axis.
In order to be able to compare more precisely with , Fig. 4 shows the shadow of the KBHSH of configuration III as observed from the equatorial plane () and for a smaller field of view (as). The left panel shows the shadow as observed from the Earth, while the right panel shows the shadow as observed by a very close observer, at the same radial coordinate as used in .
4 Accretion torus and images
We consider a simple toroidal accretion structure surrounding the various BHs of configurations I, II and III. This model is the same as presented in  and recently applied to boson stars . It is based on the magnetized torus model of . We will only present this accretion model very briefly and refer to  for details. The accretion torus model that we consider is made of a non-self-gravitating perfect polytropic fluid, circularly orbiting with constant specific angular momentum , where is the fluid 4-velocity, which is completely fixed by imposing the constancy of and the circular motion. The original work of  considers a toroidal magnetic field. However,  showed that the synchrotron images are not sensitive to the direction of the magnetic field. We thus simplify the original model of  by considering an isotropized magnetic field. The energy-momentum conservation is readily integrated by considering a polytropic-like equation of state with the gas pressure and enthalpy (equal to the sum of the gas pressure and total energy density) related through
where is the polytropic exponent, and is the polytropic constant. Conservation of energy-momentum then leads to
where is a potential (known throughout spacetime because the 4-velocity is fixed once is chosen) and is its value at the surface of the torus. This immediately gives
where is the central enthalpy, and is the potential value at the center of the torus. Thus, the enthalpy is analytically known throughout the torus. The values of gas pressure (from the polytropic relation), magnetic pressure (, where is a chosen parameter), magnetic field () and temperature (from the perfect-gas relation) immediately follow. This torus emits thermal synchrotron radiation, following the prescription given in . We note that the term accretion may be misleading given that our model is stationary. However, we consider this torus as a simple model for an instantaneous snapshot of a more realistic time-evolving accretion flow, so we keep referring to it as an accretion torus. The torus model is fully described by the choice of a particular background spacetime plus the choice of a set of astrophysical parameters. These are the torus constant angular momentum and inner radius (fixing these two parameters sets the outer radius of the torus), the inclination of the observer , the torus central electron number density and temperature , the polytropic exponent relating pressure and enthalpy, and the gas-to-magnetic pressure ratio . Among these astrophysical parameters, only will be varied. The other are fixed to the values given in Table 2 and were chosen to give reasonable values of fluxes as compared to millimeter observed data of Sgr A*, and to get a rather compact structure ( not a very extended torus but rather a structure extending over a small radial distance).
|central density ()|
|central electron temperature (K)|
|gas to magnetic pressure ratio|
For one given spacetime, the inner radius is fixed such that the (SP) radial extent of the structure (outer radius minus inner radius) is close to . This ensures to obtain an emitting structure with an angular size satisfying the constraint imposed by the first EHT data . Thus, we choose to vary from spacetime to spacetime in order to maintain an approximately constant angular size of the structure as observed from Earth. Note, however, that the range of variation of is very small and that this quantity always stays close to .
The Gyoto code allows to ray trace photons from a distant observer and integrate the radiative transfer equation through the optically thin synchrotron-emitting accretion torus, thus producing a map of specific intensity, an image. We consider an observed frequency of GHz for the ray-traced photons, corresponding to the frequency used for the early EHT science . The aim of this section is to compare each KBHsSH image to the image of the comparable Kerr configuration, i.e. that having the same ADM mass and the same total angular momentum (cf. Sec. 2).
As a first check, we have verified that Kerr BH images computed with SP coordinates were indistinguishable from Kerr images computed with the more standard Boyer-Lindquist (BL) coordinates. Fig. 5 (upper and lower-right panels) shows that this is indeed so: the comparable Kerr image of configuration I computed using the numerical solution in SP coordinates differs by only with respect to the same image computed using the analytical solution in Boyer-Lindquist coordinates.
The lower-left and upper panels of Fig. 5 show the ray-traced images of the KBHSH and comparable Kerr BH of configuration I.
The KBHSH and Kerr images are very similar, which is not surprising given that the KBHSH of configuration I was chosen to be very Kerr-like. The flux difference between the KBHSH and Kerr images is , which is vanishingly small as far as spectral observations are concerned. It may seem surprising that this flux ratio is actually smaller than the ratio between the analytical Kerr BL and numerical Kerr SP spacetimes discussed above. This fact is explained by two causes. First, the Kerr BL and Kerr SP spacetimes use different coordinates. This introduces a numerical error that is not present when comparing Kerr SP and KBHSH spacetimes that use the same coordinates. Second, the average specific intensity ratio, when comparing pixel by pixel, is of for the Kerr BL/SP comparison and for the Kerr SP/KBHSH comparison. The flux ratio is different, because most of the flux actually comes from the center of the torus, and in these regions, the numerical error due to the change of coordinates between Kerr BL and SP dominates over the difference between the Kerr SP and KBHSH spacetimes. Arguably the most interesting feature in such strong-field image is the thin ring of illuminated pixels at the center of the image, the photon ring. The KBHSH and Kerr photon rings are very similar in shape, but their sizes are differing by . The recent study of  discusses the current measured error on the ratio of the mass of Sgr A* over its distance, and shows that it is of . Moreover, this study advocates that EHT data could lead to a constraint of the photon ring angular size (the most advanced goal of the EHT) with a precision of . The size difference of of the KBHSH and Kerr photon rings that we report here is thus most probably unobservable.
Fig. 6 shows the images of the accretion torus surrounding the KBHSH and comparable Kerr solution of configuration II.
The flux difference between the KBHSH and Kerr images is of , which is still very small with respect to the error bars of spectral observations at GHz. The first difference of the KBHSH image with respect to Kerr is the fact that there are two edges in the intensity distribution: there is a brighter ring at the center of the image that looks like a distorted Kerr photon ring (the left part is visibly different from a portion of a circle), and inside this first ring, a second, fainter one. The interpretation of these two features is helped by considering the shadows in Fig. 5 of . These images show clearly that as the spacetime becomes more non-Kerr-like, a region develops around the shadow that is affected by intense and complex lensing effects (this is particularly clear in the third row from the top in Fig. 5 of ). In the following we will call this region the hyper-lensed region, and its outer boundary the lensing ring (see Fig. 7).
The brighter ring in the left panel of our Fig. 6 is this lensing ring, while the fainter ring is the photon ring. We have checked that photons forming the photon ring approach closer to the event horizon (in radial coordinate) than photons forming the lensing ring. This is the expected behavior given that the photon ring is the projection of the innermost photon orbit and marks the innermost limit a photon can visit without falling into the event horizon. It is very probable that EHT would detect only the brighter lensing ring, so we will focus on it in the following. The KBHSH lensing ring is distorted, in the sense that it does not look like the photon ring of any Kerr BH (see  for an overview of Kerr photon rings). However, this distortion is extremely tiny and would most probably be unnoticed by EHT observations. Still, the KBHSH lensing ring is smaller in angular size, as compared to the Kerr photon ring. The difference in angular size between the KBHSH lensing ring and Kerr photon ring reaches , which is bigger than the error on the ratio as discussed above, and also bigger than the foreseen precision of the EHT data for constraining the photon ring angular size of Sgr A*. Consequently, should the EHT be capable of giving a constraint on the angular size of the photon ring of Sgr A* to within and should this value be too small to be compatible with a Kerr BH of mass at distance , this would support the existence of an alternative compact object such as a KBHSH at the Galactic center. Note that the comparable Kerr BH of configuration II has a spin of . A Kerr BH with spin would have a photon ring smaller by approximately , so the difference of between the Kerr and KBHSH ring angular sizes cannot be accounted for by varying the Kerr BH spin.
Finally, Fig. 8 shows images of our accretion torus model in the KBHSH and comparable Kerr spacetimes of configuration III. The flux difference between the two images is of , still smaller than the observational error bars.
The most striking feature of the left panel in this image is the central “noisy” region, full of radiation. This is the hyper-lensed region. Its outer boundary, the lensing ring, nearly coincides with the comparable Kerr photon ring, and is bigger and less distorted than the lensing ring of the KBHSH of configuration II. This is in agreement with the findings of . Why is the hyper-lensed region of the configuration III KBHSH “noisy” while the hyper-lensed region of its configuration II counterpart is empty of radiation? Because of the very small angular size of the hyper-lensed region in configuration II. Photons from inside the hyper-lensed region will carry a non-negligible amount of radiation provided they visit the central regions of the accretion torus where most of the radiation is produced. In configuration II, no photon forming the inside of the hyper-lensed region visit these central parts of the torus. Only photons forming the lensing ring do so, because the lensing ring corresponds to extremely bent photons and a large region of spacetime (including the inner parts of the torus) is projected to this thin ring of pixels, which is thus bright. On the contrary, many photons from the much bigger hyper-lensed region of the KBHSH of configuration III visit the central parts of the torus thus leading to a lot of radiation being located in this region of the image. The BH shadow is not visible in the left panel of our Fig. 8. It is actually so small in angular size (see Fig. 3) that it is nearly completely erased by the radiation emitted by the part of the torus located in between the BH and the observer. The photon ring of this spacetime, being the outer boundary of the shadow, is invisible. In this spacetime, the important features, as far as EHT observations are concerned, are the hyper-lensed region and lensing ring. The shadow and photon ring are not interesting observation-wise.
The hyper-lensed region with a large angular scale is a very interesting feature of the KBHSH spacetime of configuration III because it may also lead to observational difference, for different reasons than for the KBHSH of configuration II. The obvious difference between the two panels of Fig. 8 is the fact that there is flux all around the hyper-lensed region in the KBHSH spacetime, while the shadow of the Kerr spacetime is free of radiation, except for the radiation emitted in the foreground by the part of the torus in between the BH and the observer. In the KBHSH spacetime, of the total flux is located in the hyper-lensed region. In the Kerr spacetime, of the total flux is located in the shadow, due to emission in the foreground. Thus, approximately of the total flux of the KBHSH image of configuration III is a “hairy flux”, due to the hyper-lensing effects specific to the KBHSH spacetime. It is probable that an algorithm (see [43, 45]) trying to detect a Kerr BH shadow on an image similar to the left panel of Fig. 8 would not converge because there is nowhere in the image a photon-ring-like structure: there is no edge ( strong and localized gradient) in the intensity distribution.
5 Final Remarks
The perspective of the near-future EHT observations of Sgr A* makes it very timely to study the observable counterparts of compact objects alternative to the Kerr BH. Among the many such objects, KBHsSH are particularly interesting because (1) they are exact solutions of Einstein field equations, (2) they only necessitate the addition of a scalar field, a rather ubiquitous object in theoretical physics, with the Higgs boson being an example of a fundamental scalar field in Nature, (3) they do not imply adding any astrophysically unclear elements (like the thin shell of gravastars ).
This article shows that KBHsSH might be observationally differentiated with respect to Kerr BHs by using EHT observations. A too Kerr-like KBHSH would obviously be impossible to differentiate, as our configuration I illustrates. However, for sufficiently non-Kerr-like KBHSH, we have highlighted two features that may allow making an observational difference. The first such feature, illustrated by our configuration II, is linked to the angular size of the photon/lensing rings. The lensing ring is a specific feature to KBHsSH spacetimes, defined in Section 4, which would observationally be interpreted as a Kerr photon ring. A sufficiently non-Kerr-like KBHSH of the same mass and spin as a Kerr BH has a lensing ring smaller in angular size than the photon ring of the comparable Kerr image. This size difference reaches in our configuration II, which is sufficient to be detectable by EHT data and to be non degenerate with a Kerr BH. The increasing non-Kerrness of a KBHSH spacetime is accompanied by the development of the central hyper-lensed region (also defined in Section 4), the outer boundary of which is the lensing ring. This hyper-lensed region is the second feature that may allow differentiating KBHsSH from Kerr BHs. Its outer boundary is rather close to the photon ring of a comparable Kerr BH for very non-Kerr-like KBHsSH. As a consequence, the shadow region of a Kerr spacetime, which is characterized by an edge in the intensity distribution, is replaced by a central “noisy” region without an edge. This edge in intensity distribution being the signal that algorithms investigating such images look for, it is likely that such algorithms would fail finding a shadow region in a KBHSH spacetime such as our configuration III.
Finally, whereas KBHsSH are certainly an interesting theoretical model for phenomenological deviations from the Kerr paradigm,
one may wonder about their realization as astrophysical objects. An obvious necessary condition is the existence of appropriate scalar fields in Nature, as discussed in the Introduction, either as fundamental fields, or, eventually, as a coarse-graining of more fundamental degrees of freedom. Another central point is the dynamical stability of these configurations. In this respect, we would like to stress that the vacuum Kerr BH (with no scalar field excited) is an unstable solution of the model (1).
We would like to thank P. Cunha and H. Rúnarsson for discussions. C. H. and E. R. acknowledge funding from the FCT-IF programme. This work was partially supported by the H2020-MSCA-RISE-2015 Grant No. StronGrHEP-690904, by the CIDMA project UID/MAT/04106/2013, and by the ANR Grant No. 12-BS01-012-01.
- Here we are taking the viewpoint that astrophysical BHs are essentially electrically neutral.
- The instability of Kerr BHs against low frequency modes of bosonic fields was first discussed by Press and Teukolsky in the setup of a Kerr BH surrounded by a mirror . Subsequently, Damour et al. found that the confining mechanism provided by the mirror is naturally present if the bosonic field is massive . See  for an overview of superradiant instabilities.
- P. V. P. Cunha, C. A. R. Herdeiro, E. Radu and H. F. Rúnarsson, Phys. Rev. Lett. 115, 211102 (2015), [1509.00021].
- R. Narayan and J. E. McClintock, Observational Evidence for Black Holes, in General Relativity and Gravitation: A Centennial Perspective, edited by A. Ashtekar, B. K. Berger, J. Isenberg and M. MacCallum, Cambridge Univ. Press, Cambridge, 2015, 1312.6698.
- Virgo, LIGO Scientific, B. P. Abbott et al., Phys. Rev. Lett. 116, 061102 (2016), [1602.03837].
- E. Berti et al., Class. Quant. Grav. 32, 243001 (2015), [1501.07274].
- B. Carter, Phys. Rev. Lett. 26, 331 (1971).
- S. W. Hawking, Commun. Math. Phys. 25, 152 (1972).
- P. T. Chrusciel, J. L. Costa and M. Heusler, Living Rev. Rel. 15, 7 (2012), [1205.6112].
- R. P. Kerr, Phys.Rev.Lett. 11, 237 (1963).
- T. Johannsen, Phys. Rev. D 87, 124017 (2013), [1304.7786].
- L. Rezzolla and A. Zhidenko, Phys. Rev. D 90, 084009 (2014), [1407.3086].
- V. Cardoso, P. Pani and J. Rico, Phys. Rev. D 89, 064007 (2014), [1401.0528].
- R. Konoplya, L. Rezzolla and A. Zhidenko, Phys. Rev. D 93, 064015 (2016), [1602.02378].
- M. Ghasemi-Nodehi and C. Bambi, 1604.07032.
- C. A. R. Herdeiro and E. Radu, Phys. Rev. Lett. 112, 221101 (2014), [1403.2757].
- C. A. R. Herdeiro and E. Radu, Int. J. Mod. Phys. D23, 1442014 (2014), [1405.3696].
- C. Herdeiro and E. Radu, Class. Quant. Grav. 32, 144001 (2015), [1501.04319].
- F. Schunck and E. Mielke, Class.Quant.Grav. 20, R301 (2003), [0801.0307].
- S. L. Liebling and C. Palenzuela, Living Rev. Rel. 15, 6 (2012), [1202.5809].
- O. Chodosh and Y. Shlapentokh-Rothman, 1510.08025.
- B. Kleihaus, J. Kunz and S. Yazadjiev, Phys. Lett. B744, 406 (2015), [1503.01672].
- C. A. R. Herdeiro, E. Radu and H. Rúnarsson, Phys. Rev. D 92, 084059 (2015), [1509.02923].
- C. Herdeiro, E. Radu and H. Rúnarsson, 1603.02687.
- A. Arvanitaki, S. Dimopoulos, S. Dubovsky, N. Kaloper and J. March-Russell, Phys.Rev. D81, 123530 (2010), [0905.4720].
- H. Falcke, F. Melia and E. Agol, Astrophys.J. 528, L13 (2000), [astro-ph/9912263].
- S. Doeleman et al., Imaging an Event Horizon: submm-VLBI of a Super Massive Black Hole, in astro2010: The Astronomy and Astrophysics Decadal Survey, , ArXiv Astrophysics e-prints Vol. 2010, p. 68, 2009, [0906.3899].
- A. E. Broderick, T. Johannsen, A. Loeb and D. Psaltis, ApJ 784, 7 (2014), [1311.5564].
- A. Bohn et al., Class. Quant. Grav. 32, 065002 (2015), [1410.7775].
- F. H. Vincent, T. Paumard, E. Gourgoulhon and G. Perrin, Class. Quant. Grav. 28, 225011 (2011), [1109.4769].
- F. H. Vincent, W. Yan, O. Straub, A. A. Zdziarski and M. A. Abramowicz, A&A 574, A48 (2015), [1406.0353].
- F. H. Vincent, Z. Meliani, P. Grandclément, E. Gourgoulhon and O. Straub, Class. Quant. Grav. 33, 105015 (2016), [1510.04170].
- P. Grandclément, C. Somé and E. Gourgoulhon, Phys. Rev. D 90, 024068 (2014).
- J. Balakrishna, E. Seidel and W.-M. Suen, Phys. Rev. D 58, 104004 (1998), [gr-qc/9712064].
- C. A. R. Herdeiro and E. Radu, Int. J. Mod. Phys. D24, 1544022 (2015), [1505.04189].
- S. Gillessen et al., ApJ 692, 1075 (2009), [0810.4674].
- M. Grould, T. Paumard and G. Perrin, A&A, in press (2016), [1605.04195].
- http://luth.obspm.fr/ luthier/grandclement/kadath.html.
- W. Schonauer and R. Weiß, J. Comput. Appl. Math. 27, 279 (1989).
- S. S. Komissarov, MNRAS 368, 993 (2006), [astro-ph/0601678].
- S. S. Doeleman et al., Nature 455, 78 (2008), [0809.2442].
- D. Psaltis, F. Özel, C.-K. Chan and D. P. Marrone, ApJ 814, 115 (2015), [1411.1454].
- C.-K. Chan, D. Psaltis and F. Özel, ApJ 777, 13 (2013), [1303.5057].
- T. Johannsen et al., Physical Review Letters 116, 031101 (2016), [1512.02640].
- P. O. Mazur and E. Mottola, Proceedings of the National Academy of Science 101, 9545 (2004), [gr-qc/0407075].
- W. H. Press and S. A. Teukolsky, Nature 238, 211 (1972).
- T. Damour, N. Deruelle and R. Ruffini, Lett.Nuovo Cim. 15, 257 (1976).
- R. Brito, V. Cardoso and P. Pani, Lect. Notes Phys. 906, pp.1 (2015), [1501.06570].
- N. Sanchis-Gual, J. C. Degollado, P. J. Montero, J. A. Font and C. Herdeiro, Phys. Rev. Lett. 116, 141101 (2016), [1512.05358].
- P. Bosch, S. R. Green and L. Lehner, Phys. Rev. Lett. 116, 141102 (2016), [1601.01384].