Remarks on the black hole shadows in Kerr-de Sitter space times

This work is geared towards analysis of shadows cast by Kerr-de Sitter (KdS) and Kerr-de Sitter Revisited (RKdS) black holes. Considering observers in the vicinity of the static radius, we derive the impact parameters defining the apparent positions of the shadows. Such observers are of interest to our work because embedding diagrams have shown that de Sitter space-time is analogous to an asymptotically flat one in the vicinity of the static radius. We also perform a comparative analysis between our result with that in Ref.[1]. Furthermore, we numerically obtain the radii of curvature, vertical diameters and horizontal diameters of the shadows. We find that for $\Lambda=1.11\times 10^{-52} m^{-2}$, M87* observations cannot distinguish a RKdS black hole shadow from that of a Kerr black hole. Additionally, for the same value of $\Lambda$, KdS and RKdS black hole shadows are identical. Previously, it has also been shown that when $\Lambda=1.11\times 10^{-52} m^{-2}$, KdS and Kerr black hole shadows are indistinguishable. Utilizing the 2017 EHT observations of M87* on the allowed range of the characteristic radius of the shadow, we obtain constraints on both black holes. When, a/M>0.812311, we observe that large angles of inclination ($\theta>30.5107^{\circ}$) do not pass the constraints for both KdS and RKdS black holes.


