Abstract

We review recent results about the modelling of gravitational collapse to black holes in higher dimensions. The models are constructed through the junction of two exact solutions of the Einstein field equations: an interior collapsing fluid solution and a vacuum exterior solution. The vacuum exterior solutions are either static or containing gravitational waves. We then review the global geometrical properties of the matched solutions which, besides black holes, may include the existence of naked singularities and wormholes. In the case of radiating exteriors, we show that the data at the boundary can be chosen to be, in some sense, arbitrarily close to the data for the Schwarzschild-Tangherlini solution.

1. Introduction

Higher dimensional black holes have recently been the subject of an increasing number of research works (see [1] for a review) and play an important role in theoretical physics, particularly in String Theory. However, the study of the actual dynamical process of collapse resulting in black hole formation has deserved little attention from the mathematical point of view.

Mathematically, black hole formation can be analysed by constructing appropriate matched spacetimes which settle down through gravitational collapse to black hole solutions. The junction (or matching) of two spacetimes requires the equality of the respective first and second fundamental forms at some matching boundary hypersurface. These conditions amount to a set of differential equations that the spacetime metric functions need to satisfy at that matching boundary. This is often not an easy problem as the metric functions are also required to satisfy the respective Einstein field equations (EFEs) on both sides of the matching surface. In what follows, we review some results about models of black hole formation with emphasis on models resulting from the spacetime matching with an exterior which is either vacuum or a cosmological background.

For a zero cosmological constant , there are well known spherically symmetric examples, one of the simplest being the Oppenheimer-Snyder model which results from the matching of a Friedman-Lemaître-Robertson-Walker (FLRW) dust metric to a Schwarzschild exterior. For spherical cases, the results are less well-known: the matching of -FLRW with the -Schwarzschild solution (also called Kottler solution) has been investigated by Balbinot et al. [2], Nakao [3], and Markovic and Shapiro [4] and its -Lemaître-Tolman-Bondi counterpart by Lake [5]. The matching of the dust -Szekeres solution with Kottler was recently studied by Debnath et al. [6] (following a work of Bonnor [7] for ), while the matching of a collapsing fluid with tangential pressure and to Kottler has been investigated by Madhav et. al. [8]. The above examples give rise to black hole solutions with in spherical symmetry.

The collapse to nonspherical black holes has been less studied. Smith and Mann [9] have shown that one can match a collapsing FLRW spacetime to an asymptotically anti-de Sitter (AdS) exterior, as a model of gravitational collapse to higher genus asymptotically AdS black holes. In a related work, Lemos [10] matched a flat FLRW metric to a radiating Vaidya exterior whose collapse results in a toroidal black hole. Subsequently, there has been considerable interest in models of toroidal and higher genus black holes (also called topological black holes), shells, and horizons, see, for example, [1114], partly due to the existence of a “landscape’’ of vacua states in String theory with positive, negative, and zero (see, e.g., [15]). More recently, [16] generalised the earlier work and constructed models for the collapse of inhomogeneous and anisotropic fluids to topological asymptotically AdS black holes in planar and hyperbolic symmetry.

There is also a vast literature about the formation of black holes in cosmological models. There are, essentially, three types of models which are geometrically different, namely, models with (i) junctions with no shells where two metrics are involved, for example, swiss-cheese-type models (see, e.g., [17] and references therein); (ii) junctions with shells (these can be thin or thick), for example, models with domain walls and bubbles mentioned in the previous paragraph; (iii) no junctions, where a single metric is involved in the problem, for example, models of collapsing over-densities resulting from perturbations on cosmological backgrounds.

