A generalization of photon sphere based on escape/capture cone

In asymptotically flat spacetimes, bearing the null geodesics reaching the future null infinity in mind, we propose new concepts, the"dark horizons"as generalizations of the photon sphere. They are defined in terms of the structure of escape/capture cones of photons with respect to a unit timelike vector field. More specifically, considering a two-sphere that represents a set of emission directions of photons, the dark horizons are located at positions where a hemisphere is marginally included in the capture and escape cones, respectively. We show that both of them are absent in the Minkowski spacetime, while they exist in spacetimes with black hole(s) under a certain condition. We derive the general properties of the dark horizons in spherically symmetric spacetimes and explicitly calculate the locations of the dark horizons in the Vaidya spacetime and the Kerr spacetime.


Introduction
Recently, Event Horizon Telescope (EHT) Collaboration succeeded in obtaining the image of the regions around black holes in M87 [1] and in Sagittarius A* [2].Future observations of the black hole shadow are expected to test the general relativistic magnetohydrodynamics (GRMHD) models [3], to constrain theories of gravity [4][5][6], and so forth.In static and spherically symmetric spacetimes, the edge of the shadow is given by the photon sphere under the assumption that light sources are distributed at infinity [7,8], which brings the photon sphere to attention in recent years.We note that the realistic observed shadow depends strongly on the distribution of the light source [9], but the photon sphere is helpful to understand the geometric structures near the black hole.Following the definition of Ref. [7], the photon sphere is provided as a SO(3) × R-invariant photon surface in a static and spherically symmetric spacetime.Here, the photon surface is defined as a nowhere-spacelike hypersurface S such that, any null geodesic tangent to S at any point on S is included in S. Interestingly, the outermost photon sphere has been shown to satisfy the areal inequality [10] (see also [11]) where A H is the area of the event horizon, A ph,out is the area of the outermost photon sphere, A sh,out is the area of the shadow observed at infinity, and M is the Arnowitt-Deser-Misner (ADM) mass.However, the definition of the photon sphere does not work in spacetimes with the rotation.In addition, since the observed intensity depends on the position and the motion of the light source, it would be more useful if we could take these into account in the definition.There have been several attempts to generalize the photon sphere [12][13][14][15], but most of them are defined with local geometrical quantities.The escape probability of photons to infinity (or sufficiently far away from black holes) is not very clear in such definitions because the global analysis for the geodesic is crucial for each photon coming all the way from light sources around black holes.In this sense, we will focus on generalizations of the photon sphere with an attention to the global properties, especially the escape probability of photons, associated with the motion of the light sources.So far, at least to our knowledge, there have been no generalized photon spheres that incorporate the effect of the position and the motion of the light source 1 .
In this paper, we introduce new generalizations of the photon sphere in terms of the escape probability of photons from the light sources incorporating the effect of its motion and the beaming effect due to their motion, which we name the outer dark horizon (ODH) and the inner dark horizon (IDH).These new concepts are defined by using the escape/capture cones, which depends on the frame associated with a unit timelike vector T µ .We take T µ as the tangent vector field of the trajectory of the light source (in the region where light sources exist).Together with this T µ , the escape/capture cones tell us how brightly the source is observed from infinity.Note that we can deal with the situation where the light sources are possibly distributed everywhere in the spacetime 2 .Note that all of the rigorous theorems/propositions/lemmas are irrelevant to how to choose T µ unless specified.
The rest of this paper is organized as follows.In Sec. 2, we provide a brief review of the asymptotic behavior of null geodesics near future null infinity as a preparation.In Sec. 3, we introduce the ODH and IDH as generalizations of the photon sphere in general asymptotically flat spacetimes and probe the properties of the ODH and IDH to see that the definition is fairly reasonable.In Sec. 4, we focus on a static and spherically symmetric case as a specific example and study basic properties of the ODH and IDH, mainly focusing on the relation to the photon sphere.In Sec. 5, we see their properties in the Vaidya spacetime as a simple example of dynamical spacetimes.In Sec. 6, we investigate the ODH and IDH in the Kerr spacetime as an example of spacetimes with angular momenta.The last section gives a summary and discussions.In App.A, we provide the derivation of the inequalities appeared in Sec. 6.We use the units in which the speed of light and the Newtonian constant of gravitation are unity, c = 1 and G = 1.We assume the metric to be C 2− functions (i.e., class C 1,1 ).The metric sign convention is (−, +, +, +).

Review of Null Asymptotics in the Bondi coordinate
In this section, as a preparation for the proof of Lem.3.13 and Lem.3.17, we review the asymptotic behavior of the metric and null geodesics near future null infinity in four-dimensional asymptotically flat spacetimes.The detailed analyses are presented in Refs.[24][25][26][27].
We begin with reviewing the Bondi coordinates based on Refs.[18,19], which describes the asymptotic behavior of the metric near future null infinity (for the generalization to higher dimensions, see also Refs.[20][21][22][23]).The non-zero components of the metric near future null infinity in the Bondi coordinates can be expanded in the power of 1/r as (2) Here, u denotes the retarded time, r is the areal radius, x I stands for the angular coordinates, and ω IJ is the metric for the unit two-sphere.The non-zero components of the inverse metric behave as where ω IJ is the inverse of ω IJ , and h (1)IJ is defined as h (1)IJ := ω IK ω JL h KL .Future null infinity is supposed to be in the limit of r → ∞ while u is kept finite.Let us impose the gauge condition where ω 2 is the volume element of the two-dimensional unit sphere.Here, h IJ − ω IJ corresponds to gravitational waves.In general relativity, the integration of m(u, x I ) over the angular coordinates gives us the Bondi mass, Next, let us look at the asymptotic behavior of null geodesics near future null infinity based on Refs.[24][25][26][27].By combining the r-component of the geodesic equation and the condition for the geodesic to be null, we have where Ω IJ is defined as Here, the prime denotes the derivative with respect to the affine parameter.We also define Ω u, x I ; dx J /du and Ω i as In Refs.[24,26,27], the careful analysis of the global behavior of null geodesics with Eq. (12) gives us a sufficient condition to reach future null infinity.In particular, the result for Ω i > 0, which we use later, is given by the following statement: Proposition 2.1.Consider a four-dimensional asymptotically flat spacetime with Ω i > 0 in which the metric near future null infinity is written as Eqs.
(2)-( 5) with the Bondi coordinates by C 2− functions.We define β crit as Take a point p with a sufficiently large radial coordinate value r = r 0 .Any null geodesic emanating from p reaches future null infinity if holds.
Here, Eq. ( 17) roughly means dr/du| p ≳ β crit .In the case with ∂m/∂u and ∂h (1) IJ /∂u being small enough, which corresponds to the situation that matter radiations and gravitational waves are both weak enough near future null infinity, the value of β crit is approximately given by β crit ≈ −(1 − 1/ √ 3) ≈ −0.423.By this proposition, we see that photons emitted with dr/du ≥ 0 and ones with dr/du = O r −1 0 reach future null infinity under the assumptions in the proposition.For higher dimensions, see Refs.[24,26].

Outer Dark Horizon and Inner Dark Horizon
In this section, we introduce two new concepts, "outer dark horizon (ODH)" and "inner dark horizon (IDH)," as generalizations of the photon sphere in general asymptotically flat spacetimes.The definitions of the ODH and IDH refer to future null infinity as the event horizon does.
This section is organized as follows.In subsection 3.1, we present the definitions of the ODH and IDH, and discuss their properties.In subsection 3.2, the relation of the ODH and IDH to the escape and capture cone is discussed.In subsections 3.3 and 3.4, we examine the properties of the ODH and IDH, respectively.

