arxiv: v3 [gr-qc] 18 Oct 2017

Transcription

1 CHARACTERIZATION OF NULL GEODESICS ON KERR SPACETIMES CLAUDIO F. PAGANINI, BLAZEJ RUBA AND MARIUS A. OANCEA arxiv: v3 [gr-qc] 18 Oct 2017 Albert Einstein Institute, Am Mühlenberg 1, D Potsdam, Germany Faculty of Physics, Astronomy and Computer Science, Jagiellonian University, Lojasiewicza 11, Cracow, Poland Faculty of Physics, Babes-Bolyai University, , Cluj-Napoca, Romania Abstract. We consider null geodesics in the exterior region of a sub-extremal Kerr spacetime. We show that most well known fundamental properties of null geodesics can be represented in one plot. In particular, one can see immediately that the ergoregion and trapping are separated in phase space. Furthermore, we show that from the point of view of any timelike observer outside of a black hole, trapping can be understood as a smooth sets of spacelike directions on the celestial sphere of the observer. Contents 1. Introduction 2 Overview of this paper 3 2. The Kerr Spacetime Symmetries and Constants of Motion 4 3. Geodesic Equations The Radial Equation The θ Equation 9 4. Special Geodesics Radially In-/Out-going Null Geodesics The Trapped Set T-Orthogonal Null Geodesics Trapping as a Set of Directions Framework The trapped sets in Schwarzschild The trapped sets in Kerr Application 21 Acknowledgements 22 Appendix A. 22 References 23 address: claudio.paganini@aei.mpg.de, blazej.ruba@student.uj.edu.pl, marius.oancea@aei.mpg.de. Date: October 19, 2017 File:Kerrnullgeod.tex. 1

2 2 C. F. PAGANINI, B. RUBA AND M. A. OANCEA 1. Introduction In recent years there has been a lot of progress in the Kerr uniqueness and the Kerr stability problem. Having a thorough understanding of geodesic motion and in particular the behavior of null geodesics in Kerr spacetimes is helpful to understand many of the harder problems related to these spacetimes. In this paper we study the properties of null geodesics in the exterior region in the sub-extremal case, where a [0, M). The geodesic structure of Kerr spacetimes has been subject of a lot of research. Our aim here is to give a unified and accessible presentation of the most important properties of geodesics in the Kerr spacetime, with regard to the open problems mentioned above. The Mathematica notebook that has been developed for this paper is intended to help the reader gain an intuition on the influence of various parameters on the geodesic motion, despite of the complexity of the underlying equations. It can be downloaded under the permanent link [1]. In Section 6 we explain where these plots give useful insights. An extensive discussion of geodesics in Kerr spacetimes and many further references can for example be found in [4, p.318] and [15]. See [19] for a nice treatment of the trapped set in Kerr including many explicit plots of trapped null geodesics at different radii. Here we focus more on global properties of the null geodesics and less on the details of motion. Analyzing the turning points for a dynamical system is a powerful tool to extract information about its global behaviours. For example, in any dimensional system stable bounded orbits only exist if there exist two disjoint turning points in the spacial direction between which the system can oscillate. For geodesics in Kerr this has been studied in detail by Wilkins [20]. The techniques used here are very close to that paper. Many of the results discussed in the first four sections of this note are in fact well known, however our focus is on proving that all these properties can be read of from one simple plot. Thereby making it easier for people to understand the general behaviour of null geodesics in Kerr. Hence we do show that the various plots provided in the notebook actually do cover the whole space of relevant parameters. A different representation of the forbidden regions in phase space can be found in [15, p.214] and also in [17]. The presentation chosen here is adapted to help understand the phase space decomposition used in the proof for the decay of the scalar wave in subextremal Kerr [7]. The novel material in this note which is presented in the section 5 is related to the notion of black hole shadows. The shadow of the black hole is defined as the innermost trajectory on which light from a background source passing a black hole can reach the observer. The first discussion of the shadow in Schwarzschild spacetimes can be found in [18], and for extremal Kerr at infinity it was later calculated in [2]. Analyzing the shadows of black holes is of direct physical interest as there is hope for the Event Horizon Telescope to be able to resolve the black hole in the center of the Milky Way well enough that one can compare it to the predictions from theoretical calculations, see for example [8]. This perspective has led to a number of advancements in the theoretical treatment of black hole shadows in recent years [6, 9, 10, 11, 13, 14]. The novel material in the present work is the first rigorous proof of Theorem 14. The significance of Theorem 14 is due to the fact that for any subextremal Kerr spacetime, including Schwarzschild, the past and the future trapped sets at any point are topologically an S 1 on the celestial sphere of any observer in the exterior region. We would like to stress that Theorem 14 therefore describes a property of trapping which does not change when going from Schwarzschild to Kerr. Further we present numerical evidence that the radial degeneracy for the shape of the shadow, as it exists in Schwarzschild and for observers on the symmetry axis of Kerr spacetimes, is broken away from the symmetry axis. For future observations

3 CHARACTERIZATION OF NULL GEODESICS ON KERR SPACETIMES 3 this means that the distance from the observed black hole can not be ignored when one tries to extract the black holes parameters from the shadow. Overview of this paper. In section 2 we collect some background on the Kerr spacetime. We discuss its symmetries and the conserved quantities for null geodesic which arise from them in section 2.1. In section 3 we discuss the geodesic equations in its separated form, focusing on the radial and the θ equation. For the radial equation we show that many properties of its solutions can be understood from one plot. In section 4 we use the analysis of section 3 to discuss the properties of a number of special solutions. In particular we discuss the radially in-/outgoing null geodesics (i.e. the principal null directions) and the trapped set which are of relevance in the black hole stability problem. Next we discuss the T-orthogonal null geodesics which are relevant in the black hole uniqueness problem. We want to stress that most material in the sections up to this point is well known. The purpose here is to collect these facts and show that these properties can actually be understood from a few plots. In section 5 we prove a Theorem on the topological structure of the past and future trapped sets. In the same section we also present numerical evidence for a breaking of the radial degeneracy of the black hole shadows in Kerr spacetimes. Finally in section 6 we discuss how the plots developed for this paper can be used. 2. The Kerr Spacetime The Kerr family of spacetimes describes axially symmetric, stationary and asymptotically flat black hole solutions to the vacuum Einstein field equations. We use Boyer-Lindquist (BL) coordinates (t, r, φ, θ), which have the property that the metric components are independent of φ and t. The metric has the form: g = where ( 1 2Mr Σ ) The zeros of (r) are given by: dt 2 2Mar sin2 θ 2dtdφ + Σ Σ dr2 + Σdθ 2 + sin2 θ Σ Adφ2, (2.1) Σ = r 2 + a 2 cos 2 θ, (2.2) (r) = r 2 2Mr + a 2 = (r r + )(r r ), (2.3) A = (r 2 + a 2 ) 2 a 2 sin 2 θ. (2.4) r ± = M ± M 2 a 2, (2.5) and correspond to the location of the event horizon at r = r + and of the Cauchy horizon at r = r. In the present work we are only interested in the exterior region hence r (r +, ). For our considerations it is useful to introduce an orthonormal tetrad. A convenient choice is: e 0 = 1 ( (r 2 + a 2 ) ) t + a φ, (2.6a) Σ e 1 = Σ r, e 2 = 1 Σ θ, e 3 = 1 Σ sin θ ( φ + a sin 2 θ t ). (2.6b) (2.6c) (2.6d) This frame is a natural choice as the principal null directions can be written in the simple form e 0 ± e 1. These generate a congruences of radially outgoing and ingoing null geodesics. We will come back to this fact in section 4.1. For the further