I. INTRODUCTION
In the General Theory of Relativity, the gravitational force is seen as the effect of the deformation of space-time. Black holes are regions of very strong gravity that are sufficient to warp space, bend light and give rise to space-time singularities. Since black holes results from the supernovae explosions that leave behind a positive angular momentum of the remaining matter, it is most probable that all black holes in nature are rotating and are therefore described by the Kerr solution [2].
In the presence of a cosmological constant Λ > 0, an expected component of the universe due to the dark energy phenomena, a generalization of the Kerr metric is given by the KdS metric, a solution that was firstly found by Carter [3]. The KdS metric belongs to the Plebanski-Demianski family of solutions [4]. These solutions describe the most general stationary, axially symmetric Petrov Type-D metrics of Einstein-Maxwell equations with a cosmological constant. The Plebanski-Demianski solutions are characterized by seven parameters that in certain instances are related to acceleration, magnetic and electric charges, mass, cosmological constant, NUT parameter and angular momentum, [5].
Recently, Ref. [1] , through gravitational decoupling for axially symmetric systems ( [6], [7], [8]), have proposed a revisited solution, the RKdS metric. RKdS solution is a rotating version of the Schwarzschild-de Sitter (SdS) solution, representing a black hole with a cosmological constant. The solution is asymptotically de Sitter and reduces to Kerr solution as a special case. Unlike the KdS solution, the RKdS solution is neither a Λ−Vacuum solution -it does not belong to the Plebanski-Demianski class of metrics -nor is it a constant curvature solution (i.e it exhibits warped curvature except on the equatorial plane where curvature remains constant).
As a consequence of the cosmological constant, there arises a cosmological horizon behind which the geometry of space-time is dynamic. For astrophysical processes, another radius associated with cosmic repulsion is relevant [9], the so-called static radius. On the static radius boundary, gravitational attraction due to the central compact object and cosmic repulsion counterbalance each other. This radius represents a natural boundary for gravitationally bound systems in an expanding universe [10]. A cosmological horizon is relevant for cosmology and not for astrophysical processes that are limited by the static radius [10].
When a black hole is in front of a luminous background, its unstable photon region, a region containing null geodesics at a constant radius, will be projected on the observer's sky to form the so-called black hole shadow. The shadow appears as a dark disc with a bright ring around it. Static black holes cast circular shadows owing to their spherical symmetry. Rotating black holes on the other hand cast elongated shadows as a consequence of frame dragging effects. The field of studying black hole shadows as a feasible and reliable observational probe is still in the first days of life but the direct detection of gravitational waves by Laser Interferometer Gravitational-Wave Observatory (LIGO) and Virgo collaboration [11] and the images of the super-massive black hole at the centre of Messier 87 galaxy by the Event Horizon Telescope (EHT) [12]- [13] have been advances of great significance towards observation of black holes.
In this work, we investigate shadows cast by black holes in KdS and RKdS space-times. Our aim is to promote a comparative analysis between various results presented in the literature. The motivation for our study is as follows. In space-times that are non asymptotically flat, the calculation of the angular radius of the shadow solely requires the use of finite distance observers. For such observers, a choice of an orthonormal tetrad is necessary, as shown by [14], [15] [16] and [17]. In asymptotically flat space-times on the other hand, the angular radius of the shadow can be analogous to that in Minkowski space-time [18]. As a result, in asymptotically flat space-times, one calculates the parameters associated to the black hole shadows via the expressions where α and β are the impact parameters defining apparent positions of the image of the black hole on the celestial   sphere [19]. Actually α is the apparent displacement of the image perpendicular to the projected axis of symmetry of the black hole while β is the apparent displacement of the image parallel to the projected axis of symmetry [19].
(r o , θ o ) are the positions of the observer. dφ dr and dθ dr are calculated from the null geodesic equations for the respective space-time. Properties of the geometry around black holes can be visualized by embedding their 2D sections of t =constant hypersurfaces onto 3D Euclidean geometry. The resulting diagrams are known as embedding diagrams. Embedding diagrams expedite in acquiring an intuitive understanding of the gravitational field rendered into curved space-time. Besides, these diagrams aid in analyzing the effect of parameters such as magnetic field, electric charge, black hole spin or cosmological constant [20]. By the use of embedding diagrams of both ordinary geometry and optical reference geometry, Refs. [21] and [22] have shown that in the vicinity of the static radius, the geometry of de Sitter space-time is analogous to an asymptotically flat space-time. Thus, eq. (1) and (2) can therefore be applied in the vicinity of the static radius where the KdS (or SdS) space-time is close enough to the asymptotically flat space-time case. This approach has been considered in Ref. [1], where the size of a shadow in KdS and RKdS space-times has been compared. We however note that there is a discrepancy on the analytical solution for the form of α and β they use with the results we have obtained in our work. In Ref. [1], attention was not given to the analysis of the black hole shadows; they were proposing the RKdS solution for the first time and basically used the shadow to support their result that a KdS black hole has a larger event horizon than the RKdS black hole. As a result, the error does not affect their result. To our knowledge, there has been no work in the literature that has given a detailed analysis of the KdS and RKdS black hole shadows for observers in the vicinity of the static radius. Moreover, RKdS solution has just been recently proposed and hence has no detailed study of it's black hole shadows alongside a comparison with those of KdS. For astrophysical purposes, a detailed analysis of the black hole shadows is necessary and this is what we intend to do in this work. In our work by considering different black hole parameters, we have shown that the behaviour of the black hole shadows, for observers in the vicinity of the static radius, is consistent with the behaviour of the corresponding spherical photon orbits in both KdS and RKdS space-time. We hope that from this work, more aspects of research can be considered in the vicinity of the static radius because a lot of study in KdS space-time has only been geared towards this space-time being asymptotically de sitter. It is worth noting that black hole shadows can also be studied through other techniques such as general-relativistic radiative transfer and general-relativistic magnetohydrodynamic simulations ( [23], [24], [25], [26]). The work is organized as follows. In sections II and III we give a summary of KdS and RKdS space-times. We provide analytic solutions to their horizons and radii of equatorial circular photon orbits. We further analyze the behaviour of the radii of the equatorial circular photon orbits in each case. In section IV, we obtain the celestial coordinates for the shadow cast by a KdS and RKdS black hole. We analyze the shadows for different values of black hole spin -the parameter a -and cosmological constant Λ. In section V, we numerically evaluate the curvature radii of the shadows at specific characteristic points. We further obtain numerical values for the horizontal and vertical diameters of the shadows. In section VI, we use the radius of curvature to constrain a KdS and RKdS black hole. Lastly, in section VII, we give a conclusion for our results. The metric signature in this work is −+++ and G = c = 1.