Basics of Outer Dark Horizon and Inner Dark Horizon
In this subsection, we introduce the outer dark horizon (ODH) and the inner dark horizon (IDH) as generalizations of the photon sphere by examining whether photons emitted from photon sources reach future null infinity or not.The ODH and IDH are defined as the boundary of "outer dark domain (ODD)" and "inner dark domain (IDD)", respectively, which are also defined below.
Let M be a four-dimensional asymptotically flat and connected spacetime, and T µ be a unit timelike vector field with T µ T µ = −13 .Definition 3.1 (outer dark domain (ODD) associated with T µ ).Suppose a set S O consists of all p ∈ M satisfying the following condition: for all spacelike vectors n µ orthogonal to T µ at p, there exists a null geodesic emanating from p, whose tangent vector k µ is orthogonal to n µ at p, such that it will not reach future null infinity.Then, we call S O the outer dark domain (ODD) associated with T µ .Definition 3.2 (outer dark horizon (ODH) associated with T µ ).Let a region S O be the ODD associated with T µ in M.Then, the boundary ∂S O is called the outer dark horizon (ODH) associated with T µ .Definition 3.3 (inner dark domain (IDD) associated with T µ ).Suppose a set S I consists of all p ∈ M satisfying the following condition: there exists a spacelike vector n µ orthogonal to T µ at p, such that all null geodesics emanating from p whose tangent vector k µ is orthogonal to n µ at p will not reach future null infinity I + .Then, we call S I the inner dark domain (IDD) associated with T µ .Definition 3.4 (inner dark horizon (IDH) associated with T µ ).Let a region S I be the IDD associated with T µ in M.Then, the boundary ∂S I is called the inner dark horizon (IDH) associated with T µ .
The outer dark domain and the inner dark domain are sometimes simply called the dark domains in short.Similarly, the outer dark horizon and the inner dark horizon are simply called the dark horizons.The ODD and IDD are complementary concepts to each other in the sense that the roles of "there exists" and "for all" are switched in the definitions.We required the vector n µ to be orthogonal to T µ in the definitions of 3.1 and 3.3 in order to consider the escape condition of photons in the frame associated with T µ .As a result, the locations of dark domains and dark horizons depend on T µ .We take T µ as tangent vectors of the light source worldlines.The escape condition of the emitted photon depends on the motion of the sources due to the effect of the relativistic beaming.Here, it is also possible to consider situations where T µ is defined only in a subregion of the spacetime and define the dark domain and dark horizon in this subregion, but we mainly focus on situations where T µ is defined globally in the whole spacetime.
In addition, we can consider a particular class of the outer/inner dark domain/horizon.Suppose that a time coordinate t is present, and a spacelike hypersurface Σ t is given by t = constant.From the time coordinate t, we can naturally introduce a timelike unit vector field by Adopting this vector field as T µ in the definitions of the ODD and IDD, we can introduce outer/inner dark domain/horizon associated with the time slice Σ t as follows:4 Definition 3.5 (outer/inner dark domain/horizon associated with Σ t ).Suppose a time coordinate t is present, and the spacelike hypersurface t =constant is denoted as Σ t .The ODD/IDD and the ODH/IDH associated with the timelike vector field T µ = −∇ µ t/ √ −∇ ν t∇ ν t are called the ODD/IDD and the ODH/IDH associated with time slice Σ t .
For the outer/inner dark domain/horizon associated with Σ t , the vector n µ in the definitions of 3.1 and 3.3 is tangent to Σ t .In Sec. 4, we will see that the ODH and IDH are generalizations of the photon sphere.We will also see that they are both absent in the Minkowski spacetime and both exist in a spacetime with black hole(s) under certain conditions.These facts imply that the existence of the ODH/IDH is expected to be an indicator for strong gravity.The words "outer" and "inner" refer to the fact that the IDD is generally included by the ODD, which will be confirmed in Cor.3.9.
Since the geodesic equations are second-order differential equations, the initial position and its first derivative provide us with the unique solution 5 .Therefore, at any point in the spacetimes, the emission angle of a photon (in the frame of the photon source) determines whether it reaches future null infinity or not.

