Spherical photon orbits in the field of Kerr naked singularities

For the Kerr naked singularity (KNS) spacetimes, we study properties of spherical photon orbits (SPOs) confined to constant Boyer-Lindquist radius r. Some new features of the SPOs are found, having no counterparts in the Kerr black hole (KBH) spacetimes, especially stable orbits that could be pure prograde/retrograde, or with turning point in the azimuthal direction. At r>1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r>1$$\end{document} (r<1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r<1$$\end{document}) the covariant photon energy E>0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal{E}> 0$$\end{document} (E<0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal{E}< 0$$\end{document}), at r=1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r=1$$\end{document} there is E=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal{E}= 0$$\end{document}. All unstable orbits must have E>0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal{E}> 0$$\end{document}. It is shown that the polar SPOs can exist only in the spacetimes with dimensionless spin a<1.7996\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a < 1.7996$$\end{document}. Existence of closed SPOs with vanishing total change of the azimuth is demonstrated. Classification of the KNS and KBH spacetimes in dependence on their dimensionless spin a is proposed, considering the properties of the SPOs. For selected types of the KNS spacetimes, typical SPOs are constructed, including the closed paths. It is shown that the stable SPOs intersect the equatorial plane in a region of stable circular orbits of test particles, depending on the spin a. Relevance of this intersection for the Keplerian accretion discs is outlined and observational effects are estimated.


Introduction
Large effort has been devoted to studies of null geodesics in the gravitational field of compact objects, because the null geodesics govern motion of photons that carry information on the physical processes in strong gravity to distant observers, giving thus direct signatures of the extraordinary character of the spacetime around compact objects that influences both the physical processes and the photon motion. The most detailed studies were devoted to the Kerr geometry that is assumed to describe the spacetime of black holes governed by the Einstein gravity [6,7,9,18,19,23,61]. Extension to the a e-mail: daniel.charbulak@fpf.slu.cz b e-mail: zdenek.stuchlik@fpf.slu.cz photon motion in the most general black hole asymptotically flat spacetimes, governed by the Kerr-Newman geometry, has been treated, e.g., in [2,26,35,36,41,62]. The case of spherical photon orbits in the Reissner-Nordström(-de Sitter) spacetimes has been intensively studied in [63]. However, the recent cosmological tests indicate presence of dark energy, probably reflected by a relict cosmological constant, that could play a significant role in astrophysical processes [4,24,43,48,51,56]. Therefore, the photon motion in the field of Kerr-de Sitter black holes, where the cosmological horizon exists along with the static radius giving the limit on the free circular motion [42,46], has been studied in [11,17,29,34,45,47]; extension to spacetimes containing a charge parameter has been discussed in [12,22,44]. There is a large number of studies related to the photon motion in generalizations of the Einstein theory, e.g., for regular black hole spacetimes of the Einstein theory combined with non-linear electrodynamics [38,55], and for black holes in alternative approaches to gravity [1,3,5,35,54,60].
A crucial role in the photon motion is played by the spherical photon orbits, i.e., photons moving along orbits of constant (Boyer-Lindquist) radial coordinate -their motion constants govern the local escape cones of photons related to any family of observers, and specially the shadow (silhouette) of black holes located in front of a radiating source [7], e.g. an orbiting accretion disk [13,20,30,31,45,50]. The properties of the spherical photon orbits outside the outer horizon of Kerr black holes were studied in [58]. Spherical photon orbits in the field of Kerr-de Sitter black holes were discussed in [17,45] -in this case, photons with negative covariant energy could be relevant, in contrast to the case of pure Kerr black holes where only spherical photons with positive covariant energy enter the play [58].
Recently, growing interest in Kerr naked singularity spacetimes is demonstrated, mainly due to the possibility of existence of Kerr superspinars proposed in the framework of String theory by Hořava and his co-workers, with interior governed by the String theory and exterior described by the Kerr naked singularity geometry [25,50,52]. The presence of Kerr superspinars in active galactic nuclei or in microquasars could give clear signatures in the ultra-relativistic collisional processes [53], in the high-frequency quasiperiodic oscillations in Keplerian disks [27,28,59], or in the Lense-Thirring precession effects [16]. The instability of test fields in the Kerr naked singularity backgrounds has been studied in [14,21], possible stabilizing effects were demonstrated for the Kerr superspinars in [33]. The classical instability of Kerr naked singularity (superspinar) spacetimes, converting them to black holes due to standard Keplerian accretion, has been shown to be slow enough in order to enable observation of primordial Kerr superspinars -at least at cosmological redshifts larger then z = 2 [49]. The Kerr naked singularity spacetimes could be applied also in description of the exterior of the superspinning quark stars with spin violating the black hole limit a = 1 [57]. The optical phenomena related to the Kerr naked singularity (superspinar) spacetimes were treated in [37,40,50,52], and in more general case including the influence of the cosmological constant in [17,45,47].
Here we focus our attention to the properties of the spherical photon orbits in the Kerr naked singularity spacetimes, generalizing thus the study of spherical photon orbits in the Kerr black hole spacetimes [58]. For the spherical photon orbits we give the motion constants in dependence on their radius and dimensionless spin, and present detailed discussion of their latitudinal and azimuthal motion. We introduce detailed classification of the Kerr spacetimes according to the properties of the spherical photon orbits, including the stability of the spherical orbits and the role of the spherical photon orbits with negative energy relative to infinity, extending thus an introductory study in more general Kerr-de Sitter spacetimes [17]. We also discuss astrophysically important interplay of the spherical photon orbits and the Keplerian accretion disks, with matter basically governed by the circular geodesic motion; we shortly discuss possible observational effects related to the irradiation of the Keplerian disks by the photons following the spherical orbits.

