Black hole shadow of a rotating polytropic black hole by the Newman--Janis algorithm without complexification

In this work, starting from a spherically symmetric polytropic black hole, a rotating solution is obtained by following the Newman--Janis algorithm without complexification. Besides studying the horizon, the static conditions and causality issues of the rotating solution, we obtain and discuss the shape of its shadow. Some other physical features as the Hawking temperature and emission rate of the rotating polytropic black hole solution are also discussed.


I. INTRODUCTION
After the historical direct detection of gravitational waves from a black hole (BH) merger by the LIGO collaboration in 2015 [1], and the first image of the supermassive BH located at the centre of the giant elliptical galaxy Messier 87 (M87) by the Event Horizon Telescope (EHT) project [2] only a few weeks back, there is currently a lot of ongoing research in investigating the properties of BHs in several different contexts. The LIGO direct detection provided us with the strongest evidence so far that BHs do exist in Nature and that they merge, while at the same time it offers us the tools to test strong gravitational fields using radio wave astronomy. However, the LIGO detection has provided us with no information about the horizon of the BH, which after all is its defining property.
The gravitational field produced by a BH is so strong that nothing, not even photons, can escape. Therefore, * econtreras@yachaytech.edu.ec † jmramirez@yachaytech.edu.ec ‡ angel.rincon@pucv.cl § grigorios.panotopoulos@tecnico.ulisboa.pt ¶ p.bargueno@uniandes.edu.co a BH cannot be seen directly. There is, however, the possibility of seeing a dark shadow of a BH via strong gravitational lensing and photon capture at the horizon, if the BH stands between the observer and nearby light sources. The photons emitted from the source, or of the radiation emitted from an accretion flow around the event horizon of the BH, are expected to create a characteristic shadow-like image, that is a darker region over a brighter background. Indeed, recently the EHT project, a global very long baseline interferometer array observing at a wavelength of 1.3 mm, announced and showed the first image of the supermassive BH located at the centre of M87 [3], while the corresponding image from the centre of the Milky Way is yet to come. For instrumentation, data processing and calibration, physical origin of the shadow, etc, see [4][5][6][7][8].
The observation of the shadow probes the spacetime geometry in the vicinity of the horizon, and therefore it tests the existence and properties of the latter [9]. It should be noted, however, that other horizonless compact objects that possess light rings also cast shadows [10][11][12][13][14][15][16][17], and therefore the presence of a shadow does not by itself imply that the object is necessarily a BH. Therefore, strong lensing images and shadows offer us an exciting opportunity not only to detect the nature of a compact object, but also to test whether or not the gravitational field around a compact object is described by a rotating or non-rotating geometry. For a recent brief review on shadows see [18]. Within the framework of Einstein's General Relativity [19] the most general BH solution is the Kerr-Newman geometry characterized by its mass, angular momentum and electric charge, see e.g. [20]. Since, however, astrophysical BHs are expected to be electrically neutral, the most interesting cases to be considered are either the Schwarzschild [21] or the Kerr geometry [22]. More rotating BH solutions may be generated starting from nonrotating seed spacetimes applying the Newman-Janis algorithm (NJA), described in [23,24]. Non-rotating solutions have been obtained in nonstandard scenarios, such as polytropic BHs [25,26] or BHs with quintessencial energy [27], to mention just a few. Over the years the shadow of the Schwarzschild geometry was considered in [28,29], while the shadow cast by the Kerr solution was studied in [30] (see also [31]). Shadows of Kerr BHs with scalar hair and BH shadows in other frameworks have been considered in [32,33] and [34][35][36][37][38][39][40][41][42][43][44][45][46][47], respectively. To explore the physics behind the so-called BH shadows, an alternative tool is provided by the well-known NJA.
As already mentioned, this method allows one to pass from a static spherically symmetric BH solution to a rotating one.
To be more precise, in the present work we will take a variation of the usual NJA, the only difference being the omission of one of the steps of the NJA, namely the complexification of coordinates [27]. Instead of this, we will follow an "alternate" coordinate transformation, which will be explained in the next section.
In the present work we propose to investigate the shadow of the rotating polytropic BH. The non-rotating, static, spherically symmetric geometry was obtained in [25]. The metric tensor is a solution to Einstein's field equations with negative cosmological constant, and the thermodynamics of the BH precisely matches that of a polytropic gas. One one hand, as stated in [26], from the point of view of a possible astrophysical tests of the nonrotating polytropic solution, the so-called static radius, which defines the equilibrium region between gravitational attraction and dark energy repulsion, would be of importance [48]. On the other hand, in the rotating case, which will be studied in the present work by using the NJA, the computation of the shadow would constitute a valuable tool in order to confirm or refute theoretical predictions regarding the intimate structure of space and time at the strong field regime. The plan of our work is the following. After this introduction, we briefly summarize how the NJA works in the next section. In section III we study the conditions leading to unstable null trajectories for a general parametrization of a rotating BH while section IV is devoted to the study of the Hawking temperature and the emission rate of a generic 3 + 1 rotating BH. In section V we construct the rotating solution starting from a static and spherically symmetric polytropic BH as the seed geometry and we study some features of the solution as for example horizon and static conditions, causality issues, BH shadow, Hawking temperature and emission rate. Finally we conclude our work in the last section. We adopt the mostly negative metric signature (+, −, −, −), and we choose natural units where c = 1 = G.