II. KERR-DE SITTER SPACE TIME
KdS metric is a limited case of the Plebanski-Demianski solution with zero acceleration, NUT parameter, electric and magnetic charges. In Boyer Lindquist coordinates, KdS space time is characterized by a metric of the form [3], where the terms appearing in the metric coefficients are defined as, Coordinates t and r range over all R while θ ∈ [0, π] and φ ∈ [0, 2π]. M is the total mass of the system, a the angular momentum per unit mass and Λ is the cosmological constant. The coefficients of metric eq. (3) are independent of t and φ thus ∂ t and ∂ φ are Killing vector fields. Any linear combination of these two Killing vector fields will also be a Killing vector, [27].
• The region where g tt > 0 is the ergoregion. In the ergoregion, ∂ t is space-like. The boundary of the ergoregion occurs at g tt = 0 which is known as the static limit. The static limit defines how close to the black hole static observers (observers for which spatial coordinates along their world lines do not change with time) can get. In the ergoregion, trajectories for which the observer remains at fixed values of Boyer-Lindquist coordinates are not possible. Observers in this region remain at fixed r and θ by rotating in φ direction.
• For g φφ < 0, ∂ φ becomes time-like. This will allow the presence of closed time-like curves whose existence violates causality. Observers moving on closed time-like curves find themselves in their own past. Closed time-like curves were first discovered by Kurt Godel [26].
Moreover, the metric has coordinate singularities at ∆ r = 0, θ = 0, π and a curvature singularity at Σ = 0 (i.e r = 0 and θ = π/2). Kerr-de Sitter horizons are determined by the roots of ∆ r , thus, Equation 8 is a quartic polynomial hence admits four roots r −− , r h− , r h+ and r c , with 0, 2 or 4 of the roots being real. r h− and r h+ are the Cauchy and Event horizons respectively while r c is the cosmological horizon and r −− is interpreted as the dual of the cosmological horizon. Applying Ferrari's solution for quartic polynomials, we obtain the relations for the horizons as, where we have defined, The condition that horizons (9) - (12) have to obey to be real is ω > 0. Thus, ω > 0 is the condition for a regular KdS black hole.

A. Photon Region
Due to the symmetries of KdS space-time, the trajectories possesses the conserved quantities, E is interpreted as the energy of the particles per unit mass and is related to the space-time being stationary while Φ is the angular momentum per unit mass in the z direction and is related to axial symmetry. There exist a third conserved quantity related to the hidden symmetry of KdS space-time which is the Carter's constant Q [28]. Q comes up as a result of separation of variables in the Hamilton-Jacobi equation. For photon orbits, only the sign of energy has a physical meaning. Hence the conserved quantities can be rescaled in terms of energy as, Thus the equations describing motion of photon orbits in KdS space-times are given by the Carter's equations [3], In the above expressions we have defined, with usual definitions p µ = dx µ /dσ, where σ is the affine parameter along the geodesic. An important class of photon orbits pivotal to the formation of a black hole shadow are the spherical photon orbits. For such orbits, the conditions, must be satisfied. Solving both conditions in eq. (24) simultaneously yields, In the case a = 0, eq. (24) reduces to SdS. Due to spherical symmetry of SdS space-time, orbits will be planar hence we can choose the plane θ = π/2. On this plane, η = 0. Thus substituting for (a = 0, η = 0) in eq. (24), we obtain the relation for λ in SdS space-time as [29], Setting Λ = 0 in eq. (25) and (26) yields parameters obtained in Kerr space-time [30], Equation 25 and (26) are the relations for the constants of motion governing photon orbits at r = constant. The radii of these orbits are bound by the equatorial circular photon orbits which are obtained by the condition η = 0. We note that solving for η = 0 is a cubic polynomial hence results in three roots. However, we will give the two roots that are relevant for this work (roots located outside the event horizon). By using Cardano's formula, the roots read, The parameter κ is defined via the following expression, Thus the photon region in KdS black hole exterior exist in the range r ∈ [r ph+,KdS , r ph−,KdS ]. Photon orbits at r =constant are unstable with respect to radial perturbations if d 2 R/dr 2 > 0 and stable if d 2 R/dr 2 < 0. Figure 2 illustrates the behaviour of the radius of equatorial circular photon orbits, eq. (30) and (31).