In the context of primordial black hole (PBH) formation, the gravitational collapse due to first-order density perturbations in FLRW models of the early universe radiative phase was first studied by Zeldovich and Novikov [18] as well as Hawking [19]. These works were then generalised [20, 21], in particular, to include black hole formation from the collapse of cosmic strings produced in phase transitions with symmetry breaking [22]. The collision of bubble walls formed at phase transitions was also considered as sources of PBHs in [2325]. Different scenarios were subsequently studied and a review is in [26]. More recently, cosmological black hole formation due to QCD and electroweak phase transitions were studied by, for example, [27] following earlier works of [28] and bubble collision has been investigated in detail (for various potentials and using planar, hyperbolic, and spherical symmetry in both dS and AdS phases) by, for example, [11, 29]. Furthermore, PBH formation in cosmological models with has very recently been analysed in [30], which contains an interesting review of numerical results up to the year 2011 (for a review on the case see [31]). In turn, the collapse to black holes in FLRW models have also been investigated using nonlinear scalar metric perturbations on radiative [32] and scalar field backgrounds [33], with those recent works establishing interesting observational bounds on the cosmological spectra nongaussianity.

In this paper, we will be mostly interested in cases which involve spacetime junctions with no shells and, in particular, we will focus on nonspherical spacetimes containing a non-zero . We start by reviewing the theory of spacetime junctions and we then summarise recent results about the gravitational collapse to black holes in higher dimensions [34] which generalise the results of [16]. We start by considering a family of solutions to the -vacuum Einstein equations in dimensions which contains black hole solutions. We find some possible interior collapsing solutions with dust as source and study the corresponding matching problem. We then analyse the global geometrical properties of the matched spacetime and find conditions for which the spacetimes have no naked singularities.

We also review the interesting case of a model of radiating gravitational collapse. In particular, we take the anisotropic Bizoń-Chmaj-Schmidt (BCS) solution in dimensions and prove that it can be matched to some collapsing dust interiors. The BCS solution is known [35] to settle down via radiation to the Schwarzschild-Tangherlini (SchT) solution. In turn, the latter is known to be stable [35, 36] and we show that, for some interiors, the data at the matching surface can be chosen to be arbitrarily close to the data for a SchT exterior [34]. So, the resulting spacetime models the gravitational collapse of a fluid in five dimensions with an exterior emitting gravitational waves which settles down to the SchT solution.

2. Geometric Theory of Spacetime Junctions in Brief

Let and be two -dimensional space times with oriented boundaries and , respectively, such that and are diffeomorphic. The matched space time is the disjoint union of with the points in identified such that the junction conditions are satisfied (Israel [37], Clarke and Dray [38] and Mars and Senovilla [39]). Since are diffeomorphic, one can then view those boundaries as diffeomorphic to a 3-dimensional oriented manifold which can be embedded in and . Let and be coordinate systems on and , respectively, where and . The two embeddings are given by the following maps: such that . The diffeomorphism from to is .

Given the basis of the tangent plane at some , the push forwards map into three linearly independent vectors at represented by : On the other hand, using the pull backs of the maps , the metrics can be mapped to given two symmetric 2-covariant tensors and whose components in the basis are The first matching conditions are given by the equality of the first fundamental forms (Israel [37]) The bases and can then be identified, as can the hypersurfaces , so henceforth we represent both by .

We now define a 1-form , normal to the hypersurface , as The vectors constitute a basis on the tangent spaces to at . Since the first junction condition allows the identification of with , we only have to ensure that both bases have the same orientation and that is satisfied in order to identify the whole 4-dimensional tangent spaces of at , .

The second fundamental forms are given by and the second matching conditions, for nonnull surfaces, are the equality of the second fundamental forms We note that these matching conditions do not depend on the choice of the normal vectors. The first matching conditions ensure the continuity of the metric across while the second conditions prevent infinite jump discontinuities in the Riemann tensor so that the Einstein field equations are well defined in the distributional sense.

The theory is also fully developed for the case where the matching boundary is null or changes character. In that case, the normal vectors are substituted by the so-called rigging vectors, which are vectors nonwhere tangent to and are used to define generalised second fundamental forms, see [39] for the details.

We are interested in the particular cases for which the matching surface inherits a certain symmetry of the two space times . Such matching is said to preserve the symmetry. In practice one demands that the matching hypersurface is tangent to the orbits of the symmetry group to be preserved, see [40].

