Abstract
We study the general motion of photons in the Kerr–de Sitter black-hole and naked singularity spacetimes. The motion is governed by the impact parameters X, related to the axial symmetry of the spacetime, and q, related to its hidden symmetry. Appropriate ‘effective potentials’ governing the latitudinal and radial motion are introduced and their behavior is examined by the ‘Chinese boxes’ technique giving regions allowed for the motion in terms of the impact parameters. Restrictions on the impact parameters X and q are established in dependence on the spacetime parameters \(M, \Lambda , a\). The motion can be of orbital type (crossing the equatorial plane, \(q>0\)) and vortical type (tied above or below the equatorial plane, \(q<0\)). It is shown that for negative values of q, the reality conditions imposed on the latitudinal motion yield stronger constraints on the parameter X than that following from the reality condition of the radial motion, excluding the existence of vortical motion of constant radius. The properties of the spherical photon orbits of the orbital type are determined and used along with the properties of the effective potentials as criteria of classification of the KdS spacetimes according to the properties of the motion of the photon.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
In the framework of inflationary paradigm [45], recent cosmological observations indicate that a very small relict vacuum energy (equivalently, repulsive cosmological constant \(\Lambda > 0\)), or, generally, a dark energy demonstrating repulsive gravitational effect, has to be introduced to explain dynamics of the recent Universe [3, 7, 16, 40, 41, 49, 92]. These conclusions are supported strongly by the observations of distant Ia-type supernova explosions indicating that starting at the cosmological redshift \(z \approx 1\) expansion of the Universe is accelerated [55]. The total energy density of the Universe is very close to the critical energy density \(\rho _{\mathrm{crit}}\) corresponding to almost flat Universe predicted by the inflationary scenario [62], and the dark energy represents about \(70\%\) of the energy content of the observable Universe [15, 63]. These conclusions have been confirmed by recent measurements of cosmic microwave background anisotropies by the space satellite observatory PLANCK [1, 52].
The dark energy equation of state is very close to those corresponding to the vacuum energy [15]. Therefore, it is relevant to study the astrophysical consequences of the effect of the observed cosmological constant implied by the cosmological tests to be \(\Lambda \approx 1.3\times 10^{-56}\,\mathrm{cm^{-2}}\), and the related vacuum energy \(\rho _{\mathrm{vac}} \sim 10^{-29}\,\mathrm{g\,cm^{-3}},\) close to the critical density of the Universe. The repulsive cosmological constant changes significantly the asymptotic structure of black-hole, naked singularity, or any compact-body backgrounds as such backgrounds become asymptotically de Sitter spacetimes, and an event horizon (cosmological horizon) always exists, behind which the geometry is dynamic.
A substantial influence of the repulsive cosmological constant has been demonstrated for astrophysical situations related to active galactic nuclei and their central supermassive black holes [69]. The black-hole spacetimes with the \(\Lambda \) term are described in the spherically symmetric case by the vacuum Schwarzschild–de Sitter (SdS) geometry [34, 72], while the internal, uniform density SdS spacetimes are given in [13, 68]. The axially symmetric, rotating black holes are determined by the Kerr–de Sitter (KdS) geometry [17, 25, 44].
In spacetimes with a repulsive cosmological term, the motion of photons was extensively investigated in many papers [8, 38, 39, 43, 47, 57, 58, 71, 73, 91]. The motion of massive test particles was studied in [2, 19, 20, 29, 31, 32, 35,36,37, 48, 51, 59, 66, 67, 72, 74, 86]. The KdS geometry can be relevant also for the so-called Kerr superspinars representing an alternative to black holes [14, 26, 27, 83], breaking the black-hole bound on the dimensionless spin and exhibiting a variety of unusual physical phenomena [21, 22, 30, 64, 76, 80, 83, 84]. It is worth to note that the SdS and KdS spacetimes are equivalent to some solutions of the f(R) gravity representing black holes and naked singularities [50, 88].
The role of the cosmological constant can be significant for both the geometrically thin Keplerian accretion discs [46, 61, 69, 72, 86] and the toroidal accretion discs [6, 42, 53, 54, 60, 87, 89] orbiting supermassive black holes (Kerr superspinars) in the central parts of giant galaxies. Both high-frequency quasiperiodic oscillations and jets originating at the accretion discs can be reflected by current carrying string loops in SdS and KdS spacetimes [28, 33, 77, 78, 93]. In the spherically symmetric spacetimes, the Keplerian and toroidal disc structures can be precisely described by a pseudo-Newtonian potential of Paczynski type [79, 88], which appears to be useful also in studies of the motion of interacting galaxies [56, 81, 82] demonstrating relation of the gravitationally bound galactic systems to the so-called static radius of the SdS or KdS spacetimes [4, 5, 23, 24, 66, 67]. This idea has been confirmed by the recent study of general relativistic static polytropic spheres in spacetimes with the repulsive cosmological constant [75, 85].
The present paper is devoted to a detailed study of the properties of the motion of the photon in the KdS black-hole and naked singularity spacetimes. We concentrate attention to the behavior of the effective potentials determining the regions allowed for the motion of the photon. Such a study is necessary for a full understanding of the optical phenomena occurring in the black-hole or naked singularity spacetimes with a repulsive cosmological constant. We generalize the previous work concentrated on the properties of the motion of the photon in the equatorial plane [73], discussing the properties of the effective potential of the latitudinal motion in terms of the constant of the motion related to the equatorial plane, and then continuing by study of the effective potential of the radial motion. We concentrate our study on the spherical photon orbits representing a natural generalization of the photon circular geodesics that enables a natural classification of the KdS spacetimes according to the properties of the null geodesics representing the motion of the photon.
2 Kerr–de Sitter spacetime and Carter’s equations of geodesic motion
2.1 Kerr–de Sitter geometry
The line element describing the Kerr–de Sitter geometry is in the standard Boyer–Lindquist coordinates, using the geometric system of units (\(c = G = 1\)), given by
where
Here, as usual, we denote by M the mass of the central gravitating body, by a its specific angular momentum (\(a = J/M\)) and by \(\Lambda \) the cosmological constant. In order to simplify the discussion of the following equations, it is convenient to introduce a new cosmological parameter \(y = \frac{1}{3}\Lambda M^2,\) and use dimensionless quantities, redefining them such that \(s/M \rightarrow s\), \(t/M \rightarrow t, r/M \rightarrow r, a/M \rightarrow a\), which is equivalent to putting \(M = 1.\) The above expressions then read
with Eq. (5) being left unchanged.
The physical singularity is located, as in the Kerr geometry, at the ring \(r=0,\) \(\theta =\pi /2.\)
The black-hole horizons are determined by the condition
and their loci can be determined by the relation
The zeros of \(y_{h}(r;\,a^2),\) determining the loci of black-hole horizons in pure Kerr spacetimes, are given by the relation
the loci of its extrema are given by the functions
where the function \(a^2_{\mathrm{{ex}}(h)}-(r)<0\) in its whole definition range and hence is irrelevant. The functions \(y_{h}(r;\,a^2),\) \(a^2_{z(h)}(r)\) and \(a^2_{\mathrm{{ex}}(h)\pm }(r)\) will be needed in the section devoted to the discussion of the radial motion.
Three event horizons, two black hole \(r_-,\) \(r_+,\) and the cosmological horizon \(r_\mathrm{c},\) (\(r_-< r_+ < r_\mathrm{c} \)) exist for \(y_{{\text {min}}(h)}(a^2)< y < y_{{\text {max}}(h)}(a^2),\) where the limits \(y_{{\text {min}}/{\text {max}}(h)}(a^2)\) correspond to a local minimum or local maximum of the function \(y_{h}(r;\,a^2),\) respectively, for given rotational parameter a. For \(0<y<y_{{\text {min}}(h)}(a^2)\) or \(y>y_{{\text {max}}(h)}(a^2)\) Kerr–de Sitter naked singularity spacetimes exist. The limit case \(y=y_{{\text {min}}(h)}(a^2)\) corresponds to an extreme black-hole spacetime, when the two black-hole horizons coalesce. If \(y=y_{{\text {max}}(h)}(a^2),\) the outer black-hole and cosmological horizon merge. There exists a critical value of the rotational parameter \(a^2_{\text {crit}}=1.212 02,\) for which the two local extrema of the function \(y_{h}(r;\,a^2)\) coalesce in an inflection point at \(r_{\text {crit}}=1.616 03\) with the critical value \(y_{\text {crit}}=0.0592.\) Thus, for \(a^2>a^2_{\text {crit}}\) only the Kerr–de Sitter naked singularity can exist for any \(y>0.\)
Properties of the event horizons for the more general case of the Kerr–Newman–de Sitter spacetimes can be found in [73].
2.2 Carter’s equations of geodesic motion
The motion of test particles and photons following geodesics in the Kerr–de Sitter spacetime is described by the well-known Carter equations [17],
where
and
Here E and \(\Phi \) are the constants of motion connected, respectively, with the time and axial symmetry of the Kerr–de Sitter geometry, and \(\mathcal {K}\) is the fourth Carter constant of motion connected with the hidden symmetry of the Kerr–de Sitter geometry. Another constant of motion is the rest mass m (energy) of the test particle; for photons \(m = 0\). Recall that E and \(\Phi \) cannot be interpreted as energy and axial component of the angular momentum of the test particle at infinity, because, due to the presence of the cosmological \(\Lambda \) term, the geometry is not asymptotically flat, but de Sitter [72].
A detailed discussion of the equatorial motion of photons in the Kerr–Newman–de Sitter spacetimes has been published in [73]. Circular motion of test particles in the Kerr–de Sitter spacetimes has been presented in [86]. Here we restrict our attention on the general motion of photons in Kerr–de Sitter spacetimes.
In fact, the motion of photons is independent of the constant of motion E and depends only on the ratio \(\Phi /E\) (\(E \ne 0\)), usually referred to as the impact parameter \(\ell \), and on the parameter \(\mathcal {K}/E^2.\) For our general discussion it is convenient to use \(Q = \mathcal {K} - I^2(\Phi - a E)^2\), which vanishes for the equatorial motion. For our purposes it is, however, following Ref. [73], convenient to introduce a new constant of motion \(X \equiv \ell - a.\) Further the constant \(q \equiv Q/I^2E^2\) will be applied. Then Eqs. (17) and (18) simplify to the form
Following Ref. [73] we study the general motion of the photon in terms of the parameter X. However, since we consider the non-equatorial motion here, it is also necessary to find the restrictions to be imposed on the parameter X that follow from the reality conditions of the latitudinal motion. The latitudinal motion in the Kerr–de Sitter spacetimes has been investigated yet [66]; however, the discussion has been related to the constant of the motion \(\mathcal {K}.\) Here we give the discussion of the effective potential of the latitudinal motion related to the constant of the motion Q, as it is convenient for the purposes of our study.
3 Latitudinal motion
Because it is more convenient to work with algebraic functions instead of trigonometric ones, we introduce a new variable
This implies replacing the Eq. (13) by
where
with the notation
Note that \(\,\mathrm{d}m/\mathrm{d}\lambda =0\) does not necessarily imply \(\,\mathrm{d}\theta /\mathrm{d}\lambda =0,\) since it can mean just transit through the equatorial plane or polar axis. Therefore, in some cases, in order to avoid any doubt, we rather discuss the behavior of the function (19).
The reality condition \(M(m;\,a,\,y,\,X,\,q)\ge 0\) can be expressed by the relations
in regions where \(\Delta _m-a^2y > 0,\) i.e., equivalently,
where
is the solution of the equation
and by the relations
in regions where \(\Delta _m-a^2y < 0,\) i.e., \(m < m_\mathrm{d},\) which requires \(y > 1/a^2\). The functions \(X^\theta _{\pm }(m;\,q,\,y,\,a),\) regarded as ‘effective potentials’ governing the latitudinal motion, are defined by
The functions \(X^\theta _{\pm }(m;\,q,\,y,\,a)\) thus determine the regions allowed for the latitudinal motion, the conditions \(X=X^\theta _{\pm }(m;\,q,\,y,\,a)\) give the turning points. In order to understand the behavior of the functions \(X^\theta _{\pm }(m;\,q,\,y,\,a),\) it is necessary to find the reality regions, and loci of its local extrema and divergencies. Following [73], we shall perform this analysis using the well-known procedure called the ’Chinese boxes technique’ and adopting labeling of the appropriate characteristic functions in a similar way. The parameters are of various significance – q is a constant of motion, whereas a, y govern the geometry. The natural choice is therefore to give the properties of the potentials \(X^\theta _\pm (m;\,q,\,y,\,a)\) by a family of functions \(q(m;\,y,\,a),\) and the properties of these functions by other families of functions of variable m with parameters lowered by one, with spacetime parameters excluded at last.
In the following analysis the relevant range of variable m is, of course, \(0\le m\le 1,\) but somewhere, in order to better understand the behavior of the characteristic functions, we formally permit \(m\in R.\)
First we shall determine the reality region of \(X_\pm (m).\) It is given by
where
Of course, this function also determines the common points of the potentials \(X^\theta _-(m;\,q,\,y,\,a)\) and \(X^\theta _+(m;\,q,\,y,\,a),\) which values are then
Of particular importance, if defined, is the value \(X_{(\pm )}(1;\,y,\,a)=-a\) (see below).
The divergency points of the functions \(q^\theta _r(m;\,y,\,a^2),\) \(X_{(\pm )}(m;\,y,\,a)\) and \(X^\theta _-(m;\,q,\,y,\,a)\) are determined by
The functions \(X^\theta _{\pm }(m;\,q,\,y,\,a)\) both can diverge, if well defined, for \(m=0,\) other divergencies are given by the function \(y^\theta _\mathrm{d}(m;\,a^2)\) for the potential \(X^\theta _{-}(m;\,q,\,y,\,a),\) but there are no other divergencies for the potential \(X^\theta _{+}(m;\,q,\,y,\,a),\) as can be seen, if we rewrite the definition (28) in the alternative form
The function \(y^\theta _\mathrm{d}(m;\,a^2)\rightarrow \infty \) for \(m\rightarrow 1\) from the left. There are no local extrema of this function and for \(0\le m<1\) it is increasing. For \(m=0\) we get \(y^\theta _\mathrm{d}(0;\,a^2)=1/a^2.\)
The point \(m_\mathrm{d}\) given by the definition (25) determines the loci where the functions \(q^\theta _r(m;\,y,\,a^2),\) \(X_{(\pm )}(m;\,y,\,a)\) and \(X^\theta _-(m;\,q,\,y,\,a)\) diverge; it occurs at relevant interval (0; 1) for \(y>1/a^2\) and \(m_\mathrm{d}\rightarrow 1\) for \(a^2y\rightarrow \infty .\) In such case, \(q^\theta _r(m;\,y,\,a^2)\rightarrow +\infty \;(-\infty )\) for \(m\rightarrow m_\mathrm{d}\) from the left (right).
From the equality
one can see that the function \(q^\theta _r(m;\,y,\,a^2)\) has no local extrema and is decreasing for \(y<1/a^2,\) or piecewise increasing with the discontinuity point \(m_\mathrm{d}\) for \(y>1/a^2,\) i.e., \(q^\theta _r(m;\,y,\,a^2) \rightarrow +\infty \;(-\infty )\) for \(m \rightarrow m_\mathrm{d}\) from the left (right). We always have \(q^\theta _r(m=0;\,y,\,a^2)=0\) and \(q^\theta _r(m=1;\,y,\,a^2)=-a^2.\)
In the special case \(y=1/a^2\) we get
with
Based on the conditions (29), (30) and the above characteristic functions, we can complete the setting of the definition range of the potentials \(X^\theta _\pm (m;\,q,\,a,\,y),\) which we leave to the end of this section.
Now we shall determine the loci of local extrema of the effective potentials \(X^\theta _{\pm }(m;\,q,\,y,\,a).\) They can be derived from the condition \(\partial X^\theta _{\pm }/\partial m =0,\) which implies the equation
It can be verified that the function \(X^\theta _+(m;\,q,\,y,\,a)\) has local extrema given by the relation
A discussion of this function is trivial, so we only note that it is independent of the cosmological parameter y and renders the loci of extrema only for \(-a^2\le q \le 0,\) while, as we shall see below, they can exist even for \(q<-a^2.\) The character of these extrema reveals inserting this expression into the second derivative, which yields
clearly they must be maxima.
From Eq. (36) we find that other extrema of the potentials \(X^\theta _{\pm }(m;\,q,\,y,\,a)\) are determined by the condition
The divergencies of the functions \(q^\theta _{\mathrm{{ex}}(\pm )}(m;\,y,\,a^2)\) are determined by the relation
The function \(y^\theta _{\mathrm{d}(\mathrm{{ex}}\pm )}(m;\,a^2)\) is positively valued at \(0\le m<0.5,\) where \(y^\theta _{\mathrm{d}(\mathrm{{ex}}\pm )}(m;\,a^2)\rightarrow +\infty \) for \(m\rightarrow 0.5\) from the left. For \(m=0,\) there is \(y^\theta _{\mathrm{d}(\mathrm{{ex}}\pm )}(0;\,a^2)=y^\theta _{\mathrm{{d}}(r)}(0;\,a^2)=1/a^2.\)
From the properties of the function \(y^\theta _{\mathrm{d}(\mathrm{{ex}}\pm )}(m;\,a^2)\) we deduce that the function \(q^\theta _{\mathrm{{ex}}(\pm )}(m;\,y,\,a^2)\) can diverge only if \(y>1/a^2,\) at
located such that \(0\le m_{\mathrm{d}(\mathrm{{ex}})} <0.5.\) Obviously
In the following we shall decide on the monotony and possible existence of local extrema of the function \(q_{ex(\pm )}.\) From
it is clear that there are no local extrema of \(q^\theta _{\mathrm{{ex}}(\pm )}(m;\,y,\,a^2)\) in the interval \(m\in (0;1).\) The derivative changes sign at the divergent point \(m_{\mathrm{d}(\mathrm{{ex}})},\) which reflects the behavior given by (41).
For \(y<1/a^2,\) there is
that is, \(q^\theta _{\mathrm{{ex}}(\pm )}(m;\,y,\,a^2)\) is decreasing.
In the limit case \(y=1/a^2,\) we get
Comparing the two functions \(q^\theta _{\mathrm{{ex}}+}(m;\,a^2)\) and \(q^\theta _{\mathrm{{ex}}(\pm )}(m;y,\,a^2),\) we find that
and have common points at \(m = 0, 1\) with
and
In the next step we shall characterize the extrema given by \(q^\theta _{\mathrm{{ex}}(\pm )}(m;\,y,\,a^2)\).
First we find that
If we now require
somewhere at \(0< m < 1,\) we obtain the condition
which ensures
Considering the previous results, we can conclude that the functions \(q^\theta _{\mathrm{{ex}}(\pm )}(m;\,y,\,a^2)\) determine local maxima of the potential \(X^\theta _+(m;\,q,\,y,\,a)\) for \(q<-a^2\) that occur on this curve in the case \(y>1/a^2\) in the interval \(m \in (m_{\mathrm{d}(\mathrm{{ex}})};1).\)
Proceeding in the same way with the function \(X^\theta _-(m;\,q,y,\,a),\) we first find that the equation
has always a solution for some \(m\in (0; 1)\) in the case \(y<1/a^2,\) but for \(y\ge 1/a^2\) this solution must fulfill \(m<m_{\mathrm{d}(\mathrm{{ex}})}.\)
Substituting \(q=q^\theta _{\mathrm{{ex}}(\pm )}\) into the second derivative of \(X^\theta _-(m;\,q,\,y,\,a)\) yields the same expression as that in (44), but now with the above conditions we have
indicating local minima. Therefore, the function \(q^\theta _{\mathrm{{ex}}(\pm )}\) gives local minima of \(X^\theta _-(m;\,q,\,y,\,a)\) for \(y<1/a^2\) at the whole interval \(m \in (0;1),\) and for \(y>1/a^2\) at \(m \in (0;m_{\mathrm{d}(\mathrm{{ex}})}).\)
In the special case \(y=1/a^2,\) the function \(q^\theta _{\mathrm{{ex}}(\pm )}\) reduces to the form
we can easily convince ourself that the function \(X^\theta _-\) has no local extremum in such case, and the extrema of \(X^\theta _+\) are given by the function \(q^\theta _{\mathrm{{ex}}+}(m;\,y,\,a).\)
The conditions (24) and (27), ensuring the allowance of the latitudinal motion, must be complemented by the case that the functions \(X^\theta _\pm \) are not defined. Their definition range is given by Eqs. (29) and (30), but it can be shown that the violation of the latter one imply \(M(m;\, X,\,q,\,y,\,a)>0.\) In such a case, the latitudinal motion is allowed for any impact parameter X (see for details the discussion below). All characteristic functions are depicted in Fig. 1, and the graphs of the potentials in Fig. 2 for selected representative values of parameters.
Now we are able to discuss the behavior of the potentials \(X^\theta _{\pm }(m;\,q,\,y,\,a)\) for various representative values of its parameters. The intersections of a line \(X = \text{ const. }\) with the curves \(X^\theta _{\pm }(m;\,q,\,y,\,a)\) represents the turning points in variable m. From the knowledge of these functions we can thus get qualitative insight into the character of the latitudinal motion. This entitles us to write the following classification of Kerr–de Sitter spacetimes and brief description of the latitudinal motion. The basic division apparently consists of the cases \(y<1/a^2,\;y=1/a^2,\;y>1/a^2\):
-
1.
Case \(y<1/a^2\)
-
\(q<-a^2\)
-
the definition range of the potentials is an empty set; the latitudinal motion is not possible;
-
-
\(q=-a^2\)
-
the potentials \(X^\theta _\pm (m;\,q,\,y,\,a)\) are defined only for \(m=1,\) where
$$\begin{aligned} X^\theta _+(1;\,q,\,y,\,a)=X^\theta _-(1;\,q,\,y,\,a)=-a; \end{aligned}$$photons with such values of parameters are the special case of the so called PNC photons ‘radially’ moving along the spin axis [10];
-
-
\(-a^2<q<0\) (Fig. 2a, b)
-
the two potentials are defined for \(m\in \langle m_l;1\rangle ,\) where the lower limit
$$\begin{aligned} m_l = \frac{q(a^2y-1)}{a^2(qy+1)}>0 \end{aligned}$$(46)is the solution of the equation \(q=q^\theta _r(m;\,y,\,a)\) (see Fig. 1b); the limits of the interval are the common points of the potentials, where
$$\begin{aligned} X^\theta _\pm (m=m_l;\,q,\,y,\,a)=X^\theta _{(\pm )}(m_l)=\frac{a(1+qy)}{a^2y-1}<0; \end{aligned}$$(47)
-
-
-
the latitudinal motion is allowed for values of the parameter X between some local minimum \(X^\theta _{{\text {min}}(-)}=X^\theta _-(m_{{\text {min}}(-)};\,q,\,y,\,a)\) and maximum \(X^\theta _{{\text {max}}(+)}=X^\theta _+(m_{{\text {max}}(+)};\,q,\,y,\,a),\) for which
$$\begin{aligned} X^\theta _{{\text {min}}(-)}<-a<X^\theta _{{\text {max}}(+)}<0; \end{aligned}$$the loci \(m_{{\text {min}}(-)}\) of minimum \(X^\theta _{{\text {min}}(-)}\) is given by the Eq. (39), the loci \(m_{{\text {max}}(+)}\) of maximum \(X^\theta _{{\text {max}}(+)}\) is determined by Eq. (37);
-
if X takes one of these extremal values, then the trajectory of such photon lies entirely on cones \(\theta =\arccos \sqrt{m_\mathrm{{ex}}},\) \(\theta =\pi -\arccos \sqrt{m_\mathrm{{ex}}},\) where \(m_\mathrm{{ex}}\in \{m_{min(-)},m_{max(+)}\};\) such photons are called PNC photons [10];
-
for \(X^\theta _{{\text {min}}(-)}<X<X^\theta _{{\text {max}}(+)}\) there are two solutions \(m_1<m_2\) of each of the two equations \(X=X^\theta _\pm (m;\,q,\,y,\,a),\) implying that the photon executes a so-called vortical motion, which is restricted between two pairs of cones, symmetrically placed relative to equatorial plane:
$$\begin{aligned} 0< \arccos \sqrt{ m_2} \le \theta \le \arccos \sqrt{m_1} < \frac{\pi }{2} \end{aligned}$$and
$$\begin{aligned} \frac{\pi }{2}< \pi - \arccos \sqrt{m_1} \le \theta \le \pi -\arccos \sqrt{m_2} < \pi ; \end{aligned}$$ -
in the special case \(X=-a\) one of the turning points is \(m_2=1,\) which represents a transition through the spin axis; such photons therefore oscillate above one of the poles in the cone which is delimited by the angle \(\theta =\arccos \sqrt{m_{1}};\)
-
from the preceding discussion it follows that we can expect that the case \(X=-a\) represents a change in azimuthal direction with respect to some privileged family of observers;
-
\(q=0\) (Fig. 2c)
-
the expression in the definition (28) can be reduced to
$$\begin{aligned} X^\theta _\pm (m;\,y,\,a)=\frac{-a(1\mp \sqrt{(1-m)\Delta _m})}{\Delta _m-a^2y}, \end{aligned}$$(48)which validity can be enlarged, without any repercussion on the correctness of the analysis, even for \(m=0;\) the definition range of the potentials is thus \(\langle 0;1\rangle ;\)
-
from the equality \(W(\theta =\pi /2;\,X,\,q,\,y,\,a)=q\) it follows that at least in the equatorial plane the (radial) motion always exists for \(q=0,\) where it can be both stable or unstable (see below); for \(q>0\) the equatorial plane is crossed, for \(q<0\) it cannot be reached;
-
there are no extrema of the potentials – \(X^\theta _+(m;\,q,y,\,a)\) is decreasing, \(X^\theta _-(m;\,q,\,y,\,a)\) is increasing; the permissible values of X for which \(\,\mathrm{d}\theta /\mathrm{d}\lambda >0\) are still confined to an interval with limits
$$\begin{aligned} X^\theta _{{\text {min}}(-)}= & {} X^\theta _-(m=0;\,q=0,\,y,\,a) \nonumber \\= & {} \frac{2a}{(a^2y-1)},\end{aligned}$$(49)$$\begin{aligned} X^\theta _{{\text {max}}(+)}= & {} X^\theta _+(m=0;\,q=0,\,y,\,a)=0, \end{aligned}$$(50)where \(X^\theta _{{\text {min}}(-)}<X^\theta _{{\text {max}}(+)};\)
-
if \(X\le X^\theta _{{\text {min}}(-)}\) or \(X\ge X^\theta _{{\text {max}}(+)}\) then the requirement \(W(\theta )\ge 0\) is fulfilled only if \(\theta =\pi /2,\) and in such case \(\,\mathrm{d}\theta /\mathrm{d}\lambda =0,\) thus the motion is stably confined to the equatorial plane;
-
for \(X^\theta _{{\text {min}}(-)}<X<X^\theta _{{\text {max}}(+)}\) photon initially released in the direction off the equatorial plane is once reflected at \(\theta =\arccos \sqrt{m_{0}}\) or \(\theta =\pi -\arccos \sqrt{m_{0}}\), respectively, where \(m_0\) denotes the only solution of \(X=X^\theta _\pm (m;\,q,\,y,\,a);\) another point where \(\,\mathrm{d}\theta /\mathrm{d}\lambda =0\) is now in the equatorial plane, however, the equality \(\,\mathrm{d}^2\theta /\mathrm{d}\lambda ^2=0\) implies halting in the latitudinal direction; the function \(W(\theta )\) has at \(\theta =\pi /2\) local minimum, which indicates, as follows from perturbation analysis, instability in the equatorial plane;
-
if specially \(X=-a\) then \(m_0=1,\) thus photon initially directed off the equatorial plane crosses the spin axis and finally is captured in the equatorial plane;
-
-
\(q>0\) (Fig. 2d)
-
the potentials are defined for \(m\in (0,1\rangle ;\) they are monotonous in the same manner as in the case \(q=0,\) but \(X^\theta _+(m;\,q,\,a,\,y)\rightarrow +\infty \) and \(X^\theta _-(m;\,q,\,y,\,a)\rightarrow -\infty \) as \(m\rightarrow 0;\)
-
from the behavior of the potentials it follows that for \(X\ne -a\) a photon is forced to oscillate in the \(\theta \) direction through the equatorial plane between two cones governed by \(\arccos \sqrt{m_{0}}\le \theta \le \pi -\arccos \sqrt{m_{0}},\) with \(m_0\) of the same meaning as above;
-
case \(X=-a\) represents the motion above both poles;
-
the foregoing conclusion is a reason to have a suspicion that cases \(X<-a\) and \(X>-a\) differ in the azimuthal direction relative to some family of stationary observers, it corresponds to \(\ell >0\) and \(\ell <0\).
-
-
2.
Case \(y=1/a^2\)
-
the potentials simplify into the form
$$\begin{aligned} X^\theta _\pm (m;\;q,\;a)=\frac{-a\pm \sqrt{(1-m^2)(q+a^2)}}{m}; \end{aligned}$$(51)
-
\(q<-a^2\)
-
the potentials are not defined, thus the latitudinal motion is not allowed;
-
-
-
\(q=-a^2\)
-
the curves \(X=X^\theta _\pm (m;\;q=-a^2,\,y=1/a^2,\,a)\) coalesce, since
$$\begin{aligned}&X^\theta _+(m;q=-a^2,y=1/a^2,a) \nonumber \\&\quad = X^\theta _-(m;q=-a^2,y=1/a^2,a)= X^\theta _{(\pm )}(m;\;a)\nonumber \\&\quad \equiv \frac{-a}{m}; \end{aligned}$$(52)
-
-
for \(X\le -a\) there is one solution of the equation \(X=X^\theta _{(\pm )}(m;\;a),\) which gives \(m=m_{(\pm )}\equiv -a/X;\) this corresponds to PNC photons moving along cones \(\theta =\arccos \sqrt{m_{(\pm )}},\) \(\theta =\pi -\arccos \sqrt{m_{(\pm )}};\)
-
for \(X\rightarrow -\infty \) the cones approach the equatorial plane;
-
if specially \(X=-a\) the cones degenerate to spin axis, therefore, such PNC photons move along the spin axis;
-
for \(X>-a\) there is no motion allowed;
-
\(-a^2<q<0\) (Fig. 2e, f)
-
the potentials are both defined for \(m\in (0;1\rangle \); there is one local maximum \(X^\theta _{{\text {max}}(+)}\) given by (37) of the function \(X^\theta _+(m;\,q,\,y,\,a)\) and no extremum of \(X^\theta _-(m;\,q,\,y,\,a);\) we have \(X^\theta _-(m;\,q,\,y,\,a)<X^\theta _+(m;\,q,\,y,\,a)<0\) and \(X^\theta _-(m;\,q,\,y,\,a),X^\theta _+(m;\,q,y,\,a)\rightarrow -\infty \) as \(m \rightarrow 0\) from the right;
-
if \(X<-a\) or \(-a<X<X^\theta _{{\text {max}}(+)},\) the vortical motion exists;
-
for \(X=-a\) the ‘inner’ cones coalesce with the spin axis, thus the vortical motion involves crossing the poles;
-
for \(X=X^\theta _{{\text {max}}(+)}\) both the ‘inner’ and ‘outer’ cones coalesce, giving thus rise to PNC photons;
-
if \(X>X^\theta _{{\text {max}}(+)},\) no motion is allowed;
-
-
\(q=0\) (Fig. 2g)
-
the same discussion holds as in the case \(y<1/a^2,\) except that the motion exists for X arbitrarily small;
-
-
\(q>0\) (Fig. 2h)
-
the same conclusions hold as in the case \(y<1/a^2;\)
-
-
3.
Case \(y>1/a^2\)
-
\(q<-a^2\) (Fig. 2i)
-
the definition range of both potentials is an interval \((0;m_u\rangle \) (see the purple curve in Fig. 1d), where the upper limit \(m_u<1\) is given as \(m_l\) in the previous case by (46);
-
there is \(X^\theta _+(m;\,q,\,y,\,a)\rightarrow -\infty \) and \(X^\theta _-(m;\,q,y,\,a)\rightarrow +\infty \) as \(m\rightarrow 0,\) moreover, \(X^\theta _-(m;\,q,y,\,a)\) now diverge at \(m=m_\mathrm{d},\) which is the solution of (33), and \(X^\theta _-(m;\,q,\,y,\,a)\rightarrow +\infty \;(-\infty )\) as \(m\rightarrow m_\mathrm{d}\) from the left (right);
-
there are thus two regions of permissible values X in the (m, X)-plane for which the motion can exist; the lower one bounded by the graph of \(X^\theta _+\) and the lower branch of \(X^\theta _-,\) which at \(m=m_u\) join into continuous curve, and the upper region given by the upper branch of \(X^\theta _-;\) the motion is therefore allowed for \(X\le X^\theta _{{\text {max}}(+)}<-a\) or \(X\ge X^\theta _{{\text {min}}(-)}>0,\) where the loci of local extrema \(X^\theta _{{\text {max}}(+)},\) \(X^\theta _{{\text {min}}(-)}\) are given by (39) (see the blue curve in Fig. 1d);
-
if \(X<X^\theta _{{\text {max}}(+)}\) or \(X>X^\theta _{{\text {min}}(-)}\) photon executes vortical motion, the cases \(X=X^\theta _{{\text {max}}(+)}\), \(X=X^\theta _{{\text {min}}(-)}\) correspond to PNC photons;
-
for \(X=X^\theta _-(m_u)=X^\theta _+(m_u)=a(1+qy)/(a^2y-1),\) the inner cones delimiting the vortical motion are the narrowest;
-
for \(X\rightarrow -\infty \) or \(X\rightarrow +\infty \) the outer cones given by angles
$$\begin{aligned}\theta =\arccos \sqrt{m_1},\quad \theta =\pi -\arccos \sqrt{m_1}\end{aligned}$$approach the equatorial plane since \(m_1\rightarrow 0;\) for the inner cones
$$\begin{aligned}\theta =\arccos \sqrt{m_2},\quad \theta =\pi -\arccos \sqrt{m_2},\end{aligned}$$we have
$$\begin{aligned}m_2\rightarrow m_\mathrm{d}=1-1/a^2y;\end{aligned}$$
-
-
\(q=-a^2\) (Fig. 2j)
-
there is no local extremum of the function \(X^\theta _+(m;\,q,\,y,\,a),\) which is now increasing; we have \(m_u=1,\) \(X^\theta _-(m_u)=X^\theta _+(m_u)=X^\theta _{+({\text {max}})}=-a,\) hence for \(X=-a\) both the inner and outer cones coalesce with the spin axis, which again corresponds to ‘axial’ PNC photon;
-
other PNC photons exist for \(X=X^\theta _{{\text {min}}(-)}>0;\)
-
there are no other qualitative differences from the case \(q<-a^2;\)
-
-
\(-a^2<q<0\) (Fig. 2k)
-
the definition range is an interval \((0;1\rangle \) and the divergencies of the potentials are the same as above;
-
the function \(X^\theta _+(m;\,q,\,y,\,a)\) has now local maximum \(X^\theta _{{\text {max}}(+)},\) \(-a<X^\theta _{{\text {max}}(+)}<0,\) \(X^\theta _{{\text {max}}(+)}\rightarrow 0\) for \(q\rightarrow 0,\) determined by Eq. (37);
-
case \(X=-a\) now corresponds to vortical motion above the poles – the inner cones have coalesced with the spin axis, the outer ones stay open;
-
the vortical motion exists as in the previous cases and above that for \(-a<X<X^\theta _{{\text {max}}(+)};\)
-
-
\(q=0\) (Fig. 2l)
-
the definition (48) holds, the functions \(X^\theta _{+(-)}(m;q,\,y,\,a)\) are defined at \(\langle 0,1\rangle \) (\( \langle 0,1\rangle \setminus \{m_\mathrm{d}\});\) the values for \(m=0\) are given by (50), but now \(X^\theta _{{\text {max}}(+)}<X^\theta _{{\text {min}}(-)};\)
-
the potential \(X^\theta _+(m;\,y,\,a)\) is decreasing in its whole definition range, \(X^\theta _-(m;\,y,\,a)\) is piecewise increasing because of the divergent point \(m_\mathrm{d};\)
-
if \(X\le X^\theta _{{\text {max}}(+)}=0\) or \(X\ge X^\theta _{{\text {min}}(-)}\) then the same conclusions can be made as in the case \(y<1/a^2\) for \(X^\theta _{{\text {min}}(-)}\le X\le X^\theta _{{\text {max}}(+)};\)
-
for \(X^\theta _{{\text {max}}(+)}< X< X^\theta _{{\text {min}}(-)}\) we have \(W(\theta =\pi /2;\,X,\,q=0,\,y,\,a) =0\) again, otherwise \(W(\theta ;\,X,\,q=0,\,y,\,a)<0,\) therefore photons can radially move in the equatorial plane;
-
-
\(q>0\) (Fig. 2m, n)
-
the function \(X^\theta _+(m;\,q,\,y,\,a)\) is defined at \(\langle m_l,1\rangle ,\) the function \(X^\theta _-(m;\,q,\,y,\,a)\) at \(\langle m_l,1\rangle \setminus \{m_\mathrm{d}\},\) where \(m_l\) is given by (46) with the difference that now \(X^\theta _{(\pm )}(m_l)>0\); the graphs of both functions now form a single open curve, which intersects a line \(X={\text {const.}}\) at a single point;
-
in the interval \(m \in \langle 0; m_l \rangle \) the latitudinal motion is allowed for arbitrarily large or small value of the motion constant X;
-
for arbitrary \(X\ne -a\) there exists oscillatory motion through the equatorial plane as described in the case \(y<1/a^2, q>0;\)
-
if \(X=X_{(\pm )}(m_l)\) the boundary cones are closest to equatorial plane, they are given by angles
$$\begin{aligned} \theta =\arccos {\sqrt{m_l}},\quad \theta =\pi -\arccos {\sqrt{m_l}}; \end{aligned}$$ -
the case \(X=-a\) corresponds to orbits above both poles crossing also the equatorial plane;
-
there is no vortical motion or PNC photons.
-
-
We finish this section with setting the allowed region in the (X, q)-plane delimiting such combinations of the constants of the motion, for which the latitudinal motion is possible, in dependence on the spacetime parameters a, y.
From the requirement that the function \(M(m;\,a,\,y,\,X,\,q)\) defined in Eq. (22) has to be non-negative somewhere in the interval \(m \in \langle 0;1\rangle ,\) one can derive that the allowed region of the \((X-q)\)-plane is determined by the condition
where \(q_{\text {min}}(X,\,y,\,a)\) is defined using functions
and
-
Case \(y<1/a^2\) (Fig. 3a)
$$\begin{aligned} q_{\text {min}}(X,\,y,\,a) \equiv \left\{ \begin{array}{ll} 0, &{} \text {for}\quad X< \frac{2 a}{a^2 y - 1}\quad \text {or}\quad X > 0;\\ q_2(X;\,y,\,a), &{} \text {for}\quad \frac{2 a}{a^2y-1} \le X < -a;\\ q_1(X), &{} \text {for}\quad -a \le X \le 0;\\ \end{array}\right. \end{aligned}$$(56) -
Case \(y=1/a^2\) (Fig. 3b)
$$\begin{aligned} q_{\text {min}}(X,\,y,\,a) \equiv \left\{ \begin{array}{ll} -a^2, &{} \text{ for }\quad X \le -a;\\ \\ q_1(X),&{} \text{ for }\quad -a \le X \le 0;\\ 0,&{}\text {for}\quad X\ge 0; \end{array}\right. \end{aligned}$$(57) -
Case \(y>1/a^2\) (Fig. 3c)
$$\begin{aligned} q_{\text {min}}(X,\,y,\,a) \equiv \left\{ \quad \begin{array}{ll} q_2(X;\,y,\,a),&{}\text {for}\quad X< -a \quad \text {or} \quad \frac{2 a}{a^2 y - 1}< X;\\ q_1(X),&{}\text {for}\quad -a \le X \le 0;\\ 0,&{}\text {when}\quad 0 < X \le \frac{2 a}{a^2 y - 1}; \end{array}\right. \end{aligned}$$(58)
The case \(y<1/a^2\) qualitatively corresponds to both black-hole and naked singularity spacetimes, the other two cases, \(y=1/a^2\) and \(y>1/a^2\), describe the naked singularity spacetimes (see Fig. 7 in the next section). The \(q={\text {const.}}\) slices of the function \(q_{\text {min}}(X,\,y,\,a)\) give for \(q<0\) extremal values \(X_{{\text {min}}(-)},\) \(X_{{\text {max}}(+)}\) of the potentials \(X^\theta _{\pm }(m;\,q,\,y,\,a)\) discussed in the text.
4 Radial motion
From Eq. (14) it is clear that the radial motion can exist if \(R(r)\ge 0,\) where the equality gives the turning points of the radial motion. This condition can be rewritten in terms of an ‘effective potential’ \(X_\pm \) in the form
if
or
where
We start the analysis by determining the reality region of the effective potential \(X_\pm .\) From Eq. (61) it follows that this function is well defined for
where we have introduced the reality function
There are thus two different types of the boundary points of the definition range of \(X_{\pm }(r;\,q,\,a,\,y).\) The points of the first type lie stably for given spacetime parameters on the borders of the static regions determined by Eq. (9), i.e. at the event horizons \((r=r_h).\) At these horizons, if they exist, for arbitrary parameter q, the functions \(X_{\pm }\) have common values,
(cf. [73]). The points of the second type, which are also common points of \(X_{\pm },\) depend on the value of parameter q and are given by the equality \(q = q_{r}(r;\,y,\,a^2).\) If we denote them \(r = r_q\) then we have
The divergencies of the function \(q_{r}(r;\,y,\,a^2)\), which are incident with the divergent points of \(X_{+}(r;\,q,\,y,\,a),\) are located at radii where
which one can express by the relation
The function \(X_{-}(r;\,q,\,y,\,a)\) cannot diverge at radii \(r_\mathrm{d}\) given by (66), since using the alternative expression
it can be shown that it has finite value
Another point where the functions \(X_{\pm }(r;\,q,\,y,\,a)\) diverge is \(r=0,\) with \(X_{\pm }(r;\,q,\,y,\,a)\rightarrow \pm \infty \) as \(r\rightarrow 0\) for \(q>0,\) but for \(q=0\) we have \(X_{\pm }(r;\,q=0,\,y,\,a)\rightarrow 0.\)
The character of the function \(y_\mathrm{d}(r;\,a^2)\) has been discussed in [70, 73], therefore we briefly repeat that the only zero of \(y_\mathrm{d}(r;\,a^2)\) is at \(r=2,\) the extrema, which for \(a^2>0\) must be maxima, yields the relation
The only zero of \(q_{r}(r;\,y,\,a^2)\) is at \(r = 0.\) For \(r \rightarrow \infty \) we have \(q_{r}(r;\,y,\,a) \rightarrow -1/y.\) Its extrema are determined by
The divergency of \(y_{\mathrm{{ex}}(r)}(r; \,a^2)\) is at \(r = 0\) and \(y_{\mathrm{{ex}}(r)}(r; \,a^2)\rightarrow -\infty \) for \(r\rightarrow 0.\) For \(r \rightarrow \infty \) it approaches the line \(1/a^2\) from below. The zero is at \(r = 3\) and its extrema do not exist, the function is purely increasing.
Now we shall specify the local extrema of the effective potential, which determine the radii of spherical photon orbits. They are given by the condition \(\partial X_{\pm }/\partial r = 0,\) which implies
or
This can be rewritten in terms of parameter q as
and
Note that the function \(q_{\mathrm{{ex}}1}(r;\,a^2)\) is independent of the cosmological parameter. The two functions \(q_{\mathrm{{ex}}1}(r;\,a^2),\) \(q_{\mathrm{{ex}}}(r;\,y,\,a^2)\) have common points determined by
and
i.e., they are located at event horizons and so-called static radius \(r_{s}=1/\root 3 \of {y}\), where the gravitational attraction is just compensated by cosmological repulsion [66, 72]. The function \(q_{\mathrm{{ex}}1}(r;\,a^2)\) is negative valued and hence, as we shall see below, the extrema of the potentials \(X_\pm \) determined by this function lie in regions forbidden by conditions for the reality of latitudinal motion.
The divergencies of \(q_{\mathrm{{ex}}}(r;\,y,\,a^2)\) are determined by the relation
its asymptotic behavior is given by \(q_\mathrm{{ex}}(r;\,y,\,a^2)\rightarrow -(I/2ay)^2\) as \(r\rightarrow \infty .\)
The function \(y_{\mathrm{d}(\mathrm{{ex}})}(r;\,a^2)\) diverges for \(r=0\) and \(y_{\mathrm{d}(\mathrm{{ex}})}(r;\,a^2)\rightarrow -\infty \) as \(r\rightarrow 0.\) For \(r\rightarrow \infty \) we have \(y_{\mathrm{d}(\mathrm{{ex}})}(r;\,a^2)\rightarrow 0.\) The zero of this function is at \(r=1\) and its local extrema are determined by the relation
where the label ‘max’ indicates that at the relevant range, \(r\ge 3/2,\) these extrema must be maxima.
The zero point of the function \(q_{\mathrm{{ex}}}(r;\,y,\,a^2)\) is at \(r=0,\) other zeros determine the loci of the circular equatorial photon orbits. They are given by the relation
Since the function \(y_{z(\mathrm{{ex}})-}(r;\,a^2) < 0\) for \(r>0,\) it is irrelevant in our discussion. The function \(y_{z(\mathrm{{ex}})+}(r;\,a^2)\) is real valued for all \(r>0\) and it diverges at \(r=0,\) with \(y_{z(\mathrm{{ex}})+}(r;\,a^2)\rightarrow \infty \) as \(r\rightarrow 0.\) For \(r\rightarrow \infty ,\) we find \(y_{z(\mathrm{{ex}})+}(r;\,a^2)\rightarrow -1/a^2.\) Its zeros represent the equatorial circular photon orbits in the Kerr spacetimes, being determined by the relation [65]
The extrema of the function \(y_{z(\mathrm{{ex}})+}(r;\,a^2)\) are determined by the equation
hence the loci of extrema of the functions \(y_{z(\mathrm{{ex}})+}(r;\,a^2)\) and \(y_{h}(r;\,a^2)\) coalesce.
The function \(a^2_{\mathrm{{ex}}(z(\mathrm{{ex}})+)-}(r)\) should be excluded from further analysis since for \(r>0\) we have \(a^2_{\mathrm{{ex}}(z(\mathrm{{ex}})+)-}(r)<0.\)
It remains to determine loci of the local extrema of the function \(q_{\mathrm{{ex}}}(r;\,y,\,a^2).\) Proceeding in the usual way we find that their occurrence is governed by the relations
and
Using Eq. (83), one can show that the extrema of the two functions \(q_{\mathrm{{ex}}}(r;\,y,\,a^2)\) \(q_{r}(r;\,y,\,a^2)\) coalesce. At this point let us add that other common points of the functions \(q_{\mathrm{{ex}}}(r;\,y,\,a^2),\) \(q_{r}(r;\,y,\,a^2)\), as well as of the functions \(q_{\mathrm{{ex}}1}(r;\,y,\,a^2),\) \(q_{r}(r;\,y,\,a^2),\) are also given by
i.e. they are located at the event horizons.
The reality conditions of the functions \(y_{\mathrm{{ex}}(\mathrm{{ex}})\pm }(r;\,a^2)\) read
and
where
and
The marginal radius \(\hat{r}\) has value
and we have
where \(a^2_{\text {crit}}\) corresponds to a local maximum of the function \(a^2_{\mathrm{{ex}}(h)+}(r)\) (see e.g. [73] for details).
The functions \(y_{\mathrm{{ex}}(\mathrm{{ex}})\pm }(r;\,a^2)\) have the divergency point at \(r=0\) and \(y_{\mathrm{{ex}}(\mathrm{{ex}})\pm }(r;\,a^2)\rightarrow \pm \infty \) for \(r\rightarrow 0.\) For \(r\rightarrow \infty \) we find that \(y_{\mathrm{{ex}}(\mathrm{{ex}})+}(r;\,a^2)\rightarrow \infty \) and \(y_{\mathrm{{ex}}(\mathrm{{ex}})-}(r;\,a^2)\rightarrow 0\) from above.
The zero point of \(y_{\mathrm{{ex}}(\mathrm{{ex}})}(r;\,a^2)\) is at \(r=3\) and the function is increasing for all \(r>0.\)
The zeros of the functions \(y_{\mathrm{{ex}}(\mathrm{{ex}})\pm }(r;\,a^2)\) are given by
The condition for the stationary points
leads to
which can be solved with respect to \(a^2\) with the same result as given by (12). However, substitution into the second derivative concurrently with the requirement \(y_{\mathrm{{ex}}(\mathrm{{ex}})\pm }(r;\,a^2)>0\) implies
therefore the function
determines the loci of the inflection points of the functions \(y_{\mathrm{{ex}}(\mathrm{{ex}})\pm }(r;\,a^2).\)
If we compare the asymptotic behavior of all characteristic functions \(y(r;\,a^2),\) we find that the following inequality is satisfied: \(1/a^2>y_{\mathrm{{ex}}(\mathrm{{ex}})}(r;\,a^2)>y_{\mathrm{{ex}}(\mathrm{{ex}})-}(r;\,a^2)>y_h(r;\,a^2)>y_\mathrm{d}(r;\,a^2)>y_{\mathrm{d}(\mathrm{{ex}})}(r;\,a^2)\) \(>0>y_{z(\mathrm{{ex}})}(r;\,a^2)>-1/a^2\) as \(r\rightarrow \infty .\)
In Fig. 5 we present all the characteristic functions related to spin parameter governing the effective potential on the lowest level:
-
\(a^2_{z(h)}(r)\);
-
\(a^2_{\mathrm{{ex}}(h)+}(r)\)=\(a^2_{\mathrm{{ex}}(z(\mathrm{{ex}}))+}(r)\)=\(a^2_{{\text {inf}}(\mathrm{{ex}}(\mathrm{{ex}})\pm )+}(r)\);
-
\(a^2_{{\text {max}}(\mathrm{d})}(r)\);
-
\(a^2_{{\text {max}}(\mathrm{d}(\mathrm{{ex}}))}(r)\);
-
\(a^2_{z(z(\mathrm{{ex}})+)}(r)\);
-
\(a^2_{r(\mathrm{{ex}}(\mathrm{{ex}})\pm )+}(r)\);
-
\(a^2_{r(\mathrm{{ex}}(\mathrm{{ex}})\pm )}(r)\);
-
\(a^2_{z(\mathrm{{ex}}(\mathrm{{ex}})\pm )}(r).\)
These functions determine the behavior of the characteristic functions related to the cosmological parameter, characterizing the functions \(q(r;y,a^2)\) and then the effective potentials on the higher level:
-
\(y_{h}(r;\,a^2)\);
-
\(y_\mathrm{d}(r;\,a^2)\);
-
\(y_{\mathrm{{ex}}(r)}(r;\,a^2)=y_{\mathrm{{ex}}(\mathrm{{ex}})}(r;\,a^2)\);
-
\(y_{\mathrm{d}(\mathrm{{ex}})}(r;\,a^2)\);
-
\(y_{z(\mathrm{{ex}})+}(r;\,a^2)\);
-
\(y_{\mathrm{{ex}}(\mathrm{{ex}})\pm }(r;\,a^2).\)
From the significance of the individual characteristic functions \(a^2(r)\) depicted in Fig. 5, one can infer that there are just two values of \(a^2\) of particular importance and leading to a qualitatively different behavior of the functions \(y(r;\,a^2)\):
-
\(a^2=1\)– the common local maximum of the functions \(a^2_{z(h)}(r)\) and \(a^2_{z(z(\mathrm{{ex}}))}(r)\) at \(r=1,\) which coincides with the inflection point of the function \(a^2_{z(\mathrm{{ex}}(\mathrm{{ex}})\pm )}(r;\,a^2)\) and with the intersection with the curve \(a^2_{\mathrm{{ex}}(h)+}(r)\);
-
\(a^2=a^2_{\text {crit}}=1.21202\)– the local maximum \(a^2_{\mathrm{{ex}}(h)+}(r)\) which is the intersection of the curves \(a^2_{r(\mathrm{{ex}}(\mathrm{{ex}})\pm )+}(r)\), \(a^2_{r(\mathrm{{ex}}(\mathrm{{ex}})\pm )}(r)\) and \(a^2_{{\text {max}}(\mathrm{d}(\mathrm{{ex}}))}(r).\)
The graphs of characteristic functions \(y(r;\,a^2)\) depicted for some values of spin parameter a representing the cases \(0<a^2<1,\) \(1<a^2<a^2_{\text {crit}}\) and \(a^2_{\text {crit}}<a^2\) are presented in Fig. 6.
In general, the behavior of the characteristic functions \(q_{r}(r;\,y,\,a^2)\) and \(q_\mathrm{{ex}}(r;\,y,\,a^2)\) will be qualitatively different, if for fixed parameter a we take the y-values from different intervals, which are limited by intersections and/or extrema of the characteristic functions \(y(r;\,a^2)\) that are demonstrated in Fig. 6. We therefore need to determine the curves \(y(a^2)\) that separate the \(a^2\)–y plane into regions that correspond to that different behavior of the characteristic functions \(q_{r}(r;\,y,\,a^2)\) and \(q_\mathrm{{ex}}(r;\,y,\,a^2)\). The number of these functions is substantially lowered by the fact that all the local extrema are multiple intersections with other curves and coincide with other extrema. Moreover, as explained below, the behavior of the characteristic functions \(q_{r}(r;\,y,\,a^2)\) and \(q_\mathrm{{ex}}(r;\,y,\,a^2)\) in their negative values we can omit as irrelevant for the character of the motion of the photon. The functions we need are the following:
-
\(y_{{\text {max}}(h)}(a^2)=y_{{\text {max}}(z(\mathrm{{ex}})+)}(a^2)=y_{{\text {inf}}(\mathrm{{ex}}(\mathrm{{ex}})-)}(a^2)= y_{\mathrm{d}(\mathrm{{ex}})\text{- }h\text{- }(z(\mathrm{{ex}})+)\text{- }\mathrm{{ex}}(\mathrm{{ex}})-}(a^2)\);
-
\(y_{{\text {min}}(h)}(a^2)=y_{{\text {min}}(z(\mathrm{{ex}})+)}(a^2)=y_{{\text {inf}}(\mathrm{{ex}}(\mathrm{{ex}})+)}(a^2)= y_{\mathrm{d}(\mathrm{{ex}})\text{- }h\text{- }(z(\mathrm{{ex}})+)\text{- }\mathrm{{ex}}(\mathrm{{ex}})+}(a^2)\);
-
\(y_{{\text {max}}(\mathrm{d})}(a^2)=y_{\mathrm{d}\text{- }\mathrm{d}(\mathrm{{ex}})\text{- }(z(\mathrm{{ex}})+)\text{- }\mathrm{{ex}}(\mathrm{{ex}})}(a^2)\);
-
\(y_{\mathrm{d}\text{- }(z(\mathrm{{ex}})+)}(a^2)\);
-
\(y_{{\text {max}}(\mathrm{d}(\mathrm{{ex}}))}(a^2)=y_{\mathrm{d}(\mathrm{{ex}})\text{- }\mathrm{{ex}}(\mathrm{{ex}})-}(a^2)\);
-
\(y_{\mathrm{{ex}}(\mathrm{{ex}})\text{- }(\mathrm{{ex}}(\mathrm{{ex}})+)}(a^2)\).
Here the dashes between two labels denote affiliation to intersection of appropriate functions (it can be proved that there are no other intersections of these functions than shown in Fig. 6). These functions are projections of the extremal values or intersections of characteristic functions \(y(r;\,a^2)\) into the \(a^2\)–y plane and they are demonstrated in Fig. 7.
The functions \(y_{\mathrm{{ex}}(h)}(a^2)\) divide the parameter plane \((a^2-y)\) into regions describing Kerr–de Sitter black-hole and naked singularity spacetimes, the curve \(y_{{\text {max}}(\mathrm{d})}(a^2)\) divides spacetimes with so-called divergent and restricted repulsive barrier of the motion of the photon. A detailed discussion of these functions has been presented e.g. in [73, 86] and will not be repeated here. The significance of the remaining functions can be understood from the depiction of the characteristic functions \(q_{r}(r;\,y,\,a^2),\) \(q_\mathrm{{ex}}(r;\,y,\,a^2)\) in Fig. 8. They are given parametrically by appropriate functions \(a^2(r),\) \(y(r;\,a^2(r))\) with r being the parameter:
-
the functions \(y_{{\text {max}}(h)}(a^2)\) and \(y_{{\text {min}}(h)}(a^2)\) are both determined by \(a^2_{\mathrm{{ex}}(h)+}(r)\) and \(y_h(r;\,a^2=a^2_{\mathrm{{ex}}(h)+}(r));\)
-
\(y_{{\text {max}}(\mathrm{d})}(a^2)\) we obtain from \(a^2_{{\text {max}}(\mathrm{d})}(r)\) with \(y_\mathrm{d}(r;\,a^2=a^2_{{\text {max}}(\mathrm{d})}(r));\)
-
the curve \(y_{{\text {max}}(\mathrm{d}(\mathrm{{ex}}))}(a^2)\) is given by functions \(a^2_{{\text {max}}(\mathrm{d}(\mathrm{{ex}}))}(r)\) and \(y_{\mathrm{d}(\mathrm{{\mathrm{{ex}}}})}(r;\,a^2=a^2_{{\text {max}}(\mathrm{d}(ex))}(r));\)
-
\(y_{\mathrm{d}\text{- }z(\mathrm{{ex}})+}(a^2)\) is determined by
$$\begin{aligned} a^2_{\mathrm{d}\text{- }z(\mathrm{{ex}})+}(r)\equiv \frac{r}{8}(1-4r+\sqrt{40r+1}) \end{aligned}$$and \(y_\mathrm{d}(r;\,a^2=a^2_{\mathrm{d}\text{- }z(\mathrm{{ex}})+}(r)),\) where the function \(a^2_{\mathrm{d}\text{- }z(\mathrm{{ex}})+}(r)\) is a solution of \( y_\mathrm{d}(r;\,a^2)=y_{z(\mathrm{{ex}})+}(r;\,a^2) \) with respect to parameter \(a^2;\) all such functions are obtained by analogous manner;
-
\(y_{\mathrm{{ex}}(\mathrm{{ex}})\text{- }\mathrm{{ex}}(\mathrm{{ex}})+}(a^2)\) are constructed from
$$\begin{aligned} a^2_{\mathrm{{ex}}(\mathrm{{ex}})\text{- } \mathrm{{ex}}(\mathrm{{ex}})\pm }(r) \end{aligned}$$$$\begin{aligned} \equiv \frac{r}{2}(4r^2{-}12r+3\pm \sqrt{16r^4-96r^3{+}156r^2-36r+9}) \end{aligned}$$and \( y_{\mathrm{{ex}}(\mathrm{{ex}})}(r;\,a^2=a^2_{\mathrm{{ex}}(\mathrm{{ex}})\text{- } \mathrm{{ex}}(\mathrm{{ex}})+}(r)).\)
There exist other functions \(y(a^2)\), corresponding to intersections of the characteristic functions \(y(r;\,a^2),\) which are not displayed in Fig. 7. The reason is that all the functions \(y(a^2)\) lie under the curve \(y=1/a^2,\) and thus we have to take into account the restriction \(q\ge -a^2\) (see Sect. 1). Therefore, the changes of the characteristic functions \(q_{r}(r;\,y,\,a^2),\) \(q_\mathrm{{ex}}(r;\,y,\,a^2)\) in the values under this limit can be omitted as irrelevant. Moreover, we can easily show that in the case \(q<0,\) the restrictions (53) imposed on the latitudinal motion yields stronger constraints on the value X than that given by Eqs. (59), (60), and (61) conditioning the reality of the radial motion. Indeed, for any triad \((q<0,y,a^2)\) there is no intersection of the curves \(X=X_\pm (r;\,q,\,y,\,a)\) with the lines \(X^\theta _{{\text {min}}(-)},\) \(X^\theta _{{\text {max}}(+)},\) where \(X^\theta _{{\text {min}}(-)},\) \(X^\theta _{{\text {max}}(+)}\) are extrema of the functions \(X^\theta (m;\,q,\,y,\,a)\) introduced in Sect. 1 (see Fig. 9e–g, \(\omega \)). To verify this, it is convenient to regard the curves \(X=X_\pm (r;\,q,\,y,\,a)\) as \(q={\text {const}}\)-slices of the surface \(q=q_{{\text {max}}}(r;\,X,\,a,\,y),\) where
is an alternative expression of the reality condition \(R(r;\,X,\,q,y,\,a)\ge 0,\), and to search instead for intersections of surfaces \(q=q_{{\text {max}}}(r;\,X,\,y,\,a)\) and \(q=q_{\text {min}}(X,\,y,\,a),\) defined by Eqs. (54)–(58).
We therefore solve the two equations, \(q_{{\text {max}}}(r;\,X,\,y,\,a^2)=q_1(X),\) with the result
and \(q_{{\text {max}}}(r;X,\,a,\,y)=q_2(X;\,a,\,y),\) which gives
The solution (96) yields \(X>0\), which, however, does not apply to the case \(q<0\) for \(y<1/a^2.\) Moreover, this solution represents touching points of the surface \(q_r(r;\,X,\,y,\,a^2)\) with the parabolic surface \(q_1(X)\) at \(X=+\sqrt{-q},\) and hence can be omitted even in the case \(y\ge 1/a^2,\) since theses values lie in the region forbidden by Eqs. (57)–(58). The solutions (97) are evidently irrelevant, since in stationary regions \(\Delta _r>0\) they are imaginary.
The above analysis shows that in the case \(q<0\) the ‘potentials’ \(X_\pm (r;\,q,\,y,\,a)\) have values in regions forbidden by reality conditions of the latitudinal motion and hence play no role at all. The limits for impact parameter X of photons with \(q<0\) are thus given by Eq. (53); photons satisfying Eq. (53) have thus no turning points of the radial motion. In the rest of this treatise we can thus focus on the behavior of the characteristic functions for \(q\ge 0.\) In Fig. 8 we present all possible variants of behavior of the characteristic functions \(q(r;\,y,\,a^2).\) These variants involve the cases:
-
I:
\(y \le y_{\mathrm{d}\text{- }z(\mathrm{{ex}})+}(a^2)\) for \(a^2\le 0.5\);
-
II:
\(y_{\mathrm{d}\text{- }z(\mathrm{{ex}})+}(a^2)\le y\le y_{{\text {max}}(\mathrm{d})}(a^2)\) for \(a^2\le 0.5,\) or \(y\le y_{{\text {max}}(\mathrm{d})}(a^2)\) for \(0.5 \le a^2\le 1,\) or \(y_{{\text {min}}(h)}(a^2)\le y\le y_{{\text {max}}(\mathrm{d})}(a^2)\) for \(1\le a^2\le 1.08316;\)
-
III:
\(y_{{\text {max}}(\mathrm{d})}(a^2)\le y \le y_{{\text {max}}(h)}(a^2)\) for \(a^2\le 1.08316,\) or \(y_{{\text {min}}(h)}(a^2)\le y \le y_{{\text {max}}(h)}(a^2)\) for \(1.08316\le a^2\le 1.21202=a^2_{\text {crit}};\)
-
IVa:
\(y\le y_{{\text {min}}(h)}(a^2)\) for \(1\le a^2\le 1.08316,\) or \(y\le y_{{\text {max}}(\mathrm{d})}(a^2)\) for \(1.08316\le a^2\le 1.28282,\) or \(y\le y_{(\mathrm{{ex}}(\mathrm{{ex}})+)\text{- }(\mathrm{{ex}}(\mathrm{{ex}})-)2}(a^2)\) for \(1.28282\le a^2\le 6\sqrt{3}-9=1.3923;\)
-
IVb:
\(y_{(\mathrm{{ex}}(\mathrm{{ex}})+)\text{- }(\mathrm{{ex}}(\mathrm{{ex}})-)2}(a^2)\le y\le y_{{\text {max}}(d)}(a^2)\) for \(1.28282\le a^2\le 1.3923,\) or \(y\le y_{{\text {max}}(\mathrm{d})}(a^2)\) for \(1.3923\le a^2\le 9,\) or \(y_{\mathrm{{ex}}(\mathrm{{ex}})\text{- }\mathrm{{ex}}(\mathrm{{ex}})+}(a^2)\le y\le y_{{\text {max}}(\mathrm{d})}(a^2)\) for \(a^2\ge 9;\)
-
V:
\(y\le y_{\mathrm{{ex}}(\mathrm{{ex}})\text{- }\mathrm{{ex}}(\mathrm{{ex}})+}(a^2)\) for \(a^2\ge 9;\)
-
VIa:
\(y_{{\text {max}}(\mathrm{d})}(a^2)\le y\le y_{{\text {min}}(h)}(a^2)\) for \(1.08316\le a^2\le 1.21202,\) or \(y_{{\text {max}}(\mathrm{d})}(a^2)\le y\le y_{(\mathrm{{ex}}(\mathrm{{ex}})+)\text{- }(\mathrm{{ex}}(\mathrm{{ex}})-)2}(a^2)\) for \(1.21202\le a^2\le 1.28282;\)
-
VIb:
\(y_{(\mathrm{{ex}}(\mathrm{{ex}})+)\text{- }(\mathrm{{ex}}(\mathrm{{ex}})-)2}(a^2)\le y\le y_{{\text {max}}(\mathrm{d}(\mathrm{{ex}}))}(a^2)\) for \(1.21202\le a^2\le 1.28282,\) or \(y_{{\text {max}}(\mathrm{d})}(a^2)\le y\le y_{{\text {max}}(\mathrm{d}(\mathrm{{ex}}))}(a^2)\) for \(a^2\ge 1.28282;\)
-
VII:
\(y_{{\text {max}}(h)}(a^2)\le y\le 1/a^2\) for \(a^2\le 1.21202,\) or \(y_{{\text {max}}(\mathrm{d}(\mathrm{{ex}}))}(a^2)\le y\le 1/a^2\) for \(a^2\ge 1.21202;\)
-
VIII:
\(y\ge 1/a^2.\)
Now it remains to assign to each region of the \(a^2\)–y plane functions \(q(y,\,a^2),\) which by themselves represent marginal values of the parameter q corresponding to some qualitative shift in the behavior of the potentials \(X_\pm (r;\,q,\,y,\,a).\)
In regions I and II, which describe black-hole spacetimes with divergent repulsive barrier, we have to compare the two local maxima \(q_{{\text {max}}(\mathrm{{ex}}+)}(y,\,a^2)\) located under the inner horizon and \(q_{{\text {max}}(\mathrm{{ex}})}(y,\,a^2)=q_{{\text {min}}(r)}(y,\,a^2)\) between the outer and cosmological horizons, respectively (see Fig. 8a, b). The function \(q_{{\text {max}}(\mathrm{{ex}}+)}(y,\,a^2)\) is given parametrically by the functions \(y_{\mathrm{{ex}}(\mathrm{{ex}})+}(r;\,y,\,a^2)\) and \(q_{\mathrm{{ex}}}(r;\,y=y_{\mathrm{{ex}}(\mathrm{{ex}})+}(r;\,y,\,a^2),\,a^2)\) with r being the parameter; similarly \(q_{{\text {max}}(\mathrm{{ex}})}(y,\,a^2)\) is given by \(y_{\mathrm{{ex}}(\mathrm{{ex}})}(r;\,y,\,a^2)\) and \(q_{\mathrm{{ex}}}(r;\,y=y_{\mathrm{{ex}}(\mathrm{{ex}})}(r;\,y,\,a^2),\,a^2).\)
The extrema function \(q_{{\text {max}}(\mathrm{{ex}})}(y,\,a^2)\) diverges at the curve \(y_{{\text {max}}(\mathrm{d})}(a^2)\), which forms the boundary between regions II–III and IV–VI, i.e. \(q_{{\text {max}}(ex)}(y=y_{{\text {max}}(\mathrm{d})}(a^2),\,a^2) \rightarrow +\infty \) (cf. Fig. 8b–f). In region III, corresponding to black-hole spacetimes with the restricted repulsive barrier, only local maximum \(q_{{\text {max}}(\mathrm{{ex}}+)}(y,\,a^2)\) located under the inner horizon remains.
As can be seen from the behavior of the characteristic functions in Fig. 6b, for \(y \rightarrow y_{{\text {min}}(h)}(a^2)\) from above, the ‘inner’ local maximum of \(q_{(\mathrm{{ex}})}(r;\,y,\,a^2),\) determined by \(y_{\mathrm{{ex}}(\mathrm{{ex}})+}(r;\,y,\,a^2),\) approaches from the left its divergency point given by \(y=y_{\mathrm{d}(\mathrm{{ex}})}(r;\,a^2),\) where \(q_{\mathrm{{ex}}}\rightarrow -\infty \) (Fig. 8b), so that \(q_{ex}(r;\, y=y_{{\text {min}}(h)}(a^2),\,a^2)\) becomes continuous. For \(y \le y_{{\text {min}}(h)}(a^2),\) the divergency of the function \(q_\mathrm{{ex}}(r;\,y,\,a^2)\) appears again with \(q_{\mathrm{{ex}}}\rightarrow +\infty \) and a local minimum has formed on the right (cf. Fig. 8b, d). Hence the curve \(y_{{\text {min}}(h)}(a^2)\) forms a boundary on which the local maxima \(q_{{\text {max}}(\mathrm{{ex}}+)}(y,\,a^2)\) convert into local minima. We denote them \(q_{{\text {min}}(\mathrm{{ex}}\pm )}(y,\,a^2),\) since, as follows from Eqs. (86)–(91), for \(1.125\le a^2 \le a^2_{\text {crit}}\) and
or \(a^2_{\text {crit}} \le a^2 \le 1.3923\) and
they are given by \(y_{\mathrm{{ex}}(\mathrm{{ex}})-}(r;\,y,\,a^2).\) The function \(y_{(\mathrm{{ex}}(\mathrm{{ex}})+)\text{- }(\mathrm{{ex}}(\mathrm{{ex}})-)1}(a^2)\) is given parametrically by \(a^2_{r(\mathrm{{ex}}(\mathrm{{ex}})\pm )}(r)\) and, e.g., by \(y_{\mathrm{{ex}}(\mathrm{{ex}})-}(r;\,y,\,a^2=a^2_{r(\mathrm{{ex}}(\mathrm{{ex}})\pm )}(r)),\) and the function \(y_{(\mathrm{{ex}}(\mathrm{{ex}})+)\text{- }(\mathrm{{ex}}(\mathrm{{ex}})-)2}(a^2)\) by \(a^2_{r(\mathrm{{ex}}(\mathrm{{ex}})\pm )+}(r)\) and \(y_{\mathrm{{ex}}(\mathrm{{ex}})-}(r;\,y,\,a^2=a^2_{r(\mathrm{{ex}}(\mathrm{{ex}})\pm )+}(r)).\) The analytical expressions in Eqs. (98) and (99) can then be derived by eliminating the radius r. Both these functions have their relevant parts entirely in regions IV and VI.
The function \(y_{(\mathrm{{ex}}(\mathrm{{ex}})+)\text{- }(\mathrm{{ex}}(\mathrm{{ex}})-)2}(a^2)\) corresponds to the local minima of the potential \(X_{-}(r;\,q,\,y,\,a)\) reaching the value \(X=-a,\) i.e., \(\ell =0.\) The photons corresponding to these minima, and having appropriate constant of the motion q, persist on ’spherical’ orbits with \(r={\text {const}}\), which are crossing the spacetime rotation axis alternately above both poles. Below we shall call them ’polar’ spherical orbits – in the following section we shall see that such polar spherical orbits form a border surface between prograde and retrograde spherical photon orbits, as related to the locally non-rotating observers. The function \(y_{(\mathrm{{ex}}(\mathrm{{ex}})+)\text{- }(\mathrm{{ex}}(\mathrm{{ex}})-)2}(a^2)\) has therefore an important meaning, since it represents a boundary between regions of qualitatively different KdS spacetimes in the \(a^2\)–y plane. From this point of view, the function \(y_{(\mathrm{{ex}}(\mathrm{{ex}})+)\text{- }(\mathrm{{ex}}(\mathrm{{ex}})-)2}(a^2)\) creates another qualitative shift in the parameter plane \((a^2-y)\) with regard to the character of the motion of the photon, however, no qualitative shift in the mathematical properties of the characteristic functions \(q_{r}(r;\,y,\,a^2)\) and \(q_{\mathrm{{ex}}}(r;\,y,\,a^2)\) in their relevant values \(q\ge 0.\) The parts of the \(a^2\)–y plane corresponding to different behavior of the characteristic functions are in Fig. 7 distinguished by Roman numerals, the curve \(y_{(\mathrm{{ex}}(\mathrm{{ex}})+)\text{- }(\mathrm{{ex}}(\mathrm{{ex}})-)2}(a^2)\) then induces an additional division a / b.
Further we have to relate the minima \(q_{{\text {min}}(\mathrm{{ex}}\pm )}(y,\,a^2)\) with the maxima \(q_{{\text {max}}(\mathrm{{ex}})}(y,\,a^2)=q_{{\text {min}}(r)}(y,\,a^2).\) In the region V, the minima of \(q_{\mathrm{{ex}}}(r;\,y,\,a^2)\) coalesce with the minima of \(q_{r}(r;\,y,\,a^2)\) (Fig. 8e). We therefore have to compare the minima function \(q_{{\text {min}}(\mathrm{{ex}})}(y,\,a^2)= q_{{\text {min}}(r)}(y,\,a^2)\) determined by \(y_{\mathrm{{ex}}(\mathrm{{ex}})}(r;\,y,\,a^2)\) and, e.g., \(q_{\mathrm{{ex}}}(r;\,y=y_{\mathrm{{ex}}(\mathrm{{ex}})}(r;\,y,\,a^2),\,a^2),\) with the maxima function \(q_{{\text {max}}(\mathrm{{ex}}+)}(y,\,a^2)\) parametrized by \(y_{\mathrm{{ex}}(\mathrm{{ex}})+}(r;\,y,\,a^2)\) and \(q_{\mathrm{{ex}}}(r;\,y=y_{\mathrm{{ex}}(\mathrm{{ex}})+}(r;\,y,\,a^2),\,a^2).\) The boundary curve \(y_{\mathrm{{ex}}(\mathrm{{ex}})\text{- }\mathrm{{ex}}(\mathrm{{ex}})+}(a^2)\) then represents such combinations of parameters \(a^2,y\) for which the local extrema of \(q_{\mathrm{{ex}}}(r;\,y,\,a^2)\) have coalesced into an inflection point. For parameters from region VI, corresponding to naked singularity spacetimes with restricted repulsive barrier (as well as from the remaining regions VII, VIII), the function \(q_{\mathrm{{ex}}}(r;\,y,\,a^2)\) has one local minimum (Fig. 8f), and we therefore construct a function \(q_{{\text {min}}(\mathrm{{ex}}\pm )}(y,\,a^2)\) determined by functions \(y_{\mathrm{{ex}}(\mathrm{{ex}})\pm }(r;\,a^2)\) and \(q_{\mathrm{{ex}}}(r;\,y=y_{\mathrm{{ex}}(\mathrm{{ex}})\pm }(r;\,a^2),\,a^2),\) where the minus sign has to be chosen for \(1.17007 \le a^2 \le a^2_{\text {crit}}\) and \(y_{{\text {max}}(\mathrm{d})}(a^2)\le y \le y_{(\mathrm{{ex}}(\mathrm{{ex}})+)\text{- }(\mathrm{{ex}}(\mathrm{{ex}})-)1}(a^2),\) or \(a^2_{\text {crit}} \le a^2 \le 1.2828\) and \(y_{{\text {max}}(\mathrm{d})}(a^2) \le y \le y_{(\mathrm{{ex}}(\mathrm{{ex}})+)\text{- }(\mathrm{{ex}}(\mathrm{{ex}})-)2}(a^2)\) (see Fig. 7). For \(y\rightarrow y_{{\text {max}}(\mathrm{d}(\mathrm{{ex}}))}(a^2)\) and \(a^2\ge a^2_{\text {crit}}\) we have \(q_{{\text {min}}(\mathrm{{ex}}+)}(y,\,a^2)\rightarrow +\infty ,\), and for \(y> y_{{\text {max}}(\mathrm{d}(\mathrm{{ex}}))}(a^2),\) i.e. in region VII, it converts into the local maximum (cf. Fig. 8f, g). The transition into region VII from region III can be inferred from a comparison of Fig. 8c with g. Therefore, in region VII we have to follow up the values of the function \(q_{{\text {max}}(\mathrm{{ex}}+)}(y,\,a^2)\) determined by the functions \(y_{\mathrm{{ex}}(\mathrm{{ex}})+}(r;\,a^2)\) and \(q_{\mathrm{{ex}}}(r;\,y=y_{\mathrm{{ex}}(\mathrm{{ex}})+}(r;\,a^2),\,a^2).\) The functions \(q(y,\,a^2)\) are demonstrated in Fig. 9.
With the knowledge of the behavior of the extremal values \(q_{{\text {min}}/{\text {max}}(\mathrm{{ex}})}(y,\,a^2)\) at each region of the \(a^2\)–y plane, we can finally construct all qualitatively different types of the behavior of the effective potentials \(X_{\pm }(r;\,q,\,y,\,a).\) They are presented in Fig. 3 for appropriately chosen representative combinations \((q,y,a^2).\)
5 Spherical photon orbits and classification of the Kerr–de Sitter spacetimes due to properties of the motion of the photon
We shall demonstrate by using the behavior of the effective potentials \(X_{\pm }(r;\,q,\,a,\,y)\) that the null geodesics create qualitatively different structures in the various cases of Kerr–de Sitter spacetimes with the spacetime parameters chosen from different parts of the \(a^2\)–y plane labeled by numerals I–VIII. Hence the regions of the spacetime parameter space of these labels can be considered as representatives of the classification of the Kerr–de Sitter spacetimes due to the motion of the photon (null geodesics). Similarly to [73], there are three (four) criteria used – the main criterion for the classification is the existence (number) of the event horizons. The other differentiating factors follow from the nature of the motion of the photon. First, there is some kind of repulsive barrier preventing light from reaching the ring singularity, which is always created in its vicinity for photons with \(q>0\). However, a similar barrier can emerge between the outer black-hole horizon and the cosmological horizon in black-hole and naked singularity spacetimes, repelling photons towards one of these horizons. In the naked singularity spacetimes, the occurrence of an additional barrier, which reflects photons towards the ring singularity, leads to the occurrence of the phenomenon of bound photon orbits. Such bound photon orbits are not present in the case of the black-hole spacetimes. The presence and character of this barrier we take as another criterion in the following classification. The other aspect that authorizes us to make such a distinction between the KdS spacetimes will be the existence and character of the spherical photon orbits. In the KdS naked singularity spacetimes the bound orbits are concentrated around the stable spherical photon orbits.
5.1 Spherical photon orbits
The spherical photon orbits are determined by the conditions \(R(r)=0\) and \(\mathrm{d}R/\mathrm{d}r=0\), which have to be solved simultaneously. The physically acceptable solution is governed by the relations for the constants X and q for the motion of the photon that are expressed as functions of the radius r and the spacetime parameters a, y, and they take the form
where the function \(q_{\mathrm{{ex}}}(r;\,y,\,a^2)\) is defined by Eq. (75). These solutions governing the spherical photon orbits are allowed in the interval of radii limited by the equatorial photon circular orbits. The stability of the spherical photon orbits relative to radial perturbations is determined by the sign of the expression
evaluated at appropriate radii. It can be shown that the local maxima of the potential \(X_{+}\) and local minima of \(X_{-}\) correspond to unstable orbits in the black-hole spacetimes. However, the local minima of \(X_{+}\) and local maxima of \(X_{-}\) represent stable orbits, which occur in the naked singularity spacetimes.
It is useful to relate the parameters (wave-vector components) of photons orbiting along the spherical null geodesics to the locally non-rotating frames (LNRFs), which are the most convenient frames for a description of physical processes in the Kerr (dS) spacetimes [9]. The LNRF tetrad of differential one-forms is given by the relations
the corresponding tetrad of the dual vectors reads
The wave-vector components related to the LNRFs are then determined by the relations
hence
where
and
is the angular velocity of the LNRFs related to distant static observers.
In order to determine the orientation of the spherical orbits, we have chosen as the azimuthal direction indicator the sign of the ratio \(k^{(\phi )}/k^{(t)}.\) If we define the directional angle \(\Psi \) in such a way that \(\Psi =0\) for the motion in the direction of the latitudinal tetrad vector \(e_{(\theta )},\) while \(\Psi =\pi /2\) for motion in the direction of the azimuthal tetrad vector \(e_{(\phi )},\) then \(k^{(\phi )}/k^{(t)}=\sin \Psi \) and we find the relation
If the sign of \(\mathrm{{sin}}\Psi \) is positive, we call the spherical orbit prograde, if it is negative, we call the spherical orbit retrograde. The special case of limiting spherical orbits corresponds to the equatorial circular orbits that are again co-rotating (prograde), respectively, counter-rotating (retrograde). It can be shown that the sign of the directional angle remains fixed at any latitude of any particular spherical orbit, i.e., the locally non-rotating observers see the motion of the photon in a fixed azimuthal direction.Footnote 1 Since the functions \(A, \Omega _{\text {LNRF}}\) are positive [73], it is clear that all photons with \(X<-a\) (\(\ell <0\)) are retrograde. However, photons with
can be retrograde as well. Considering in such a case the relation for the tetrad LNRF component (112), we can see that in order to keep for \(k^{(t)}\) the standard physical meaning, i.e., \(k^{(t)}>0,\) we have to put \(E<0.\) Footnote 2 In order to find the conditions under which such a situation occurs, it is convenient to express by the alternate relation
the impact parameter \(\ell \) in the form
and reverse the problem by searching for conditions under which \(\ell >0.\) Such a relation is evidently fulfilled if \(\sin \psi >0,\) i.e., the positive impact parameters pertain to prograde photons. However, there is another possibility, namely to have \(\sin \psi <0\) together with
from which it follows that
However, the last inequality can be written in the form
which implies
Hence, such a situation can occur only in the ergosphere. Of course, the impact parameter of such photons must fulfill the condition (119). The function
has for \(\Delta _{r}=0\) common points with the potentials \(X_{\pm }\) given by Eq. (64). There are no other intersections with the potentials, hence, the reality condition of the radial motion together with (119) imply \(X>X_{+}>0.\) Therefore, the motion of photons with negative energy E, which appear to be retrograde in the LNRFs, is governed by effective potentials \(X_{+}(r;\,q,\,a,\,y)\) with positive values.
Using the properties of the effective potentials \(X_{\pm }(r;\,q,a,\,y)\), we can identify the radii \(r=r_{0}\) of the spherical photon orbits as loci of the local extrema of the effective potentials and determine their stability and orientation as described above. At each allowed radius \(r_{0}\), located between the radii of the equatorial photon circular orbits, we can assign corresponding limits \(\theta _{\text {min}},\) \(\theta _{{\text {max}}}\) on the latitudinal motion by solving the equation
which due to the results of Sect. 3 has one real positive root \(m_{0},\) since \(q_{\text {sph}}(r_{0})\ge 0\) (\(q_{\text {sph}}(r_{0})=0\) for \(r_{0}=r_{{\text {ph}}\pm },\) i.e., equatorial circular co-rotating or counter-rotating photon orbit). The marginal latitudes (turning points of the latitudinal motion) then read
for details see the discussion of the latitudinal motion in Sect. 3. We can thus easily determine for a spherical orbit at an allowed radius \(r_0\) the impact parameters of the orbit, the extension of the latitudinal motion, and the orientation of the azimuthal motion.
5.2 Classification
In the following classification we introduce ten classes of the Kerr–de Sitter spacetimes and demonstrate the properties of the motion of the photon using the spherical photon orbits that serve as crucial characteristic for the classification. We give the loci of the spherical photon orbits and their extension in latitude, stability against radial perturbations, and orientation of their azimuthal motion. The classification is represented by a family of characteristic figures corresponding to the separated classes of the KdS spacetimes. For easy interpretation of the family of the figures representing the classification, we introduce an auxiliary figure, Fig. 11, with detailed explanatory notes. In order to fully and clearly characterize the KdS spacetimes and their horizon and ergosphere structure, and to demonstrate the spheroidal character of the applied coordinate system, we now use the so-called Kerr–Schild coordinates x, y, z, which are connected to the Boyer–Lindquist coordinates \(r,\theta \) by the relations
In the figures we, of course, use the meridional sections of \(y=0\). The characteristics of the classes of the KdS spacetime according to the photon orbits are presented as follows.
-
Class I:
Black-hole spacetimes with the divergent repulsive barrier of the radial motion of the photon, having one equatorial counter-rotating circular unstable orbit with negative energy located under the inner black-hole horizon (\(0<r<r_{-}\)), which is limiting the range of the spherical photon orbits with negative energy. There exist stable orbits, corresponding to local minima \(X_{{\text {min}}(+)}\) of the effective potential \(X_{+}\) at \(0<r<r_{{\text {max}}(\mathrm{{ex}})1},\) and unstable orbits, corresponding to local maxima \(X_{{\text {max}}(+)}\) of \(X_{+}\) at \(r_{{\text {max}}(\mathrm{{ex}})1}<r<r_{z(\mathrm{{ex}})1}\) for \(0<q<q_{{\text {max}}(\mathrm{{ex}}+)}(y,\,a^2)\) (Fig. 12a). We denote by \(r_{{\text {min}}/{\text {max}}(\mathrm{{ex}})}\) the local extrema, and by \(r_{z(\mathrm{{ex}})}\) the zero point, of the function \(q_{ex}(r;\,y,\,a^2)\) hereafter. Such a structure is present under the inner horizon of any KdS black-hole spacetime. Outside the ergosphere, one unstable co-rotating equatorial circular orbit, located at \(r=r_{{\text {ph}}+}=r_{z(\mathrm{{ex}})2}\), and a polar spherical orbit with \(r=r_{\text {pol}},\) \(r_{{\text {ph}}+}<r_{\text {pol}},\) limit the range of unstable prograde spherical orbits given by the local minima \(X_{{\text {min}}(-)}\) of the effective potential \(X_{-},\) for which \(X_{{\text {min}}(-)}>-a\). The radius of the polar spherical orbit is found by solving \(X_{-}(r_{\text {pol}};q_{\mathrm{{ex}}}(r_{\text {pol}}))=-a.\) The counter-rotating equatorial circular orbit at \(r=r_{\mathrm{{ph}}-}=r_{z(\mathrm{{ex}})3}\) gives the limit of the region of unstable retrograde spherical orbits, given by the local minima \(X_{{\text {min}}(-)}<-a\), and maxima \(X_{{\text {max}}(+)}<-a\) for \(0<q<q_{{\text {max}}(\mathrm{{ex}})}(y,\,a^2),\) such that \(r_{\text {pol}}<r_{{\text {ph}}-}.\)
-
Class II:
Black-hole spacetimes with the same features as in the class I, but now the ergosphere enters the region of the spherical photon orbits (Fig. 12b). No spherical orbit is fully immersed in the ergosphere and photons at all the spherical orbits have positive energy. The presence of the ergosphere in region of the spherical photon orbits influences character of the light escape cones [18].
-
Class III:
Black-hole spacetimes with the restricted repulsive barrier of the radial motion of the photon. The ergosphere spreads over all radii. The prograde spherical orbits are given by the local minima \(X_{{\text {min}}(-)}>-a\) at \(r_{{\text {ph}}+}<r<r_{\text {pol}},\) while the retrograde spherical orbits with \(E>0\) are given by the minima \(X_{{\text {min}}(-)}<-a\) at \(r_{\text {pol}}<r<r_{\mathrm{d}(\mathrm{{ex}})}.\) The spherical orbits given by the local maxima \(X_{{\text {max}}(+)}\) (see Fig. 10k–n) at \(r_{\mathrm{d}(\mathrm{{ex}})}<r<r_{{\text {ph}}-},\) where \(r_{\mathrm{d}(\mathrm{{ex}})}\) denotes the divergence point of \(q_{\mathrm{{ex}}}(r;\,y,\,a^2)\) (see Fig. 8c), are fully immersed in the ergosphere. Such areas are drawn in Fig. 12 in green and the spheres with \(r=r_{\mathrm{d}(\mathrm{{ex}})}\) as full/dashed green ellipses. Photons in such regions have \(E<0\).
-
Class IVa:
Naked singularity spacetimes with divergent repulsive barrier of the motion of the photon. At radii \(0<r<r_{\mathrm{d}(\mathrm{{ex}})}\) (Fig. 8d can be used for illustration), there are local minima of the potential \(X_{+}\) (for illustration use Fig. 10p–r) corresponding to the stable retrograde spherical orbits with negative energy (\(E<0\)) (Fig. 12d). The stable retrograde orbits with positive energy corresponding to the local maxima of \(X_{-}\) are at \(r_{\mathrm{d}(\mathrm{{ex}})}<r<r_{{\text {pol}}1}\). These maxima exceed the value \(X_{{\text {max}}(-)}=-a\) at \(r_{{\text {pol}}1}<r<r_{{\text {min}}(\mathrm{{ex}})},\) where they yield stable prograde spherical orbits. At radii \(r_{{\text {min}}(\mathrm{{ex}})}<r<r_{pol2}\), we have the local minima of \(X_{-}\) with values \(X_{{\text {min}}(-)}>-a,\) – these radii are thus occupied by the unstable prograde orbits. The local minima of \(X_{-}\) and the local maxima of \(X_{+}\) at \(r>r_{{\text {pol}}2}\) correspond to the unstable retrograde orbits. There are thus two polar spherical orbits enclosing the region of prograde orbits – the inner at the radius \(r=r_{{\text {pol}}1}\) being stable, the outer at the radius \(r=r_{{\text {pol}}2}\) being unstable.
-
Class IVb:
Naked singularity spacetimes with the same features as in the class IVa, but the two polar orbits have coalesced, therefore, there are no prograde spherical orbits (Fig. 12e).
-
Class V:
Naked singularity spacetimes having the structure of the spherical orbits corresponding to the previous case (Fig. 12f), but with is a small region of bound orbits for photons with constants of the motion \(q_{{\text {min}}(\mathrm{{ex}})}(y,\,a^2)<q<q_{{\text {max}}(\mathrm{{ex}}+)}(y,\,a^2)\) and X between the appropriate local extrema of \(X_{+}\) (cf. Fig. 10q, u), which is not contained in the other cases.
-
Class VIa:
Naked singularity spacetimes with the restricted repulsive barrier of the radial motion of the photon. For \(0<r<r_{\mathrm{d}(\mathrm{{ex}})1}\) (the function \(q_\mathrm{{ex}}(r;\,y,\,a^2)\) has two divergence points \(r_{\mathrm{d}(\mathrm{{ex}})1}, r_{\mathrm{d}(\mathrm{{ex}})2}\) – see Fig. 8f) the minima of \(X_{+}\) correspond to stable retrograde orbits with \(E<0\); for \(r_{\mathrm{d}(\mathrm{{ex}})1}<r<r_{{\text {pol}}1}\), there are the local maxima of \(X_{-}\) with values \(X_{{\text {max}}(-)}<-a\) giving retrograde orbits with \(E>0.\) The local minima of \(X_{-}\) at \(r_{{\text {pol}}1}<r<r_{{\text {min}}(\mathrm{{ex}})}\) give the stable prograde orbits. At radius \(r=r_{{\text {pol}}1}\) the stable polar orbit is located. For \(r_{{\text {min}}(\mathrm{{ex}})<r<r_{{\text {pol}}2}}\), the function \(X_{-}\) has minima with values \(X_{{\text {min}}(-)}>-a\), giving the unstable prograde orbits. For \(r_{{\text {pol}}2}<r<r_{\mathrm{d}(\mathrm{{ex}})2}\), they correspond to the unstable retrograde orbits. At radius \(r=r_{{\text {pol}}2}\), the unstable polar orbit exists. The local maxima of the function \(X_{+}\) at \(r_{\mathrm{d}(\mathrm{{ex}})2<r<r_{{\text {ph}}-}}\) correspond to the retrograde unstable spherical orbits with \(E<0.\)
-
Class VIb:
The structure of the spherical orbits corresponds to the class VIa with the exception that the local extrema of the potential \(X_{-}\) have values \(X<-a\), implying that there are no polar spherical orbits, nor the prograde spherical orbits (Fig. 12h).
-
Class VII:
Naked singularity spacetimes with the restricted repulsive barrier of the radial motion of the photon having stable retrograde spherical orbits at \(0<r_{{\text {max}}(\mathrm{{ex}})}\) corresponding to local minima of \(X_{+}\) (Fig. 10 \(\alpha \)), and unstable retrograde spherical orbits at \(r_{{\text {max}}(\mathrm{{ex}})}<r<r_{{\text {ph}}-}\) corresponding to local maxima of \(X_{+}.\) All these spherical orbits, including the counter-rotating equatorial circular orbit at \(r=r_{{\text {ph}}-}\), correspond to photons with \(E<0.\)
-
Class VIII:
Special class of the naked singularity spacetimes demonstrating the same features of the radial motion of photons with \(q\ge 0\) as the class VII, but differing from all previous cases by the existence of null geodesics for arbitrary \(q<0.\) The allowed values of the impact parameter X are for \(q<0\) confined to the intervals \(X<X^\theta _{{\text {max}}(+)}<0\) or \(X>X^\theta _{{\text {min}}(-)}>0\) (see Sect. 3). The potentials governing the radial motion of the photon are fully immersed in the forbidden region (Fig. 10 (\(\gamma \)) and (\(\delta \))), thus in the radial direction the photons with such parameters move freely in the whole range between the ring singularity and the cosmological horizon.
6 Conclusions
We can summarize our results by the following concluding remarks.
-
1.
In any kind of the black-hole spacetimes, there are no radially bound null geodesics in the stationary region, i.e., the trajectory of a photon has at most one turning point in radial direction between the outer and cosmological horizon, or the photons can move freely between the outer black-hole and the cosmological horizons. However, such bounded photon orbits exist in each naked singularity spacetime for photons with parameters \(q>0\) and X chosen appropriately.
-
2.
No photons with \(q>0\) can reach the ring singularity at \(r=0\) in any of the Kerr–de Sitter spacetimes.
-
3.
In the Kerr–de Sitter spacetimes of classes I–VII, i.e., with the spacetime parameters satisfying the condition \(y<1/a^2,\) there is a lower limit \(q=-a^2\) of the parameter \(q<0,\) for which the motion of the photon is allowed. The range of the allowed values of the impact parameter X is then an interval given by Eqs. (53)–(56). Photons with such tuned parameters have no turning point in radial direction, since the effective potential lies entirely in the forbidden region (Fig. 10e–g). Furthermore, by the results of Sect. 3, only such photons execute the vortical motion, or their trajectory lies completely on the cones of \(\theta = {\text {constant}}\). We can therefore reject the possibility of the existence of a vortical motion of the photon of constant radius, or off-equatorial circular photon orbits.
-
4.
In the Kerr–de Sitter spacetimes of class VIII (\(y>1/a^2\)), the motion of the photon is allowed for any \(q<0.\) The permissible values of the parameter X are then two disjunct unlimited intervals determined by Eq. (58). In the extreme case \(y=1/a^2,\) we must have \(q\ge -a^2\) again and for negative q, the parameter X can take a smaller value than a certain negative value, given by (57). The consequences for the motion of the photon are then the same as in previous note (Fig. 10 (\(\gamma \)) and (\(\delta \))).
-
5.
In the Kerr–de Sitter spacetimes with the divergent repulsive barrier of the radial motion of the photon, there exists a critical value \(q_{{\text {max}}(\mathrm{{ex}})}(y,\,a^2),\) for which this barrier becomes impermeable between the outer black-hole horizon and cosmological horizon, or, in naked singularity spacetimes, between the ring singularity and cosmological horizon, for photons with any impact parameter X. In spacetimes with the restricted repulsive barrier of the radial motion of the photon, the height of this barrier slowly grows with increasing parameter q, but it stays finite for any \(q>0\) (Fig. 10n, y).
-
6.
In the Kerr–de Sitter spacetimes of classes I–III, IVa and VIa, there exist spherical photon orbits, which can be both prograde or retrograde as seen by the family of locally non-rotating observers. Additionally, each of the two types can be stable or unstable with respect to radial perturbations. The regions of spherical orbits of different orientations are separated by the so-called polar spherical orbit, at which photons cross the spacetime rotation axis alternately above both poles. In the naked singularity spacetimes of class IVa, VIa, there are two polar spherical orbits, the inner one being stable, the outer one being unstable.
-
7.
In the Kerr–de Sitter spacetimes of classes IVb, V, VIb, VII, VIII there are no prograde or polar spherical orbits.
-
8.
In each class of the Kerr–de Sitter spacetimes, there exists a region where the effective potential \(X_{+}\) have positive values. Photons with impact parameter X exceeding these values appear to move in retrograde direction as seen in LNRFs. This region must be located inside the ergosphere, and photons with such impact parameters must have negative energy, \(E<0\). In the black-hole spacetimes with the divergent barrier of the radial motion of the photon, the ergosphere has two parts above the black-hole outer event horizon – the inner one, which is limited to the outer vicinity of the outer event horizon, and the outer one, limited to the inner vicinity of the cosmological horizon. In the black-hole spacetimes with restricted repulsive barrier of the radial motion of the photon, the two regions of the ergosphere merge in the equatorial plane, and they spread at any radii except for a certain region in the vicinity of the rotation axis.
-
9.
In the LNRFs, trajectories of photons moving along any spherical orbit have no turning point of the azimuthal motion.
We have thus demonstrated a variety of very extraordinary phenomena related to the motion of the photon in the KdS spacetimes, of both black-hole and naked singularity types. Especially relevant effects are found in the case of spherical photon orbits that can be directly related to the observational phenomena. It is quite interesting that we could expect another interesting phenomenon related to the charged Kerr–Newman or Kerr–Newman–de Sitter naked singularity spacetimes (with both the standard electric charge and the tidal charge of the braneworld models), especially in the case of the so-called mining Kerr–Newman spacetimes [12], containing a special type of equatorial stable photon orbits.
Notes
However, we have to note that, similarly to the case of Kerr black holes [90], the sign of the variation of the azimuthal coordinate can be changed at some latitude, if related to distant observers.
Keeping \(E>0\) means \(k^{(t)}<0,\) i.e., a photon in negative-root state with time evolution directed to the past – for details see [11].
References
C. Adami, F. Durret, L. Guennou, Rocha C. Da, Diffuse light in the young cluster of galaxies CL J1449\(+\)0856 at \(z=2.07\). Astron. Astrophys. 551, A20 (2013)
A.N. Aliev, Electromagnetic properties of Kerr–anti-de Sitter black holes. Phys. Rev. D 75(8), 084041 (2007)
C. Armendariz-Picon, V. Mukhanov, P.J. Steinhardt, Dynamical solution to the problem of a small cosmological constant and late-time cosmic acceleration. Phys. Rev. Lett. 85(21), 4438 (2000)
I. Arraut, Komar mass function in the de Rham-Gabadadze-Tolley nonlinear theory of massive gravity. Phys. Rev. D 90, 124082 (2014)
I. Arraut, The astrophysical scales set by the cosmological constant, black-hole thermodynamics and non-linear massive gravity. Universe. 3(2), 45 (2017)
B. Aschenbach, Measurement of mass and spin of black holes with QPOs. Chin. J. Astron. Astrophys. Suppl. 8, 291–296 (2008)
N. Bahcall, J.P. Ostriker, S. Perlmutter, P.J. Steinhardt, The cosmic triangle: revealing the state of the universe. Science 284, 1481–1488 (1999)
P. Bakala, P. Čermák, S. Hledík, Z. Stuchlík, K. Truparová, Extreme gravitational lensing in vicinity of Schwarzschild–de Sitter black holes. Cent. Eur. J. Phys. 5(4), 599–610 (2007)
J.M. Bardeen, Timelike and null geodesics in the Kerr metric, in Black Holes (Les Astres Occlus), ed. by C. Dewitt, B.S. Dewitt (1973), pp. 215–239
J. Bičák, Z. Stuchlík, On the latitudinal and radial motion in the field of a rotating black hole. Bull. Astronom. Inst. Czechoslov. 27(3), 129–133 (1976)
J. Bičák, Z. Stuchlík, V. Balek, The motion of charged particles in the field of rotating charged black holes and naked singularities. Bull. Astron. Instit. Czechoslov. 40, 65–92 (1989)
M. Blaschke, Z. Stuchlík, Efficiency of the Keplerian accretion in braneworld Kerr–Newman spacetimes and mining instability of some naked singularity spacetimes. Phys. Rev. D 94, 086006 (2016)
C.G. Böhmer, Eleven spherically symmetric constant density solutions with cosmological constant. Gen. Relativ. Gravit. 36, 1039–1054 (2004)
E.K. Boyda, S. Ganguli, P. Hořava, U. Varadarajan, Holographic protection of chronology in universes of the Gödel type. Phys. Rev. D 67, 106003 (2003)
R. Caldwell, M. Kamionkowski, Cosmology: dark matter and dark energy. Nature 458(7238), 587–589 (2009)
R.R. Caldwell, R. Dave, P.J. Steinhardt, Cosmological imprint of an energy component with general equation of state. Phys. Rev. Lett. 80(8), 1582 (1998)
Carter, B: Black hole equilibrium states, in Black Holes (Les Astres Occlus), ed. by C. Dewitt, B.S. Dewitt (Gordon @ Breach Science Publishers, New York, 1973), pp. 57–214
D. Charbulák, Z. Stuchlík, Light escape cones in Kerr–de Sitter black hole spacetimes in special local reference frames Europ. Phy. J. C. (under review)
J.-H. Chen, Y.-J. Wang, Influence of dark energy on time-like geodesic motion in Schwarzschild spacetime. Chin. Phys. B 17(4), 1184 (2008)
N. Cruz, M. Olivares, J.R. Villanueva, The geodesic structure of the Schwarzschild anti-de Sitter black hole. Class. Quantum Gravity 22(6), 1167–1190 (2005)
F. de Felice, Repulsive phenomena and energy emission in the field of a naked singularity. Astron. Astrophys. 34, 15–19 (1974)
F. de Felice, Classical instability of a naked singularity. Nature 273, 429–431 (1978)
V. Faraoni, Turnaround radius in modified gravity. Phys. Dark Universe 11, 11–15 (2016)
V. Faraoni, M. Lapierre-Léonard, A. Prain, Turnaround radius in an accelerated universe with quasi-local mass. J. Cosmol. Astropart. Phys. 2015(10), 013–013 (2015)
G.W. Gibbons, S.W. Hawking, Cosmological event horizons, thermodynamics, and particle creation. Phys. Rev. D 15, 2738–2751 (1977)
E.G. Gimon, P. Hořava, Over-rotating black holes, Godel holography and the hypertube (2004). arXiv:hep-th/0405019
E.G. Gimon, P. Hořava, Astrophysical violations of the Kerr bound as a possible signature of string theory. Phys. Lett. B 672, 299 (2009)
Z. Gu, H. Cheng, The circular loop equation of a cosmic string in Kerr–de Sitter spacetimes. Gen. Relativ. Gravit. 39(1), 1–7 (2007)
E. Hackmann, B. Hartmann, C. Lämmerzahl, P. Sirimachan, Test particle motion in the space-time of a Kerr black hole pierced by a cosmic string. Phys. Rev. D 82(4), 044024 (2010)
K. Hioki, K.-J. Maeda, Measurement of the Kerr spin parameter by observation of a compact object’s shadow. Phys. Rev. D 80(2), 024042 (2009)
L. Iorio, Constraining the cosmological constant and the DGP gravity with the double pulsar PSR J0737–3039. New Astron. 14(2), 196–199 (2009)
V. Kagramanova, J. Kunz, C. Lammerzahl, Solar system effects in Schwarzschild–de Sitter space–time. Phys. Lett. B 634(5–6), 465–470 (2006)
M. Kološ, Z. Stuchlík, Current-carrying string loops in black-hole spacetimes with a repulsive cosmological constant. Phys. Rev. D 82(12), 125012 (2010)
F. Kottler, Über die physikalischen Grundlagen der Einsteinschen Gravitationstheorie. Annalen der Physik 361(14), 401–462 (1918)
G.V Kraniotis, Precise theory of orbits in general relativity, the cosmological constant and the perihelion precession of Mercury, eds. by H.V. Klapdor-Kleingrothaus, R. Arnowitt (Springer, Berlin, Heidelberg, 2006). pp. 469–479
G.V. Kraniotis, Precise relativistic orbits in Kerr and Kerr–(anti-)de Sitter spacetimes. Class. Quantum Gravity 21, 4743–4769 (2004)
G.V. Kraniotis, Periapsis and gravitomagnetic precessions of stellar orbits in Kerr and Kerr–de Sitter black hole spacetimes. Class. Quantum Gravity 24, 1775–1808 (2007)
G.V. Kraniotis, Precise analytic treatment of Kerr and Kerr–(anti) de Sitter black holes as gravitational lenses. Class. Quantum Gravity 28, 085021 (2011)
G.V. Kraniotis, Gravitational lensing and frame dragging of light in the Kerr–Newman and the Kerr–Newman-(anti) de Sitter black hole spacetimes. Gen. Relativ. Gravit. 46(11), 1818 (2014)
L.M. Krauss, The end of the age problem, and the case for a cosmological constant revisited. Astrophys. J. 501(2), 461–466 (1998)
L.M. Krauss, M.S. Turner, The cosmological constant is back. Gen. Relativ. Gravit. 27(11), 1137–1144 (1995)
H. Kučáková, P. Slaný, Z. Stuchlík, Toroidal configurations of perfect fluid in the Reissner-Nordstr-(anti-)de Sitter spacetimes. J. Cosmol. Astropart. Phys. 2011(01), 033 (2011)
K. Lake, Bending of light and the cosmological constant. Phys. Rev. D 65(8, B), 087301 (2002)
K. Lake, T. Zannias, Global structure of Kerr–de Sitter spacetimes. Phys. Rev. D 92, 084003 (2015)
A.D. Linde, Particle Physics and Inflationary Cosmology (Gordon and Breach, New York, 1990)
A. Müller, B. Aschenbach, Non-monotonic orbital velocity profiles around rapidly rotating Kerr (anti-)de Sitter black holes. Class. Quantum Gravity 24, 2637–2644 (2007)
T. Müller, Falling into a Schwarzschild black hole. Gen. Relativ. Gravit. 40(10), 2185–2199 (2008)
M. Olivares, J. Saavedra, C. Leiva, J.R. Villanueva, Motion of charged particles on the Reissner–Nordström (anti)-de Sitter black hole spacetime. Mod. Phys. Lett. A 26(39), 2923–2950 (2011)
J.P. Ostriker, P.J. Steinhardt, The observational case for a low-density universe with a nonzero cosmological constant. Nature 377(6550), 600–602 (1995)
D. Pérez, G.E. Romero, S.E. Bergliaffa, Accretion discs around black holes in modified strong gravity. Astron. Astrophys. 551, A4 (2013)
C. Pierre-Henri, T. Harko, Bose–Einstein Condensate general relativistic star. Phys. Rev. D 86(6), 064011 (2012)
Planck Collaboration, P.A.R Ade, N. Aghanim, C. Armitage-Caplan, M. Arnaud, M. Ashdown, F. Atrio-Barandela, J. Aumont, C. Baccigalupi, A.J. Banday, et al., Planck 2013 results. XII. Diffuse component separation. Astron. Astrophys. 571, A12 (2014)
D. Pugliese, Z. Stuchlík, Ringed accretion disks: equilibrium configurations. Astrophys. J. Suppl. 221(2), 25 (2015)
L. Rezzolla, O. Zanotti, J.A. Font, Dynamics of thick discs around Schwarzschild–de Sitter black holes. Astron. Astrophy. 412(3), 603–613 (2003)
A.G. Riess et al., Type Ia supernova discoveries at \(z>1\) from the hubble space telescope: evidence for past deceleration and constraints on dark energy evolution. Astrophys. J. 123, 145 (2004)
J. Schee, Z. Stuchlík, M. Petrásek, Influence of the cosmic repulsion on the MOND model of the Magellanic Cloud motion in the field of Milky Way. J. Cosmol. Astropart. Phys. 12, 026 (2013)
T. Schücker, N. Zaimen, Cosmological constant and time delay. Astron. Astrophys. 484(1), 103–106 (2008)
M. Sereno, On the influence of the cosmological constant on gravitational lensing in small systems. Phys. Rev. D 77(4), 043004 (2008)
Z. Sheng, C. Ju-Hua, W. Yong-Jiu, Time-like geodesic structure of a spherically symmetric black hole in the brane-world. Chin. Phys. B 20(10), 100401 (2011)
P. Slaný, Z. Stuchlík, Relativistic thick discs in the Kerr–de Sitter backgrounds. Class. Quantum Gravity 22(17), 3623–3651 (2005)
P. Slaný, Z. Stuchlík, Comment on ‘non-monotonic orbital velocity profiles around rapidly rotating Kerr (anti-)de Sitter black holes’. Class. Quantum Gravity 25(3), 038001 (2008)
D.N. Spergel, R. Bean, O. Dore, M.R. Nolta, C.L. Bennett, J. Dunkley, G. Hinshaw, N. Jarosik, E. Komatsu, L. Page, H.V. Peiris, L. Verde, M. Halpern, R.S. Hill, A. Kogut, M. Limon, S.S. Meyer, N. Odegard, G.S. Tucker, J.L. Weiland, E. Wollack, E.L. Wright, Three-year Wilkinson microwave anisotropy probe (WMAP) observations: implications for cosmology. Astrophys. J. Suppl. 170(2), 377–408 (2007)
D.N. Spergel, R. Bean, O. Dore, M.R. Nolta, C.L. Bennett, J. Dunkley, G. Hinshaw, N. Jarosik, E. Komatsu, L. Page, H.V. Peiris, L. Verde, M. Halpern, R.S. Hill, A. Kogut, M. Limon, S.S. Meyer, N. Odegard, G.S. Tucker, J.L. Weiland, E. Wollack, E.L. Wright, Three year Wilkinson microwave anisotropy probe ( WMAP ) observations: implications for cosmology. Astrophys. J. Suppl. 170(2), 377–408 (2007)
Z. Stuchlík, Equatorial circular orbits and the motion of the shell of dust in the field of a rotating naked singularity. Bull. Astron. Inst. Czechoslov. 31, 129–144 (1980)
Z. Stuchlik, The radial motion of photons in Kerr metric. Bull. Astron. Inst. Czechoslov. 32, 40–52 (1981)
Z. Stuchlík, The motion of test particles in black-hole backgrounds with non-zero cosmological constant. Bull. Astron. Inst. Czechoslov. 34(3), 129–149 (1983)
Z. Stuchlík, An Einstein–Strauss–de Sitter model of the universe. Bull. Astron. Inst. Czechoslov. 35(4), 205–215 (1984)
Z. Stuchlík, Spherically symmetric static configurations of uniform density in spacetimes with a non-zero cosmological constant. Acta Phys. Slovaca 50(2), 219–228 (2000)
Z. Stuchlík, Influence of the relict cosmological constant on accretion discs. Mod. Phys. Lett. A 20(8), 561–575 (2005)
Z. Stuchlík, G. Bao, E. Østgaard, S. Hledík, Kerr–Newman–de Sitter black holes with a restricted repulsive barrier of equatorial photon motion. Phys. Rev. D 58(8), 084003 (1998)
Z. Stuchlík, M. Calvani, Null geodesics in black-hole metrics with nonzero cosmological constant. Gen. Relativ. Gravit. 23(5), 507–519 (1991)
Z. Stuchlík, S. Hledík, Some properties of the Schwarzschild-de Sitter and Schwarzschild-anti-de Sitter spacetimes. Phys. Rev. D 60(4), 044006 (1999)
Z. Stuchlík, S. Hledík, Equatorial photon motion in the Kerr–Newman spacetimes with a non-zero cosmological constant. Class. Quantum Gravity 17(21), 4541–4576 (2000)
Z. Stuchlík, S. Hledík, Properties of the Reissner–Nordström spacetimes with a nonzero cosmological constant. Acta Phys. Slovaca 52(5), 363–407 (2002)
Z. Stuchlík, S. Hledík, J. Novotný, General relativistic polytropes with a repulsive cosmological constant. Phys. Rev. D 94, 103513 (2016)
Z. Stuchlík, S. Hledík, K. Truparová, Evolution of Kerr superspinars due to accretion counterrotating thin discs. Class. Quantum Gravity 28(15), 155017 (2011)
Z. Stuchlík, M. Kološ, Acceleration of string loops in the Schwarzschild–de Sitter geometry. Phys. Rev. D 85(6), 065022 (2012)
Z. Stuchlík, M. Kološ, String loops in the field of braneworld spherically symmetric black holes and naked singularities. J. Cosmol. Astropart. Phys. 2012, 008 (2012)
Z. Stuchlík, J. Kovář, Pseudo-Newtonian gravitational potential for Schwarzschild–de Sitter spacetimes. INTJMD 17(11), 2089–2105 (2008)
Z. Stuchlík, J. Schee, Appearance of Keplerian discs orbiting Kerr superspinars. Class. Quantum Gravity 27(21), 215017 (2010)
Z. Stuchlík, J. Schee, Influence of the cosmological constant on the motion of Magellanic Clouds in the gravitational field of Milky Way. J. Cosmol. Astropart. Phys. 9, 018–018 (2011)
Z. Stuchlík, J. Schee, Comparison of general relativistic and pseudo-Newtonian description of Magellanic-clouds motion in the field of Milky Way. Int. J. Mod. Phys. D 21(4), 1250031 (2012)
Z. Stuchlík, J. Schee, Observational phenomena related to primordial Kerr superspinars. Class. Quantum Gravity 29(6), 065002 (2012)
Z. Stuchlík, J. Schee, Ultra-high-energy collisions in the superspinning Kerr geometry. Class. Quantum Gravity 30(7), 075012 (2013)
Z. Stuchlík, J. Schee, B. Toshmatov, J. Hladík, J. Novotný, Gravitational instability of polytropic spheres containing region of trapped null geodesics: a possible explanation of central supermassive black holes in galactic halos. J. Cosmol. Astropart. Phys. 2017(06), 056 (2017)
Z. Stuchlík, P. Slaný, Equatorial circular orbits in the Kerr–de Sitter spacetimes. Phys. Rev. D 69, 064001 (2004)
Z. Stuchlík, P. Slaný, S. Hledík, Equilibrium configurations of perfect fluid orbiting Schwarzschild–de Sitter black holes. Astron. Astrophys. 363(2), 425–439 (2000)
Z. Stuchlík, P. Slaný, J. Kovář, Pseudo-Newtonian and general relativistic barotropic tori in Schwarzschild–de Sitter spacetimes. Class. Quantum Gravity 26(21), 215013 (2009)
Z. Stuchlík, P. Slaný, G. Török, M.A. Abramowicz, Aschenbach effect: unexpected topology changes in the motion of particles and fluids orbiting rapidly rotating Kerr black holes. Phys. Rev. D 71(2), 024037 (2005)
E. Teo, Spherical photon orbits around a Kerr black hole. Gen. Relativ. Gravit. 35, 1909–1926 (2003)
J.R. Villanueva, J. Saavedra, M. Olivares, N. Cruz, Photons motion in charged anti-de Sitter black holes. Astrophys. Space Sci. 344(2), 437–446 (2012)
L. Wang, R.R. Caldwell, J.P. Ostriker, P.J. Steinhardt, Cosmic concordance and quintessence. Astrophys. J. 530(1), 17–35 (2000)
L. Wang, H. Cheng, The evolution of circular loops of a cosmic string with periodic tension. Phys. Lett. B 713(1), 59–62 (2012)
Acknowledgements
Z.S. acknowledges the Albert Einstein Centre for Gravitation and Astrophysics supported by the Czech Science Foundation Grant no. 14-37086G. D.Ch. acknowledges the Silesian University in Opava Grant no. SGS/14/2016.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Funded by SCOAP3
About this article
Cite this article
Charbulák, D., Stuchlík, Z. Photon motion in Kerr–de Sitter spacetimes. Eur. Phys. J. C 77, 897 (2017). https://doi.org/10.1140/epjc/s10052-017-5401-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-017-5401-9