II. NEWMANN-JANIS ALGORITHM WITHOUT COMPLEXIFICATION
In this section we review the main aspects on the NJA to generate rotating solutions introduced by M. Azreg-Aïnou in Ref. [49]. In this work, the author performed a modification in the algorithm with the purpose to avoid the complexification process in the original protocol as follows.
As usual, the starting point is a static spherically symmetric metric parametrized as The next step consists in introducing the advanced null coordinates (u, r, θ, φ) defined by from where, the non-zero components of the inverse metric can be written as and l µ l µ = m µ m µ = n µ n µ = l µ m µ = n µ m µ = 0 and l µ n µ = −m µm µ = 1. Now, introducing the complex transformation and assuming that we obtain Using the above transformations, the line element, in the so-called rotating Eddington-Finkelstein coordinates reads In order to write the metric (14) in the Boyer-Lindquist coordinates, we proceed to perform the global coordinate transformation where λ and χ must depend on r only to ensure the integrability of Eq. (15). As it is well known, in the original NJA the next step in the construction of the rotating metric, consists in complexifying r. However, in order to circumvent the complexification, Azreg-Aïnou (see Ref. [49]) proposed an ansatz for the unknown functions involved. Namely, taking where and the metric (14) takes the Kerr-like form with where a = J/M , with M, J being the mass and the rotation speed, respectively, of the black hole. At this point some comments are in order. First, note that the function Ψ(r, θ, a) remains unknown but it must satisfy the following differential equation which corresponds to imposing that the Einstein tensor satisfies G rθ = 0. Second, it can be shown (see appendix A in Ref. [49]) that the metric (22) satisfies Einstein's field equations G µν = 8πT µν with the source given by with Even more, in order to ensure the consistency of Einstein's field equations, the unknown Ψ must satisfy another constraint, namely Finally, the above expressions can be simplified in the particular case G = F and H = r 2 . Indeed, in this case it can be shown that one solution of Eq. (24) is given by and the metric (22) takes the form with ρ 2 = r 2 + a 2 cos 2 θ, It is easy to verify that when a = 0 we recover the non rotating black hole solution.