4 4 C. F. PAGANINI, B. RUBA AND M. A. OANCEA considerations we will use e 0 as the local time direction. Furthermore we define the local rotation frequency of the black hole to be: a ω(r) = r 2 + a 2, (2.7) which has the rotation frequency of the horizon as a limit for r r + : ω H = ω(r + ). (2.8) The name choice for ω(r) is motivated by noting that a particle at rest in the local inertial frame given by the tetrad (2.6) will move in the φ direction in Boyer- Lindquist coordinates with dφ dt = ω(r) with respect to an observer at rest in this frame at infinity Symmetries and Constants of Motion. The independence of the metric components of the coordinates t and φ is a manifestation of the presence of two Killing vector fields ( t ) ν, ( φ ) ν, which both satisfy the Killing equation: (µ K ν) = 0 (2.9) and they commute. Furthermore the Kerr spacetime features an irreducible Killing tensor σ µν, cf. [4, p.320]. It is symmetric and satisfies the Killing tensor equation: (α σ βγ) = 0 (2.10) but it can not be written as a linear combination of tensor products of Killing vectors. In terms of the tetrad (2.6) the Killing tensor can be written as: σ µν = a 2 cos 2 θg µν + Σ ((e 2 ) µ (e 2 ) ν + (e 3 ) µ (e 3 ) ν ). (2.11) In the limiting case a = 0 it takes the form: σ (a=0) = r 4 (dθ 2 + sin 2 dφ 2 ). (2.12) This is the tensor field obtained by taking the sum of the second tensor product of the three generators of spherical symmetry. From the geodesic equation: γ µ µ γ ν = 0 (2.13) and the Killing equations it follows that a contraction of the tangent vector γ µ of an affinely parametrized geodesic γ with a Killing tensor gives a constant of motion. In Kerr spacetimes those are the mass 1, the energy, the angular momentum with respect to the rotation axis of the black hole and Carter s constant [3]: m 2 = g µν γ µ γ ν, E = ( t ) ν γ ν, L z = ( φ ) ν γ ν, K = σ µν γ µ γ ν. (2.14a) (2.14b) (2.14c) (2.14d) It follows from equation (2.12) that in a Schwarzschild spacetime K is the square of the total angular momentum of the particle. Carter s constant is non-negative for all time like or null geodesics, which can be seen immediately from equation (2.11) and the fact that g µν γ µ γ ν 0 for any future directed causal geodesic. In the case of a 0 it is even strictly positive for any time like geodesic. It turns out that some combinations of these conserved quantities are more convenient to work with, so we give them their own names: Q = K (ae L z ) 2, (2.15) L 2 = L 2 z + Q. (2.16) One can think of L 2 as the total angular momentum of the particle, in the sense that it is this quantity that is replaced with the spheroidal eigenvalue in the potential 1 The metric satisfies the Killing tensor equation (2.10) trivially and therefore it gives us a conserved quantity as well.

5 CHARACTERIZATION OF NULL GEODESICS ON KERR SPACETIMES 5 of the wave-equation, as we will show in Section 6. One can then think of Q as the component of the angular momentum in direction perpendicular to the rotation axis of the black hole. 2 It is important though that these interpretations should not be taken to strictly, because geodesics in Kerr spacetimes do not feature a conserved total angular momentum vector. Remark 1. In contrast to K, Q is not positive anymore but from the equation of motion that we will later define in (3.1d) we get the condition that Θ 0 for any geodesic to exist at a point and hence that: holds. Q a 2 E 2 cos 2 θ (2.17) In the case of null-geodesics we can rescale the tangent vector without changing the properties of the geodesic. We will use this to reduce the number of parameters. We define the conserved quotients to be: An alternative set of conserved quotients is given by: Q = Q L 2, E = E. (2.18) z L z K = K E 2, L = L z E. (2.19) These are more commonly used in the literature as they are more suited for certain calculations. However for the present work we will mostly use the first set. Properties of the Killing Fields and Tensor. The vector field ( φ ) ν is spacelike for all r > 0. The vector field ( t ) ν is timelike in the asymptotically flat region. It becomes spacelike in the interior of the ergoregion, which is defined by the inequality g( t, t ) 0 or in terms of BL-coordinates by + a 2 sin 2 θ 0. The case of equality determines the boundary of the ergoregion which is often referred to as the ergosphere. Solving for the case of equality we get the radius of the ergosphere to be: r ergo (θ) = M + M 2 a 2 cos 2 (θ). (2.20) At the equator the ergosphere lies at r = 2M while it corresponds to the horizon r = r + on the rotation axis. As mentioned above the two Killing vector fields generate one-parameter groups of isometries. It is natural to ask if the Killing tensor present in the Kerr spacetimes can also be related to some sort of symmetry. This question can be answered using Hamiltonian formalism. For a Hamiltonian flow parametrized by λ with Hamiltonian H the derivative of any function f(x, p) is given by the Poisson bracket: df H f = {H, f} dλ p µ x µ H f x µ. (2.21) p µ Each smooth function on phase space can be taken as a Hamiltonian and therefore gives rise to a local flow. Geodesic motion is generated by the function m 2. E and L z generate translations in t and φ. In the case a = 0 the function K generates rotations in the plane orthogonal to the particle s total angular momentum vector with angular velocity equal to the angular momentum of the particle. For rotating black holes it involves a change of all spatial coordinates, but it leaves quantities conserved under geodesic flow invariant. This flow explicitly depends on fiber coordinates p µ and can not be projected to a symmetry of the spacetime manifold itself. 2 The three quantities K, Q and L 2 are often labeled differently by different authors.