We note that if the exterior is a vacuum spacetime and the interior contains a fluid then the normal pressure of the fluid has to vanish at the boundary . This is a consequence of the fact that the matching conditions imply We also note that if the interior contains dust then the boundary is ruled by geodesics, and we will explore this fact to solve the matching problem in what follows.

3. Black Holes in Higher Dimensions

An -dimensional Lorentzian manifold is a solution of the Einstein equations with cosmological constant and energy-momentum tensor if its Ricci tensor satisfies where is a constant, and the indices . We will consider the family of higher-dimensional black holes given in the following proposition (see also [41, 42]).

Proposition 3.1. Let be an -dimensional Riemannian Einstein manifold with Ricci scalar , and let where , are constants. If is an open interval where is well defined and does not vanish then the -dimensional Lorentzian manifold given by and is a solution of the vacuum Einstein equations with cosmological constant .

It is easy to see that the above metric generalises the Kottler metric to arbitrary dimensions. Nevertheless, black holes with the above metric do not immediately integrate into the usual intuition of a black hole in four dimensions. For instance, since the metric on the sections of null infinity does not need to be a metric of constant curvature, the spacetimes are not asymptotically flat (Nor asymptotically de Sitter (dS), resp., anti de Sitter (AdS), in the case of , resp., .); however, by analysing the equations for null geodesics in and , one can show that the spacetimes are weakly asymptotically simple [34].

Proposition 3.2. The metrics (3.3) are conformally compactifiable at infinity and are weakly asymptotically-simple.

Therefore, by comparison with the four-dimensional case [16], the following properties hold.(i)With , and , has a single zero, corresponding to an event horizon, and asymptotes to a positive constant at large . This is a black-hole solution, which can be thought of as generalising the Schwarzschild metric. The (degenerate) metric on the horizon is , which is also the conformal metric on future null infinity .(ii)With , and large enough positive , is positive between two zeroes, corresponding to a black-hole event horizon and a cosmological event horizon. The solution generalises the asymptotically dS Kottler solution.(iii)With and , again has a single zero, corresponding to an event horizon, and the solution generalises the asymptotically AdS Kottler solution.(iv)The solutions in the previous class with may have no global symmetries except the staticity Killing vector. This is because compact, negative scalar curvature Einstein manifolds have no global symmetries, nor does, for example, the Ricci-flat metric on (an example with and ).

3.1. Generalised Friedman-Lemaître-Robertson-Walker Interiors

We will fill in the solutions of the previous section with (interior) collapsing dust solutions, so that the resulting global spacetime is a model of gravitational collapse with a black hole as the end state. In order to do that, we take the following classes of generalised FLRW spacetimes.

Proposition 3.3. The -dimensional Lorentzian metric, is a solution of the Einstein equations with cosmological constant and energy-momentum tensor , corresponding to a dust fluid with density and velocity , if and only if and satisfy the conservation equation: for constant , and the Friedman-like equation:

We consider the matching of an interior metric (3.4) to a static exterior (3.3) at a hypersurface ruled by radial time-like geodesics in (3.3), that is, a surface of constant (say ) in (3.4). We then solve the matching conditions (2.5) and (2.9) together with the EFEs (3.6) to get the following result [34].

Theorem 3.4. The metric (3.3) can be matched to the FLRW-like metric (3.4) at provided that and .

The matching boundary is comoving with the collapsing fluid whose dynamics is given by (3.6). We call the matched spacetime of the previous theorem FLRW-Kottler spacetime. We now summarize some properties of the FLRW-Kottler spacetime in the three cases , , and : when the Einstein manifold is not cobordant to a point (e.g., ) the solutions we find cannot have a regular origin, though they can be regular with spacetime wormholes or a “cusp” at the origin. When there is a singularity at the origin, it may or may not be visible from infinity. If the singularity is visible it is called naked singularity as it is not hidden by an horizon. The next three propositions exploit these aspects, and Figures 1, 2, and 3 show Penrose diagrams which represent the global structure of the matched spacetimes in the cases , , and [34]. In particular, by comparing the conformal lifetime of the FLRW universe with the supremum of the possible values of we show [34] the following.