The Kerr spacetimes
The line element of the Kerr spacetime is in the standard Boyer-Lindquist spheroidal coordinates t, r, θ, φ, with geometric system of units (c = G = 1), described by the well known formula ds 2 = − dt − a sin 2 θ dφ 2 + sin 2 θ a dt − r 2 + a 2 dφ where = r 2 + a 2 − 2Mr, (2) = r 2 + a 2 cos 2 θ. ( Here M is the gravitational mass parameter of the Kerr spacetime, and a is its angular momentum per unit gravitational mass. Without any loss of generality we can assume the parameter a to be positive. If a < M, the Kerr spacetime describes black holes that posses two pseudosingularities (horizons) determined by the condition located at radii If a = M, the horizons coincide at r = M and the spacetime describes an extreme Kerr black hole [7]. In the present paper, we focus ourselves to the case a > M corresponding to the Kerr naked singularity spacetimes where no horizons exist. In all the Kerr spacetimes, the physical singularity is located at r = 0, θ = π/2. The stationary limits of the Kerr spacetimes, determined by the condition g tt = 0, give the boundary of the so called ergosphere where extraction of rotational energy is possible due to the Penrose process [32]. At r < 0, allowed in the Kerr spacetimes, the so called causality-violation region exists [15] -the Kerr superspinars are constructed in such a way that both the causality-violation region and the physical singularity are removed, and it is assumed that they are substituted by an interior regular solution governed by the String theory [25,50]. In order to demonstrate clearly behaviour of the Kerr spacetimes at region close to the physical singularity, it is convenient to use the ("flat") Kerr-Schild coordinates that are connected to the (spheroidal) Boyer-Lindquist coordinate r by the relations x 2 + y 2 = (r 2 + a 2 ) sin 2 θ, z 2 = r 2 cos 2 θ (6) and enable proper visualization of the Kerr spacetime in the innermost regions.

Carter equations of geodesic motion
The symmetries of the Kerr spacetimes imply existence of the time Killing vector ξ (t) = ∂/∂t and the axial Killing vector ξ (φ) = ∂/∂φ. Then the projections E = −ξ (t) · p (energy) and = ξ (φ) · p (angular momentum about the φ-axis) of a particle four-momentum p = dx μ dλ are constants of the particle motion. Here, x μ are the coordinate components of the four-momentum and λ is an affine parameter. Two another motion constants are the rest energy m of the particle (m = 0 for photons) and Q, Carter's constant connected with total angular momentum of the particle. The motion of photons can be described by the Carter equations where the functions P, R, W are defined by the relations It is convenient to use following rescaling λ → λ = λE and put M = 1, i.e., express the radius and time (line element) in units of gravitational mass M. Further, it is usual to introduce the substitution which enables one to replace dealing with the trigonometric functions by the algebraic ones. The Carter equations then take the form Here the dot indicates differentiation with respect to the rescaled affine parameter λ . We have introduced new param- where is the impact parameter. 1 Introduction of the Carter constant Q enables simple classification of the latitudinal 1 In case of more complex spacetimes with non-zero cosmological constant and possibly non-zero electric charge of the gravitating source it is more convenient to introduce the modified impact parameter X = −a, which simplifies the radial equation of motion 15 and its discussion, as shown in [17], where the more general case of the Kerr-de Sitter spacetimes was considered. However, its introduction is unnecessary in this paper. motion in all the Kerr spacetimes, as Q = 0 governs the equatorial motion with θ = π/2, Q > 0 governs the so called orbital motion crossing the equatorial plane 2 , while Q < 0 governs so called vortical motion when crossing of the equatorial plane is forbidden [9,19,42]. It is crucial for the purposes of the present paper that in both Kerr black hole [58] and Kerr (-de Sitter) naked singularity [17] spacetimes the spherical orbits with r = const must be of the orbital type, having thus necessarily Q > 0.

Spherical photon orbits
The SPOs have a crucial role in characterization of the Kerr spacetimes as they govern the shadows of Kerr black holes or Kerr superspinars, and could have influence on the behaviour of the Keplerian disks or more complex accretion structures due to effect of self-illumination [6,41,50].

Covariant energy of photons following spherical orbits
The photons (test particles) whose motion is governed by the Carter equations of motion can be in the classically allowed positive-root states where they have positive energy as measured by local observers, and, equivalently, they evolve to future ( dt/ dλ > 0), or in the classically forbidden negativeroot states with negative energy measured by local observers, evolving into past ( dt/ dλ < 0) [32]. Above the outer horizon of a Kerr black hole, the situation is simple and all photons on the SPOs have positive covariant energy E > 0 and are in the positive-root states. However, the situation is more complex in the case of Kerr naked singularities. For distinguishing of the positive-root and negative-root states in the case of Kerr naked singularities, and under the inner horizon of Kerr black holes, we can use the time equation of motion (18); alternatively, the projection of the photon 4-momentum on the time-like tetrad vector of physical observers can be used for this purpose [10,17,45]. The time equation can be written in the form Now it is clear that in the regions of spherical orbits where the bracket on the r.h.s. is positive, we have positive-root states for E > 0, while at regions where the bracket takes negative values, the positive-root states must have the energy relative to infinity E < 0. Notice that in the case of E < 0 the photons with < 0 have positively-valued impact parameter = E . We present the results of the determination of the covariant energy of the SPOs later, and then we use them in the classification of the Kerr spacetimes where this property is considered as one of the criteria of the classification.

Motion constants of spherical photon orbits
The simultaneous solution of the equations where R(r ) denotes the r.h.s. of the equation (15), yields the motion constants of the spherical orbits 3 The functions (22), (23) determine the constants of motion of the SPOs in dependence on their radius, hence, they represent the basic characteristics of the spherical orbits and require careful analysis.

Existence of spherical orbits
First, we shall devote attention to the function in Eq. (22) giving limits on the parameter q. Discussion of this function was already performed in [58] for the case of the Kerr black holes, here we extend the discussion to the case of the Kerr naked singularities, including also the properties of the SPOs under the inner horizon of Kerr black holes that were not discussed in [58]. The function q sph (r ; a) is well defined for any radius r = 1, while it diverges at r = 1 for a = 1, with lim r →1 q sph = ∓∞ as a ≶ 1.
In the case of extreme Kerr black holes, one can find by substituting a = 1 into the conditions (21) that the function in (22) should be replaced by which is now continuous even for r = 1, c. f. [58]. As we shall see bellow, only non-negative values of the parameter q permit the motion of constant radius. 4 The zeros of (22) determining the photon equatorial circular orbits are given by the condition For a = 0 this function gives the SPO around the Schwarzschild black hole, located at r = 3; if 0 < a < 1, 3 There exists a second family of solution of these equations, but it is not physically relevant -see [17,58]. 4 The same applies to a more general case of the Kerr-de Sitter spacetimes, see [17].
it determines one circular orbit under the inner black hole horizon with radius and two circular orbits located above the outer black hole horizon, which are the main subject of astrophysical interest. The inner one being co-rotating, the outer being counterrotating, both being unstable with respect to radial perturbations. Their radii r ph+ , r ph− , where r ph+ < r h− , can be expressed by the relation [7] r ph± = 2 1 + cos 2 3 arccos(∓a) .
Above the outer horizon of the KBH spacetimes, the SPOs are thus located at radii In the extreme KBH case a = 1, the counter-rotating orbit is located at r = 4, while the co-rotating orbit shares the same (Boyer-Lindquist) radius r = 1 with the black hole horizon, although they are in fact separated by a non-zero proper radial distance [7]. In the Kerr naked singularity spacetimes, only the counterrotating equatorial photon orbit exists. For a > 1, we express the photon orbit radii r ph− by the expression [40] The co-rotating circular photon orbit should be (formally) located at r = 0, θ = π/2, representing thus the limit of unstable circular equatorial orbits of test particles; however, such an orbit has no physical meaning as it coincides with the physical singularity of the Kerr spacetimes. Therefore, in the KNS spacetimes, the SPOs can be located at all radii

Stability of spherical photon orbits
Stability of the spherical null geodesics against radial perturbations is governed by the condition d 2 R/dr 2 < 0 considered in the loci of the spherical orbits, i.e., by the extrema of the function q sph (r ; a). The function q sph (r ; a) has one local extreme q ex = 27 located at r = 3, independently of the rotational parameter a. This extreme is a local maximum for 0 < a < 3, while for a > 3 it becomes a minimum. The significance of this extreme, as follows from (17), is that the photon orbit at r = 3 intersects the equatorial plane perpendicularly, i.e.,φ(r = 3, θ = π/2) = 0 (c.f. [58]). For the stability criterion, another extreme of the function q sph (r ; a) is relevant that is determined by the condition a = a stab (r ) ≡ (r − 1) 3 + 1, The function a stab (r ) governs the stability of the SPOs -the orbits above/bellow the curve a stab (r ) are stable/unstable which is implied also directly by the condition d 2 R/dr 2 = 0 determining the marginal stability of the spherical orbits. The stability condition d 2 R/dr 2 < 0 can be thus written in the form a > a stab (r ). The function a stab (r ) governing the marginally stable SPOs is increasing with increasing r , having an inflexion point at the special radius r = 1; its behaviour is demonstrated in Fig. 1. It is immediately clear that the marginally stable spherical orbits are located at r > 1 in KNS spacetimes, while in the KBH spacetimes it must be located under the inner horizon r < r − < 1. The radius of the marginally stable SPO, r ms± , can be directly expressed by the simple formulas or r ms+ = 1 + (a 2 − 1) 1/3 for a 2 ≥ 1.
The stable (unstable) spherical orbits are located for given parameter a at r < r ms± (r > r ms± ).

Polar spherical photon orbits
Now we consider behaviour of the function sph (r ; a). In the case of non-extreme Kerr BH spacetimes, this function is monotonically decreasing in the stationary region above the outer horizon and its point of discontinuity at r = 1 is hidden between the black hole horizons. In the extreme KBH case a = 1, the equation (21) reduces to the form having no discontinuity For the KNS spacetimes, a > 1, there is lim r →1 ∓ sph = ±∞ and the discontinuity occurs at r = 1.
In order to find the special case of spherical orbits covering whole the range of the latitudinal coordinate, reaching thus the symmetry axis at θ = 0, we have to find when the function sph (r ; a) takes the significant value of sph (r ; a) = 0 corresponding to photons with zero angular momentum, since only such photons can reach the symmetry axis. We adhere notation introduced in [17] and denote the radii of the polar SPOs crossing the symmetry axis by r pol . These radii can be found by solving the equation The function a pol(r ) has zeros at r = 0 and at r = 3, and a local maximum at r pol = √ 3 for The solution of Eq. (35) can be written in whole the relevant region of the rotation spacetime parameter a ∈ (0, a pol(max) ) in the form where r pol+ < r pol− . In the case of KBH spacetimes, 0 < a < 1, only the formula for r pol− is relevant and it gives the only polar SPO in the stationary region, where r pol− < 3 [58]. In the KNS spacetimes with 1 < a < a pol(max) = 1.17996, two polar SPOs exist at the radii r pol± given by Eq. (37). For a = a pol(max) these radii coalesce at r pol = √ 3. For a > a pol(max) no polar SPOs exist. Combining Eq. (22) and Eq. (35) one can find the values of parameter q of the polar orbits to be given by The graph of the function (38) is depicted in Fig. 3. Its minimum value is q pol(min) = 19.3923 for r = r pol and in the limit points of its definition range it is We shall still mention the polar SPOs in section devoted to latitudinal motion. The local extrema of the function sph (r ; a) are determined by the relation (31), i.e., they are located at the same radii as the local extrema of the function q sph (r ; a). The functions (25), (31), (35) are, together with the function

Spherical photon orbits with negative energy
For further discussion it is now necessary to determine the loci of the SPOs with negative energy E < 0 (Fig. 3). Substituting = sph into the time equation (20), we arrive to equation where Clearly, the equation (40) describes the motion of photons with positive energy E > 0 directed to the future ( dt/ dλ > 0) on spherical orbits at radii r > 1. 5 On the contrary, for 0 < r < 1 the condition dt/ dλ > 0 demands E < 0 and m < m zE (r ; a), or E > 0 and m > m zE (r ; a), where In the region of interest, 0 < r < 1, the function m zE (r ; a) is increasing, and m zE (0; a) = 0 and m zE (r ; a) → ∞ as r → 1. As we shall see bellow, at the interval 0 < r < 1, the inequality m > m zE (r ; a) is inconsistent with the reality condition for the latitudinal motion, hence the latter alternative is irrelevant and the range 0 < r < 1 corresponds to region of the SPOs with E < 0. Of course, this region must be hidden under the inner black hole horizon in case of the Kerr black holes. From discussion of the stability of the SPOs it follows that all the SPOs with E < 0 must be stable against 5 Here we restrict our discussion on the stationary region > 0. radial perturbations. Of course, photons from the stable SPOs cannot escape to infinity due to a perturbation.

Spherical photon orbits with E = 0
As we have shown above, the radii 0 < r < 1 correspond to range of SPOs with negative energy, which suggests that there are SPOs with zero energy at radii r = 1. It is supported by the fact that the functions (22,23,41) diverge at this radii. For this reason, let us consider the Carter equations of motion for the photons (7-9) with explicitly expressed energy E and revise the discussion assuming E = 0. The conditions (21) then imply r = 1 as expected, and the ratio of the motion constants of the SPO with E = 0 in the KNS spacetimes reads Notice that the impact parameters and q are not defined in this special case. In the extreme KBH spacetimes, = 0 is required for arbitrary Q > 0 for the spherical orbits at r = 1.

Trajectories of photons on the spherical null geodesics
In order to construct trajectories of the photons following the spherical null geodesics, we have to discuss in detail the latitudinal and azimuthal motion at the r = const surfaces.
In the context of the spherical motion of photons the natural question arises, what is the range of the latitudinal coordinate in dependence on the allowed motion constants and the dimensionless spin of the Kerr spacetime. Simultaneously, the important question is on the possible existence and number of the turning points of the azimuthal motion.

Latitudinal motion
The latitudinal motion can be of the so called orbital type, where the photons oscillate between two latitudes θ 0 , π −θ 0 , crossing repeatedly the equatorial plane or even being con- Fig. 3 Dependence of the parameter q pol of the polar spherical orbits on its radii r pol depicted with corresponding dimensionless spin a pol (left) and its projection q pol (r pol ) onto the (r − q) plane (right) fined to the equatorial plane 6 , or of the so called vortical type, where the photons oscillate 'above' or 'bellow' the equatorial plane between two pairs of cones coaxial with the symmetry axis of the spacetime, with latitudes θ 1 , θ 2 , θ 1 < θ 2 and π − θ 1 , π − θ 2 . The special case is the vortical motion along the symmetry axis θ = 0, or the motion at any constant latitude -such photons are called PNC photons and have a generic role in the Kerr spacetimes [9]. Now one can ask, which of these types is possible for the spherical photon motion. First, we demonstrate that for the spherical orbits there is q ≥ 0 necessarily, i.e., the motion is of the orbital type. This can be shown easily using the Carter equation of the latitudinal motion 16), from which it can be immediately seen that for q < 0 the condition M(m) ≥ 0 ensuring the existence of the latitudinal motion is fulfilled only if On the other hand, we can show that a 2 − q sph (r ; a) + l 2 sph (r ; a) < 0 so that the possibility q < 0 must be rejected [17,58]. In fact, there is even more restrictive condition for the motion constants than that given by (44). To show this, let us express the reality condition M(m) ≥ 0 using linearity in parameter q by The equality gives the turning points in the latitudinal coordinate. Notice that using the coordinate m = cos 2 θ , we have to restrict the range of the solutions to m ∈ 0; 1 . Of course, there is obvious zero of (16) given by m = 0, emerging due to the used substitution, which indicates just a transit through the equatorial plane. Behaviour of the function q m (m; a, ) in the limit points of the interval 0; 1 is as follows: Another zero of q m (m; a, ) is at The local extrema can be expressed from the condition dq/ dm = 0 in an implicit form or, equivalently Clearly, the extrema exist for | | ≤ a. The subscript 'min' in (47) indicates that these extrema must be minima, as follows from the inequality The values of these minima read q min = −(a −l) 2 for 0 ≤ ≤ a, or q min = −(a + l) 2 for −a ≤ ≤ 0. The above results can be written in a compact form This is the relation giving stronger limitation on the impact parameter in case q < 0 than that given by (44).
From the behaviour of the function q m (m; a, ), which is shown in Fig. 4, it can be seen that the vortical motion exists for negative values of the parameter q. According to the relation (48), the negative values of the parameter q allowing the latitudinal motion exist for | | ≤ a, where the lowest value is q = −a 2 and occurs for = 0. However, as can be verified by calculation, at radii where |l sph | ≤ a, there is q sph > 0, which confirms that the vortical motion of constant radius is impossible for the spherical orbits.
The latitudinal turning points of the SPOs are given by zeros of the r. h. s. of Eq. (16) with substitution of = sph . The maximum latitude reached by a photon following a spherical null geodesic at a particular radius can be inferred from the relation where the latitudinal turning function (relevant in the range m ∈ (0, 1)) is defined by The function m θ (r ; a) is real everywhere in the stationary region. Its zeros are determined by the function a co (r ). For KBHs (0 < a < 1), the function m θ (r ; a) has two local maxima (Fig. 5a); one under the inner black hole horizon, located at radius determined by the function a stab (r ), i.e., at r = r ms− . The second maximum of m θ(max) = 1 is determined by the function a pol (r ). In extreme KBHs case where a = 1 (Fig. 5b), the function m θ (r ; a) has one local maximum m θ(max) = 1 given by a pol (r ). For fixed rotation parameter a, the radii r pol± are given implicitly by relation (35) and explicitly by Eq. (37). The corresponding impact parameter q pol = q sph (r pol (a), a) is represented in Fig. 3; the other motion constant, l pol = 0 by definition. For the KNS spacetimes with a < a pol(max) (Fig. 5c), two maxima exist at m θ(max) = 1 that are given by the function a pol (r ), and one local minimum located at r = r ms+ . If a = a pol(max) , the three extrema coalesce into maximum m θ(max) = 1 (Fig. 5d). In these KNS spacetimes thus polar SPOs can exist for properly chosen motion constants at the properly chosen radii.
For the KNS spacetimes with a > a pol(max) (Fig. 5e, f), the function m θ (r ; a) has a local maximum at m θ(max) < 1 at r = r ms+ . In such KNS spacetimes the polar SPOs cannot exist.