6 6 C. F. PAGANINI, B. RUBA AND M. A. OANCEA 3. Geodesic Equations We now focus our attention on null geodesics. The constants of motion introduced in (2.14) can be used to decouple the geodesic equation to a set of four first order ODEs, cf. [4, p. 242 ]: Σ ṫ = AE 2MarL z, Σ φ = 2MarE + (Σ 2Mr) L z sin 2 θ, Σ 2 ṙ 2 = R(r, E, L z, K) = ((r 2 + a 2 )E al z ) 2 K, ( ) L Σ 2 2 θ2 = Θ(θ, E, L z, Q) = Q z sin 2 θ E2 a 2 cos 2 θ, (3.1a) (3.1b) (3.1c) (3.1d) where the dot denotes differentiation with respect to the affine parameter λ. 3 For L z 0 the four equations are homogeneous in L z, when written in terms of the conserved quotients. For the radial and the angular equations we have: R(r, E, L z, Q) = L 2 zr(r, E, 1, Q), (3.2) Θ(θ, E, L z, Q) = L 2 zθ(r, E, 1, Q). (3.3) Remark 2. To avoid introducing new functions whenever we change between different sets of conserved quantities we use R(r, E, L z, Q) = R(r, E, L z, K(Q, L z, E)). From the homogeneity of the equations of motion (3.1) we get that the only conserved quantities which affect the dynamics are conserved quotients like E, Q or Q E in the case of L 2 z = 0. This is due to the fact that an affine reparametrization λ αλ, γ α 1 γ changes the values of E, L z and Q while leaving the trajectories and the aforementioned quotients unchanged. The case L z = 0 can be seen as the limit of Q and E tending to infinity. In these notes we will omit a separate discussion of this case as it is essentially equivalent to the Schwarzschild case and it is not needed for the understanding of the phase space decomposition in [7]. It is sufficient to consider future directed geodesics as the past directed case follows from the symmetry of the metric when replacing (t, φ) with ( t, φ). In Schwarzschild the condition ṫ > 0 guarantees that the geodesic is future directed. For Kerr a suitable condition is to require that g( γ, e 0 ) 0 is satisfied. From that we obtain the condition for causal geodesic in the exterior region to be future directed to be: E ω(r)l z. (3.4) In terms of the conserved quotients this condition takes the form: +1 if E > ω(r) sgn(l z ) = 1 if E < ω(r) undet. if E = ω(r) (3.5) which eventually allows us to represent the pseudo potential, which will be introduced in the next subsection, for the co-rotating and the counter-rotating geodesics in the same plot The Radial Equation. In this section we characterize the radial motion by locating the turning points of a geodesic in r direction. Turning points are characterized by the fact that the component of the tangent vector γ in the radial direction satisfies ṙ = 0. From equation (3.1c) we see immediately that the radial turning points are given by the zeros of the radial function R. In the following we will investigate the existence and location of these zeros. In this section we will be working 3 The radial and the angular equation can be entirely decoupled by introducing a new non-affine parameter κ for the geodesics. It is defined by dκ dλ = 1 Σ.

7 CHARACTERIZATION OF NULL GEODESICS ON KERR SPACETIMES 7 with the conserved quantity Q as it is well suited to describe the phenomena we are interested here, namely the trapping and the null geodesics with negative energy. Lemma 3. R(r, E, L z, Q) is strictly positive in the exterior region for Q < 0. Proof. The radial function can be written as: R(r, E, L z, Q) = E 2 r 4 + (a 2 E 2 Q L 2 z)r 2 + 2MKr a 2 Q, (3.6) which is clearly positive for large r. For the proof we make use of the Descartes rule, which states that if the terms of a polynomial with real coefficients are ordered by descending powers, then the number of positive roots is either equal to the number of sign differences between consecutive nonzero coefficients, or is less than it by an even number. Powers with zero coefficient are omitted from the series. For a proof of Descartes rule see for example [5, p.172]. Applied to (3.6) with Q < 0 we get that for two zeros of R to exist in r (0, ) the conserved quantities of the geodesic have to satisfy the inequality: a 2 E 2 Q L 2 z < 0. (3.7) Otherwise there are no zeros at all and the proposition is true. Assume the contrary. Then for geodesics with certain parameters to exist at a given point, additionally to R 0 we also need to have that Θ 0. Applying this condition to equation (3.1d) and combining it with inequality (3.7) we obtain the following estimate: This is clearly a contradiction. cos 4 θa 2 E 2 Q > 0. (3.8) Lemma 3 tells us that geodesics with Q < 0 either come from past null infinity I and cross the future event horizon H + or come out of the past event horizon H and go to future null infinity I +. We will discuss the property of these geodesics in section 4.1. For the rest of this section we will restrict to the case of Q 0. Note that even though some of the discussions and proofs might be simpler when working with K we are going to work with Q for all discussions that feed into the plot in section 6. This seems to be the natural choice to describe trapping and the ergoregion in phase space. To find the essential properties of the radial motion, we use pseudo potentials. The pseudo potential V (r, Q) is defined as the value of E such that the radius r is a turning point. In other words it is a solution to the equation: R(r, V (r, Q), 1, Q) = Σ 2 ṙ 2 = 0. (3.9) This equation is quadratic in V (r, Q) and for non-negative Q there exist two real solutions at every radius, denoted by V ±. They are given by: V ± (r, Q) = 2Mar ± r ((1 + Q)r 3 + a 2 Q(r + 2M)) r[r(r 2 + a 2 ) + 2Ma 2. (3.10) ] Remark 4. The pseudo potentials should not be mistaken for potentials known from classical mechanics, where the equation of motion is given by 1 2ẋ2 +V (x) = E. However the potentials of classical mechanics can always be considered as pseudo potentials in the above sense. The radial function can be rewritten as: R(r, E, 1, Q) = L 2 zr[r(r 2 + a 2 ) + 2Ma 2 ] (E V + (r, Q)) (E V (r, Q)). (3.11) This form of R reveals the significance of the pseudo potentials: The only turning points that can exist for fixed Q > 0 are those where either E = V + (r, Q) or E = V (r, Q). Analyzing the properties of the V ± allows us to extract all the information we are interested in.

