arXiv:gr-qc/0604015 v1 4 Apr 2006 ste motn olo o etrwy oetatmore extract to It ways 7]. better 6, for 5, look ringdown 4, to 3, the important 2, to then [1, year is up past right the in inspiral possible the become the have from of starting orbits process, few Sim- merger last years. entire few the last of the ulations in black improved of greatly simulations have Numerical holes horizon. hole of dynamics black the the namely, ingredient regime, important this an understanding fea- discuss for we paper, non-perturbative this new, In tures. qualitatively regime contain This equations. may Einstein non-linearities full the beyond the of go confront dynamics and to to us and require grav- theory will the star, perturbation or a holes of black collapse two itational of coalescence the merger the the understanding of say, phase to, However, due arising waveforms gravitational holes. ade- processes black usually astrophysical are involving many particle, understanding point which a for as quate calculations hole black post-Newtonian the treat with along thereof. perturbations black This, and axisymmet- regarding stationary, solutions, Kerr-Newman the intuition ric studying our from of comes Most holes gen- of relativity. them regime eral about dynamical known non-perturbative, is fully little the relatively in holes, black of ing de/ ¶ § ‡ † ∗ URL: fl[email protected] address: Electronic [email protected] address: Electronic [email protected] address: Electronic URL: nsieo udmna dacsi u understand- our in advances fundamental of spite In http://www.cct.lsu.edu/about/focus/numerical/ http://numrel.aei.mpg.de/ 2 a-lnkIsiu ¨rGaiainpyi,Albert-Ei f¨ur Gravitationsphysik, Max-Planck-Institut h ue feeg rsigtehrzn hc ecie ho describes do which radiation. settles horizon, and/or hole the black crossing energy the of that fluxes verify the to used be can calcul which We p hole (iii) were which situations. momentum, dynamical to angular horizons and isolated mass hole tim black spacelike, the import either of be An tion can which evolution. horizon, the time of their signature study In and situations. surfaces time-dependent marginal to horizons hori isolated dynamical on the within work done is This simulations. merical ASnmes 42.m 47.w 53.f 97.60.Lf, 95.30.Sf, 04.70.Bw, 04.25.Dm, numbers: PACS w collision hole black non-axisymmetric mentum. bl a two and of star, collision neutron head-on axisymmetric the namely, cases, hspprpeet us-oa ehdo tdigtephy the studying of method quasi-local a presents paper This edsrb u ueia mlmnaino hs concepts these of implementation numerical our describe We .INTRODUCTION I. nrdcint yaia oiosi ueia relativit numerical in horizons dynamical to Introduction rkSchnetter, Erik URL: ; 1 oiin tt nvriy ao og,L 00,USA 70803, LA Rouge, Baton University, State Louisiana etrfrCmuainadTcnlg,32Jhso Hall, Johnston 302 Technology, and Computation for Center http://www.aei.mpg. ,2, 1, ∗ ar Krishnan, Badri Dtd pi ,2006) 4, April (Dated: senIsiu,A ¨hebr ,D146Gl,Germany Golm, D-14476 M¨uhlenberg 1, Am nstein-Institut, usiuefreethrzn,d aemn useful many have do not horizons, while event surfaces, for trapped substitute it marginally a above, that mentioned shown horizons theoreti- is quasi-local the in on Similarly, works time. cal real in on surface holes Cauchy black to locate a common to instead surfaces is trapped marginally it use simulations, sur- in- numerical Cauchy the In given a locate face. with to horizon event used the be of tersection can known no that is from condition there hori- quasi-local and, perspective, event experimentally relativity that observed numerical implies space- the be This the never of it. can history locate zons entire to order the in know time to have we cept, without quasi-locally holes horizons. black event study The to reference reviews. to for is 26] 25, aim [24, See basic to 23]. horizons related 22, trapping [21, on closely Hayward work by very earlier the turn, by stationary motivated in and Both isolated are, spacetime. models frameworks dynamical which otherwise these an 20], sig- framework in 19, holes a horizon black 18, are isolated 17, horizons the Dynamical [16, of extension holes. nificant black study to methods. non-gauge-dependent extend analysis emphasise to to want want We space- we and non-vacuum 13]). times, and 12, non-axisymmetric 11, to 10, work 9, this [8, axisym- in e.g. detail (see in metry past the in studied numerically been ingredient. necessary a dy- is understanding holes and black itself, namical in task non-trivial a This be etc. can conditions gauge using systems, coordinate performed different simulations different compare two to from and results simulations from information physical ic h vn oio sagoa,tlooia con- teleological global, a is horizon event the Since 15] [14, horizons dynamical of formalism the use We have horizons event and apparent of dynamics The t h oremlioemmnso h black the of moments multipole source the ate ,† 2, lk,o ul i)W eeaietecalcula- the generalize We (ii) null. or elike, o rmwr,wihetnsteearlier the extends which framework, zon c oe,teaiymti olpeo a of collapse axisymmetric the holes, ack lc oegosa tacee matter accretes it as grows hole black a w of kinds various locate We (i) particular: nt erslto.(v eas study also We (iv) solution. Kerr a to wn n lra Beyer Florian and t o-eoiiilobtlaglrmo- angular orbital initial non-zero ith n nrdeti h aclto fthe of calculation the is ingredient ant n pl hmt he pcfi test specific three to them apply and iso yaia lc oe nnu- in holes black dynamical of sics eiul endfraxisymmetric for defined reviously § ,‡ 2, E-0608 LSU-REL-033006 AEI-2006-018, y ¶ 2 properties and can be used fruitfully to study black hole II. BASIC NOTIONS AND DEFINITIONS physics. Thus, it is likely that ideas and results from dynamical/trapping horizons can be very useful for nu- A. Trapped surfaces and apparent horizons merical relativity. Information obtained from the quasi- local horizons complements the information obtained Let S be a closed, orientable spacelike 2-surface in from the event horizon. Once a simulation is complete a 4-dimensional spacetime ( , gab). The expansion of and ready for post-processing, event horizons are use- any such surface can be definedM invariantly without any ful for studying global properties and the causal struc- reference to a time slicing of the spacetime. Since S is ture of the spacetime, and also phenomena such as the smooth, spacelike, and 2-dimensional, the set of vec- topology change of the horizon during a black hole co- tors orthogonal to it at any point form a 2-dimensional alescence. Reliable and computationally efficient codes Minkowskian vector space. Thus, we can define two lin- are now available for locating event horizons (see e.g. early independent, future-directed, null vectors ℓa and [27]). Such information cannot be obtained at the quasi- na orthogonal to S such that local level, which is instead better for tracking the physi- cal parameters and geometry of a black hole in real time. g ℓanb = 1. (2.1) ab − In particular, we consider the following applications: Note that this convention is different from that used in (i) We study the behavior of various marginally trapped [15]. We shall assume that we know a priori what the surfaces under time evolution. This leads to greater outgoing and ingoing directions on are. By conven- insights about the trapped region of a spacetime. An tion, ℓa will denote an outgoing nullM normal and na an important ingredient here is the signature of the world ingoing one. The null normals are specified only up to a tube of marginally trapped surfaces. This world tube is boost transformation known to be null for isolated horizons, and more gener- a a a 1 a ally, it can be either spacelike or timelike; we show that ℓ f ℓ , n f − n (2.2) → → both types occur frequently in numerical simulations. (ii) We give meaningful definitions for the angular mo- where f is a, positive definite, smooth function on S. All mentum, mass, and higher multipole moments for the physical quantities must be invariant under this gauge dynamical black hole. The multipole moments capture transformation. gauge invariant geometrical information regarding the The Riemannian 2-metric q˜ab on S induced by the horizon geometry, and should be useful for understand- spacetime metric gab is ing fundamental issues such as the final state of black q˜ = g + ℓ n + n ℓ . (2.3) hole collapse. For example, we would expect that after ab ab a b a b a black hole has formed and settled down, its multipole The tensor q˜b can be viewed as a projection operator on moments should be identical to the source multipoles of a to S. The null expansions are a Kerr black hole. We show that it is, in principle, pos- sible to verify this conjecture and to calculate the rate at ab ab Θ ℓ = q˜ ℓ , Θ = q˜ n . (2.4) which a black hole approaches equilibrium. (iii) We also ( ) ∇a b (n) ∇a b describe and implement methods for calculating the en- These expansions tell us how the area element of S ergy flux falling into the horizon. This gives us detailed changes as it is deformed along ℓa and na respectively. information on how black holes grow as they swallow ℓa The shear of , σ(ℓ)ab, is the symmetric trace-free part matter and radiation. of the projection of ℓ : ∇a b
This paper is organized as follows. Section II sets c d 1 σ ℓ = q˜ q˜ ℓ Θ ℓ q˜ . (2.5) up notation, and summarizes the basic definitions and ( )ab a b∇(c d) − 2 ( ) ab properties of trapped surfaces and dynamical horizons. a Section III describes the various physical quantities that Similarly, the shear of n is we calculate using dynamical horizons, and also their 1 c d Θ numerical implementation. Section IV presents three σ(n)ab = q˜aq˜b (cnd) (n)q˜ab . (2.6) concrete, well known numerical examples where these ∇ − 2 concepts are applied and finally, section V discusses Note that these definitions only involve derivatives tan- some open issues and directions for further work. Un- gential to S. Thus ℓa and na can, if necessary, be ex- less mentioned otherwise, we use geometrical units tended arbitrarily away from S while computing these with G = c = 1, the spacetime signature is ( , +, +, +), − quantities. all manifolds and fields are assumed to be smooth, and The closed 2-surface S is said to be a trapped surface the Penrose abstract index notation is used throughout. if both expansions Θ(ℓ) and Θ(n) are strictly negative. The derivative operator compatible with the spacetime This is very different from a sphere in normal flat space metric g is and, following Wald [28], the Riemann ab ∇a which has positive outgoing expansion and negative in- tensor is defined via ( )ω = R dω . going expansion. This definition was first introduced ∇a∇b − ∇b∇a c abc d 3 by Penrose [29], who recognized its importance in the a smooth 3-surface foliated by MOTSs. formation of singularities. On a marginal surface, one The existence of MTTs: Numerically, it has been ob- of the two null expansions vanish. Of particular in- served that marginal surfaces (though not apparent terest are the marginally outer trapped surfaces (MOTSs), horizons – see below) usually behave smoothly under for which the outgoing null rays along ℓa have zero ex- time evolution and produce a smooth MTT. This obser- pansion. In addition, we shall mostly deal with future vation is placed on a more rigorous footing by the re- marginally outer trapped surfaces (FMOTSs), i.e., MOTSs cent result of Andersson et al. [33], which proves the lo- with Θ < 0. (n) cal existence of MTTs for a large class of MOTSs. Their There are three main reasons why closed trapped sur- results require the MOTS to be strictly-stably-outermost. faces are important for studying black holes. First, the An MOTS S on Σ is said to be strictly-stably-outermost existence of a trapped surface implies the existence of if there exists an infinitesimal first order outward defor- a singularity in the future [29, 30]. Secondly, they are mation which makes S strictly untrapped. Working with guaranteed to always lie within the event horizon. Fi- a radial coordinate r on Σ such that S is a level set of nally, in stationary spacetimes, the null generators of the r, and r increases in the outward direction, a sufficient event horizon have zero expansion. Thus for stationary (but not necessary) condition for S to be strictly-stably- spacetimes, the cross-section of the event horizon is a 1 ∂ Θ ℓ (r) > S MOTS. outermost is r ( ) 0 everywhere on . Here it is understood that we obtain Θ ℓ as a function of r by While trapped and marginally outer trapped surfaces ( ) calculating Θ ℓ for the constant-r surfaces in the vicin- are defined in the full four dimensional spacetime, in ( ) ity of S. In principle, for an unfortunate choice of r, it numerical relativity, one usually considers trapped sur- might happen that ∂ Θ ℓ < 0 even though there is a faces in conjunction with a foliation of (partial) Cauchy r ( ) different choice for which this condition is satisfied. In surfaces containing S; it is numerically much easier to any case, this is sufficient for verifying that S is strictly- look for closed surfaces on the Cauchy surface rather 2 than in the full spacetime manifold. For concreteness, stably-outermost. This condition, unlike the outermost we shall work in the ADM formalism where the rele- condition for an AH, is a quasi-local condition. We have vant portion of spacetime is foliated by spacelike sur- found in our simulations that most physically interest- faces, and Σ shall denote one of the leaves of this folia- ing MOTSs, such as ones which asymptote to the event tion. However, it will be obvious that the formalism is horizon, and also AHs, satisfy this condition quite gen- applicable no matter how Einstein’s equations are im- erally. However, as we shall see, there exist also MOTSs plemented. which are not strictly-stably-outermost. In practice, in- stead of checking ∂ Θ ℓ > 0 directly, we look for a sur- The trapped region Σ on Σ is defined to be the set of r ( ) points in Σ through whichT there passes a trapped sur- face with a small positive (or negative) non-vanishing face contained entirely in Σ. Note that there could be expansion, and check that it lies completely outside (or inside) the MOTS. points in Σ not contained in Σ, but through which there T passes a trapped surface not contained in Σ. Thus, Σ is It is shown in [33] that if a MOTS S is strictly-stably- T a subset of the intersection of Σ with the 4-dimensional outermost, then at least locally in time, S is a cross- trapped region in the full spacetime. A connected com- section of a smooth MTT. More explicitly, this result ponent of the boundary of Σ is called an apparent hori- shows that given a foliation of the spacetime by Cauchy T zon (AH). Under suitable regularity conditions, the AH surfaces Σt, if there is a MOTS S0 on Σ0 which is strictly- can be shown to be a MOTS [31, 32]. Thus, an appar- stably-outermost, then MOTSs S exist on Σ for ǫ < t t − ent horizon is the outermost MOTS on Σ. Due to this t < ǫ (for sufficiently small ǫ) such that the union St “outermost” property, an AH is not a quasi-local object is a smooth MTT. The MTT will exist for at least as longS on Σ. The behavior of AHs under time evolution can be as the MOTS remains strictly-stably-outermost. This is quite irregular. For example, they can “jump” discontin- a conceptually important result for numerical relativity uously. On the other hand, as we shall soon see, MOTSs because it shows that a large class of MOTSs behave reg- are more regular. ularly under time evolution. How is this to be reconciled with the known fact that AHs can “jump” during a time evolution? The reason is simply because of the outer- B. Dynamical horizons most property. It is possible that a new MOTS can ap- pear on the outside of a given MOTS. The “old” MOTS is 1. Definition and examples then no longer the globally outermost one even though
We can use marginal surfaces to extract physically in- teresting information about the black hole. The key idea 1 More precisely, ∂ Θ ℓ (r) 0 with ∂ Θ ℓ (r) > 0 somewhere on S. is to look not at a single MOTS by itself, but rather a r ( ) ≥ r ( ) 2 It is harder to show that a MOTS is not strictly-stably-outermost. world tube H of MOTSs constructed by stacking up the This can be done by calculating the signature of the horizon (see MOTSs obtained by time evolution. Such a world tube is below) or by calculating the principle eigenvalue of the stability op- called a Marginally Trapped Tube (MTT). An MTT is thus erator defined in [33]. 4 it is locally outermost, and it continues to evolve in a for this spacetime is shown in figure 1. The 4-metric is perfectly regular manner, but it is no longer an AH. There are, as yet, no similar existence proofs for 2M(v) g = 1 ∂ v∂ v + 2∂ v∂ r MOTSs which are not strictly-stably-outermost. How- ab − − r a b (a b) ever, as we shall see later, we find in all the examples + r2(∂ θ∂ θ + sin2 θ∂ φ∂ φ) , (2.7) we have looked at, that MOTSs evolve smoothly even a b a b in this case, forming a regular world tube. and the stress energy tensor is Isolated and dynamical horizons: An MTT is null in equi- librium situations when no matter or radiation is falling ˙ into it; the rest of the spacetime is still allowed to be M(v) Tab = ∂av∂bv (2.8) highly dynamical. This situation is formalized by the 4πr2 notion of an isolated horizon [16, 17, 18, 19, 20]. Us- ing isolated horizons, it has been possible to derive the where M˙ (v) dM/dv. The coordinates (v, r, θ, φ) are ≡ laws of black hole mechanics, use it as a basis for the analogous to the ingoing Eddington-Finkelstein coordi- quantum black hole entropy calculations and find unex- nates for Schwarzschild spacetime. The prescribed mass pected properties of hairy black holes in Einstein-Yang- function M(v) is a positive, non-decreasing function of Mills theory; see [24] and references therein. Most im- the retarded time coordinate v; the Schwarzschild space- portantly for our purposes, isolated horizons have also time is recovered when M(v) is a positive constant. Just proved to be useful in numerical relativity. For exam- as in Schwarzschild, the 2-surfaces r = 2M(v) (for con- ple, isolated horizons provide a coordinate invariant stant v) are FMOTSs. Unlike in the Schwarzschild space- method of calculating the angular momentum and mass time where the 3-surface r = 2M is null and coincides of a black hole [34]. They can be used to obtain bound- with the event horizon, in this case the r = 2M(v) sur- ary conditions for constructing quasi-equilibrium initial face is spacelike if M˙ > 0, and it lies strictly inside the data sets [35, 36]. They might have a role in waveform event horizon. While H is the only spherically symmet- extraction [20]. A pedagogical review of isolated hori- ric dynamical horizon in this spacetime, and there are zons from the numerical relativity perspective can be no spherically symmetric MOTSs outside H, it is shown found in [26]. in [44] that H is not the boundary of the trapped re- In this paper, we are more interested in the dynamical gion. There exist many more non-spherically-symmetric MOTSs and MTTS which come arbitrarily close to the regime when the MTT is not null. A spacelike MTT con- sisting of future-marginally trapped surfaces is called event horizon. Thus the event horizon is the most likely candidate for the boundary of the trapped region [45]. a Dynamical Horizon (DH). Thus, a dynamical horizon is a spacelike 3-surface equipped with a given foliation More generally, figure 2(a) depicts a dynamical hori- by FMOTSs. The properties of a dynamical horizon are zon H bounded by two MOTSs S1 and S2. S is a typi- a studied in detail in [14, 15, 37]. The case when the hori- cal member of the foliation. The vector τˆ is the future a zon is very close to being isolated but still evolving dy- directed unit timelike normal to H, rˆ is tangent to H namically has been studied in [38, 39] and its Hamilto- and is the unit outward pointing spacelike normal to the nian treatment is considered in [40]. Note that the lo- cross-sections. A fiducial set of null normals is cal existence of DHs follows from the local existence of < 1 MTTs because if Θ(n) 0 at any given time, it will con- ℓa = (τˆ a + rˆa) , (2.9) tinue to be strictly negative for at least a short duration. √2 We elaborate on the spacelike property below. 1 na = (τˆ a rˆa) . (2.10) A timelike MTT will be called a timelike membrane √2 − (TLM). A TLM cannot be considered to represent the < surface of a black hole since a time-like surface is not As before, Θ(ℓ) = 0 and Θ(n) 0. The area of a cross- a one-way membrane, and both ingoing and outgoing section S will be denoted by AS and its radius by RS := causal curves can pass through it. In some instances, we AS/4π. A radial coordinate on H will be denoted by shall use the term “horizon” loosely to refer to a generic rp; the cross sections of H are the constant r surfaces. The marginal surface or a MTT without any further quali- 3-metric and extrinsic curvature of H will be denoted fiers. The exact meaning should hopefully be clear from respectively by qab and Kab, and q˜ab is the 2-metric on S. the context. An explicit example of a dynamical horizon is pro- Figure 2(b) shows a Cauchy surface Σ intersecting a vided by the Vaidya spacetime which describes the dynamical horizon H. This intersection S will always gravitational collapse of null dust [41, 42] (see [43] for be assumed to be one of the given cross-sections of H. further examples). This example is not meant to be an The unit timelike normal to the horizon is Ta and the astrophysically realistic model of gravitational collapse, unit outward pointing spacelike normal to S within Σ but it nevertheless provides a good illustration of the is Ra. The three metric and extrinsic curvature of Σ are properties of a dynamical horizon. The Penrose diagram denoted by q¯ab and K¯ ab respectively. The fiducial set of 5
i+ S2
+ H I ℓa a E τˆ io a rˆ na H S v = 0 1 0 = r (a) Dynamical horizon v bounded by S1 and S2.
S2 i− Ta FIG. 1: The Vaidya spacetime. The dashed line indicates the S singularity. This diagram is valid for a strictly increasing mass Σ Ra function M(v) which vanishes for v 0, and asymptotes to a finite value for v ∞. Furthermore,≤ in drawing this diagram, → ˙ < M(v) is assumed to satisfy M(0) 1/16 so that the singular S1 point r = 0, v = 0 is locally naked, and for large v, M(v) is such that the singularity is not globally naked. The Penrose (b) A dynamical horizon intersecting a Cauchy diagram may change qualitatively for other mass functions; surface. see [42] for details. The dynamical horizon is the spacelike surface r = 2M(v) denoted by H, and the event horizon is FIG. 2: The top panel shows a dynamical horizon H. S1 and denoted by E. The shaded portion of the spacetime is flat. a S2 are the initial and final FMOTSs, ℓ is the outgoing null normal, na is the ingoing null normal, rˆa is the unit space- like normal to the cross-sections, and τˆa is the unit timelike Σ normal to H. The bottom panel shows a dynamical horizon null normals to S arising naturally from are and a Cauchy surface Σ intersecting in a 2-sphere S. Ta is the Σ a 1 unit timelike normal to and R is the unit space-like outward ℓ¯a = (Ta + Ra) , (2.11) pointing vector normal to S and tangent to Σ. √2 1 n¯ a = (Ta Ra) . (2.12) √ − 2 Second Law: The area of the cross-sections of a DH in- a A boost transformation of the form of equation (2.2) con- creases along rˆ . This is simply a consequence of < nects (ℓa, na) and (ℓ¯a, n¯ a): Θ(ℓ) = 0 and Θ(n) 0: ℓa ℓa a 1 a 1 = f ¯ , n = f − n¯ . (2.13) Darˆ = q˜ab (ℓ n ) a 2 ∇a b − b When the horizon settles down and becomes null, an 1 a a = (Θ ℓ Θ ) > 0. (2.14) infinite boost ( f ∞) is required to go from (ℓ¯ , n¯ ) to 2 ( ) − (n) (ℓa, na). → If we choose a time evolution vector field ta for which t rˆ > 0, then the area of the dynamical · 2. Summary of basic properties horizon will increase in time, and this result can be called the second law for dynamical horizons. An analogous calculation for TLMs shows that the Topology: The cross-sections of a DH can be either area decreases if Θ < 0, and increases if Θ > 0. spherical or toroidal [14, 15, 21, 33]. Toroidal topol- (n) (n) ogy is possible only in exceptional cases when Foliation: Any given spacelike surface cannot have ℓb σ(ℓ)ab, the scalar curvature R˜ of S, ℓΘ(ℓ), Rab , more than one distinct dynamical horizon struc- and ζa (defined in section III) all vanishL on S [15]. ture on it [37]. This means that a DH can have We shall therefore always take the cross-sections to one, and only one foliation by FMOTSs. This im- be spherical. There are no similar results for cross- plies that if a Cauchy surface Σ does not inter- sections of TLMs. However, we use an apparent sect a given DH in one of the preferred cross- horizon tracker which can only locate spherical sections, then the intersection cannot be a MOTS AHs [46] and therefore all observed MOTSs have at all. Thus, different choices of Cauchy surfaces spherical topology. in general lead to different dynamical horizons. 6
Uniqueness: As mentioned above, the dynamical hori- Σ zon depends on the choice of foliation t, i.e., Singularity choosing a different time slicing would lead to a different dynamical horizon (if it exists). There are however some constraints on the location of dy- a H a ℓa a ℓ 2 namical horizons and trapped surfaces as proved n n H1 by Ashtekar and Galloway [37]. For example, they Σ show that given a dynamical horizon H (along with a mild genericity assumption), there cannot be any trapped surfaces (and therefore no DHs) EH contained entirely in the past domain of depen- Surface of star dence of H. See also [44] for further discussion.
III. APPLICATIONS FIG. 3: Portion of Oppenheimer-Snyder collapse spacetime de- scribing spherically symmetric collapse of pressure-less dust. This section discusses some possible applications of Trapped surfaces are formed the point where the surface of the dynamical horizons. These ideas are illustrated using star intersects the event horizon. There are two sets of trapped concrete numerical examples later in Section IV. surfaces. The ones lying on the timelike surface H1 are the in- ner marginally trapped surfaces while the ones lying on the null surface H2, which is part of the event horizon and also a A. The signatureofa MTT an isolated horizon, are the outermost ones. Since n points away from the trapped region, deforming the inner MOTSs on a > H1 along n makes them untrapped. Therefore nΘ(ℓ) 0 for 1. Background L H1.
As discussed above, MTTs have been shown to exist for a large and physically interesting class of MOTSs, by H1 and H2 in figure 3; H2 is an isolated horizon, i.e., and this is borne out in a large number of numeri- it is null, is part of the event horizon E, and has constant cal simulations where MOTSs are located and evolved area. H1 has decreasing area, is timelike, and goes into smoothly. How many of these MTTs are actually dy- the singularity. A Cauchy surface such as Σ would con- namical horizons? In other words, when is a MTT space- tain two FMOTSs. As expected, Θ ℓ < 0 for H , but Ln ( ) 2 like? The first result in this direction was obtained by nΘ(ℓ) > 0 for H1. There are no spherically symmetric Hayward [21] (see also [34]). Using the Raychaudhuri dynamicalL horizons in this dynamical black hole space- equation for ℓa, it can be shown that an MTT is space- time. like if α < 0, null if α = 0 and timelike if α > 0, where The issue of the signature has been considered in [33]. There it is shown that if a MOTS S is strictly stably out- ab ℓaℓb ab ℓaℓb ermost, and if the quantity σ(ℓ)abσ(ℓ) + Rab is non- σ(ℓ)abσ(ℓ) + Rab α . (3.1) zero somewhere on S (and assuming the null energy con- ≡ Θ ℓ Ln ( ) dition), then the MTT containing S is spacelike in a neighborhood of S. This result is stronger than Hay- In writing this expression, it is assumed that ℓa and a ward’s result (Eq. (3.1)) and it shows clearly that the n are extended off H geodetically, so that Θ ℓ is Ln ( ) spacelike case is physically the most interesting because meaningful. The term in the numerator is strictly pos- ab a b σ ℓ σ + R ℓ ℓ will not vanish in a non-stationary itive in the case of dynamical horizons if the matter ( )ab (ℓ) ab fields satisfy, say, the null energy condition. It vanishes situation. It also shows, somewhat surprisingly, that for isolated horizons. The denominator is negative for even if matter or radiation is falling into a black hole the Vaidya spacetime and also for the stationary Kerr- only in the form of say, a single narrow beam from a par- Newman family. This captures the notion that as we ticular direction, the entire MTT is spacelike. One might go inside the black hole, the outgoing null rays become naively have thought that the MTT would be spacelike more and more converging. Assuming that the numer- only on portions where the energy flux is non-zero, and ator of Eq. (3.1) is nowhere vanishing on H, the hypoth- null otherwise. This is not the case because of the elliptic nature of the equations governing the deformations of a esis that H is spacelike is equivalent to nΘ(ℓ) < 0. As shown by Ben-Dov [47], this lastL condition is not MOTS. satisfied for all MTTs; in Oppenheimer-Snyder collapse In figure 3, the inner MOTS on Σ is not strictly sta- [48], there exists a timelike world tube of FMOTSs with bly outermost; an outward deformation takes it into the trapped region. The results of [33] do not place any re- nΘ(ℓ) > 0. This is illustrated in figure 3 showing a por- tionL of the Oppenheimer-Snyder spacetime. There are strictions on the signature of H1. While the outer MOTS ab ℓaℓb two sets of FMOTSs in this spacetime which are denoted is strictly-stably-outermost, σ(ℓ)abσ(ℓ) + Rab vanishes 7