Azimuthal motion
Finally, let us consider the relations of the function m θ (r ; a) with a function m φ (r ; a) determining the latitude at which a turning point of the azimuthal motion occurs. The equation (17), after performing substitution = sph , can be written in the form where The condition dφ/ dλ ≥ 0 implies where the azimuthal turning function is defined as The change of sign in denominator of (51), while crossing the divergence point r = 1, now plays no role, since, as follows from the preceding discussion, in order to have dt/ dλ > 0, the relation E/(1 − r ) < 0 holds at any radius r since E < 0(E > 0) at r < 1(r > 1). Therefore, the motion in the φ-direction is fully governed by the function (52) and the inequalities presented above. If we compare the expressions in (35) and (52), we see that The turning point of the azimuthal motion thus exist for all KBHs and KNSs. The functions m θ (r ; a)and m φ (r ; a) have common points at r = 1, where they have value m θ (1; a) = m φ (1; a) = 1/a 2 , and at the local extrema determined by the curve a pol (r ). As is evident from Eq. (53), the function m φ (r ; a) has a local maximum at r = √ 3 = r pol with m φ(max) = (6 √ 3 − 9)/a 2 = a 2 pol(max) /a 2 . The position of the SPOs with turning point of the azimuthal motion is represented in Fig. 14. The latitudinal angle where the azimuthal turning point occurs is given by the relation The existence of the azimuthal turning points enables existence of "oscillatory" orbits with change of the azimuthal angle for half period in the latitudinal motion φ(a) = 0.