8 8 C. F. PAGANINI, B. RUBA AND M. A. OANCEA First we note that for r big enough we have that V + > 0 and V < 0 for all Q 0. However in the limit we have that: lim V ± = 0. (3.12) r At the horizon the limit of the pseudo potential and its derivative are given by: lim V ± (r) = ω H (3.13) r r + dv ± lim (r) = ±. (3.14) r r + dr Lemma 5. For a fixed value of Q the pseudo potentials V ± (r, Q) have exactly one extremum as a function of r in the interval (r +, ). Proof. It is clear from the above properties that V + (V ) has at least one maximum (minimum) in the DOC. From the fact that the two pseudo potentials have the same limiting value at and at r + together with (3.11) we get that in both limits we have that R(r, E, 1, Q) 0. Therefore R(r, E, 1, Q) has to have an even number of zeros in the interval (r +, ). Given the fact that R(r, E, 1, Q) is a fourth order polynomial it can have at most 4 zeros. From the asymptotic behaviour of the potentials V ± we get that they need to have an odd number of extrema. Therefore if for some value of Q one of the potentials has more than one extremum there exists E such that R(r, E, 1, Q) has three zeros in r (r +, ). Applying Descartes rule to (3.1c) we infer that R(r, E, 1, Q) can have at most three zeros in r [0, ). But R(0, E, 1, Q) 0 and R(r +, E, 1, Q) 0. Hence there is at least one zero of R(r, E, 1, Q) in the interval [0, r + ], so it is impossible for R(r, E, 1, Q) to have three zeros in (r +, ). From that it follows that V ± can both only have one extremum in that interval. Stationary points occur at the extrema of the pseudo potentials. So Lemma 5 tells us that for every fixed value of Q 0 there exist exactly two spherical geodesics with radii r = r ± trap and energies E trap± = V ± (r ± trap). They will be studied in depth in section 4.2. Bounded geodesics with non-constant r would only be possible between two extrema of one of the pseudo potentials. These are excluded by the Lemma. From (3.1c) we have that for any geodesics to exist we have to have R(r, E, 1, Q) 0. This condition is satisfied except if V (r, Q) < E < V + (r, Q). This set is therefore a forbidden region in the (r, E) plane. Furthermore it follows that R(r, E, 1, Q) 0 for E = ω(r) with equality only in the limits r r + and r. Therefore we have that: V (r) ω(r) V + (r). (3.15) Again, the equality can occur only at the horizon and in the limit r. This fact combined with (3.5) shows that for future pointing null geodesics { +1 if E V + (r), sgn(l z ) = (3.16) 1 if E V (r). Therefore the pseudopotential V + determines the behaviour of co-rotating null geodesics and V that of counter-rotating ones. Furthermore it is worth noting that E V + (r) implies E > 0. Finally we observe that for every fixed radius r r + we get from inspection of (3.10) that: V Q 0 V + (3.17) Q holds. The equality in the relation occurs again only in the limits r r + and r. This means that for every radius r > r + the range of forbidden values of E is strictly expanding as Q increases. This fact will be quite useful for the considerations in section 4.2.

9 CHARACTERIZATION OF NULL GEODESICS ON KERR SPACETIMES 9 Figure 1. This figure shows shapes of function Θ L 2 z of values of conserved quotients. for four choices 3.2. The θ Equation. In Schwarzschild spacetimes, due to spherical symmetry the motion of any geodesics is contained in a plane. This means that for every geodesic there exists a spherical coordinate system in which it is constrained to the equatorial plane θ = π 2. This is no longer true in Kerr spacetimes, but most geodesics are still constrained in θ direction. The allowed range of θ is obtained by solving the inequality Θ(θ, E, L z, Q) 0. After multiplication with sin 2 (θ), Θ(θ, E, L z, Q) can be expressed as a quadratic polynomial in the variable cos 2 (θ). Hence Θ(θ, E, L z, Q) = 0 has two solutions given by: cos 2 θ turn = a2 E 2 L 2 z Q ± (a 2 E 2 L 2 z Q) 2 + 4a 2 E 2 Q 2a 2 E 2. (3.18) For Q > 0 only the solution with the plus sign is relevant and the motion will always be contained in a band θ min < θ eq < θ max symmetric about the equator θ eq = π 2. As L z increases, this band shrinks. In fact only in the case L z = 0 it is possible for a geodesic to reach the poles θ = 0, θ = π. Otherwise Θ(θ, E, L z, Q) blows up to there. If Q < 0 both solutions are positive and the inclination of the geodesic with respect to the equator is also constrained away from the equator, so either θ eq < θ min < θ max or θ eq > θ max > θ min. These trajectories will be contained in a disjoint band which intersects neither the equator nor the pole. This band can degenerate to a point, i.e. there exist null geodesics which stay at θ = const. The relevance of these trajectories and how they are connected to the Schwarzschild case will be discussed in the next section. All possibilities for the potentials that constrain the motion in θ direction are summarized in the Figure Special Geodesics We will now apply the discussions of the last section to describe a number of special geodesics in Kerr geometries. All of these are in some way related to either the black hole stability problem or the black hole uniqueness problem Radially In-/Out-going Null Geodesics. In this section we find geodesics which generalize the radially ingoing and outgoing congruences in Schwarzschild spacetimes. In section 3.1 we saw that the geodesics with Q < 0 extend from the horizon to infinity. In section 3.2 we saw that Q < 0 is again a special case, as these null geodesics can never intersect the equator and in the extreme case are even constrained to a fixed value of θ. At first this behaviour seems odd, but a similar situation can be observed in Schwarzschild. If we look at geodesics which move in a plane with inclination θ 0 about the equatorial plane we see that there exists a set of

10 10 C. F. PAGANINI, B. RUBA AND M. A. OANCEA null geodesics with similar properties as the ones with Q < 0 in Kerr. It is clear that the radially ingoing geodesic which moves orthogonally to the axis around which the plane of motion was rotated, moves at fixed θ value, namely that at which the plane is inclined with respect to the equatorial plane, hence θ = π/2 ± θ 0. Furthermore some null geodesics reach the horizon before intersecting the equatorial plane. They don t necessarily move at fixed θ but their motion in θ direction is still constrained away form the equatorial plane and away from the poles of the coordinate system. Now we want to investigate the null geodesics which move at fixed θ in Kerr. Demanding that θ = const. is equivalent to requiring Θ = d dθ Θ = 0. From these conditions we obtain: L z = ±ae sin 2 θ, Q = a 2 E 2 cos 4 θ, K = 0. (4.1a) (4.1b) (4.1c) By choosing the plus sign for L z in the above equation, it follows from the remaining equations of motion that: φ ṫ = dφ dt = ω(r), (4.2a) ṙ ṫ = dr dt = ± r 2 + a 2. (4.2b) This congruence is generated by the principal null directions e 0 ± e 1. When choosing the minus sign for L z in (4.1) the remaining equations of motions can not be simplified in a similar manner. In the case a = 0 these are the radially in-/outgoing null geodesics. An interesting observation is that along these geodesics the inner product of the ( t ) µ vector field is monotone. A simple calculation shows that: γ µ µ (( t ) ν ( t ) ν ) = ṙ 2M(r2 a 2 cos 2 θ) Σ 2 + θ 2Ma2 r sin 2θ Σ 2. (4.3) For the principal null congruence we have θ = 0, the coefficient of ṙ is positive and there is no turning point in r. This property might be interesting in the context of the black hole uniqueness problem. If one could show a similar monotonicity statement for a congruence of null geodesics in general stationary black hole spacetimes, one could conclude that the ergosphere in such spacetimes has only one connected component enclosing the horizon. This is a necessary condition if one wants to show that no trapped T-orthogonal null geodesics can exist in that case The Trapped Set. One of the most important features of geodesic motion in black hole spacetimes is the possibility of trapping. A geodesic is called trapped if its motion is bounded in a spatially compact region away from the horizon. In Kerr this corresponds to the geodesics motion being bounded in r direction. This is only possible if r = const. or if the motion is between two turning points of the radial motion. For null geodesics in Kerr we ruled out the second option in Lemma 5. We will now discuss orbits of constant radius. 4 These null geodesics are stationary points of the radial motion, hence null geodesics with ṙ = r = 0. Dividing equation (3.1c) by Σ and taking the derivative with respect to λ we see that this condition is equivalent to R(r) = d dr R(r) = 0. The solutions to these equations can be parametrized explicitly by, cf. [19]: ( a(r M) E trap (r) = = ω(r) 1 2r ) A(r) A(r) (4.4) 4 Null geodesics of constant radius are often referred to as "spherical null geodesics" but it is important to note that r = const. does not imply that the whole sphere is accessible for such geodesics.

11 CHARACTERIZATION OF NULL GEODESICS ON KERR SPACETIMES 11 Figure 2. Plot of the relation between the radius of the equatorial trapped null geodesics at r 1 and r 2, the trapped null geodesic with L z = 0 at r 3 and the horizon at r + for all values of a. Q trap (r) = B(r) A 2 (r) (4.5) A(r) = r 3 3Mr 2 + a 2 r + a 2 M = (r r 3 )P 2 (r) (4.6) B(r) = r 3 (r 3 6Mr 2 + 9M 2 r 4a 2 M) = (r r 1 )(r r 2 )P 4 (r) (4.7) where P 2 and P 4 are polynomials in r, quadratic and quartic respectively, which are strictly positive in the DOC. The following three radii are particularly important: ( ( 2 r 1 = 2M 1 + cos ( 3 arccos a ) )) (4.8) M ( ( 2 ) )) r 2 = 2M 1 + cos (4.9) r 3 = M + 2 satisfying the inequalities: M 2 a2 3 cos ( a 3 arccos M ( ( )) 1 3 arccos M(M 2 a 2 ) (M 2 a2 3 ) 3 2 (4.10) M < r + < r 1 < r 3 < r 2 < 4M (4.11) for a (0, M). Orbits of constant radius are allowed only inside the interval [r 1, r 2 ], because outside of it Q would have to be negative. This possibility has already been excluded in Lemma 3. The boundaries of the interval at r = r 1 and r = r 2 correspond to circular geodesics constrained to the equatorial plane with Q = 0. The trapped null geodesics at r = r 3 have L z = 0 which is the reason why the functions E trap and Q trap blow up there. From the second representation in (4.4) we see that E trap (r) ω(r) is positive in [r 1, r 3 ) and negative in (r 3, r 2 ]. Combined with (3.16) this implies that the stationary points in [r 1, r 3 ) correspond to extrema of V + and the stationary points in (r 3, r 2 ] correspond to extrema of V. In Lemma 5 we showed that V + and V both have exactly one extremum. Since extrema of the pseudo potentials always correspond to orbits of constant radius, we get that the extrema of V + (r, Q) and V (r, Q) have to be within the intervals [r 1, r 3 ) and (r 3, r 2 ] respectively for any value of Q. In Figure 2 we plot the behaviour of these intervals as a function of a/m.

12 12 C. F. PAGANINI, B. RUBA AND M. A. OANCEA We now know that the maps given by: [0, ) Q r + trap [r 1, r 3 ) [0, ) Q r trap (r 3, r 2 ] which take Q into radii of trapped geodesics corresponding to the unique maximum of V + (r, Q) and minimum of V (r, Q) respectively are one-to-one and therefore monotone. By using (4.5) the sign of their derivatives can be easily evaluated in some ɛ-neighbourhood of r = r 3 where the term of highest order in 1 r r 3 dominates: r trap Q < 0 < r+ trap Q. (4.12) By the equation (3.17) and the fact that radii of trapping always correspond to global extrema of the pseudo potentials we get that: Q E trap(rtrap(q)) < 0 < Q E trap(r trap(q)). + (4.13) Using the chain rule and combining these two facts we obtain: E trap > 0. (4.14) r These inequalities provide an important piece of the picture of the influence of Q on the trapped geodesics. We have Q = 0 for the outermost circular geodesics and as we increase it, the radii of trapping converge towards r = r 3 while E blows up to ±, with the sign depending on the direction from which we approach r 3. We can also desciribe the function E trap (r): it starts with some finite positive value at r = r 1 and increases monotonically to + as r approaches r 3. There it jumps to and increases again to a finite negative value at r = r 2. It is interesting to ask what region in physical space is accessible for trapped geodesics. By plugging (4.5) and (4.4) into the equation Θ = 0 we get that for a geodesic with r = const.: cos 2 θ turn = 2 Mr 2 (2r 3 3Mr 2 + a 2 M) r(r 3 3M 2 r + 2a 2 M) a 2 (r M) 2 (4.15) holds. This gives two turning points in θ direction which are symmetric about the equatorial plane. The whole region of trapping in the (r, θ) plane is bounded by curves defined implicitly by (4.15) and r 1 r r 2. Figure (3) presents this set for a particular value of a. Remark 6. Two warnings: (1) One has to be careful when interpreting Figure 3 (and the plots in the Mathematica notebook). Despite the fact that the region in physical space occupies finite range of r values, every individual trapped null geodesic is still constrained to a fixed radius. For an insight on what those trajectories look like in detail we recommend the study of [19]. (2) When taking a M in the Mathematica notebook the ergosphere appears to develop a kink on the rotation axis. This is an artifact of the coordinate system, as the ergosphere coincides with the horizon there and is thus orthogonal to itself T-Orthogonal Null Geodesics. In the ergoregion there exist null geodesics with negative values of E. In physical space they are constrained to the region defined by equation (2.20). From Lemma 3 we know that geodesics with Q < 0 reach either I + or come from I and can therefore not have negative values of E. This allows us to use the pseudo potentials to give a more precise characterization of the ergoregion in phase space. It is located in the region where V (Q) > 0,

13 CHARACTERIZATION OF NULL GEODESICS ON KERR SPACETIMES 13 Figure 3. The region accessible for trapped null geodesics for a = The shaded region represents the black hole, r r +. The only qualitative change in this picture occurs at a = 1 2 because at this value the region of trapping starts intersecting with the ergoregion. between E = 0 and V (Q). An immediate consequence of that is, that all future pointing null geodesics with negative E begin at the past event horizon and end at the future event horizon. Furthermore they must have L z < 0. Those null geodesics with E = 0 can reach the boundary of the ergoregion. In this case equation (3.1d) gives us, that: Q = cos2 (θ max ) sin 2 (θ max ). (4.16) When calculating the turning points from equation (3.1c) we get that: ( ) sin 2 (θ max )R r, 0, 1, cos2 (θ max ) sin 2 = r 2 + 2Mr a 2 cos 2 (θ max ) = 0. (4.17) (θ max ) The only solution to this equation in the exterior region is: r turn (θ max ) = M + M 2 a 2 cos 2 (θ max ) (4.18) which is exactly the location of the ergo sphere (2.20). So V (Q) > 0 can be considered as the boundary of the ergoregion in phase space. From this considerations we see immediately that T-orthogonal null geodesics are clearly non-trapped in Kerr. In fact there do not even exist any trapped null geodesics orthogonal to: where ɛ min = min[ V + (0, r 1 ), V (0, r 2 ) ]. K ν = ( t ) ν + ɛ min ( φ ) ν (4.19) 5. Trapping as a Set of Directions In this section we will link the previous discussion to the black hole shadows. We introduce a more formal framework for the discussion. This allows us to give a more technical discussion of the trapped sets in Schwarzschild and Kerr black holes.

14 14 C. F. PAGANINI, B. RUBA AND M. A. OANCEA 5.1. Framework. First we have to introduce the basic framework and notations. Let M be a smooth manifold with Lorenzian metric g. At any point p in M you can find an orthonormal basis (e 0, e 1, e 2, e 3 ) for the tangent space, with e 0 being the timelike direction. It is sufficient to treat only future directed null geodesics as the past directed ones are identical up to a sign flip in the parametrization. The tangent vector to any future pointing null geodesic can be written as: γ(k p ) p = α(e 0 + k 1 e 1 + k 2 e 2 + k 3 e 3 ) (5.1) where α = g( γ, e 0 ) and k = (k 1, k 2, k 3 ) satisfies k 2 = 1, hence k S 2. The geodesic is independent of the scaling of the tangent vector as this corresponds to an affine reparametrization for the null geodesic. We will therefore set α = 1 in the following discussion. The S 2 is often referred to as the celestial sphere of a timelike observer at p whose tangent vector is given by e 0, cf. [16, p.8]. For the further discussions we fix the tetrads. We can make the following definition: Definition 7. Let γ(k p ) denote a null geodesic through p for which the tangent vector at p is given by equation (5.1). It is clear that γ(k a p ) and γ(k b p ) are equivalent up to parametrization if k a = k b. Suppose now that M is the exterior region of a black hole spacetime with complete future and past null infinity I ± and a boundary given by the future and past event horizon H + H. We can then define the following sets on S 2 at every point p. Definition 8. The future infalling set: Ω H +(p) := {k S 2 γ(k p ) H + }. The future escaping set: Ω I +(p) := {k S 2 γ(k p ) I + }. The future trapped set: T + (p) := {k S 2 γ(k p ) (H + I + ) = }. The past infalling set: Ω H (p) := {k S 2 γ(k p ) H }. The past escaping set: Ω I (p) := {k S 2 γ(k p ) I }. The past trapped set: T (p) := {k S 2 γ(k p ) (H I ) = } We finish the section by defining the trapped set to be: Definition 9. The trapped set: T(p) := T + (p) T (p). The region of trapping in the manifold M is then given by: Definition 10. Region of trapping: A := {p M T(p) } The trapped sets in Schwarzschild. The discussion of Schwarzschild serves as an easy example for the various concepts. Definition 11. We refer to the set Ω H (p) T (p) as the shadow of the black hole. However for any practical purposes information about its boundary which is given by T (p) is good enough. An observer in the exterior region can only observe light from background sources along directions in the set Ω I (p). From the radial equation we get immediately that if k = (k 1, k 2, k 3 ) T + (p) then k = ( k 1, k 2, k 3 ) T (p). Hence the properties of the past and the future sets are equivalent. This is true both in Schwarzschild and Kerr. An explicit formula for the shadow of a Schwarzschild black hole was first given in [18]. In Schwarzschild the orthonormal tetrad (2.6) reduces to: e 0 = 1 1 2M/r t, e 1 = 1 2M/r r, (5.2) e 2 = 1 r θ, e 3 = 1 r sin θ φ. To determine the structure of T ± (p) in Schwarzschild it is sufficient to consider p in the equatorial plane and k = (cos Ψ, 0, sin Ψ) with Ψ [0, π]. The entire sets T ± (p) are then obtained by rotating around the e 1 direction. Note that from

15 CHARACTERIZATION OF NULL GEODESICS ON KERR SPACETIMES (a) (b) 15 (c) Figure 4. The trapped set on the celestial sphere of a standard observer at different radial location in a Schwarzschild DOC. Observer (a) is located outside the photon sphere at r = 3.9M, observer (b) is located on the photon sphere at r = 3M and finally observer(c) is located between the horizon and the photonsphere at r = 2.5M. One can see that the future trapped set moves from the ingoing hemisphere in (a) to the outgoing hemisphere in (c) as one varies the location of the observer. The future and past trapped set coincide on the r = 0 line when the observer is located on the photon sphere at r = 3M in (b) the tetrad it is obvious that E(k) = E(r) is independent of Ψ. On the other hand Lz (k) is zero for Ψ = {0, π} and maximal for Ψ = π/2. Away from that maximum, Lz is a monotone function of Ψ. Note that the geodesic that corresponds to Ψ = π/2 has k1 = 0 and thus a radial turning point. Thus the E/Lz value of this geodesic corresponds to the minimum value any geodesic can have at this point in the manifold. For r 6= 3M this is smaller then the value of trapping and thus there exist two k with the property that E/Lz (k) = 1/ 27M 2. One of them has k1 > 0 and therefore r > 0 and one has r < 0. For r > 3M the first corresponds to T (p) and the second corresponds to T+ (p). For 2M < r < 3M the roles are reversed. For r = 3M we have T+ (p) = T (p) = (0, k2, k3 ). In Figure 4 we depict these three cases for some fixed radii. To conclude we see that T+ (p) and T (p) are circles on the celestial sphere independent of the value of r(p). In [16, p.14] it is shown that Lorentz transformations on the observer correspond to conformal transformations on the celestial sphere. They are equivalent up to Möbius transformations on the complex plane when considering the S2 to be the Riemann sphere. Hence circles on the celestial sphere stay circles under Lorentz transformations. As a consequence if r(pa ) 6= r(pb ) then there exists a Lorentz transformation (LT) such that T (pa ) = LT[T (pb )]. This concept is sufficiently important that it deserves a proper definition. Definition 12. The shadows at two points p1, p2 are called degenerate if, upon identification of the two celestial spheres by the orthonormal basis, there exists an element of the conformal group on S 2 that transforms T (p1 ) into T (p2 ). Remark 13. The shadow at two points p1, p2 being degenerate implies that for every observer at p1 there exists an observer at p2 for which the shadow on S 2 is identical. Because this notion compares structures on S 2, the two points do not have to be in the same manifold for their shadows to be degenerate. Just from the shadow alone an observer can not distinguish between these two configurations. The only reliably information an observer can obtain is thus, that T (p) is a proper circle on his celestial sphere The trapped sets in Kerr. We will now discuss the properties of the sets T± (p) in Kerr. Note that the equations of motion for r (3.1c) and θ (3.1d) have two solutions that differ only by a sign for a fixed combination of E, Lz, K. Therefore

16 16 C. F. PAGANINI, B. RUBA AND M. A. OANCEA we know that the trapped sets will have a reflection symmetry across the k 1 = 0 and the k 2 = 0 planes. A sign change in k 2 maps both sets T ± (p) to themselves, while a sign flip in k 1 maps T + (p) to T (p) and vice versa. In the following we will use the parametrization of [10]. We introduce the coordinates ρ [0, π] and ψ [0, 2π) on the celestial sphere. Thus (5.1) can be written as: γ(ρ, ψ) p = α(e 0 + cos(ρ)e 1 + sin(ρ) cos(ψ)e 2 + sin(ρ) sin(ψ)e 3 ) (5.3) The direction towards the black hole is given by ρ = π. Following [10] one finds the following parametrization of the celestial sphere in terms of constants of motion: sin(ψ) = (L a) + a cos2 (θ) K sin(θ) K sin(ρ) = r 2, a(l a) p, (5.4a) p (5.4b) Analog to the functions (4.4) and (4.5) which give the value of the conserved quantities in terms of the radius of trapping, we can give such relations for the conserved quotients K and L in this formulation they can be found for example in [10]. To differ between the trapped radius and the observers radius we use x [r min (θ), r max (θ)] to parametrize the conserved quantities of the trapped set. Here r min (θ) and r max (θ) are given as the intersection of a cone of constant θ with the boundary of the area of trapping. Note that it is part of the following proof to show that the given interval is in fact the correct domain for the parameter x. The parameter x here corresponds to the radius of the trapped null geodesic which a particular future trapped direction is asymptoting to. We have then: K = 16x2 (x) ( (x)) 2 (5.5a) ( a(l a) = x 2 4x (x) ) (5.5b) (x) where (x) = 2x 2M is the partial derivative of (x) with respect to x.plugging (5.5) into (5.4) we obtain: sin(ψ) = (x){x 2 + a 2 cos 2 (θ(p))} 4x (x) 4x (x)a sin(θ(p)) sin(ρ) = 4x (r(p)) (x) (x)(r(p) 2 x 2 := h(x) ) + 4x (x) (5.6b) We are now ready to prove the following Theorem. (5.6a) Theorem 14. The sets T + (p) and T (p) are smooth curves on the celestial sphere of any timelike observer at any point in the exterior region of any subextremal Kerr spacetime. Proof. We start by anayzing the right hand side of (5.6a): ( ) d (x){x 2 + a 2 cos 2 (θ(p))} 4x (x) dx 4x (x)a sin(θ(p)) = {x2 + a 2 cos 2 (θ(p))}((m x) 3 M(M 2 a 2 )), 2x 2 (x) 3/2 a sin(θ(p)) (5.7) which is strictly negative for x (r +, ). The limit of the right hand side of (5.6a) is given by for x r + and for x. Therefore the right hand side is invertible and x(sin(ψ)) is a smooth function of ψ with extrema at the extremal points of sin(ψ). As was shown in [10] the minimum x min (θ(p)) at ψ = π/2 and the maximum of x max (θ(p)) at ψ = 3π/2 correspond exactly to the intersections

17 CHARACTERIZATION OF NULL GEODESICS ON KERR SPACETIMES 17 of a cone with constant θ with the boundary of the region of trapping. This can be seen by setting the left hand side of (5.6a) equal to ±1, taking the square of the equation, solving for cos 2 (θ) and comparing to (4.15). Important here is that [x min (θ(p)), x max (θ(p))] [r 1, r 2 ] for all values of θ(p). Now we take a look at the right hand side of equation (5.6b): d dx (h(x)) = 8(r(p)2 x 2 ) (r(p))((x M) 3 + M(M 2 a 2 )) (r(p)) (x)(4x (x) + (r(p)2 x 2 ) (x)). (5.8) 2 This is positive when x < r(p) and negative when x > r(p). The denominator never vanishes for x (r +, ) because: (4x (x) + (r(p) 2 x 2 ) (x)) {r(p)=r+,x=r +} = 0 (5.9) and d dx (4x (x) + (r(p)2 x 2 ) (x)) = 2(3x 2 6Mx + 2a 2 + r(p) 2 ) > 0 (5.10) d dr(p) (4x (x) + (r(p)2 x 2 ) (x)) = 2r(p) (x) > 0, (5.11) where we used r(p) > r + > M > a in (5.10). If we set x = r(p) in (5.6b) then the right hand side is equal to 1. Furthermore in any of the limits r(p) r +, r(p), x r +, and x it goes to zero. Case 1. If p / A hence if r(p) / [x min (θ(p)), x max (θ(p))] then the two functions: ( ρ 1 (ψ) = arcsin(h(x(ψ))) : [0, 2π) [ρ 1min, ρ 1max ] 0, π ) (5.12) ( 2 π ) ρ 2 (ψ) = π arcsin(h(x(ψ))) : [0, 2π) [ρ 2min, ρ 2max ] 2, π (5.13) are both smooth with ρ 1 (0) = ρ 1 (2π) and ρ 2 (0) = ρ 2 (2π). If p is between the region of trapping and the asymptotically flat end, the function ρ 2 (ψ) corresponds to T + (p) and ρ 1 (ψ) corresponds to T (p). Because (π/2, π] corresponds to the geodesic with ṙ < 0. If p is between the region of trapping and the horizon then the role of ρ 1 (ψ) and ρ 2 (ψ) are switched. Case 2. If p A we need to do some extra work. For simplicity we only consider the interval ψ [π/2, 3π/2] as the rest follows by symmetry of sin(ψ) in [0, π] across π/2 and in [π, 2π] across 3π/2. Let arcsin(x) map into this interval, then we define: ( ) (r(p)){r(p) 2 + a 2 cos 2 (θ(p))} 4r(p) (r(p)) ψ 0 (r(p)) = arcsin 4r(p). (5.14) (r(p))a sin(θ(p)) The two functions: { arcsin(h(x(ψ))) ρ 3 (ψ) = π arcsin(h(x(ψ))) { π arcsin(h(x(ψ))) ρ 4 (ψ) = arcsin(h(x(ψ))) if ψ [π/2, ψ 0 (r(p))] if ψ (ψ 0 (r(p)), 3π/2] if ψ [π/2, ψ 0 (r(p))] if ψ (ψ 0 (r(p)), 3π/2] (5.15) (5.16) are then smooth on [π/2, 3π/2]. For a proof see Appendix A and note that at ψ 0, h(x(ψ)) satisfies the conditions required in the appendix. Since p A we have that x min (θ(p)) < r(p) < x max (θ(p)). Therefore the geodesic on the celestial sphere parametrized by x max (θ(p)) has to have ṙ > 0 and thus has to be in [0, π/2). On the other hand the geodesic on the celestial sphere parametrized by x min (θ(p)) has to have ṙ < 0 and thus has to be in (π/2, π]. In fact by the monotonicity of the right hand side of (5.6a) and the fact that x(ψ 0 ) = r(p) we know that for ψ [π/2, ψ 0 ) we have x(ψ) < r(p) and for ψ (ψ 0, 3π/2] we have x(ψ) > r(p). Thus we can conclude that for p A, ρ 4 corresponds to T + (p) and ρ 3 corresponds to T (p) and thus both sets are smooth.

18 18 C. F. PAGANINI, B. RUBA AND M. A. OANCEA Case 3. In the special case when r(p) = x max (θ(p)) or r(p) = x min (θ(p)) the functions ρ 1 and ρ 2 describing T ± (p) do reach ρ = π/2 at ψ = 3π/2 or ψ = π/2 respectively. However since in these cases we have that: d 2 (h(x(ψ))) = 0, (5.17) dψ2 the two sets meet at this point tangentially and do not cross over into the other hemisphere. This concludes the proof. Remark 15. In [10] it was observed that ρ max of T + (p) always corresponds to the trapped geodesic with x min (θ(p)) and ρ min of T + (p) always corresponds to the trapped geodesic with x max (θ(p)). When p is outside the region of trapping h(x) xmax is a local maximum of h(x(ψ)) (as a function of ψ) and h(x) xmin is a local minimum of h(x(ψ)). When p is between the region of trapping and the horizon h(x) xmax is a local minimum of h(x(ψ)) and h(x) xmin is a local maximum of h(x(ψ)). Since outside T + (p) is always described by ρ 2 (ψ) and inside by ρ 1 (ψ), ρ min then always corresponds to x min and ρ max always corresponds to x max. This also holds for p A. This observation means that the null geodesic approaching the innermost photon orbit has the smallest impact parameter (deviation from the radially ingoing null geodesic) and the null geodesic approaching the outermost photon orbit has the largest impact parameter. For T (p) the correspondence is switched. Remark 16. The parametrization for sin(ψ) breaks down on the rotation axis, the one for sin(ρ) however remains valid with only one possible value for x = r 3. Due to the symmetry at these points we know that T ± (p) are described by proper circles on the celestial sphere. The PND is aligned with the axis of the rotation symmetry and hence the sets are symmetric under rotation along ψ. The situation is therefore equivalent to Schwarzschild and an observer can not distinguish whether he observes a Schwarzschild black hole or a Kerr black hole in the direction of the rotation axis. Remark 17. We have only proved Theorem 14 for one standard observer at any particular point. However since any other observer at this point is related to the standard observer by a Lorentz transformation and the Lorentz transformations act as conformal transformations on the celestial sphere [16, p.14], the Theorem indeed holds for any observer. In [9] the quantitative effect on the shape of the shadow of boosts in different directions are discussed. Remark 18. The parametrization for the trapped set on the celestial sphere of any standard observer in [10, 11] was derived for a much more general class of spacetimes. Therefore Theorem 14 might actually hold for these cases as well. However this is beyond the scope of this paper. From Theorem 14 we immediately get the following Corollary: Corollary 19. For any observer at any point p in the exterior region of a subextremal Kerr spacetime we have that for any k T + (p) and any ɛ > 0 B ɛ (k) Ω H +(p) B ɛ (k) Ω I +(p). So if we interpret the celestial sphere as initial data space for null geodesics starting at p, the Corollary is a coordinate independent formulation of the fact that trapping in the exterior region of subextremal Kerr black holes is unstable. See Figure 5 as an example on how the trapped sets change under a variation of the radial location of the observer.

19 CHARACTERIZATION OF NULL GEODESICS ON KERR SPACETIMES (a) (b) (d) 19 (c) (e) Figure 5. The trapped set on the celestial sphere of a standard observer at different radial locations in the equatorial plane of the exterior region of a Kerr black hole with a = 0.9. Observer (a) is located outside the region of trapping at r = 5M. Observer (b) is located on the outer boundary of the region of trapping at r = r2 = 3.535M. Observer (c) is located inside the region of trapping at r = r3 = 2.56M. Observer (d) is located on the inner boundary of the area of trapping at r = r1 = 1.73M Finally observer (e) is located between the horizon and the area of trapping at r = 1.59M. Again one can observe how the two trapped sets move in opposite directions on the celestial sphere as the observer approaches the black hole. In (a) the future trapped set is on the ingoing hemisphere and the past trapped set is on the outgoing hemispere. In (b) they meet in one point tangentially but are still entirely in one hemisphere except for that one point. In (a) the trapped sets intersect in two points and both have parts in both hemispheres. In (d) they only meet at one point tangentially again (now on the "other" side of the celestial sphere) and finally in (e) the future trapped set is entirely in the outgoing hemisphere and the past trapped set is entirely in the ingoing hemisphere. Even though the qualitative features of T± (p) do not change under a change of parameters, the quantitative features do. In [13, 14] it is discussed what information can be read of from the shadow at infinity. For points inside the manifold as considered in this work, it is unclear to the authors what the maximal amount of information is that an observer can obtain from the shape of the shadow. In the following we will present numerical calculations that suggest that the radial degeneracy of the shadow is broken in Kerr spacetimes away from the symmetry axis. We take the stereographic projection: [16, p.10] x + iy ζ= (5.18a) 1 z x = sin(ρ) sin(ψ) y = sin(ρ) cos(ψ) z = cos(ρ) (5.18b) of T (p) in order to compare the exact shape of the shadow for different observers. Here the values of ρ and ψ are given by the parametrization in (5.5). The above projection for the celestial sphere of a standard observer in the exterior region of