III. KERR-DE SITTER REVISITED SPACE-TIME
The RKdS solution has been recently proposed by Ovalle, Contreras & Stuchlik in Ref. [1] and is defined by the metric, with the following definitions, In this solution, Λ is interpreted as a vacuum energy that suffers distortions in the vicinity of the black hole due to rotation. Coordinates r, t ∈ R while θ ∈ [0, π] and φ ∈ [0, 2π). The metric is stationary and axially symmetric (i.e ∂ t and ∂ φ are Killing vector fields). The geometric interpretation for metric components g tt and g φφ revised in section II can be immediately applied now. From the condition, making use of Ferrari's solution we obtain the horizons as of the RKdS metric such that, The above results made use of the definitions, Whereasr h+ and r h− are the event horizon and Cauchy horizon respectively,r c is the cosmological horizon andr −− is interpreted as the dual of the cosmological horizon. These roots will be real provided the discriminant of eq. (37) is positive, Thus eq. (43) is the condition for the existence of a regular RKdS black hole.

A. Photon Region
Due to RKdS solution being stationary and axially symmetric, we have the conserved quantities, Where E Λ is the energy per unit mass and Φ Λ is the angular momentum in the z direction. Despite RKdS solution being a non-vacuum solution, we find that the Hamilton-Jacobi equation is separable and thus this space-time possesses the Carter's constant Q Λ , For photons, we can rescale the conserved quantities in terms of energy as, Using the Hamilton-Jacobi formalism and making use of these symmetries, we obtain the equations for null geodesics as, with the following definitions, Also, it is worth noting that we proceed as before and by simultaneously solving R Λ (r) = R Λ (r) = 0, which results in the following relations, The RKdS space-time reduces to SdS in the limit a = 0. Substituting for (a = 0 η Λ = 0) in R Λ (r) = R Λ (r) = 0 and solving for λ then gives [29], η Λ = 0 since we chose the plane θ = π/2 as a result of SdS space-time being spherically symmetric. For Λ = 0, eq. (53) and (54) yields equations same as those obtained in Kerr space-time [30], Equation 53 and (54) are the relations for constants governing motion of photon orbits at constant r in RKdS space-time. The zeros of eq. (53), yields the radii of equatorial circular photon orbits, κ is defined via the expression, Thus the radii of photon orbits with r =constant is bound by r ∈ [r ph+,RKdS , r ph−,RKdS ]. These orbits are unstable to radial perturbations when R Λ > 0 and stable if R Λ < 0. Figure 3 illustrates the behaviour of the radius of equatorial circular photon orbits, eq. (58) and (59). Up to this point we have extensively compared the KdS and RKdS space times from the qualitative point of view. Henceforth we investigate quantitatively the predictions for the black hole shadows in both space-times and compare them.