Relation to escape and capture cones
The escape cone is a useful concept to describe the condition for photons to reach future null infinity.In order to explain this concept, we introduce the tetrad basis, T µ and (e i ) µ with i = 1, 2, and 3, and the projection tensor γ µν = g µν + T µ T ν onto the three-dimensional spacelike tangent subspace spanned by (e i ) µ .Then, the emission direction in the frame of the photon source is expressed by the tangent null vector k µ of a photon Since E µ is a unit vector, there is the apparent one-to-one correspondence between E µ and a point of a unit two-sphere in the tangent subspace, which we call the "twosphere of emission directions."For each point on the two-sphere, we can judge whether the corresponding photon escapes to future null infinity or not.As a result, the twosphere is divided into the escape region and the capture region.If we connect the central point of the two-sphere and the points on the boundary of these two regions with straight lines, the escape and capture cones appear.Below, we use the words the "escape cone" and the "capture cone" with the same meanings as the escape region and the capture region on the two-sphere of emission directions, respectively.Null geodesics that correspond to the boundary of the escape and capture cones coincide with the future wandering null geodesics defined by Siino [16,17].Since the wandering null geodesics never arrive at future null infinity, the boundary is regarded to belong to the capture cone.Therefore, the capture cone is a closed set, while the escape cone is an open set.
Fig. 1 Representation of the two-sphere of emission directions at r = 3M in the Schwarzschild spacetime, at which the photon sphere is located.The escape region and the capture region are shown by the grey region and the white region, respectively.The boundary of the two regions is given by an orthodrome.
In a static and spherically symmetric spacetime, the boundary of an escape region associated with a constant-time slice is given by an exact circle.In particular, on the photon sphere in the Schwarzschild spacetime, the boundary of the escape region is given by an orthodrome of the two-sphere (see Fig. 1).One might think that it would be reasonable to define a generalization of the photon sphere as a set of points at which the boundary of the escape region is an orthodrome associated with a constant-time slice, at least for stationary spacetimes.However, in general spacetimes, the boundary of the escape region is not necessarily an exact circle.This leads us to the physical interpretation for our new concepts as described below.
To characterize emission directions E µ orthogonal to a spacelike vector n µ in (negations of) Defs.3.1 and 3.3, the notion of orthodrome is useful.In the tangent subspace orthogonal to a unit timelike vector field T µ , an orthodrome is defined as an intersection between a plane containing the origin of the tangent subspace and the two-sphere of emission directions (see circles on the two-spheres in Fig. 2).Note that an orthodrome orthogonal to n µ represents a set of emission directions E µ orthogonal to n µ .
The conditions for the ODD and IDD associated with T µ are now rephrased using the escape and capture cones as follows.On the one hand, the negation of Def.3.1 is equivalent to that at the spacetime point p outside the ODD, there exists an orthodrome that is entirely included in the escape cone.On the other hand, Def.3.3 is equivalent to that at a spacetime point p in an IDD, there exists an orthodrome that is entirely included in the capture cone.Note that each of the escape and capture cones can have multiple components.For example, at a point between two distant black holes, there would be two capture cones and one "escape belt".In the case that each of the escape and capture cones has a single component, a hemisphere is included in the capture cone at a point in the IDD, and a hemisphere is included in the escape cone at a point outside the ODD.
Figure 2 illustrates the situation where each of the escape and capture cones has a single component.Spacetime points can be divided into three classes.In the first class (the left sphere in Fig. 2), there exists an orthodrome which is totally included by the capture cone, and the corresponding point in the spacetime is in the IDD.This point is also in the ODD because the escape cone is included in a hemisphere and it cannot include an orthodrome (a strict proof will be given in Cor.3.9).In the second class (the center sphere in Fig. 2), there is no orthodrome that is entirely included by the escape region or by the capture region.The corresponding points in the spacetime are in the ODD, but not in the IDD.In the last class (the right sphere in Fig. 2), there exists a hemisphere which is totally included in the escape region.In this case, the corresponding point is not in the ODD as mentioned above.In addition, since the capture cone is included in a hemisphere, this point is not in the IDD (see a strict proof for Cor.3.9).
In short, any point in the ODD has a small escape cone, and any point not in the ODD has a large escape cone.As for the IDD, any point in the IDD has a large capture cone, and any point not in the IDD has a small capture cone.There is no point such that it is in the ODD and not in the IDD, which will be proved more strictly in Cor.3.9.
Fig. 2 Depiction of two-spheres of emission directions, which represent the escape regions (the grey regions) and the capture regions (the white regions) on three points in a spacetime without symmetry.A vector orthogonal to the plane including the orthodrome corresponds to the spacelike vector n µ in (negations of) Defs.3.1 and 3.3.Left sphere: An illustration of the escape and capture regions at a certain point.In this case, this point is both in the ODD and in the IDD.Center sphere: A picture of the escape and capture regions at another point.In other words, all orthodromes intersect the boundary of the two regions.In this case, the point is in the ODD but not in the IDD.Right sphere: An illustration of the escape and capture regions at another point.In this case, the point is not in the ODD nor IDD.
We now present several useful propositions in order to judge whether a given spacetime point is in an ODD/IDD or not.For this purpose, we prove the following Lemma: Lemma 3.6.Consider a unit two-sphere and a pair of antipodal points P and Q on that sphere.Let γ be any continuous curve that connects P and Q, and γ includes the two endpoints P and Q.Then, any orthodrome on that sphere intersects γ.
Proof.In the case that the orthodrome includes the points P and Q, it trivially intersects γ.In the case that the orthodrome does not include P and Q, the orthodrome divides the two-sphere into two regions, and the points P and Q are included in different regions because they are antipodal points.Since the curve γ continuously connects P and Q, it must cross the boundary of the two regions.That is, the orthodrome intersects γ.
Suppose that a pair of antipodal points are included in a connected component of the capture cone at some spacetime point.Then, there exists a curve that connects the two points in the capture cone.From the above Lem.3.6, any orthodrome crosses that curve, and hence, it cannot be included entirely in the escape cone.This means that such a spacetime point is in an ODD.For a spacetime point at which each of the escape and capture cones is path-connected, we easily see that the reverse is also true and obtain the following: Theorem 3.7.If a pair of antipodal points are included in a connected component of the capture cone, the point in the spacetime is in an ODD.In particular, suppose that each of the escape and capture cones is path-connected.Then, such a spacetime point is in an ODD if and only if the capture cone contains a pair of antipodal points.
Similarly, when a pair of antipodal points are included in a connected component of the escape cone, there exists a curve that connects the two points in the escape cone.Since any orthodrome crosses that curve, it cannot be entirely included in the capture cone.This means that such a spacetime point is not in an IDD.For a point at which each of the escape and capture cones is path-connected, we see that the reverse also holds.Theorem 3.8.If a pair of antipodal points are included in a connected component of the escape cone, such a spacetime point is not in an IDD.In particular, suppose that each of the escape and capture cones is path-connected.Then, a point is not in an IDD if and only if the escape cone contains a pair of antipodal points.
At a spacetime point p in an IDD, an orthodrome is included in the capture cone.This immediately implies the existence of an infinite number of pairs of antipodal points in a connected component of the capture cone.Then, Thm.3.7 implies that the point p is also in an ODD, and hence, we have the following: Corollary 3.9.An IDD is always included in an ODD, as long as they are associated with the same T µ .
We now present the useful criteria to find an ODH and an IDH.Assuming the continuity of the mappings from each spacetime point to a configuration of escape and capture cones, we can find the following: Proposition 3.10.Suppose that each of the escape and capture cones has a single component at each spacetime point, and there are continuous mappings from each spacetime point to a configuration of escape and capture cones.Then, at the ODH, there is an orthodrome that is included in the escape cone except at two or more points, which are on the boundary of the escape and capture cones.Similarly, at the IDH, there is an orthodrome that is entirely included in the capture cone, and furthermore, the orthodrome intersects the boundary of the escape and capture cones at two or more points.
There are several remarks on this Proposition.First, the continuity of the escape and capture cones is guaranteed if the spacetime metric has C 2− continuity, 6 and the timelike vector field T µ is continuous.Here, other than the ODH or the IDH, there might be a surface having the same property for the ODH/IDH mentioned in Prop.3.10.However, the outermost surface among them can be identified as the ODH/IDH.

Properties of Outer Dark Domain and Outer Dark Horizon
In this subsection, we provide several properties on the ODD and ODH.First, since all null geodesics in the Minkowski spacetime reach future null infinity, we have the following proposition.Proposition 3.11.In the Minkowski spacetime, the ODD and ODH are both absent for any associated timelike vector field T µ .Proposition 3.11 means that the curved spacetime is essential for the existence of the ODH.
Next, we focus on the spacetimes with black hole(s).For any point p in a black hole, all null geodesics emanating from p do not reach future null infinity by the definition of the black hole.Therefore, p is an element of the ODD.This is summarized as the following theorem: Theorem 3.12.In an asymptotically flat spacetime with black hole(s), points in the black hole(s) are also in the ODD for every associated timelike vector field.
Let us discuss the existence of the ODH in spacetimes with black hole(s) for a particular class of T µ .To see this, we first show the following lemma.Lemma 3.13.Consider a four-dimensional asymptotically flat spacetime in which the metric near future null infinity behaves as Eqs.(2)-( 5) with C 2− functions using Bondi coordinates.Let a timelike vector T µ satisfy T µ ∝ (du + dr) µ .We assume that Ω i , which is defined in Eq. ( 15), is positive.Then, all points p with a sufficiently large r are not in the ODD associated with T µ .
Proof.We will show that, for all point p with sufficiently large r, there exists a tangent vector n µ orthogonal to T µ , such that all null geodesics which are orthogonal to n µ at p reach future null infinity.Taking µ to be an orthogonal vector to T µ , we will show that any null vector orthogonal to n µ 1 satisfies k r = O r −1 .By using Eqs.( 2), (3), and (5), the orthogonal condition between n µ 1 and k µ , that is, n µ 1 k µ = 0, gives and we have Then, a simple calculation gives us and we see With Eq. ( 21), this shows us As a consequence, we have k r = O r −1 = O r −1 k u .By Prop.2.1, null geodesics with k r = O r −1 k u at sufficiently large r region reach future null infinity if Ω i > 0 holds.Thus, all null geodesics which are orthogonal to n µ 1 at p reach future null infinity.Therefore, p is not in the ODD associated with T µ .
Here, Lem.3.13 does not depend on whether there exist black hole(s) or not.We can also show a higher-dimensional version of Lem.3.13 in a similar way.Note that in higher dimensions, Ω i > 0 is not necessary to show the null geodesics to reach future null infinity [24,26].Since the same fact also holds for all other theorems, we will not mention higher-dimensional cases, hereinafter.Now, let us show the existence of the ODH in spacetimes with black hole(s) under the same conditions with Lem.3.13.Proposition 3.14.Consider an asymptotically flat spacetime with black hole(s) under the assumptions in Lem.3.13.The ODH associated with T µ is not empty, and not connected to future null infinity.
Proof.Let the region of the ODD be denoted by S O .Then, S O is not empty by Thm.3.12.Due to Lem. 3.13, a region sufficiently close to future null infinity is a subset of M\S O .Therefore, the ODH is not empty by Def.3.2, and not connected to future null infinity.This proposition is supporting evidence for the ODH to serve as an indicator for a strong gravity region.