III. NULL GEODESICS AROUND THE ROTATING BLACK HOLE
In this section we implement the standard Hamilton-Jacobi formalism to separate the null geodesic equations in the rotating space-time. Our main goal here is to obtain the celestial coordinates parametrized with the radius of the unstable null orbits and to study the shadow of a generic rotating solution. Let us start with the Hamilton-Jacobi equations [50] ∂S ∂τ where τ is the proper time and S is the Jacobi action. As usual, if we assume that Eq. (37) has separable solutions, the action takes the form where E and Φ are the conserved energy and angular momentum respectively. Now, replacing (38) in (37) we obtain that where where Q is the so-called Carter constant. As it is wellknown, the unstable photon orbits in the rotating spacetime must satisfy the constraints R = 0 and R = 0 which lead to 4 a 2 − aξ + r 2 − (a − ξ) 2 + η (rF + 2F ) = 0, (44) where ξ = Φ/E and η = Q/E 2 correspond to the impact parameters. Hence, It is worth mentioning that, in the above expressions, r corresponds to the radius of the unstable null orbits. Now, the apparent shape of the shadow is obtained by using the celestial coordinates which are defined as [51] α = lim r0→∞ −r 2 0 sin θ 0 dφ dr (r 0,θ 0 ) , where (r 0 , θ 0 ) correspond to the coordinates of the observer. Finally, the after calculating the limit in the above expressions, the celestial coordinates read and the shadow corresponds to the parametric curve of α and β with r as a parameter.