IV. KERR-DE SITTER AND KERR-DE SITTER REVISITED BLACK HOLE SHADOWS
Utilizing the null geodesics of KdS space-time, we obtain the corresponding impact parameters as, while from the null geodesics of RKdS space-time we obtain, We however note that our results eq. (61) and (62), eq. (63) and (64), are different from those used for in Ref. [1], where in our observation, for both KdS and RKdS, they consider impact parameters of the form [31], Thus, for the apparent displacement of the image perpendicular to the projected axis of symmetry of the black hole, we obtain the impact parameter as in eq. (61) and eq. (63) for KdS and RKdS respectively. These equations are clearly different from eq. (65). Moreover, for the apparent displacement of the image parallel to the projected axis of symmetry of the black hole, we obtain the impact parameter as in eq. (62) and eq. (64) for KdS and RKdS respectively. These equations are also different from eq. (66). This difference can also be observed as illustrated in fig. 4. In fig. 4, for the same values of black hole parameters, our results do not coincide with Ref. [1] in both KdS and RKdS black holes. The error in the impact parameters considered in Ref. [1] is due to the fact that they directly considered the impact parameters of Kerr space-time for which they cited Ref. [31]. Physically, one should utilize the null geodesic equations of the space-time under consideration in order to obtain the impact parameters. Thus, the impact parameters in Ref. [1] have been obtained using null geodesic equations of Kerr space-time and this leads to inconsistent black hole shadow results when applied to KdS and RKdS space-times. We however note that this error is not relevant to Ref. [1] as they show a shadow to complement their work, without this being the relevant part of their paper. For the purposes of our work and any other future studies this error is relevant. Hence, our impact parameters yield consistent black hole shadow results, as will be seen in the next sections, because we derived them using the null geodesic equations of KdS and RKdS space-time respectively. The boundary of the shadow is formed by spherical photon orbits that are unstable to radial perturbations. For a KdS black hole we found that such orbits lie in the range r ph+,KdS ≤ r ≤ r ph−,KdS while r ph+,RKdS ≤ r ≤ r ph−,RKdS for a RKdS black hole. As a result, the impact parameters will be evaluated for values of r in these ranges. Figure 1b demonstrates a general contour of a black hole shadow. The contour is a parametric curve (α, β). Every point along this curve is associated to a particular radial point in the photon region. When the observer is inclined on the equatorial plane r D = r ph+,KdS (r ph+,RKdS ) and r R = r ph−,KdS (r ph−,RKdS ), which are the radii of equatorial circular photon orbits. α increases from a negative value at r D to a positive value at r R . β on the other hand vanishes at D and R. Thus when an observer is inclined on the equatorial plane, the behaviour of the shadows at point D and R should correspond to the behaviour of the radii of equatorial circular photon orbits that we obtained in fig. 2 and  fig. 3. In other words, • For fixed Λ, r ph+,KdS (r ph+,RKdS ) both decrease with increase in black hole spin. This implies that as black hole spin increases, point D on the curve of the shadow will decrease by moving to the right. r ph−,KdS (r ph−,RKdS ) increases with increase in black hole spin, thus point R increases by moving further to the right. As a consequence, when the black hole spin increases, the shadows will appear to shift to the right.
• For fixed black hole spin, r ph+,KdS (r ph+,RKdS ) increases with increase in Λ. In this case, point D increases by moving further to the left. r ph−,KdS (r ph−,RKdS ) decreases with increasing black hole spin. Point R consequently decreases by moving to the left. Hence, with increasing Λ, the shadows undergo a shift to the left.
This behaviour maintains even when the observer is inclined away from the equatorial plane. In the subsequent section, we will analyze the qualitative behaviour of the shadows using this concept so as to identify consistencies or inconsistencies.

A. Plots of Black Hole shadows
In fig. 4, (5) and (7), we have plotted the shadows using a Λ = 0.06 and M = 1.06 as considered in Ref. [1] (figure 2 of their work). For different values of the black hole spin, a, we notice that our results qualitatively agree with Ref. [1] i.e, as the black hole spin increases, the left side of the shadows moves to the right which corresponds to the radius of equatorial circular prograde photon orbit decreasing with increase in black hole spin. The right side moves further to the right which is consistent with the radius of equatorial circular retrograde photon orbit increasing with an increase in black hole spin. Nonetheless, it is clear that there exists a quantitative discrepancy. The celestial coordinates used in Ref. [1] yield shadows that appear elliptic. On the other hand, our results yield shadows that do not deviate so much from circularity.
From fig. 6, we observe that for increasing Λ, our obtained impact parameters yield shadows that undergo a shift to the left. This is consistent with the discussion that we have given in the previous section. The impact parameters used in Ref. [1], gives rise to shadows inconsistent with the behaviour of equatorial circular photon orbits, i.e the right side of the shadows moves further to the right as Λ increases. This would imply that r ph−,KdS (r ph−,RKdS ) increases with increase in Λ which is not the case as observed in the previous section. In-addition, there is still a quantitative discrepancy. Figure 7 shows that a KdS black hole will cast a larger shadow compared to a RKdS black hole. This observation, as in Ref. [1], has been attributed to the fact that a KdS black hole has a larger event horizon compared to that of a RKdS. When an observer is inclined away from the equatorial plane as illustrated in fig. 8, the behaviour of the shadows is still consistent with that observed in fig. 5. For high values of black hole spin, the shadows will have a near circular shape when the observer is inclined away from the equatorial plane. However, as the angle of inclination moves towards π/2, the left side of the shadow begin to appear more flattened as seen in fig. 9. Thus, a change in the angle of inclination will change the quantitative aspect of the shadows i.e, the horizontal and vertical diameters of the shadows will change with the change in angle of inclination. The analytic behaviour with respect to changes in a and Λ will remain the same.