Properties of Inner Dark Domain and Inner Dark Horizon
In this subsection, we provide several properties on the IDD and IDH.In the same way to the proofs for Prop.3.11 and Thm.3.12, we can easily show the following statements.Proposition 3.15.In the Minkowski spacetime, the IDD and IDH are both absent for any associated timelike vector field T µ .Theorem 3. 16.In an asymptotically flat spacetime with black hole(s), points in the black hole are elements of the IDD for any associated timelike vector field T µ .Proposition 3.15 means that curved spacetimes are essential for the existence of the IDD and IDH.Meanwhile, Thm.3.16 guarantees the existence of the IDD in spacetimes with black hole(s).
To show the existence of the IDH in spacetimes with black hole(s), the following lemma is essential, which can easily be verified by using Cor.3.9 and Lem.3.13.Lemma 3.17.Under the assumptions in Lem.3.13, all points p with sufficiently large r are not in the IDD associated with T µ .
One may expect that the IDH would serve as an indicator for the existence of a strong gravity region.This is indeed correct as the following proposition.Proposition 3.18.Consider an asymptotically flat spacetime with black hole(s) under the assumptions in Lem.3.13.The IDH associated with T µ is not empty, and not connected to future null infinity.
Proof.The IDD S I is not empty by Thm.3.16.In addition, a region sufficiently close to future null infinity is a subset of M\S I due to Lem. 3.17.Therefore, the IDH is not empty, and not connected to future null infinity by Def.3.4.
Finally, let us prove that the ODH coincides with the IDH in a spherically symmetric spacetime for the spherically symmetric unit timelike vector field T µ : Theorem 3.19.Consider an asymptotically flat, spherically symmetric spacetime.Let x I be the angular coordinates.We adopt the timelike vector field T µ to be the spherically symmetric one that satisfies T µ (dx I ) µ = 0.Then, the ODH associated with T µ coincides with the IDH associated with T µ .
Proof.If the ODD and IDD coincide with each other, the ODH and IDH also coincide with each other.Here, the IDD is a subset of the ODD by Cor.3.9.Therefore, the task here is to show that the ODD is a subset of the IDD.
Let a point p be in the ODD.By the definition of the ODD, for all spacelike vectors n µ orthogonal to T µ , there exists a null geodesic which emanates from p, whose tangent vector k µ is orthogonal to n µ at p, such that it will not reach future null infinity.We now adopt n µ as the one satisfying n I = 0. Due to the spherical symmetry of the spacetime, all null geodesics with the tangent vector k µ satisfying k µ n µ = 0 emanating from p will not reach future null infinity.Then, p satisfies the condition of being in the IDD.

Explicit example I: static and spherically symmetric spacetimes
To understand properties of the dark horizon in a simple example as a first step, let us focus on a static and spherically symmetric spacetime with the metric; with f (r) = 1+O r −1 and g(r) = 1+O r −1 .To clarify the dependence of associating timelike unit vector T µ on the position of the dark horizon, we divide our analyses into two parts.In subsection 4.1, we assume T µ to be proportional to timelike Killing vector, and explain the similarity between the dark horizon and the photon sphere.In subsection 4.2, we generalize T µ to include the radial component, and see its dependence on the dark horizon.Subsections 4.1 and 4.2 assume static and dynamical light source flow setups, respectively.

Static Flow
In this subsection, we take an associating timelike unit vector field T µ to be T µ = (∂ t ) µ / f (r) in Defs.3.1-3.4,and consider the ODH and IDH associated with t−constant hypersurfaces.This setup corresponds to a simple model assuming that light sources are static, i.e., fixed at their own constant radii.For simplicity, we omit writing "associated with T µ " and "associated with t−constant hypersurfaces" here.
Since the ODD and IDD coincide due to Thm. 3.19, we will simply refer to the ODD and IDD as the dark domain, and the ODH and IDH as the dark horizon, respectively.Along with the above T µ , a certain point is in the dark domain if and only if null geodesics passing that point with zero radial velocity do not reach future null infinity.
Let us review a condition for a null geodesic to reach future null infinity in a static and spherically symmetric spacetime [29][30][31].Due to these symmetries, the impact parameter b := −k ϕ /k t is conserved along the geodesic.Without loss of generality, we can restrict our considerations on the equatorial plane with b ≥ 0.Then, the geodesic equation for r is where the prime is the derivative with respect to the affine parameter.Then, due to the non-negativity of the right-hand side of Eq. ( 26), b is restricted to |b| ≤ B(r) := r/ f (r).Taking the derivative of Eq. ( 26) gives us which tells us that the signs of r ′′ | r ′ =0 and dB/dr coincide.In particular, dB/dr = 0 is a necessary and sufficient condition for r ′′ | r ′ =0 = 0, which determines the locus of the photon sphere.Since b is conserved along a null geodesic, the behavior of a null geodesic is determined by a b−constant line, constrained with the potential barrier of B(r) on the (r, b)-plane (see Fig. 3), similarly to the fundamental discussion on the effective potential in the classical mechanics.For a function B(r) drawn in Fig. 3, for instance, the fate of a null geodesic emitted with r ′ = 0 (i.e., starting from a point on b = B(r) in the figure) is determined as follows.If starting from r II > r or r min < r < r I with r ′ = 0, the radial coordinate of the null geodesic increases to infinity because there is no potential barrier outside of the initial point.If starting from the other region, the radial coordinate on the null geodesic has an upper bound, and hence that null geodesic does not reach future null infinity.Then, the dark horizon is located at the boundary of these regions, that is r = r min , r I and r II .
It is easy to see from Fig. 3 that a "visible" photon sphere is a dark horizon.Here, we call a photon sphere "visible" if null geodesics passing slightly outside of the photon sphere with r ′ = 0 reach future null infinity.In particular, the outermost photon sphere, which is visible due to the asymptotic fall-off condition f (r) = 1 + O r −1 and g(r) = 1 + O r −1 , coincides with the outermost dark horizon.Indeed, the visible photon sphere captures the essential properties of the photon sphere-the properties such that it represents the boundary of the observable region, satisfies the areal inequality Eq. ( 1) for the outermost visible photon sphere [10], and serves as an indicator for a strong gravity region.
Note that there are photon spheres that are not dark horizons.For example, stable photon spheres, at which B(r) is maximal, are not dark horizons.It should be also noted that there exist dark horizons that are not photon spheres.For example, in a spacetime with B(r) of Fig. 3, there are two dark domains.The inner boundary of the outer component of the dark domain is a dark horizon, but not a photon sphere.Let us discuss the physical meaning of this dark horizon.As shown in Fig. 3, we express the radial coordinate of this dark horizon as r I , and that of the outermost dark horizon as r II .Then, we consider the observation of photons emitted from this region.For a distant observer, the image of the observed region can be regarded as a series of concentric circles, and each circle is specified by the angle ϑ from the center of the image.Here, the angle ϑ is a monotonically increasing function of the impact parameter b, and ϑ = 0 corresponds to b = 0.When one regards the minimum radial coordinate for each null geodesic as a function of b, it has a discontinuity at b c := B(r I ) = B(r II ), and hence, the observable region changes discontinuously at a certain angle ϑ c corresponding to the impact parameter b c .Namely, we observe the region r > r II for ϑ = ϑ c + ϵ for small ϵ, while we catch signals from the region r ≥ r I for ϑ = ϑ c − ϵ.Due to the above discontinuity of the observable region, the brightness would change discontinuously at ϑ = ϑ c .