Proposition 3.5. If (hence ) and is not an -sphere then the locally naked singularity of the FLRW-Kottler spacetime at is always visible from future null infinity for , but can be hidden if and (space-time dimension ).

Interestingly, [43] investigated higher dimensional dust collapsing spacetimes with and also found that for dimension the final central singularity is always hidden by an horizon, while for naked singularities may appear and thus the cosmic censorship conjecture is violated. By construction, their collapsing interiors are spherical, while here we allow interior spacetimes with other topologies.

Proposition 3.6. If (hence ) and is not an -sphere then the locally naked singularity of the FLRW-Kottler spacetime at can be always hidden except if the FLRW universe is recollapsing (hence ) and .

The case has wormhole solutions that may have causal curves which cross from one to the other and may violate causality. The next proposition, in particular, states conditions under which this can be avoided (see also a result of Galloway [44]). Those conditions are obtained by showing that the future (resp., past) horizons hit the matching surfaces at marginally outer trapped (resp., antitrapped) surfaces which are spacelike and, in turn, that these two curves touch at .

Proposition 3.7. For the FLRW-Kottler spacetime satisfies the following.(1) If and is not an -sphere then the locally naked singularity of the FLRW-Kottler spacetime at can always be hidden.(2)If then the cusp singularity is not locally naked.(3)If then no causal curve can cross the wormhole from one to the other.

3.2. Generalised Lemaître-Tolman-Bondi Interiors

We consider now higher dimensional versions of the inhomogeneous Lemaître-Tolman-Bondi (LTB) solutions, generalising those of [45].

Proposition 3.8. The -dimensional Lorentzian metric, is a solution of the Einstein equations with cosmological constant and energy-momentum tensor , corresponding to a dust fluid with density and velocity , if and only if , and satisfy for some functions and , and (where dot and prime denote differentiation with respect to and ).

This metric has three free functions of , namely, , , and , one of which can be removed by coordinate freedom. We seek to match an interior represented by the metric (3.7) to a static exterior represented by the metric (3.3) at a surface ruled by radial time-like geodesics in (3.3), that is, a surface of constant (say ) in (3.7). The solutions to the matching conditions (2.5) and (2.9) allow us to prove [34].

Theorem 3.9. The metric (3.3) can be matched to the LTB-like metric (3.7) at provided that and .

Note that the matching boundary is again comoving with the collapsing inhomogeneous fluid and its dynamics is now given by (3.9). We will use the term LTB-Kottler space-times for these matched solutions. The global properties of the LTB-Kottler spacetime obtained in Theorem 3.9 are much more diverse than what we find in the FLRW-Kottler spacetimes. For instance, one can easily find examples of black hole formation with wormholes inside the matter with positive and , see [34] (similar results in 4 dimensions are in [46]).

4. Collapse with Radiating Exteriors

In this section, we consider models of gravitational collapse with a gravitational wave exterior, so that the exterior metrics will be time-dependent generalisations of (3.3). This is a specially interesting problem for observational and experimental purposes [36]. We take the Bizoń-Chmaj-Schmidt (BCS) metric in dimensions [36], although similar ansetz can be made in other dimensions (see [47]). We then consider three different interiors with this exterior, which are anisotropic generalisations of FLRW models [34]. Consider the metric: where , and are functions of and . The one forms are left invariant for the standard Lie group structure on and can be taken as where are Euler angles on with , , and . The space-time with is interpreted as containing pure gravitational waves with radial symmetry [36] and the Schwarzschild-Tangherlini (SchT) limit is obtained for . There is a residual coordinate freedom ; in (4.1), which one can use to set arbitrarily along a timelike curve. The dimensional vacuum EFEs give [36] together with the quasi-linear wave equation for In [36], the authors solve this system numerically by giving and at with . Within the BCS class, they provide numerical evidence that linear perturbations of SchT decay in time so that SchT is linearly stable, in this sense. The nonlinear stability of SchT within the BCS class was later established in [35].