V. INTRINSIC CURVATURE OF THE SHADOWS
The black hole shadow is a one dimensional closed curve parametrized by (α(r), β(r)). Making use of differential geometry concepts of closed curves, the intrinsic properties of the black hole shadow can be obtained. In this section, we will employ differential geometry to determine the curvature radius of the shadows. The concept of curvature radius was introduced in Ref. [32] together with its applications in constraining the black hole's parameters. Moreover, by utilizing the symmetry of a black hole shadow, Ref. [33] have obtained characteristic points along the one dimensional curve of the shadow. Using the characteristic points and curvature radius, the Kerr black hole parameters have been obtained. An illustration of the characteristic points is as shown in fig. 1b. We will obtain the curvature radius and characteristic points of a KdS and RKdS black hole shadow. Using these properties, we will do a quantitative comparison of these shadows. We will further compare this quantitative behaviour to that of a Kerr black hole [33].
Utilizing the concept of parametric curves, the radius of curvature can be obtained by, We then compute the curvature radius at points T ,D,R and B as shown in fig. 1b and denote them as R T , R D , R R , R B respectively. To avoid writing the cumbersome equations that we obtain, we give the general approach by which we calculate the points, 1. To measure the curvature radius at point T and B, we work with a value of r obtained through ∂ r β = 0.
2. In evaluating curvature radius at D and R, we use r for which β = 0.