Dynamical Flow
In this subsection, let us generalize T µ to have the radial component, and study the dependence of T µ on the position of the dark horizon.This situation assumes a distribution of light sources moving in the radial direction.
We parameterize the direction of T µ with the light source's outward velocity β (|β| < 1) as If β is positive, the light source is moving in the outward direction.We do not require β to be constant but impose continuity in the spacetime.Due to the spherical symmetry of the spacetime, we set n µ in Def.3.3 to be Consider a null vector k µ with k t = −1 which is orthogonal to n µ .Without loss of generality, we assume k θ = 0 and consider a point on the equatorial plane in the spacetime.The ϕ component of k µ of corresponding to this null vector is Then, the condition for a spacetime point to be on the dark horizon is that this null vector is on the boundary of the escape cone at that point.As a concrete example, let us consider Schwarzschild spacetime, f = 1/g = 1 − 2M/r.Then Eq. ( 30) reduces to Here, we recall the fact that the impact parameter b is given by the angular momentum k ϕ divided by the energy −k t , that is, b = −k ϕ /k t = k ϕ .Since the absolute value of the impact parameter corresponding to the boundary of the escape cone is 3 √ 3M , we consider the condition which gives us the radius of the dark horizon associated with T µ .Figure 4 represents the depedence of β on the radial coordinate of the dark horizon r DH .In the case of negative β, the light source accretes onto the black hole and the dark horizon tends to be larger.Conversely, in the case of positive β, the light source explodes to distant region, and the dark horizon tends to be smaller.This corresponds to a situation where photons emitted isotropically in the comoving frame are beamed toward the direction of the source motion.In other words, the definition of the dark horizon generalizes the photon sphere to include the relativistic beaming effect.
Fig. 4 A plot showing the relation between the radial coordinate of the dark horizon (r DH ) and the value of β at that point.The unit of the vertical axis is M .The dark horizon is located at r = 3M if β = 0 on r = 3M .A point sufficiently close to the horizon can be on the dark horizon if a light source at that point is moving outward nearly at the speed of light (β ≃ 1).In addition, a point with a sufficiently large radial coordinate is on the dark horizon if a light source is moving inward nearly at the speed of light (β ≃ −1).

Explicit example II: Vaidya spacetime
In this section, as an example of dynamical spacetimes, we investigate the dark horizon in the Vaidya spacetime.The Vaidya spacetime has two classes depending on the directions of the flow-ingoing and outgoing Vaidya spacetimes.Their metrics are written as and respectively, where v and u are the advanced and retarded times, respectively.We focus our attention to the domain where r > 2M (v) and r > 2M (u) hold, and assume that the timelike vector field T µ is proportional to (∂ v ) µ and (∂ u ) µ in the ingoing and outgoing Vaidya spacetimes, respectively.Here, we have assumed that the light sources are different from the null fluid matter, and the worldlines of the light sources correspond to (r, θ, ϕ)-constant line.Due to Thm. 3.19, the ODH and IDH associated with the same unit time vector field T µ coincide with each other.In what follows, we just call the ODH and IDH the dark horizon, and do not write "associated with T µ " in this section.In Refs.[32,33], the shadow edge worldlines in the Vaidya spacetime were discussed for an idealized situation in which the observer and the light source are supposed to be located at future null infinity and past null infinity, respectively.In this sense, the authors of Refs.[32,33] defined the dynamical photon sphere as the boundary of the set of null geodesics connecting future and past null infinity, and explicitly calculated the time-dependence of the dynamical photon sphere.To compare the dark horizon with the dynamical photon sphere, let us examine the dark horizon in the same setup.
In the case with the ingoing Vaidya spacetime, we use the same mass function as that in Ref. [32]: with v f = 100.In order to specify the position of the dark horizon at each v, we numerically solved the null geodesic equation with the initial condition that corresponds to r = r 0 and k r = 0, where k µ is a tangent vector of the null geodesic.By changing the initial values of r 0 , we specified the radius of the dark horizon r DH such that the null geodesic enters the black hole region for r 0 < r DH and escapes to infinity for r DH < r 0 .
The dynamical photon sphere is also obtained numerically in the same way as Ref. [32].
Figure 5 shows the radii of the dark horizon, the dynamical photon sphere and the graph r = 3M (v).To see the difference between the dark horizon and the dynamical photon sphere, consider two sets of null geodesics-ones corresponding to the dark horizon and ones corresponding to the dynamical photon sphere.The null geodesics emitted from the dark horizon initially have k r = 0.By contrast, from the behavior of the dynamical photon sphere depicted in Fig. 5, the null geodesics that determine the locus of the dynamical photon sphere has k r > 0 because the null geodesics are included in the dynamical photon sphere, whose area increases with time.This is the reason why the dark horizon is located outside the dynamical photon sphere, shown in Fig. 5.This is a consequence of the assumption that the light sources are located at past null infinity in this case.Here, both sets of null geodesics describe the boundary of null geodesics which reach future null infinity or not, but the locations of the light sources are different.
Similarly, in the case of the outgoing Vaidya spacetime, we use the mass function as follows: 5 Comparison between the dark horizon (gray dots) and the dynamical photon sphere (red curve) in the ingoing Vaidya spacetime with the metric Eqs. ( 33) and (35).The blue dashed curve denotes r = 3M (v).
with u f = 100.Figure 6 shows the dark horizon, the dynamical photon sphere and the graph r = 3M (u).One can see that the dark horizon is located inside of the dynamical photon sphere.The reason is interpreted as follows.The null geodesics that correspond to the edge of the escape cone at the dark horizon has the vanishing radial component of the tangent vector, k r = 0.By contrast, the plot of the dynamical photon sphere in Fig. 6 tells us that the null geodesics corresponding to the dynamical photon sphere direct inward because the null geodesics are included in the dynamical photon sphere, whose area decreases with time.This gives us the reason why the dark horizon is located inside of the dynamical photon sphere, shown in Fig. 6.This is a consequence of the assumption that the light sources are located at past null infinity in this case.Here, both sets of null geodesics describe the boundary of null geodesics which reach future null infinity or not, but the locations of the light sources are different.

Explicit example III: Kerr spacetime
In this section, we study the dark horizons in the Kerr spacetime as a simple model of typical black holes in our universe, assuming for simplicity that the motion of light sources is the same as that of the zero-angular-momentum observers (ZAMOs).To be specific, we consider the ODH and IDH associated with a timeslice Σ t of the Boyer-Lindquist coordinates (t, r, θ, ϕ).This section is organized as follows.In subsection 6.1, we show basic properties of null geodesics in the Kerr spacetime.In subsection 6.2, the ODH and IDH on the rotation axis is considered.In subsections 6.3 and 6.4, we investigate the ODH and IDH in the Kerr spacetime, respectively.