We will give data , and the normal derivative at the timelike boundary of the collapsing interior with the gauge choice . Uniqueness and local existence follow as standard. From [35], one knows that if data close to that for SchT is given on an asymptotically flat hypersurface then the solution will exist forever and stay close to the SchT solution. As we will see, in some of our cases, data on can be chosen to be arbitrarily close to data for SchT.

4.1. Generalised FLRW Interiors

As interior metrics, we will consider three classes of FLRW-like solutions based on Riemannian Bianchi-IX spatial metrics which are, respectively, the Eguchi-Hanson metric (with ), the -Eguchi-Hanson metric (with excluding the case ), and the -Taub-NUT metric (with , including, as a particular case).

We summarize now our main result and leave the details of the proof to the next three sections (see also [34]).

Theorem 4.1. In each case, the interior metric gives consistent data for the metric (4.1) at a comoving time-like hypersurface. Local existence of the radiating exterior in the neighbourhood of the matching surface is then guaranteed. In the case of Eguchi-Hanson and -Taub-NUT with , the data can be chosen to be close to the data for the Schwarzschild-Tangherlini solution.

4.1.1. The Eguchi-Hanson Metric

Eguchi and Hanson found a class of self-dual solutions to the Euclidean Einstein equations with metric given by [48] with given by (4.2) and is a real constant. The generalised FLRW metric built on this Riemannian metric is with the EFEs for a dust source reducing to We match at so that is parameterised by . The equality of the first and second fundamental forms at then gives from which we calculate , , and on in terms of quantities from the interior. At this point, we have , , and on , and we use the gauge freedom to set on . By (4.10), we have and, to say that the data is close to SchT data, we want these to be small. The first term is small if . The second term will increase without bound as decreases to zero. If we restrict by its value when a marginally outer trapped surface forms on then, from the Friedman equation and with , this happens when so that we control on by controlling . So, by choice of the location of , at , and choice of , we can choose data close to SchT.

4.1.2. The -Eguchi-Hanson Metric

The -Eguchi-Hanson metric is given by [49] where . This metric is complete for and , , where in an integer. Since is related to for a complete solution, we cannot obtain the previous case from this case by taking . However, the matching formulae do formally allow this limit. The generalised FLRW metric built on this metric is and the respective EFEs for a dust source reduce to The matching conditions on give We can calculate in terms of interior data, as before, and check that this is consistent with the EFEs (4.3) and (4.4). However, in this case, it is not clear that we may choose data close to SchT data as the normal derivative of is and for this to be small, we would require to be large on outside the marginally trapped surface. It is hard to see how to arrange this and so although the solution in the exterior exists locally, we do not have a good reason to think that it will settle down to SchT.

4.1.3. -Taub-NUT

We consider the Riemannian Taub-NUT metric with a cosmological constant ( rather than with our conventions) [50, 51] where , and use it to construct the interior: The EFEs for a dust source are again (4.13). From matching the first fundamental forms we get and the second matching conditions read As before, the expression for on is consistent with calculated from (4.3) and (4.4). Now note that if then the metric (4.16) is precisely the 4-dimensional hyperbolic metric. In this case, and vanish on whatever is the value of , so that the exterior metric is precisely SchT: this is a case from the previous section as the interior is now a standard FLRW cosmology. Consequently, if we take close to we get data close to SchT data. To see this, set . Then , so that and the data can be chosen as small as desired by choosing large and small .

Acknowledgments

The work summarised in this paper was done in collaboration with José Natário and Paul Tod. The author is supported by FCT Projects PTDC/MAT/108921/2008, CERN/FP/116377/2010, and Est-C/MAT/UI0013/2011 and by CMAT, Universidade do Minho, through FCT plurianual FEDER Funds “Programa Operacional Factores de Competitividade COMPETE.”