Thermal behavior of a radially deformed black hole spacetime

In the present article, we study the Hawking effect and the bounds on greybody factor in a spacetime with radial deformation. This deformation is expected to carry the imprint of a non-Einsteinian theory of gravity, but shares some of the important characteristics of general relativity (GR). In particular, this radial deformation will restore the asymptotic behavior, and also allows for the separation of the scalar field equation in terms of the angular and radial coordinates -- making it suitable to study the Hawking effect and greybody factors. However, the radial deformation would introduce a change in the locations of the horizon, and therefore, the temperature of the Hawking effect naturally alters. In fact, we observe that the deformation parameter has an enhancing effect on both temperature and bounds on the greybody factor, which introduces a useful distinction with the Kerr spacetime. We discuss these effects elaborately, and broadly study the thermal behavior of a radially deformed spacetime.


I. INTRODUCTION
The Kerr metric is one of the remarkable findings of GR [1]. From the weak field to strong field regime, Kerr solution has passed all tests with flying colors [2][3][4][5]. All of these success stories make a strong case of GR, and even constrained some of the alternative theories of gravity [6,7]. Moreover, with the gravitational wave (GW) astronomy coming to the fore, these studies emerge with brighter possibilities [8][9][10][11]. Besides these successes of GR, there are also strong limitations which motivates to seek for alternatives. It is also known that GR fails to explain both the small and large scale structure of the nature, and therefore, a modification of the theory is required [12,13].
If we attempt to modify GR, it is also bound to happen that the resultant spacetime may loose some of its useful properties, such as axis-symmetry or separability condition. Nonetheless, if we assure that these conditions are built-in restored, and still aim to modify spacetime structure, we may end up with constraining its metric functions [14]. The deviation from GR may be coded within these constraints. One of such possibilities comes into play if we modify ∆, which in Kerr case is ∆ = r 2 − rr s + a 2 , and r s = 2M with M being the mass of the black hole (BH). In the present article, we will be concerned with this specific example, where we modify ∆ by adding a r-dependent term to it [15]. Note that this deformation is only radial, and does not effect angular distribution of the spacetime.
The motivation to study a spacetime which mimics a radially deformed Kerr spacetime is two folded. First, it provides a simple yet useful extension of GR, which is well-grounded with GW data [16]. Therefore, it can be a potential candidate of alternative theories of gravity, and studying along this line can be beneficial. The metric corresponding to these deformed BH spacetimes has the same asymptotic features as the original ones from Einstein gravity. Secondly, the radial deformation makes it possible for the field equation to be separable in terms of the radial and angular coordinates, which paves the way for the formulation of semi-classical analysis in these spacetimes. Furthermore, the horizon structure differs, as the position of the horizon is now changed due to the introduction of the deformation. Then the effects of these deformations will also be felt through the predictions of semi-classical gravity, which concerns the horizon structure. In this regard, the Hawking effect [17] is a major arena to venture in, which states that an asymptotic observer in a BH spacetime will realize a Planckian thermal distribution of particles with temperature proportional to the surface gravity of the BH's event horizon. In the deformed BH spacetime, the distortion in the horizon structure is expected to change its surface gravity, which naturally affects the spectrum of the perceived Hawking radiation.
Another important thing to note, is that the spectrum of the Hawking effect as should practically be seen by an asymptotic observer is not an absolute blackbody distribution, rather it is a greybody distribution. This greybody distribution is characterized by the transmission coefficient through the effective potential of the considered field. Greybody factor also contains the information regarding different BH parameters. Here also one can expect prominent effects of the deformation parameter. However, a straight forward exact estimation of these greybody factors is an insuperable job analytically, though one can seek the help of numerical methods [18][19][20][21][22][23][24][25]. Analytically these estimations can be performed in asymptotic frequency regimes [19,[26][27][28][29][30][31][32][33], i.e., for very high or low frequencies of the field wave modes. On the other hand, there are methods that deal with taking extremal limit to evaluate these quantities, see [34][35][36], or analytically estimating the bounds on these greybody factors, see [37][38][39][40][41][42][43]. These bounds have the advantage of being predicted in all frequency regimes, including the intermediate frequency regimes, and also for all values of the angular momentum quantum number. In particular, we are going to consider a massless minimally coupled scalar field in the radially deformed BH spacetime, and estimate these bounds to study the spectrum of the Hawking effect with the greybody factors. Especially our motivation is understanding the changes caused by the inclusion of the radial deformation parameter in a stationary and rotating BH spacetime.
In Section (II), we begin by providing a detailed investigation of the horizon structure in the radially deformed BH spacetime. In Section (III), we consider a massless minimally coupled scalar field in this radially deformed BH spacetime, and obtain the scalar field equation of motion. Decomposition of the scalar field in terms of the spheriodal harmonics provides one with a Schrödinger wave like equation, namely the Teukolsky equation for stationary Kerr BHs, in terms of the radial tortiose coordinate. In particular, from this equation the structure of the effective potential can be perceived. Subsequently, in Section (IV), a study of the Hawking effect and the corresponding temperature and spectrum of the Hawking quanta are provided. Furthermore, in Section (V), we study the bounds on the greybody factors in these radially deformed BH spacetimes considering the effective potential from Section (III). We conclude our analysis with a discussion in Section (VI).

II. HORIZON STRUCTURE OF A RADIALLY DEFORMED SPACETIME
We should mention that there are a few well-known deformations of the Kerr metric which serve specific purposes [44][45][46][47][48]. For example, the metric provided by Johannsen and Psaltis is a Kerr-like metric that provides a regular spacetime everywhere outside of the event horizon [44]. In this case, the deformation is a function of both the radial coordinate r and the angular coordinate θ, and the spacetime is asymptotically the same as the Kerr spacetime. The field equation is not separable in this background in terms of the radial and angular coordinates. However, it is imperative to understand quantum field theory in the background black hole spacetime to realize the Hawking effect properly. That becomes possible in the spacetime provided by Konoplya and Zhidenko [15], which only has radial coordinate in deformation, enabling one to separate the field equation. This motivated us to consider this second type of description for the deformed Kerr black hole spacetime [15] to study the Hawking effect.
We start with the following Kerr metric written in a more generic form [15]: with, In the above expressions, r s = 2M , where M denotes the mass of the BH, and a is the angular momentum per unit mass. In order to inject the radial deformation, we use the substitution r s → r s + η/(r 2 ). This substitution would not change any of the built-in properties of the spacetime including the separability condition of the Klein-Gordon equation. The only difference that distinguishes the deformed spacetime from the Kerr metric of Eq. (1), is in ∆, which now becomes, ∆ = r 2 − rr s + a 2 − η/r, and also Σ changes toΣ with ∆ is replaced by∆ in the expressions of Eq. (2). Due to the presence of deformation parameter η, the locations of the horizons, given as N 2 (r, θ) = 0, would differ from the usual Kerr case. Moreover, as η is clearly coupled with r, new solutions may also emerge. To be specific, the locations of the horizons are obtained from, and we will apply Descartes' sign rule to estimate the number of solution(s). Note that with η < 0, there is no positive solution for the above case, and the naked singularity always exists. With η > 0, Eq. (3) can either have one or three positive solution(s). In the later case, the event horizon, r H and these inner horizons, r 1 and r 2 (assume r 1 < r 2 ), can be expressed in terms of the other BH parameters as where, A = β 1 + 3 √ 3 β 2 , with the expression of β 1 and β 2 given by β 1 = −9a 2 r s + 2r 3 s + 27η and β 2 = 4a 6 − a 4 r 2 s − 18a 2 r s η + 4r 3 s η + 27η 2 1/2 . We may now employ the above expressions to have a deeper understanding about the horizon structure in presence of η. Based on the properties of A, whether it is positive, negative or complex, we encounter different outcomes. In case of (a) A > 0, r H is always positive and describe the event horizon, while r 1 and r 2 are complex conjugate to each others. For (b) A < 0, say A = −α where α > 0, Eq. (4) becomes: which now readily gives r H and r 2 are complex conjugate to each other, and r 1 now becomes the outer horizon. Interestingly, if we use the η = 0 limit, and the extremality condition a = r s /2, we gather α = −r 3 s /4. The above equations now gives r H = r 2 = r s /2 and r 1 = 0, which is the usual Kerr case. This serves as an useful validation of our solution. Now we consider the last case, i.e., (c) In this case, we further simplify A as, .
(6) As it can be shown that in this case, the imaginary part of all of the above expressions would identically vanish, and r H continue to be the event horizon.
To summarize, the presence of η manifests an additional horizon other than the event and Cauchy horizon. Except for the A < 0 (note that the extremal Kerr is a special case in this limit where the imaginary parts of r H and r 2 become zero), r H continue to be the event horizon, where r 1 and r 2 are the inner horizons. For a clear exposition of this horizon structure in presence of η, we illustrate the horizon structure in Fig. (1). From this figure, we observe that for a specific η, at the same a, there can be different numbers of horizons or real roots to the solution of∆. For example, for η = 0, denoted by the solid red line and signifying zero deformation, there are two roots, and these two roots merge to a single one at the point a = 0.5 = r s /2 in the extremal case. On the other hand, when η = 0.01 or η = 0.05, denoted by the dot-dashed blue and dotted green lines respectively, there can be three, two, or even one real roots to the solution for a fixed value of a.
The other intrinsic property that is associated with a rotating spacetime is known as the ergoregion, where no observer can be kept stand still. This is given by the condition g tt = 0, which according to Eq. (1) becomes N 2 = W 2 sin 2 θ. Ergoregion is closely related with frame dragging and zero angular momentum observer (ZAMO), both are well known relativistic effects [49]. Given that ZAMO is relevant for the present purpose, we shall briefly discuss the same as follows. For the present spacetime, the angular velocity of ZAMO at a given (r, θ) is If one considers the maximally extended case of θ = 0 or θ = π for the expression of this angular velocity, then which becomes Ω H = a/(r 2 H + a 2 ), on the event horizon. Now we move on to understand the nature of Ω m as a function of r, which we will see to be appearing in the scalar field equation of motion in the next section. In order to prove that Ω m is a monotonic function in r outside the event horizon, which is often required in the context of greybody factor, we may recall the case with η = 0 first. It can be shown that in the Kerr case, Ω m has a peak at r = a/ √ 3, which always lies inside the outer event horizon and henceforth, the function is monotonically decreasing outside the horizon. Unfortunately, with the η being present, the analysis becomes more involved and a straightforward solution to r turns out to be unlikely. However, we carry out an approximate analysis to incorporate the effects of η up to the linear order terms. We note that the peak of Ω m as r varies, located at r peak , has now become which shifts close to the singularity compared to the Kerr case. Similarly, we can find that the expression of the location of outer horizon changes as: The above equation clearly states that the outer horizon shifts away from the singularity. Therefore, we conclude that the addition of η will shift the outer horizon away from singularity, and move the angular velocity's peak close to the singularity. This essentially assures that outside the outer horizon, Ω m is montonic in r. For an illustration of this incident, we plot Ω m for various values of η in Fig. (2). As it can be realized, with η = 0, For different values of η, it can be shown that the peaks of Ωm always reside inside the event horizon. Therefore, outside the horizon, angular velocity of ZAMO is a monotonically decreasing function of r. We take, a = rs/2 = 0.5. say η = 10 −2 , the maxima in Ω m is covered within the outer event horizon. Larger the value of η, the peak gets shifted close to the singularity, and the outer horizon shifts away from singularity. This indicates that for any nonzero value of η, the monotonicity of Ω m is retained outside the event horizon.

III. THE SCALAR FIELD EQUATION OF MOTION
In this deformed spacetime we consider a massless minimally coupled free scalar field Φ(x) described by the action The variation of this action with respect to the scalar field Φ provides one with the scalar filed equation of motion By substituting the metric components from Eq. (1) and Eq. (2) with the expressions of∆ andΣ corresponding to the deformed spacetime, the equation of motion becomes Like the Kerr BH here also the metric components are independent of time t and azimuthal angle φ, which suggests a field decomposition of the form Φ(t, r, θ, φ) = exp(−iωt + imφ)R(r)S(θ)/ √ r 2 + a 2 , where S(θ) denotes spheroidal harmonics. Then using this field decomposition and the tortoise coordinate r , defined from the expression the previous scalar field equation of motion from Eq. (13) can be expressed as In this deformed geometry, one may express the effective potential V(r) as (16) where, A ω lm denotes the eigenvalue corresponding to the spheroidal harmonics equation. In slow rotation limit this eigenvalue can be expressed as This Schrödinger wave like equation of Eq. (15) in deformed BH spacetime resembles the Teukolsky equation [50] from general Kerr BH spacetime. It should be noted that the from this equation the Regge-Wheeler [40,51] like equation corresponding to a deformed Schwarzschild spacetime can be obtained quite easily making a = 0. In this non-rotating limit, i.e., a = 0, the expression of potential V(r) from Eq. (16) becomes Here we have used the fact that when a = 0 the eigenvalue A ω lm of spheroidal harmonics equation signifies spherical harmonics with expression A lm = l(l + 1). In the geometrical optics limit, i.e., l 1, we gather that the above equation encounters an maxima at the photon orbit, r ph , which is a solution of the equation: With η = 0, we arrive ar r ph = 3r s /2, which is the location of photon orbit in Schwarzschild BH [52,53]. In Fig. (3), we depict the potential V(r) from Eq. (16) for different values of the parameter η. As it can be evident that there exist a maxima, r peak , in each of these figures, which shifts based on the values of η. In passing, we should also touch upon an important consequence of the above potential. It can be sensed from Fig. (3) that in either side of the peak (at r = r peak ), the function V(r) is monotonic. In particular, from the event horizon r H to r peak , it is monotonically increasing, and from r peak to ∞, it is monotonically decreasing. It is similar to saying that there will be a single peak outside the event horizon. In the presence of both a and η, the concern may arise weather this feature remains intact or not. Consequently, even with η set to zero, it is unlikely to guess from Eq. (16) that the potential will have a single peak outside the event horizon. Therefore, an analytical proof is beyond expectation. What we found that within the weak rotation approximation (aω 1) [41], this property is always valid for a wide range of various BH parameters, and different modes.

IV. HAWKING EFFECT IN RADIALLY DEFORMED SPACETIME
In the original work [17], the thermal nature of the Hawking effect is realized through the usage of Bogoliubov transformation between the ingoing and outgoing field modes described in terms of the null coordinates. The Hawking effect can be realized through other various means, like using tunnelling formalism [54][55][56][57][58][59][60], path integral approach [61], conformal symmetry [62], via anomalies [63][64][65], canonical formulation [66][67][68][69], and as an effect of near horizon local instability [70][71][72]. However, the conclusion remains the same, i.e., an asymptotic observer will perceive the BH horizon with some temperature proportional to its surface gravity. In particular, in the case of a deformed BH spacetime, the number density of the Hawking quanta perceived by an asymptotic observer with frequency ω and angular momentum quan-tum number m will be given by where, κ H , Ω H , and Γ(ω) denote the surface gravity, angular velocity at the event horizon, and the greybody factor respectively. For modes with mΩ H > ω this expression gives rise to the so called super-radiance phenomenon [73]. From the Planckian distribution of Eq. (20) the characteristic temperature corresponding to the Hawking effect is T H = κ H /2π. In a radially deformed BH spacetime (with r s → r s + η/r 2 deformation) the surface gravity at the outer horizon r = r H , which gives the temperature of the horizon, can be found out to be the evaluation of which is given in Appendix A. The inner horizons, r 1 and r 2 can exist if they are real or cease to exist if they become imaginary for certain values of η. From Eq. (4), Eq. (5) it is clear that the position of the horizon now has a signature of the deformation which will also be apparent in the spectrum of the Hawking effect. Another quantity is the angular velocity at the outer horizon from Eq. (7), which will also carry the signature of the deformation in the spectrum of Eq. (20). In our following discussion we study the nature of surface gravity and angular velocity at the outer horizon, in particular, we observe how it changes with varying η. In Fig. (4) we have plotted the surface gravity corresponding to the outer horizon r = r H in a deformed spacetime with respect to varying η. The figure is obtained considering the other BH parameters to be r s = 1 and a = 0.45. In the same figure the zero η situation is also depicted by a dash-dotted line. From this figure it can be observed that as the value of the deformation η increases the value of the horizon's surface gravity also increases, signifying an increase in the characteristic temperature of the Hawking effect. In Fig. (5), we have plotted the event horizon's surface gravity, which gives the Hawking temperature, with respect to varying angular momentum parameter a for different values of the deformation parameter η. From this figure, we note that as the value of η increases, surface gravity departs further from the Kerr case and also increases corresponding to different fixed values of a.
Similarly, in Fig. (6), we demonstrate how the angular velocity of the event horizon changes with deformation parameter. Note that Ω H is implicitly affected by η through the value of r H . We retrieve the kerr case for η = 0 limit, which is shown as the curve constant along x-axis.

V. THE BOUNDS ON THE GREYBODY FACTOR
The spectrum of the Hawking effect Eq. (20) is given by a greybody distribution as perceived by an asymptotic future observer. The greybody factor arises from the transmission amplitude of the field modes through the effective potential outside the horizon, for the modes nearly escaping the formation of the horizon and travelling from the near-horizon region to an asymptotic observer. In this section we estimate the bounds on the greybody factor, see [37,39,41], which can be analytically expressed in all frequency range, for certain general field momenta and general spacetime dimensionalities. These bounds on the greybody factor can be expressed as where Here, V(r) denotes the effective potential corresponding to a massless minimally coupled scalar field and ω corresponds to the frequency of field mode. Furthermore, h ≡ h(r ) is some positive function satisfying the condition h(−∞) = h(∞) = ω. The tortoise coordinate r is obtained from the expression of Eq. (14). From Eq. (16) one can obtain the effective potential correspond to a massless minimally coupled free scalar field in a deformed BH spacetime and express ω 2 − V(r) in another convenient form as Now one can consider the simplest choice of the positive function h(r ) = ω, as also done in [37] for Schwarzschild BHs. However, it will only be fruitful for the case of m = 0, as the first quantity from Eq. (24) contributes to a diverging term in the integration of Eq. (23) for the calculation of the bound for m = 0. We then consider the evaluation of these bounds on the greybody factor in a case by case manner as done in [41]. In particular, in separate situations for the angular momentum quantum number m = 0 and m = 0.

A. The case of m=0 :
This particular case of m = 0 is the simplest and provides an overall picture of the bounds on the greybody factor. In this case we are going to consider the positive function to be h = ω, i.e., h = 0. Now as one makes the choice of m = 0 in the expression of the scalar field effective potential from Eq. (24), one gets Then with this expression of the potential and the above choice of the positive function h one can represent the bound on the greybody factor to be where, I ω l0 = I ω lm (m = 0) and the quantity I ω lm in general is defined by the integral where the expressions for I 1 , I 2 , I 3 , and I 0 4 are given in Appendix B. With this expression of the integral I ω l0 from Eq. (28) and putting it in Eq. (26) one can find out the bound on the greybody factor in a deformed spacetime for m = 0. We have further plotted this bound in Fig. (7) for varying η with different fixed ω, and in Fig. (8) for varying ω with different fixed η. From Fig. (7) it can be observed that as the value of the deformation parameter increases the bound on the greybody factor also increases. However, it never goes beyond the upper limit 1. The bound also increases with increasing frequency ω of the wave mode, which can be observed from both Fig. (7) and Fig. (8).
B. The case of m = 0 As already observed from Fig. (2) and the discussions related to Eq. (9) and Eq. (10), we consider the angular velocity of ZAMO to be monotonic, in particular monotonic decreasing, in the entire range between the outer horizon and the asymptotic infinity. We are going to consider this property and find out the bounds on the greybody factor. Because of this feature then one can also consider the form of the positive function h in terms of Ω m (r), as done in [41], which can be made monotonic in the region between the event horizon and the asymptotic infinity. With this consideration from Eq. (22) the bound on the greybody factor becomes where, we have sought the help of the triangle inequality to express the numerator of Eq. (23). This expression further simplifies to We consider Ω H ≡ Ω(r H ) to be the angular velocity at the event horizon. In regard to the monotonicity of Ω m (r), it should be noted that Ω m (r) is always smaller than Ω H in the region outside the outer horizon, which can be used for the monotonicity of h. Furthermore, this monotonicity can be achieved in two different regions, namely in the non super-radiant regimes ω > mΩ H or m < m with m = ω/Ω H , and for the case of the superradiant modes of m > m . In the following discussion we shall consider the case with the non super-radiant and super-radiant modes separately.

Non super-radiant modes m < m
First we consider the non super-radiant modes, i.e., m < m , where we can have another two possibilities either m < 0 or m ∈ (0, m ). The function h in both cases are considered to be h = ω − mΩ m (r), which is positive. Let us check for the first quantity of Eq. (30) in these two scenarios. In particular we have On the other hand, using the expression of (ω 2 − V(r)) from Eq. (24) the second quantity of Eq. (30) becomes where, the expression U(r) = U(r)(r 2 + a 2 )/∆, and U(r) is obtained from Eq. (24). Then for the scenario m < 0 of the non super-radiant modes m < m , the bound on the greybody factor from Eq. (30) can be obtained using the outcomes of Eq. (31) and Eq. (32) as On the other hand, for the scenario m ∈ (0, m ) of the non super-radiant modes the bound on the greybody factor from Eq. (30) is where, the expression of I ω lm is given in Eq. (27). The explicit evaluation of this quantity I ω lm and the other quantities I 2 , I 3 3 , and I 6 can be obtained from the evaluated integrals of Appendix B(in particular from Eq. (B.1)). Now it can be noticed that one can get the bound on the greybody factor for m = 0 case(the expression obtained by putting the result of Eq. (28) in Eq. (26)) from these bounds of Eq. (33) and Eq. (34) by simply making m = 0 in these expressions. Here the contribution of non zero m comes through specific two quantities in the bound, and they also contain the effect of the deformation parameter η for m = 0. In Fig. (9) we have plotted this bound with respect to varying η for fixed parameters l = 2, m = −1, a = 0.05 and different values of the frequency ω. Note here we have kept the quantity aω 1 so that the approximation of Eq. (17) can be made. From this figure also one can observe that the bound on the greybody factor increases with an increasing value of the deformation parameter and mode frequency.

super-radiant modes m ≥ m
The super-radiant modes are those for which mΩ H > ω or m ≥ m . In this scenario the expression of the bound on the greybody factor from Eq. (30) is further simplified using the triangle inequality and is expressed as where, U(r) is obtained from Eq. (24). One can carry out this integral by considering two regions m ∈ [m , 2m ) and m ∈ [2m , ∞), where the value of h can be taken to be positive and with a reasonable boundary condition, see [41]. Let us discuss these two situations in a case by case manner.
Case I (m ∈ [m , 2m )): In this case we consider the the function to be h(r) = max {ω − mΩ m , mΩ H − ω}, which satisfies the requirement for h(r) to be positive and also its asymptotic behaviors, see [41]. Then the first quantity of the right hand side of Eq. (35) can be evaluated to be The second quantity in the same equation becomes and the third quantity where, r 0 is obtained from equation ω − mΩ m (r 0 ) = mΩ H − ω.
Case II (m ∈ [2m , ∞)): In this case the function can be chosen to be h(r) = max {mΩ m − ω, ω}, which satisfies the requirement for h(r) to be positive and also its asymptotic behaviors, see [41]. Then the first quantity of the right hand side of Eq. (35) can be evaluated to be The second quantity in the same equation becomes and the third quantity where, r 0 is obtained from equation mΩ m (r 0 ) − ω = ω.
In these two cases also major contributions from Eq. (37) and Eq. (40), in the bound on the greybody factor, are analogous to the results obtained from non super-radiant and m = 0 cases. However, here the effects of deformation also come from the first and third quantities of the integral in Eq. (35). One may choose a suitable range of parameter values and depict this case also in a figure by plotting the lower bound on Γ(ω) with respect to η. However, that plot does not impart any new information.
C. The bound on the greybody factor when a = 0 The case of a = 0 is particularly significant as it denotes the radially deformed static Schwarzschild BH spacetime. In this case one does not need to consider the intricacies coming from the super-radiant frequency, and the contribution from different angular momentum quantum number m becomes irrelevant in the calculation of the bound on the greybody factor. Then one can consider the expressions provided for m = 0 in Eq. (26) with Eq. (27) for the evaluation of these bounds.
In particular, the expression of Eq. (27) in a = 0 case can be evaluated to provide where, we have used the fact that in the a = 0 case A ω lm denotes the eigenvalue corresponding to the Spherical harmonics equation, i.e., A ω lm | a=0 = l(l + 1), which can also be seen from Eq. (17). Then the bound on the greybody factor can be estimated using Eq. (26) as This ensures that at very high frequency the bound has the form in agreement with the result from Born approximation [74]. In Fig. (10) we have plotted the bound of Eq. (43) and observed that in this radially deformed Schwarzschild BH spacetime also the value of the bound increases with increasing value of the deformation parameter and the frequency of the wave mode.

VI. CONCLUSION
In this article, we have discussed the thermal behavior of a BH spacetime which is deformed from the Kerr solution. We have mentioned that this deformation is only radial and therefore, the separability of scalar field wave equation and asymptotic properties remain the same as the Kerr spacetime. However, the position of the horizons and near horizon geometry differs from the Kerr case. This, in fact, motivated us to pursue the analysis concerning the Hawking radiation and bounds on the greybody factors in this spacetime.
The addition of the deformation parameter η makes the horizon equation a cubic one, and there can be either one or three real positive solution(s), as seen in Section (II). In this deformed BH spacetime, we observed that other than the event horizon, i.e., the outer horizon provided by the largest positive real root, there are other two inner horizons. This would make sure that the singularity is never naked, as it is covered by at least one horizon all the time. This is in stark contrast with Kerr BH case, where the extremality condition exists, and no real positive root of the horizon equation can be found for a > r s /2.
The main objective of this article was to highlight how a radial deformation η could change the thermal behavior of a BH solution. In this regard, from Fig. (3), we observe that the height of the effective scalar field potential decreases with increasing η. It indicates more transmission of the field modes through the effective potential with increasing η, thus leading to the anticipation of an enhancement in the greybody factor. In Fig. (4) and Fig. (5), the event horizon's surface gravity of this deformed BH spacetime is depicted for varying η, which imply an increasing Hawking temperature with increasing deformation parameter. In order to describe the bounds on the greybody factors in presence of η, which we discussed in Section (V), we closely followed the analytical treatment as given in [37,39,41]. Note that these bounds are the lower bounds, while the upper bound is always 1.
In the case where bound becomes 1, all the modes can escape to infinity unaltered, while the case of the lower bound nearing zero signify that much smaller amount of the wave modes transmit to the infinity. In Section (A) of Section (V), we considered the m = 0 case, and the bounds shown in Fig. (7) and Fig. (8) imply an increase in the greybody factor for increasing deformation parameter η, which is consistent with the results obtained from Fig. (3). From these figures one also perceives that the greybody factor increases with increasing wave mode frequency ω. We observed in m = 0 case, a significant contribution in the bound coming from an integral Eq. (27) similar to the m = 0 case, which keeps the dependence of the bounds on η also same here, see Fig. (9). Therefore, it is apparent that the radial deformation enhances the transmission probability for the modes to travel to infinity, which is also evident from Fig. (3). From this figure, we observed the height of the effective potential V(r) decreases with increasing the deformation parameter η, signifying a greater possibility for the field modes to transmit through the potential barrier. Moreover, from Eq. (20), it can be observed that, at least for the m = 0 case when the horizon's angular velocity does not contribute, the number density of Hawking quanta corresponding to a certain wave mode of frequency ω in a deformed BH spacetime will be higher compared to the Kerr case. Finally, in the last part of Section (V), we found that for the Schwarzschild case also the bound on the greybody factor increases with increasing η.
In passing, we would like to remind that a more general deformation can be found if we use r s → r s +η/r n+1 with n ≥ 1, which also provides a spacetime where the field equation of motion is separable with respect to the radial and angular coordinates. Note that based on the value of n, different horizon structure emerges, and we end up with different number of horizons. However, as far as the thermal behavior is concerned, we are doubtful how much impact does n can cause on the overall numbers. This is because as n increases, the spacetime becomes more Kerr-like, and η loses its essence. Therefore, the dominant contribution from η would only come in the n = 1 case, while for all n > 1, the contribution from the deformation becomes dimmer. However, it remains an interesting arena to venture further.