Preparation
In this subsection, we conduct several calculations as a preparation for the depiction of the shapes of the ODH and IDH in the Kerr spacetime.The metric of the Kerr Fig. 6 Comparison between the dark horizon (gray dots) and the dynamical photon sphere (red curve) in the outgoing Vaidya spacetime with the metric Eqs. ( 34) and ( 36).The blue dashed curve corresponds to r = 3M (u).
spacetime is given by where Σ, ∆, and A are defined as Here, M is the ADM mass and a is the rotation parameter that is related to the ADM angular momentum as J = M a.We also sometimes use the dimensionless rotation parameter a * = a/M to specify the solution.Below, the black hole spacetime with a nondegenerate event horizon, 0 < |a * | < 1, is considered (see Sec. 4 for the case of a * = 0).The event and Cauchy horizons are located at r = r + H and r = r − H , respectively, where and we shall limit our attention to the region outside the event horizon, r > r + H . Let us introduce the tetrad in the Kerr spacetime as follows: (e 3 ) µ = 0, 0, 0, Σ A Here, (e 0 ) µ is orthogonal to the t-constant slice Σ t of the Boyer-Lindquist coordinates.
The ODD, ODH, IDD, and IDH associated with Σ t in the Def.Here, E indicates the emission direction vector given by Eq. ( 19) [the spatial part of Eq. ( 44)].
At each point p of the Kerr spacetime, we consider the (initial) null tangent vector k µ of a null geodesic that corresponds to a photon emitted from the point p.Adopting the normalization k µ (e 0 ) µ = −1 , the initial null tangent vector k µ is expressed in the form of where 0 ≤ ξ ≤ π and −π ≤ η < π.See the illustration for the definition of ξ and η in the right panel of Fig. 7.Then, at each point p, the covariant components k t and k ϕ are written down as where the values of r and θ at the point p must be used in these equations.After the emission, the motion of the photon follows the null geodesic.In a Kerr spacetime, each null geodesic possesses the conserved quantities, the energy E, the angular momentum L, and Carter's constant Q [30] (we follow the notation of Ref. [29] for Carter's constant).Their values are determined by the initial condition (i.e., the information at the point p), and for the energy and angular momentum, we have E = −k t | p and L = k ϕ | p , respectively.The impact parameter b is defined by b := L/E, and its value becomes In addition, we introduce the dimensionless Carter's constant q = Q/E 2 , and the values of q and b satisfy Eliminating η from Eqs. (47) and (48), one has where From Eq. (49), it is found that the value of F (b, q) must be nonnegative, that is, the impact parameter b and the dimensionless Carter's constant q can only take values such that F (b, q) ≥ 0 is satisfied.Taking the square root of Eq. ( 49) and substituting it into Eq.(47), we have Let us discuss the relation between ξ and η on the boundary of the escape/capture cones.It is known that the null geodesics corresponding to the boundary of the escape cone neither reach future null infinity nor fall into the black hole.These null geodesics asymptote to the spherical orbits around the black hole with specific values of b and q, which are called the "spherical photon orbits" [35].To discuss the null geodesics which correspond to the boundary of the escape cone, let us look at the relation between b and q for the spherical photon orbits.
Along the null geodesic, the following first-order differential equation holds [30]: where A spherical photon orbit with the radius r s satisfies R| r=rs = 0 and dR/dr| r=rs = 0, and these conditions are rewritten as Equations ( 54) and (55) give the relation between b and q via a parameter r s that is satisfied by various spherical photon orbits.Next, consider null geodesics asymptoting to spherical photon orbits in the future (t → ∞).These null geodesics correspond to the boundary of the escape cone at a point on these null geodesics.Since b and q of these null geodesics are the same as those of corresponding spherical photon orbits, the same relation of Eqs.(54) and (55) holds.By substituting Eqs.(54) and (55) into Eqs.(51a) and (51b), we obtain the relation between sin ξ and cos η as parametrized by r s , and this relation gives the boundary of the escape cone.However, since two values of ξ give the same sin ξ in the range 0 ≤ ξ ≤ π, the obtained relation between ξ and η includes both the emission directions of photons that asymptote to the spherical photon orbits in the future (t → ∞) and in the past (t → −∞).Therefore, we must select an appropriate boundary to specify the real escape cone.This can be done by several criteria: (i) The structures of the escape cones for a Schwarzschild spacetime are well known, and the escape cones change smoothly as the value of a * is increased; (ii) There exist continuous mappings from a spacetime point to a configuration of escape and capture cones; (iii) The boundary of the escape/capture cones must be smooth on the two-sphere of emission directions.
We have found that the formula of cos ξ that satisfies the above requirements is given by (56) In fact, the function in the square root of Eq. ( 56) has the factor (r s − r) 2 , and hence, Eq. ( 56) gives an analytic formula with respect to r s for any fixed r and θ.The explicit formula is shown as Eq.(A5b) with Eq. (A6) in App. A. Here, we present examples of the boundaries of the escape and capture cones.For this purpose, we consider the projection of the two-sphere of emission directions onto the plane spanned by e 1 and e 3 .In other words, we consider a "side view" of the two-sphere as seen from a distant position with ξ = η = π/2 in the tangent subspace.Then, the two-sphere is projected to a unit disk, and the boundary of the escape and capture cones becomes a curve that has the endpoints on the unit circle due to the symmetry in the transformation η → −η. Figure 8 presents the projected boundary of the escape and capture cones for various r values on the equatorial plane θ = π/2 for a * = 0.99.The upper region of each curve is the escape cone, while the lower region is the capture cone.Due to the effect of dragging into rotation of the Kerr black hole, the escape cone tends to be directed in the e 3 direction.
We point out important properties of the curve of the projected boundary of the escape and capture cones.Let us introduce the Cartesian coordinates (x, y) on the two-dimensional plane defined by ∂ x = e 3 and ∂ y = e 1 .The origin of the coordinates is set so that the projected two-sphere of emission direction is given by x 2 + y 2 ≤ 1.Note that (ξ, η) = (π/2, ±π/2) are both projected to the origin.The curve of the projected boundary can be expressed as y = y(x).Then, as proved in App.A, it is possible to show the inequalities dy dx ≤ 0, (57) and where equalities hold on and only on the rotation axis.Denoting the two endpoints of the projected boundary of the escape and capture cones as (ξ, η) = (ξ + , 0) and (ξ − , π), the first inequality of Eq. ( 57) indicates that ξ + ≥ ξ − always holds for a * > 0.
The second inequality of Eq. ( 58) means that the curve of the projected boundary is convex upward (i.e.concave downward).This property will be used in figuring out the ODH and IDH.For a later convenience, let us examine the behavior of a null geodesic whose tangent vector k µ of Eq. ( 44) has the parameter values (ξ, η) = (π/2, ±π/2), corresponding to (x, y) = (0, 0), at the emission point p of the corresponding photon.The motion of such a photon is initially directed toward ±(e 2 ) µ in the ZAMO frame and tangent to the r-constant surface, r ′ = 0.The photon possesses zero impact parameter, b = 0, by Eq. (47).Due to Eq. ( 52), r ′ = 0 is equivalent to R = 0, and the value of dR/dr under the condition R = 0 is computed as We can easily see that dR/dr| R=0 = 0 has three real solutions, and only one solution outside the event horizon, r > r + H , is given by This gives the radius of the spherical photon orbit, and the trajectory of a photon in this orbit crosses the rotation axis due to the property b = 0.It is easy to see that dR/dr| R=0 < 0 holds in the range r + H ≤ r < r (0) s , and the emitted photon falls into the black hole.We also see dR/dr| R=0 > 0 in the range r > r (0) s , which means that r ′′ > 0 initially.For such situations, R is positive in the region r > r (0) s , and hence, the emitted photon escapes to future null infinity.
We would like to study one more case: the case sin ξ = 0, corresponding to (x, y) = (0, ±1), at the emission point p.In this case, the photon is emitted in the ±(e 1 ) µ direction in the ZAMO frame.Equations ( 47) and ( 48) indicate b = 0 and q = −a 2 cos 2 θ p .Then, the function R of Eq. ( 52) is positive and there is no turning point for the radial motion.Therefore, the photon emitted in the +(e 1 ) µ direction escapes to infinity, while the photon emitted in the −(e 1 ) µ direction falls into the black hole.This means that the points ξ = 0 and ξ = π on the sphere of emission directions belong to the escape cone and the capture cone, respectively, at an arbitrary spacetime point outside of the horizon.