IV. HAWKING TEMPERATURE AND EMISSION RATE OF A ROTATING BLACK HOLE
In Ref. [52], the author derived the Hawking temperature of a general four-dimensional rotating BH using the null-geodesic tunneling method developed by Parikh and Wilczek [53][54][55]. In this section we shall summarize the main result obtained in [52]. First of all, the radiated particles are considered as swaves because for an observer at infinity the radiation is spherically symmetric whether the BH is rotating or not. Now, the tunnelling rate of a s-wave from the inside to the outside of the BH is given by where I is the action of the tunnelling particle and Γ 0 the normalization factor. Moreover, as it is well-known, the emission rate satisfies where E is the energy of the emitted particle, β = 2π/κ and κ is the surface gravity of the horizon. Now, from (51) and (52) we arrive to from where the Hawking temperature can be derived form the standard relation In order to obtain an expression of the Hawking temperature in terms of the components of the metric we proceed to consider a generic rotating line element of the form ds 2 = −g tt dt 2 + g rr dr 2 + g θθ dr 2 + g θθ dθ 2 +g φφ dφ 2 − 2g tφ dtdφ.
It can be demonstrated that (see section II of Ref. [52] for details) where the prime stands for derivative respect to the radial coordinate, r + is the horizon radius and with Ω + the angular velocity of the event horizon, which is a constant defined by Finally, replacing (56) in (53) an using (54), the Hawing temperature takes the form which can be written alternatively as Please, notice that the Eq. (60) is a generalization of the standard formula for the Hawking temperature. In particular, when the angular velocity is taken to be zero, we recover the well-known formula, namely: In what follows, we shall study the emission rate of a rotating BH. As it is claimed in references [56,57], the BH shadow corresponds to its high energy absorption cross section for the observer located at infinity. On the other hand, it is well known that for a BH endowed with a photon sphere, the cross section has a limiting value, namely, σ lim , which coincides with the geometrical cross section of this photon sphere [58]. Now, for a spherically symmetric BH it can be demonstrated that [59], σ lim ≈ πR 2 s with R s the radius of the photon sphere defined by (see [56]) In the previous expression, the quantities α t , α r , β t correspond to particular values of the celestial coordinates (see Eqs. (49) and (50)). In particular, α r corresponds to the most right value of α and the pair (α t , β t ) stands for coordinates of the top of the shadow [60]. With all the quantities defined in this section, the emission rate [56,[61][62][63][64][65][66][67] can be calculated by where ω represents the frequency of photons. In the next section we obtain the Hawing temperature and the emission rate for a rotating polytropic BH.

V. ROTATING POLYTROPIC BLACK HOLE SOLUTION
In this section we construct the rotating BH solution from a static and spherically symmetric polytropic BH. Then we shall study some of its physical properties as well as the BH shadow and its emission. As a starting point, we consider the polytropic BH solution obtained in [25] where F = r 2 L 2 − 2M r , L 2 = −3/Λ, with Λ being the cosmological constant, and dΩ 2 = dθ 2 + sin 2 θdφ 2 is the line element of the unit two-dimensional sphere. It is worth mentioning that F in (64) has the same form of the black string solution found in [68,69] in the context of plane symmetric solutions of Einstein's equations [20]. Now, using F = G = r 2 L 2 − 2M r , the line element (32) reads with ∆ = a 2 + r 2 r 2 At this point some comments are in order. First, the horizons are solutions of ∆(r ± ) = 0 which, in this particular case, corresponds to It is worth mentioning that from the horizon condition we can obtain a bound on the spin parameter, namely a/M . In this case, the allowed values are constrained by Note that above result differs from the obtained in the case of the Kerr solution, where the constraint is given by a/M < 1. Second, the static limit, namely, the surface from where observers can remain static, corresponds to g tt = 0. More precisely, Note that, the event horizon, r + , coincides with the static radius r st at the poles θ = 0 and θ = π as in the Kerr solution. Third, note that causality violation and closed time-like curves are possible if g φφ > 0. To be more precise, the condition is from where, given that sin 2 (θ)/ρ 2 is positive, the sign of Σ plays a crucial role in the analysis. What is more, the condition to avoid causality violation and closed time-like curves is to impose Σ > 0. In the particular case of the rotating polytropic BH, the condition on Σ reads a 2 sin 2 (θ) −a 2 − r 4 L 2 + 2M r + a 2 + r 2 2 > 0, (71) It is worth mentioning that, in contrast to Kerr solution, the causality issues can be avoided for particular choices of the parameters. For example, taking a = L, the above condition reduces to which is trivially satisfied given that all the coefficients involved in the left hand side are positive. Now, we shall focus our attention in the construction of the BH shadow. Replacing the metric function F (r) = − 2M r + r 2 L 2 in Eqs. (45) and (46), the impact parameters (ξ, η) associated to the unstable null geodesics around the rotating BH are given by from where, assuming θ 0 = π/2, the celestial coordinates read In figure 1 the shadow of the rotating BH is shown for different values of a/M and L/M . Note that the shadow undergoes a deformation as a/M increases. Indeed, for a/M = 1.5 the silhouette losses its symmetry and looks like an oval instead of an ellipse. The hawking temperature for this rotating BH is given by where r + is the event horizon radius obtained from the condition (67), which reads   (77) and (62) in (63), we can study the behaviour of the emission rate as a function of the frequency ω. The emission rate profile is depicted in figure 2, where it can be shown that it decreases with an increasing value of L/M .

VI. CONCLUSIONS
In this work we have reviewed the main aspects related to the Newman-Janis algorithm without complexification, the construction of unstable null orbits and the computation of the Hawking temperature and emission rate for general rotating black holes and we have used these tools to construct a rotating polytropic black hole solution. It is worth mentioning that, to the best of our knowledge, the construction of such a rotating solution has not been considered before. Besides, we have studied some physical properties as the position of the horizons and the static limit which defines the ergosphere as well as the causality condition. We have obtained that, in contrast to the Kerr solution, the causality issues can be avoided for certain choice of the black hole parameters, namely a, M and L. Additionally, we have demonstrated that, in contrast to the Kerr black hole, the bound on the spin parameter a/M can be grater that one because the condition of the appearance of horizons entails that this bound is proportional to L/M . We have also analysed the shadow of the rotating polytropic solution finding that its shape can change while increasing L/M . As a final result, we have studied the emission rate, showing that it increases for large values of L/M . Finally, it would be interesting to investigate how the properties of the solution here obtained are modified in light of the so called scale-dependent scenario. In such cases, the coupling constants acquire a dependence on the scale with a free parameter which encodes the quantum features (which is called scale-dependent parameter). Thus, given that the black hole shadows give observational evidence for black holes, we could establish bounds on the aforementioned parameter (see [70][71][72][73][74][75][76][77][78][79][80][81][82][83][84][85]and references therein). We hope to be able to address this issue in a future work.