Horizontal diameter, ∆α = (α
We have numerically computed these points for a KdS and RKdS black hole shadow. Additionally, with the aid of the value of gravitational radius obtained from Event Horizon Telescope (EHT) observations [34], [33], we model M87* to these shadows. The results are in table I and (II) for two different values of Λ. We have chosen Λ = 1.11 × 10 −52 m −2 because cosmological tests and CMB measurements imply that this is the relevant value of the cosmological constant [35], [36], [37] In table I, the second column contains values that we obtained by modelling M87* to a Kerr black hole shadow. The values agree with those obtained by [33]. We however note that, in [33], the values of R T and R R given in table 2 seem to have been interchanged. In the third and fourth column, we have the KdS and RKdS values using the celestial coordinates obtained in our calculation, eq. (61)-(62) and eq. (63)-(64) respectively. The fifth and sixth column has   second and third column with the corresponding column in table I (third and forth column), we observe that the radii of curvature and diameters decrease when a larger value of Λ is used. Thus, in both KdS and RKdS space-times the shadow will decrease in size for a larger value of Λ, in accordance with the celestial coordinates we obtained, eq. (61)-(62) and eq. (63)-(64). By using orthonormal tetrads approach, Ref. [14] results also show that the size of a KdS shadow would decrease when Λ is increased. However, comparing the forth and fifth column with the corresponding column in table I (fifth and sixth column), the radii of curvature and diameters increase with increase in the value of Λ. From table II, the values of curvature radii, horizontal and vertical diameters of a KdS black hole shadow are larger than those of a RKdS shadow despite the difference being small. This is in agreement with fig. 7 where we observed that a KdS black hole casts a larger shadow than a RKdS black hole. Thus, the RKdS solution shows a gravitational field less intense than the standard KdS solution.
It is worth to note that besides the approach of intrinsic curvature, there exist other alternative ways of comparing black hole shadows as illustrated in Ref. [38] and [39].

VI. CONSTRAINTS ON KERR-DE SITTER AND KERR-DE SITTER REVISITED BLACK HOLE.
Recently, in Ref. [40], the 2017 EHT observations of M87* were utilized and a constraint on the characteristic radius of the shadow was obtained at 68% confidence level. It was shown that irrespective of the underlying solution being spherically or axially symmetric, the radius of the shadow ought to lie in the range, 4.31M ≈ r sh,EHT −min ≤r sh , r sh,A ≤ r sh,EHT −max ≈ 6.08M, wherer sh and r sh,A denotes radius of the shadow in spherically symmetric and axially symmetric solution respectively. From this shadow size, highly charged dilaton black holes were ruled out of M87* observations. Additionally, in Ref. [41] this observation was used to constrain the tidal charge of a Reissner-Nordstrom black hole. In this section we will implement the concept of curvature radius that we have discussed in the previous section to this observation. In so doing we will constrain the parameters of a KdS and RKdS black hole. In Ref. [42], the concept of curvature radius has been used to constrain a Kerr and a Kerr-Newman black hole. While excluding Z 2 symmetry, the radius of curvature will have one maximum and one minimum [32]. Figure 10 illustrates the maximum and minimum points of the curvature radius of a KdS and RKdS black hole shadow. From this figure, the radius of curvature has a local maximum at r = r A and r = r B . Furthermore, as the black hole spin and θ increases, these two local maxima also increase. The local maximum at r = r A is greater than that at r = r B . Hence R max = R(r A ). Additionally, the minimum of the radius of curvature (R min ) forms at the well of these curves. This minimum point decreases with increase in spin and θ. R min and R max then gives a lower and an upper bound (a) Kerr Practically, there are no sensible differences in constraints between Λ = 0 ( [43], [42]) and Λ ∼ 10 −52 , this section thus confirms the findings of [37]. The shaded parts indicate regions where the minimum radius of curvature is less than 4.31M and the maximum curvature is greater than 6.08M . Such regions according to eq. (70) are excluded from observations by M87* hence imposing constraints on θ and a. The unshaded part is the allowed region.

VII. CONCLUSION
We have analyzed black hole shadows in KdS and RKdS space-times for observers located in the vicinity of the static radius. Embedding diagrams of both ordinary geometry and optical reference geometry have shown that the space-time in the proximity of the static radius in de Sitter space-time is analogous to an asymptotically flat spacetime. This makes such observers to be of great interest to our work. In view of the fact that a black hole shadow forms as a result of the unstable photon region being projected on the observer's sky, we have first investigated the behaviour of photon orbits in the respective photon regions.
Particularly, we have investigated the behaviour of the radius of equatorial circular photon orbits which form the lower and upper bound of the photon region. In both space-times, the radius of equatorial prograde (retrograde) circular photon orbits decreases (increases) with increase in black hole spin. For increasing value of Λ, the radius of equatorial prograde (retrograde) circular photon orbits increases (decreases). We have done a qualitative analysis of the shadows with respect to this behaviour. We find that the shadows reflect this observed behaviour of the radius of equatorial circular photon orbits. Besides, we have compared our result with Ref. [1]. For varying black hole spin, our results agree qualitatively. However, for varying Λ, the impact parameters considered in Ref. [1] yield shadows whose behaviour contradicts the behaviour of the radius of equatorial circular photon orbits.
Furthermore, we have numerically computed the radii of curvature at specific characteristic points along the curve of a KdS and RKdS black hole shadow. We have further obtained numerical values for the horizontal and vertical diameters of the shadows. These values have then been modelled to M87* observations. For Λ = 1.11 × 10 −52 m −2 , we observe that KdS and RKdS black hole shadows have radii of curvature, horizontal and vertical diameters that are approximately equal. Further, these values are approximately equal to those of a Kerr black hole shadow. Thus, astrophysically relevant observations (observations for which Λ = 1.11 × 10 −52 m −2 ) cannot distinguish between a RKdS, KdS and Kerr black hole shadow. For Λ = 0.06m −2 , see table II, the quantitative discrepancy between our result and Ref. [1] becomes obvious. But such immense value for Λ is not realistic. The values also indicate that a RKdS black hole casts a smaller shadow than that of a KdS black hole despite the difference being small. Generally, considering our results, an increase in the cosmological constant decreases the size of a KdS and RKdS black hole shadow.
Finally, utilizing the constraint on the characteristic areal radius of the shadow obtained by the EHT collaboration, we have constrained a KdS and RKdS black hole. We find that, for a/M > 0.812311, large angles of inclination θ > 30.5107 • are excluded from M87* observations in both KdS and RKdS black holes.