Dark horizons on the rotation axis
Here, we consider the ODH and IDH on the rotation axis.Strictly speaking, since the tetrad frame of Eqs. ( 40)-( 43) becomes singular at the rotation axis, θ = 0 and π, the above analysis must be reconsidered if the initial emission point is on the rotation axis.However, this case can be handled with a minor modification to the case θ ̸ = 0 or π.Let us focus on the upper axis, θ = 0, for simplicity.The most convenient way would be to consider the limit θ → 0 of the tetrad frame along ϕ = 0 with t, r = constant.In this case, (e 1 ) µ is the radial unit vector, and (e 2 ) µ and (e 3 ) µ are the unit vector directed to +θ directions along ϕ = 0 and ϕ = π/2, respectively.Then, the results of Sec.6.1 hold also on the rotation axis.
We can also obtain the radial coordinate of the ODH and IDH, analytically.Because of the axial symmetry, the escape cone at a point on the rotation axis also becomes axially symmetric, and hence, the boundary of the escape cone is given by ξ = const.Because of this property, the ODH and IDH become degenerate, similarly to the spherically symmetric case.At the ODH and IDH, the boundary of the escape cone becomes ξ = π/2.From the results of Sec.6.1, ξ = π/2 belongs to the capture cone for r + H < r < r s , belongs to the escape cone for r > r

Outer Dark Horizon in Kerr spacetime
In this subsection, we depict the ODH of the Kerr black hole by applying the results in subsection 6.1.In order to find the ODD, we can use Thm.3.7 which tells us that if a pair of antipodal points are included in the capture cone, that position is in the ODD.As a pair of antipodal points, we focus on the points (ξ, η) = (π/2, ±π/2), that correspond to the emission directions ±(e 2 ) µ .The motion of the corresponding photons are discussed in Sec.6.1:In the regions r + H < r < r s , and r > r s , they fall into the black hole, propagate along the circular orbits eternally, and escape to infinity, respectively.This means that the region r + H < r ≤ r (0) s is the ODD.We now examine whether r = r (0) s is the ODH or not by applying Prop.3.10.There exists an orthodrome that is tangent to the boundary of the escape and capture cones at (ξ, η) = (π/2, ±π/2) for the following reason.Consider the tangent subspace in which the two-sphere of emission directions exists.Since the boundary of the cones has the symmetry in the transformation η → −η, the tangent lines of the boundary at (ξ, η) = (π/2, ±π/2) are parallel to each other.Then, there exists a plane on which the two tangent lines exist, and the intersection between the plane and the two-sphere of emission directions becomes the desired orthodrome.If that orthodrome is included in the escape cone except at (ξ, η) = (π/2, ±π/2), the surface r = r (0) s is confirmed to be the ODH. Figure 9 shows an example of the projected boundary of the escape and capture cones and the projected orthodrome that is tangent to the boundary at (ξ, η) = (π/2, ±π/2) on the two-dimensional plane spanned by e 3 and e 1 for r = r (0) s on the equatorial plane in the case a * = 0.99.The projected boundary of the cones becomes a curve that passes through the origin, and the projected orthodrome becomes a tangent line to the curve of the projected boundary at the origin.The projected orthodrome divides the unit disk into the upper and lower half-disks, and then, the condition that the orthodrome is included in the escape cone except for (ξ, η) = (π/2, ±π/2) is equivalent to that the projected boundary of the cones is included in the lower halfdisk.This is guaranteed by the fact that the projected boundary is convex upward, as indicated by the inequality of Eq. ( 58).Therefore, we can safely declare that r = r (0) s is the ODH.
We depicted the ODH in the Kerr spacetime, i.e., r = r (0) s in Fig. 10 for the cases a * = 0.0, 0.5, 0.9, and 0.999.In Fig. 10, the IDHs, which will be discussed in the next subsection, are shown by blue closed curves.
At the distant place r ≫ M , the value of ξ + + ξ − is close to 2π.As r is decreased, the value of ξ + + ξ − becomes smaller.Suppose that a point p satisfying ξ + + ξ − = π is found.At that point, we can introduce the orthodrome that is tangent to the boundary of the escape cone at (ξ, η) = (ξ + , 0) and (ξ − , π) since these two points are antipodal to each other and the tangent vectors of the boundary of cones at these two points both have vanishing ξ components due to the symmetries in the transformation η → −η and in the transformation η − π → −(η − π).If that orthodrome is included in the capture cone, the point p is on the IDH due to Prop.3.10.Therefore, our strategy is firstly to solve for the point where ξ + + ξ − = π, and secondly to confirm that the obtained point is actually the point of the IDH.
Let us start from solving for a position that satisfies ξ + +ξ − = π.We focus on each θ-constant surface in the Kerr spacetime.For each r, the formulas to parametrically specify the boundary of the escape cone are given by Eqs.(51a) and (51b) with Eqs. ( 54) and (55).These formulas are symbolically written as sin ξ = sin ξ(r, r s ) and cos η = cos η(r, r s ).Writing the values of r s that satisfy η = 0 and π (i.e., cos η = ±1) as r ± s , we have ±1 = cos η(r, r ± s ).Then, ξ ± are determined by sin ξ ± = sin ξ(r, r ± s ).Since ξ + + ξ − = π leads to the condition sin ξ + = sin ξ − , we obtain the set of equations to determine the radius r IDH of the IDH for a given θ: This set of equations can be solved numerically easily.Except for θ = 0 and π, the value of r IDH is smaller than r  Let us move on to the next step at which we check the obtained surface, r = r IDH (θ), is actually the IDH. Figure 11 shows an example of the boundary of the escape and capture cones and the orthodrome tangent to it that are projected onto the twodimensional plane spanned by e 3 and e 1 for r = r IDH on the equatorial plane in the case of a * = 0.99.On the two-dimensional plane, the curve of the projected boundary and the line of the projected orthodrome have the same endpoints on the unit circle.By the projected orthodrome, the projected two-sphere of emission directions is divided into the upper and lower half-disks, and the condition that the orthodrome is included in the capture cone is equivalent to that the projected boundary of the cones is included in the upper half-disk.This is guaranteed by the fact that the curve of the projected boundary of the cones is convex upward, as indicated by the inequality of Eq. (58).Therefore, the obtained numerical solution r = r IDH (θ) is confirmed to be the IDH.
In Fig. 10, we plot the IDH with blue closed curves in the Kerr spacetime with several Kerr parameters.The ODH and IDH both exist in the Kerr spacetime due to Props.3.14 and 3.18.In the case of a * = 0, the ODH and IDH coincide with each other by Thm.3.19, and their radii are given by 3M .In the case of a * > 0, the IDH coincides with the ODH only on the rotation axis, and at other positions, the ODH is located in the outside region of the IDH due to Cor. 3.9.In each case, the ODH is given by r = r (0) s , which coincides with the edge of the shadow observed from infinity on the axis.