Shift of nodes
Now we examine the dragging of the nodes, i.e., we thus determine the azimuthal angle between the two points, where  (r ; a). Therefore, all the SPOs with negative energy are prograde the ascending, i.e., the 'northwards' directed parts of the photon track intersects the equatorial plane. For this purpose we need to evaluate the change in the azimuthal coordinate φ by integrating dφ/dm expressed from the Carter's equations 16, 17 with = sph (r ; a), q = q sph (r ; a) inserted. 7 The exact expression for the nodal shift φ was found in [58]with our notification it reads 7 The more general case of the function M(m) with non-zero cosmological parameter was studied in detail in [17].
Here m ± are the positive and negative roots of the function M(m), while Fig. 6 Dependence of the node shift φ on the radius of the spherical orbit. A discontinuity occurs at r = r pol± the jump change is 2π between the single point and each of the end points of the curve, c. f. [58]. For a = a pol(max) = 1.18 the points of the discontinuity coalesce at r = r pol = √ 3. For a > a pol(max) the function is continuous. The zero of φ first occurs for parameter a = a φ=0(min) = 1.179857 slightly lower than a pol(max) ; it indicates an existence of an oscillating orbit, which forms a boundary between prograde ( φ > 0) and retrograde orbits ( φ < 0.) For a < a pol(max) the boundary is at r = r pol− .