Conclusion & Discussion
In this paper, we have proposed new concepts dark horizons (the outer dark horizon, ODH for short, and the inner dark horizon, IDH for short) as generalizations of the photon sphere, from static and spherically symmetric spacetimes to general asymptotically flat spacetimes.In short, the dark horizons are defined with some timelike vector field T µ or the timeslice Σ t which can be chosen along the motion of light sources, and the dark horizons are defined as the boundary of the region where escape (capture) cone is large (small) enough in the frame associated with T µ or Σ t .This definition captures the escape probability of photons emitted from the light sources, and thus, the dependence of the position and the motion of light sources on the brightness observed from infinity.We have focused on four-dimensional asymptotically flat spacetimes.We can also apply the same definitions in higher dimensions as the dark horizons.In higher-dimensional case, the positive condition for Ω i , which is defined in Eq. (15), is not needed in Lems.3.13, 3.17, Props.3.14, 3.18.
We have strictly shown several statements for the dark horizons.First, the ODD and IDD are both absent in the Minkowski spacetime (see Props.3.11 and 3.15) and are both present in the spacetime with black hole(s) under a natural assumption (see Thms. 3.12 and 3.16).Second, two dark horizons both exist in black hole spacetime under a natural condition (see Props.3.14 and 3.18).Third, the IDH is typically located inside of or coinsides with the ODH (see Cor. 3.9).Fourth, in spherical symmetric spacetimes, the ODH and IDH coincide with each other.In addition, we have clarified explicit shapes of the dark horizons in the Vaidya and Kerr spacetimes with simply chosen T µ .However, we left plots with T µ determined by more reasonable motions of light sources for future work.
Through analyses in static and spherically symmetric spacetimes with T µ proportional to the timelike Killing field, we have seen that any photon sphere is a dark horizon if null geodesics passing slightly outside of that photon sphere with r ′ = 0 reach future null infinity.In this sense, the dark horizon is a generalization of a "visible" photon sphere.In addition, if we adopt T µ with nonzero radial component that represents the light source moving in radial directions, the dark horizon is located at a different position from the photon sphere.In this sense, the dark horizon is a further generalization of the visible photon spheres to include the relativistic beaming effect.
The dark horizons are defined by referring to null geodesics reaching future null infinity, at which ideal observers are supposed to be distributed.In this sense, the location of dark horizons are affected not only by the geometry around the dark horizons but also by that in the distant region.This distinguishes the concepts of the dark horizon from some of other generalizations of the photon sphere defined in terms of the local geometrical quantities, such as the extrinsic curvature (e.g.[12,15]).
The dark horizons rather resemble the concept of the black room of Ref. [17] in the sense that these are defined in terms of the global behavior of the null geodesics.Here, the black room is defined as a spacetime region such that any photon entering the black room never goes from it.According to this definition, photons emitted at a point on the boundary of the black room in the tangential directions to the boundary must either propagate on the boundary of the black room or fall into the black room.Then, if we prepare a timelike vector field T µ which is tangent to the boundary of the maximal black room, the IDD associated with it must include the boundary of the black room.
The main difference between the black room and the dark horizons is as follows.On the one hand, the condition for a specific point to be in the black room refers to the global behavior of the null geodesics both in the future and past directions of this point, and therefore, it cannot have past endpoints.Then, the black room does not exist in spacetimes in which a black hole exists only after the gravitational collapse.In addition, the black room does not exist for spacetimes without a black hole [17].On the other hand, each point of the dark horizons refers to the global behavior of null geodesics only in the future light cone of this point.As a result, it is possible for dark horizons to exist in black hole spacetimes with the gravitational collapse.In addition, it is also possible that the dark horizons exist in spacetimes without a black hole and they coincide with the visible photon sphere.These differences would originate from the difference of the supposed positions of the light sources in these concepts.The black room concept would correspond to the situation where light sources are located at past null infinity (as the simple theoretical models of the black hole shadows), whereas the dark horizons correspond to the situation where light sources are located in and around the ODH and IDH, which are typically near the black hole.We expect that the dark horizons serve as important concepts to describe the properties of the shadows whose light sources distribute in strong gravity regions in varieties of asymptotically flat spacetimes.

Fig. 3
Fig. 3 An example of a spacetime having five photon spheres.The shaded region b > B(r) is the forbidden region for a null geodesic.An orbit of null geodesic, described as a b−constant line, is reflected if reaching b = B(r) (see red arrows).The intervals [r H , r min ] and [r I , r II ] correspond to the dark domain, which are indicated by blue lines.Here, r H , r I and r II denote the horizon radius, the second outermost dark horizon and the outermost dark horizon, respectively.This spacetime has three dark horizons.

Fig. 7
Fig. 7 The spatial vectors of the tetrad basis and the emission direction vector.Left panel: The configuration of e 1 , e 2 , and e 3 at some spatial point on a t-constant hypersurface, where e 1 , e 2 , and e 3 correspond to Eqs. (41), (42), and (43), respectively.The black sphere is the event horizon, and the larger sphere indicates an r-constant surface.At each point, e 1 , e 2 , and e 3 are parallel to ∂r, ∂ θ , and ∂ ϕ , respectively.Right panel: The definition of ξ and η to parametrize the emission direction given in Eq. (44).Here, E indicates the emission direction vector given by Eq. (19) [the spatial part of Eq. (44)].

Fig. 8
Fig.8The boundaries of the escape and capture cones projected onto the plane spanned by e 3 and e 1 for r/M = 1.142, 1.145, 1.15, 1.16, 1.2, 1.3, 1.7, 2.3, 3.0, 4.0, 5.0, 6.5, and 10.0 in the equatorial plane θ = π/2 for a Kerr black hole with a * = 0.99.The upper region of each curve is the escape cone, while the lower region is the capture cone.

s
, and becomes the boundary of the escape and capture cones for r = r (0) s .Therefore, the location of the ODH and IDH on the rotation axis is r = r (0) s , i.e. the location of the spherical photon orbit with b = 0.

Fig. 9
Fig. 9 The boundary of the escape and capture cones projected onto the plane spanned by e 3 and e 1 for r/M = r (0) s /M ≈ 2.431 in the equatorial plane of a Kerr spacetime with a * = 0.99.The projected orthodrome that is tangent to the boundary at (ξ, η) = (π/2, ±π/2) is shown by the dashed line.Because the boundary of the cones is convex upward, the orthodrome is included in the escape cone except at (ξ, η) = (π/2, ±π/2).
can see blue curves of Fig. 10 in advance).

Fig. 11
Fig. 11 The boundary of the escape and capture cones projected onto the plane spanned by e 3 and e 1 for r/M = r IDH /M ≈ 2.092 in the equatorial plane of a Kerr spacetime with a * = 0.99.The orthodrome tangent to the boundary of the cones at (ξ, η) = (ξ + , 0) and (ξ − , π) is shown by the dashed line.Because the boundary of the cones is convex upward, the orthodrome is included in the capture cone.