and
(v, are the complete elliptic integrals of the first and third kind, respectively. The dependence of the nodal shift φ on the radial coordinate qualitatively differs in dependence on the presence of the polar orbits. If they are present, discontinuities occur at the radii r pol± -see Fig. 6a-d. The jump drop is always 4π , and the single points, which are inherent values in the discontinuity points r pol± , are halfway between the endpoints of the interrupted curve. This is a simple consequence of the singular behaviour of the Boyer-Lindquist coordinates at the poles, as explained in [58]. For a pol(max) these discontinuities coalesce at r = r pol , for which the nodal shift function φ < 0 (Fig. 6e).

Periodic orbits
The nodal shift function becomes to be continuous for the KNS spacetimes with a > a pol (max) . Therefore, in KNS spacetimes with a > a pol(max) ∼ 1.18, there exist "oscillatory" trajectories with φ = 0, i. e., the photons in such trajectories are following a closed path with finite extent in azimuth, ending at the starting point. They are of octal-like shape (see, e. g., Fig. 9 case k = 0). However, a detailed cal- Radii r φ=0 of orbits with zero nodal shift φ in dependence on the spin parameter a and comparison with the loci of marginally stable orbits r ms . For given dimensionless spin a the radii of the stable orbits are r < r ms (a). Hence all such orbits are stable. Note that the detail is approx. 10 4 -times horizontally stretched culation using numerical procedure reveals that the radius r φ=0 of the zero nodal shift φ = 0 first occurs for a = a φ=0(min) ≡ 1.179857 at r = r φ=0(max) = 1.71473 (see Fig. 6d). In Fig. 7 we show that the function r φ=0 (a) is descending, hence the subscript ' φ = 0(max)'. Notice that the value of a φ=0(min) is only very slightly lower than a pol(max) . Since the point r φ=0(max) and corresponding point of the 4π -discontinuity r pol+ (a φ=0(min) ) are infinitesimally separated, we can claim that the SPO at r pol+ (a φ=0(min) ) is a special case of polar oscillatory orbit. We illustrate its trajectory explicitly in Fig. 10 and for comparison we give illustration of the orbit at r = r pol for a = a pol(max) .  Fig. 6 in independence on the dimensionless spin a.
In the Fig. 8, the nodal shift is shown for the polar spherical orbits r pol± , i. e., for the isolated points in Fig. 6, in dependence on the spin parameter a. As follows from the above, varying a, the first occurrence of the zero nodal shift appears when there is φ[r pol+ (a)] = 2π . According to Figs. 6, 7, 8, we can summarize that for the very tiny interval a φ=0(min) ≤ a ≤ a pol(max) there exist orbits with φ < 0, i. e., globally retrograde, with radii r φ=0(max) < r < r pol+ (a), followed by orbits with φ > 0, i. e., globally prograde, at r pol+ (a) < r < r pol− (a), and again orbits globally retrograde with r > r pol− . Otherwise, the SPOs of this class of the KNS spacetimes possesses no new essential features in comparison with the other cases, therefore, its character is explained sufficiently (Figs. 1, 10).
The radius of the oscillatory orbits r φ=0 is the only point dividing the trajectories which are globally prograde ( φ > 0), and those that are globally retrograde ( φ < 0) in case a > a pol (max) . For the KNS spacetimes with a φ=0(min) < a < a pol(max) the role of such dividing points play the radii r φ=0 and the discontinuity points r pol± . For a < a φ=0(min) only the radius r pol− is the divider, however, the oscillatory orbits are not possible (c.f. orbits in the BH background with a = 0.9 at r r pol− , r = r pol− and r r pol− in Fig. 11). The function r φ=0 (a) is compared in Fig. 21 with another relevant functions.
Of course, another closed paths corresponding to the general periodic orbits occur whenever In Fig. 9, we present basic types of the closed SPOs for several ratios k = m/n, giving the cases of prograde (k > 0), oscillate (k = 0), and retrograde (k < 0) orbits. Of course, of special astrophysical relevance could be the oscillate orbits with φ = 0, as they immediately return to the starting point at the azimuthal coordinate related to distant static observers. Further, we consider a possibility of orbits having a plurality of revolutions about the spacetime symmetry axis per one latitudinal oscillation (|k| > 1), which are of a helix-like shape, and that of greater number of latitudinal oscillations per one revolution about the axis (|k| < 1), both with, or without a turning point in the φ-direction. Note that the cases k ≤ −1 cannot be realized, since from Eq. (55) it follows that the nodal shift φ > −2π . This is demonstrated in Fig.  12 on the left, which depicts the shift of nodes at r → r ph− .

Spherical photon orbits with zero energy and their nodal shift
For the SPOs with zero energy, located at radii r = 1, the constants of motion have to fulfil the relation 2 /Q = a 2 −1.
The latitudinal equation (8) with the substitution m = cos 2 θ then reads and the azimuthal equation (9) reads Therefore, the turning points of the latitudinal and the azimuthal motion coalesce at m = 1/a 2 , in concordance with the result for intersection of the functions m θ (r ; a), m φ (r ; a), as we have found earlier. In this case, the point with dφ/dλ = 0 is not a real azimuthal turning point, it is only the point of vanishing of the azimuthal velocity, as shown in Fig. 13. The complete change of the azimuthal coordinate per one latitudinal oscillation, denoted φ, is now given by the relation which corresponds to values of the local maxima of the function (55). Behaviour of this function is shown in Fig. 12.

Summary of properties of spherical photon orbits
We can summarize properties of the SPOs in the following way. The covariant energy E > 0 (E < 0) have the SPOs at r > 1 (r < 1) for all KBHs and KNSs; there is SPO with E = 0 at r = 1. The other properties depend on the dimensionless spin a.
In the KBH spacetimes, the stable SPOs are located under the inner event horizon at 0 < r < r ms− , the unstable ones at r ms− < r < r ph0 , all being co-rotating (prograde), none of these orbits can be polar. Above the outer event horizon only Details, including the minimum attained latitude θ min , and latitude θ φ at which the turning point in the φ-direction occurs, sign of the photon's energy E and the time period of the orbit (to be discussed bellow), characterizing the orbits are presented in the attached Table 1 unstable SPOs spread in the interval r ph+ ≤ r ≤ r ph− . They are purely co-rotating for r ph+ < r < r pol− ; with turning point in the azimuthal direction but globally counter-rotating (retrograde) for r pol− < r < 3; and purely counter-rotating at 3 < r < r ph− . One polar orbit exist above the outer horizon of the KBH spacetimes.
In the KNS spacetimes the SPOs exist at 0 < r < r ph− , being stable/unstable for r ≶ r ms+ . For the KNS spacetimes with 1 < a < a φ=0(min) the SPOs at r < 1 are purely co-rotating; for 1 < r < r pol+ they have turning point in the azimuthal direction but are co-rotating globally; between the inner stable and the outer unstable polar orbit, i. e., at r pol+ < r < r pol− , they are purely co-rotating; at r pol− < r < 3 they have turning point in the azimuthal direction and are counter-rotating globally; for 3 < r < r ph− they are purely counter-rotating. Two polar orbits can exist in such KNS spacetimes.
There exists extremely small interval of the spin parameter a φ=0(min) < a < a pol(max) , for which the KNS spacetimes possess two polar orbits and one orbit of zero nodal shift r φ=0 (a) -for 1 < r < r φ=0 (a) there are orbits with turning point in the azimuthal motion, being co-rotating in global; such orbits occur also for r φ=0 < r < r pol+ , being globally counter-rotating; at the radii r pol+ < r , the same behaviour occurs as in the previous case.
For the KNS spacetimes with a > a pol (max) , the SPOs at r < 1 are purely corotating; at 1 < r < 3 there are orbits with turning points of the azimuthal motion, which can be both co-rotating or counter-rotating globally, in dependence on the spin a; spherical orbits at r > 3 are purely counter-rotating. For the KNS spacetimes with a > 3 there is 3 < r ms+ , hence stable purely counter-rotating spherical orbits exist at 3 < r < r ms+ . All the above results are clearly illustrated in Fig. 14.   (61) in dependence on the rotational parameter a, depicted as the full curve, which represents the shift of nodes φ of orbits with zero energy. The dot-dashed curve is a depiction of the lowest negative shift of nodes of a retrogressive spherical orbits approaching the position of the equatorial circular counter-rotating orbit at r = r ph− . It can be seen that the increase of the rotational parameter a weakens the negative shift due to strong dragging of the spacetime. The shift is always φ > −2π

Classification of the Kerr spacetimes due to properties of the spherical photon orbits
Using the knowledge of the behaviour of the above described functions, we are now able to summarize the properties of the SPOs in dependence on the value of the Kerr spacetime spin parameter a, giving thus the corresponding classification of these spacetimes. The following description of the individual classes is supplemented by an explicit illustrations of the spatial structures of the SPOs, which can regarded as a spatial representation of the Fig. 14. In the figures we use the so called Kerr-Schild coordinates x, y, z that are connected to the Boyer-Lindquist coordinates r, θ by the relations We present a meridional sections of the SPOs with the ycoordinate being suppressed, hence, the surfaces of constant Boyer-Lindquist radius are depicted as an oblate ellipses.
Class I: Kerr black holes 0 < a < 1 (see Fig. 15) endowed with the SPOs of two families. First family is limited by radii 0 < r < r ph0 < r − , where is the radius of co-rotating equatorial circular photon orbit located under the inner black hole horizon r − (green dot). These are the orbits with negative energy E < 0; at 0 < r < r ms− they are stable (rich green area), for r ms− < r < r − they are unstable with respect to radial perturbations (light green area). The radius r ms− , given by (33), denotes marginally stable spherical orbit with negative energy. All the first family orbits are prograde and span small extent in latitude in vicinity of the equatorial plane -its maximum at r = r ms− is approaching the value θ min(zE) (a = 1) = arccos 2 √ 3 − 3 = 47.1 • as r → 1 when a → 1. Here and in the following, we restrict our discussion on the 'northern' hemisphere (0 ≤ θ ≤ π/2), the situation in the 'southern' hemisphere is symmetric with respect to the equatorial plane. Second family of the SPOs spreads between the inner corotating (red dot) and outer counter-rotating (blue dot) equatorial circular orbits with radii given by (27). There is no turning point of the azimuthal motion for orbits with radii r ph+ < r < r pol− , where r pol− is the radius of the polar spherical orbit (purple ellipse) given by (37), and all such photons are prograde. At radii r pol− < r < 3 (black dotted incomplete ellipse), there exist orbits with one turning point of the azimuthal motion in each hemisphere, such that the photons become retrograde as they approach the symmetry axis. At radii 3 < r < r ph− , all spherical orbits are occupied by retrograde photons with no turning point of the azimuthal motion.
Class II: Extreme KBH with a = 1 (Fig. 16). A family of stable prograde spherical orbits with negative energy occurs at radii 0 < r < 1 near equatorial plane, 8 where also non-spherical bound photon orbits with two turning points of the radial motion exist. Such orbits are not present in any kind of the most general case of the Kerr-Newmann-(anti) de Sitter black hole spacetimes. At radii 1 < r < r ph− = 4, there are orbits with positive energy and properties similar to the second family orbits in the previous KBH case. The special case of r = 1 apparently corresponds to the SPO with zero energy, which is marginally stable and prograde. This zero energy orbit has minimum latitude θ min(zE) (a = 1) = 47.1 • ; note that for a > 1, there is θ min(zE) (a) = arccos (1/a). The Fig. 13 The SPOs with zero energy, located solely at radius r = 1, are presented for KNS spacetimes with increasing value of the rotational parameter a. The turning points of the latitudinal and the azimuthal motion coalesce in these cases, hence the turning loops have shrinked into peaks. In the attached Table 2 we demonstrate the descending total change in azimuth per one latitudinal oscillation, and the descending range in the latitudinal motion, with increasing spin parameter a. In some cases closed spherical orbits emerge due to properly chosen parameters of the SPOs a = 1.0001 a = 1.001

Fig. 14 Characteristic functions a(r ) with supplementary informations about stability and latitudinal and azimuthal motion for wider
interval of the parameter a (above) and in greater detail (bellow). Significance of the characteristic functions is given in Fig. 1. Gray shading demarcates the dynamic region, colouring highlights region of the SPOs. Orbits with E < 0, all being purely prograde, i. e., without turning points in φ-direction, are highlighted in green; purely prograde orbits with E > 0 in red; orbits with turning points in φ-direction and with E > 0 in purple and purely retrograde with E > 0 in blue. Full/faint hues correspond to stable/unstable orbits. Thin black curve is the zero nodal shift function a zs (r ) dividing regions of globally prograde/retrograde orbits marked by +/− sign SPOs at radii 0 < r ≤ 1 have the same properties as those that occur in all KNS spacetimes, and we shall not repeat them in the following cases. Similarly, for all KNS spacetimes the orbits at r > 1 have positive energy.
Class III: KNS spacetimes with 1 < a < a φ=0(min) = 1.17986 (see Fig. 17). Two polar SPOs appear at radii r pol+ , r pol− (inner and outer purple ellipse, respectively) given by Eq. (37). Photons at 1 < r < r pol+ are stable, they have E > 0 and one turning point of the azimuthal motion in each hemisphere. The latitudinal coordinate where the SPO has the azimuthal turning point is given by Eq. (54), the minimum allowed latitude reads θ min = arccos √ m θ (r ; a). At the region θ min ≤ θ < θ φ , the photon motion is in negative φdirection (deep blue area), while for θ φ ≤ θ ≤ π/2 it is in positive φ-direction (deep red area). The break point dividing the globally prograde orbits from the globally retrograde ones is at r = r pol+ . The motion constants sph → −∞, q sph → ∞, as r → 1 from the right, and sph = 0, q sph = 27 at r = r pol+ . The orbits at the radii r pol+ < r < r pol− are prograde, with no change in the azimuthal direction, for r pol+ < r < r ms+ being stable (deep red area), for r ms+ ≤ r ≤ r pol− being unstable (light red area). The function m θ (r ; a) has a local minimum at r = r ms+ , hence, the marginally stable SPO is of the least extent in the latitude (black dashed curve), contrary to the case of the marginally stable SPOs with negative energy at r ms− at the KBH spacetimes, where they have the widest extent (see detail in Fig. 15). The motion constant q sph of photons on this orbit corresponds to the local minimum 0 < q sph(min) < 27 of the function defined in (22), and the local maximum 0 < sph(max) of the function defined by (23). The turning point of the azimuthal motion appears for the SPOs at the radii r pol− < r < 3, and such SPOs appear to be retrograde as whole. The motion constants for r = r pol− are q sph(min) < q sph < 27, and sph = 0. At r = 3, there is q sph = 27, corresponding to the local maximum of (22), i.e., to the photons crossing the equatorial plane with zero velocity component in the φ-direction ([58]), and sph < 0. The SPOs in the region 3 < r < r ph− are purely retrograde with sph < 0 and q sph → 0 as r → r ph− .
Class IV: KNS spacetimes with a φ=0(min) ≤ a < a pol(max) = 1.17996. In the limit case a = a φ=0(min) , there exist radius r = r φ=0(max) = 1.7147 of spherical orbit with zero nodal shift φ = 0, infinitesimally distant from the 4π -discontinuity point r pol+ (see Fig. 6d). The radius r pol+ therefore corresponds to special case of polar oscillatory orbit. For a φ=0(min) < a < a pol(max) , the radius which separates the globally prograde SPOs from the globally retrograde ones is at r = r φ=0 r φ=0(max) < r pol+ (see the detail of Fig. 14). For  Fig. 14 -in green it is depicted the region of the spherical orbits with negative energy E, in red the region of photon prograde motion in the azimuthal direction, and in blue the region of the retrograde azimuthal motion. Note that not all the orbits are purely prograde, nor purely retrograde, but can have a turning point of the azimuthal motion, hence the colouring of a particular curve representing some fixed Boyer-Lindquist radius r can be change. Full/light hues correspond to stable/unstable orbits. Grey shading demarcates the area where g tt ≥ 0 (ergosphere), grey curve is its boundary (ergosurface). Purple curve is the polar orbit. Loci of turning points of the azimuthal motion are designated by red curves, loci of turning points of the latitudinal motion are distinguished by blue curves. In all figures some outstanding radii are highlighted, namely r = 1 (full black ellipse), r = 3 (black dotted incomplete ellipse). The ring singularity is depicted by the horizontal abscissa, the spin axis by the vertical dashed line. The bold dots represent the photon equatorial circular orbits, namely, the co-rotating orbits with negative energy E < 0 (green), the co-rotating orbits with E > 0 at r = r ph+ (red), and the counter-rotating orbits at r = r ph− (blue). In addition, we included images of a possible Keplerian accretion discs, represented by the first family of the equatorial circular orbits of the test particles [40] (bold black horizontal abscissas). Their inner edge is the marginally stable circular orbit r φ=0 < r < r pol+ there appear globally retrograde orbits. The other properties remain the same as in previous case and the SPOs structure is represented by the Fig. 17.
For KNS spacetimes with a > a pol (max) , there are no spherical polar orbits. As a consequence, no purely prograde spherical orbits, neither stable nor unstable, are possible. The radius which separates the globally prograde spherical orbits from the globally retrograde ones is at r < r φ=0(max) (black dot-dashed curve in Fig. 18 below) and it slowly decreases as a → ∞ (see Fig. 7). In addition to this, the discussion is qualitatively same as in the previous case.
Class VI: Kerr naked singularity spacetimes with a ≥ 3. For a = 3 the local extrema of the function q sph (r ; a) coalesce at the inflex point at r = 3 = r ms+ with q in f = Fig. 18 Structure of the SPOs in the KNS spacetime with a = a pol(max) (above) and with a = 1.7 (below) corresponding to the Class V 27, where it becomes local minimum q sph,min = 27 for a > 3. The local maximum is then at r ms+ > 3 (see Fig. 2), which is also locus of the local maximum of latitudinal turning function m θ (r ; a) ( Fig. 5f) at the radius of the marginally stable spherical orbit. The orbits at 3 < r < r ms+ are stable and purely retrograde; in the range r ms+ < r < r ph− , there are unstable retrograde orbits (see Fig. 19).
We present a systematic construction of the spherical photon trajectories for the outstanding radii and other appropri- ately chosen representative radii for the KNS spacetime of the Class III with dimensionless spin parameter a = 1.1 in Fig. 20. Basic characteristics of these orbits, i. e., their radii r , impact parameters , q, total nodal shift φ , minimum attained latitude θ min , latitude of change in the φ-direction θ φ , sign of energy E and type, are presented in Table 3.

Spherical photon orbits related to the Keplerian disks and possible observational consequences
It is well known that in the KBH spacetimes the SPOs define the light escape cones in the position of an emission, i.e., they represent a boundary between the photons captured by the black hole and the photons escaping to infinity. In case of the KNS spacetimes, there is in addition a possibility of existence of trapped photons, which remain imprisoned in the vicinity of the ring singularity [39,50]. Such a "trapping" region spreads in the neighbourhood of the stable SPOs, where small radial perturbations cause that the photons oscillate in radial direction between some pericentre and apocentre, but remain trapped in the gravitational field -efficiency of the trapping process was for the Kerr superspinars (naked singularity spacetimes) studied in detail [50] where also possible selfirradiation of accreting matter was briefly discussed. 9 Note  that the self-irradiation of the disk is possible also due to the unstable SPOs, but they could be send away from the sphere (to infinity) due to any small perturbative influencefor this reason we focus our attention on the influence of the stable spherical photons on the Keplerian disk in a special and observationally interesting case of the oscillatory photon orbits that return to a fixed azimuthal position, and the self-irradiation thus occurs repeatedly at the fixed position relative to distant observers.

Connection of the spherical photon orbits and the stable circular geodesics
The existence of the trapped photon orbits motivates us to examine the spread of the SPOs in relation to a possible distribution of radiating matter related to accretion disks. Of special interest are the stable equatorial circular orbits of test particles, which are assumed to be governing the structure of the thin accretion (Keplerian) disks. Particularly important are the marginally stable circular orbits, representing the inner boundary of the Keplerian disks -if a mass element of a Keplerian disk reaches the marginally stable orbit after loss of the energy and angular momentum due to viscous friction, any additional loss of its energy causes its direct free fall onto the naked singularity (black hole). The thorough discussion of the equatorial circular orbits in the field of the Kerr naked singularities was done in [40]. Since the SPOs can exist only at radii r < r ph− , we can restrict our attention to the equatorial circular orbits of the so called first family, which are co-rotating at large distance from the KNS, but they could become counter-rotating at close vicinity of the ring singularity; the unstable orbits of the first family could extend down to the ring singularity [40]. On the contrary, the circular orbits of the second family, which are all counterrotating, are located at r > r ph− and hence are irrelevant for the study of the interaction of the SPOs with the Keplerian disks. The marginally stable circular orbits of the first family are determined by [8,40] where We introduce the subscript 'msc' in order to distinguish the marginally stable SPOs. In the case of thick accretion disks governed by axi-symmetric toroidal structures of perfect fluid, the radiation matter could reach the radius of the marginally bound orbit at [8,40] In the present paper, we focus our attention on the Keplerian disks, because the case of toroidal accretion configurations requires more detailed study of the relation of the SPOs and the orbiting matter, as in some Kerr spacetimes the existence of the SPOs could be excluded by the presence of the accretion torus. The stable circular orbits of the first family are located at r > r msc . It is relevant from the astrophysical point of view that these stable circular orbits interfere with the area of the stable SPOs at r < r ms for KNS spacetimes with rotation parameter 1 ≤ a ≤ a i = 4.468, since just in these spacetimes there is r ms > r msc . The limits a = 1, a i correspond to equality r ms = r msc = 1 for a = 1, and r ms = r msc = 3.667 for a = a i . Moreover, due to the results of [40], it follows that even the circular orbits with negative energy with respect to infinity (E < 0) lie at the region of stable spherical orbits. Their radii satisfy inequality r z E1 < r < r z E2 , where are radii of the zero energy orbits. We can expect observationally significant optical and astrophysical effects, particularly in the region of negative energy orbits connected with an interplay between the possible extraction of the rotational energy of the KNS (Kerr superspinar) and a subsequent trapping of the radiated energy, both for the radiated electromagnetic and gravitational waves. Both forms could substantially violate the structure and stability of the accretion disk, generating thus an observationally relevant feedback, which we presume to be the subject of our further study. Of course, similar effects can be expected in the whole overlap region of the stable circular orbits of matter and the stable SPOs, where repeated re-absorption of the radiated heat, and self-illumination and self-reflection phenomena take place, as indicated in the preliminary study of the Kerr superspinars in [50]. The illustration of the overlap region of the stable circular orbits with the stable SPOs, and the other relevant orbits in dependence on the Kerr spacetime rotation parameter a is given in Fig. 21.
Now we concentrate attention on the special case of the repeated irradiation of a Keplerian disk at an azimuthal position fixed to the distant observers that could be observationally relevant, enabling even an efficient determination of spacetime parameters of the KNS (Kerr superspinar) spacetimes allowing existence of this effect.

Time periods of the latitudinal oscillations and time sequences related to the oscillatory orbits
In order to have a clear observational signature of the effects of the SPOs, we have to calculate dependence of the time period of the nodal motion, i.e., the time interval measured by a distant static observer between the two subsequent crossing of the equatorial plane by the photons following the spherical orbits. We first give a general formula and then we discuss the case of the closed spherical orbits finishing their azimuthal motion at the starting point in the equatorial plane. Such orbits could be observationally relevant, as the time interval of the nodal motion could give exact information on the Kerr naked singularity (black hole) spin, if its mass is determined by an independent method. The time interval that elapses during one latitudinal oscillation can be computed by combining equations (16) and (18). The searched period of the latitudinal motion at a fixed radius r can be expressed by the integral where sph and q sph are taken at the fixed radius r . Of course, the nodal period, i.e. the period between two subsequent crossing of the equatorial plane, is half of the time interval given by Eq. (68). The nodal period can be expressed, using the standard procedures (see Gradshtein, S., Ryzhik,M., Tables of integrals, series and products), in terms of the elliptic integrals by the formula Here m ± are the positive/negative roots of M(m), K (s) denotes the complete elliptic integral of the first kind given by (56) and is the complete elliptic integral of the second kind. Using the general formula for the time period of the nodal motion, we are able to calculate the time period of the special nodal motion of closed oscillatory orbits with vanishing change of the azimuthal coordinate for one node. We give the time periods t nod=0 in dependence on the dimensionless spin parameter a of the spacetime in Fig. 22.
Notice that there is relatively strong decline of the time period t nod=0 with increase of the spin parameter a, although the dependence of the radius of these orbits, r φ=0 , on the spin a demonstrates only very slow descent, as shown in Figs We have found that the condition r φ=0 > r msc is satisfied for the KNS spacetimes with dimensionless spin in the interval a ∈ (a pol(max) = 1.7996, 2.47812), therefore, just for this interval of the KNS spacetimes we can consider the effect of the irradiation at a fixed azimuth. Observation of such an effect enables to estimate the spin parameter a, if the mass parameter of the spacetime can be fixed by an independent method. 10 In this case we have considered only the spherical orbits with positive covariant energy E > 0, as the oscillatory orbits with azimuthal turning point are necessarily located at r > 1. Now we have to distinguish the case of the SPOs demonstrating φ = 0 with the turning point of the azimuthal motion, and the analogous cases of the periodic SPOs where φ = k2π , with k being an integer, giving also the return to the fixed azimuthal coordinate in the equatorial plane where the Keplerian disk is located.
The case of the SPOs with t nod=2π is given in Fig. 23. We found by numerical methods that the SPOs with E > 0 and φ = 2π are limited by the radii r = 1.22854 for a = 1.09533 and, by solving Eq. 61, by the radii r = 1 for a = √ 4/3 = 1.1547. A representative of such orbits is depicted in Fig. 9 (see case k = 1).
The case of the orbits with E > 0 demonstrating t nod=4π is illustrated in Fig. 24. These orbits exist for the spin parameter 0 < a ≤ a φ=0(min) = 1.17986 with r → 2 as a → 0, and r = r φ=0(max) = 1.7147 as a = a φ=0(min) . This limit case is depicted in Fig. 10a, for the values from the inside of the interval the Fig. 9, case k = 2, is representative. In the figures we depict together with the graphs the radii corresponding to the limits of their definition range. 10 We could slightly extend the range of the considered spin of the KNS spacetime, down to the value of a φ=0(min) = 1.79857 when the SPOs with φ = 0 start to appear; however, in the case of KNS spacetimes with a ∈ (a φ=0(min) , a pol(max) ) there is an extended region of spherical orbits allowing for φ ∼ 2π that invalidates applicability of the orbits with φ = 0. Case φ = 4π. The shape of such orbits represents second picture in Fig. 9 Now we can consider also the SPOs at r < 1 (with E < 0), in the KNS spacetimes with spin 1 < a < a pol(max) , since we are not limited by the necessity of having the azimuthal turning point. Of course, for the KNS spacetimes allowing for existence of closed spherical orbits of all three types (or two of them), we can use the time sequences (the ratio of the time delays) corresponding to the relevant types of the closed orbits to obtain relevant restrictions on the allowed values of the dimensionless spin a independently of their mass parameter. For comparison, we also compute the time intervals for the closed spherical orbits in Fig. 9 in associated Table 1.
To estimate the astrophysical relevance of the time delay effect, let us present the time interval for the circular photon orbit in case of the Schwarzschild black hole: in the limit a → 0 the formula (69) for r = 3 gives the result t Schwar z = 32.6484, in accordance with the exact solution t Schwar z = 6 √ 3π obtained by integration from the Schwarzschild line element after inserting ds = dr = dθ = 0 and θ = π/2. The nodal time delay is governed by half of the quantity: t nod(Schwar z) = 16.3242. In order to have dimensional estimates governing the time delays, we have to express the resulting time delays in the dimensional form that takes in the standard units the form Therefore, in the systems with stellar mass compact object (KNS) we can estimate the time delay of the level of t nod(dim) ∼ 0.01s, but in the case of the supermassive object in the Galaxy centre we can estimate t nod(dim) ∼ 10 3 s ∼ 1/4 hour, while in the case of the galaxy M87 the central object can demonstrate time delay by three orders higher ∼ 10 6 s ∼ 12 days. Clearly, from the point of view of the observational abilities of the recent observational techniques, the most convenient candidate for testing the time delay effect seems to be the Galaxy centre SgrA* object.

Concluding remarks
We have shown that in the field of the KNS spacetimes there exist a variety of the SPOs, which are not present in the KBH spacetimes. Existence of spherical photon orbits stable relative to the radial perturbations has been demonstrated. The photons at spherical orbits located at r > 1 are of standard kind, having covariant energy E > 0, as outside the black hole horizon, but at r < 1 they must have E < 0; for the special position at r = 1, the photons at spherical orbits have E = 0. 11 The character of the spherical motion in the field of Kerr naked singularities is more complex in comparison with those of spherical orbits above the outer horizon of Kerr black holes -along with the orbits purely co-rotating or counter-rotating relative to distant observers, also orbits changing orientation of the azimuthal motion occur in the field of Kerr naked singularities. Existence of polar spherical orbits reaching the symmetry axis of the Kerr spacetime is limited to the naked singularity spacetimes with spin a < a pol(max) = 1.17996. On the other hand, the KNS spacetimes with a > a pol(max) = 1.17996 allow for existence of oscillatory orbits with azimuthal turning point that return in the equatorial plane to the fixed original azimuthal coordinate as related to distant observers.
From the astrophysical point of view it seems that the most interesting and relevant is the existence of the closed spherical photon orbits with zero total change in azimuth, which intersect themselves in the equatorial plane, where the stable circular orbits of massive particles can take place. This effect is potentially of high astrophysical relevance as it enables a relatively precise estimation of the dimensionless spin of the KNS spacetimes allowing their existence, if their mass parameter is known due to other phenomena. 11 The photons at the spherical orbits located under the inner black hole horizon have also E < 0. This is the case of the KNS spacetimes with the spin parameter a ∈ (a pol(max) = 1.7996, 2.47812) (see Fig. 21). These phenomena could be extended for the case of the closed, periodic orbits demonstrating the azimuthal angle changes φ = 2kπ with integer k, when the orbits are again closed at fixed azimuth as related to distant observers.
In the case of the effects related to a fixed azimuth as related to distant observers it is in principle possible to admit some event (e.g. a collision) in the Keplerian accretion disc causing a release of energy in the form of electromagnetic radiation that would be partly and repeatedly returned back to the same radius as the original event, and at the same azimuth as observed by a distant observer, initializing repetition of the effects under slightly modified internal conditions. It can be expected that such periodic light echo could be characteristic for particular KNS spacetime. Detailed study of the related phenomena will be the task of our future work.