Charged radiation collapse in Einstein–Gauss–Bonnet gravity

We generalise the continual gravitational collapse of a spherically symmetric radiation shell of matter in five dimensional Einstein–Gauss–Bonnet gravity to include the electromagnetic field. The presence of charge has a significant effect in the collapse dynamics. We note that there exists a maximal charge contribution for which the metric functions in Einstein–Gauss–Bonnet gravity remain real, which is not the case in general relativity. Beyond this maximal charge the spacetime metric is complex. The final fate of collapse for the uncharged matter field, with positive mass, is an extended, weak and initially naked central conical singularity. With the presence of an electromagnetic field, collapse terminates with the emergence of a branch singularity separating the physical spacetime from the complex region. We show that this marked difference in singularity formation is only prevalent in five dimensions. We extend our analysis to higher dimensions and show that for all dimensions N≥5\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N\ge 5$$\end{document}, charged collapse ceases with the above mentioned branch singularity. This is significantly different than the uncharged scenario where a strong curvature singularity forms post collapse for all N≥6\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N\ge 6$$\end{document} and a weak conical singularity forms when N=5\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N=5$$\end{document}. A comparison with charged radiation collapse in general relativity is also given.


Introduction
The allure of gravitational collapse is undeniable in gravitational physics, cosmology and astrophysics. This is due to it being a fundamental mechanism in the structural formation of the entire universe. Over a time period, initially smooth matter fields will contract to form small localised dense distributions which form the building blocks for galaxies, stellar groups, stars and planets. Galaxy formation is thought to be a e-mail: drbrasselint@gmail.com b e-mail: maharaj@ukzn.ac.za (corresponding author) c e-mail: Goswami@ukzn.ac.za the result of monolithic collapse, i.e. every galaxy formed as the result of the contraction of a single, turbulent gas cloud [1]. Gravitational collapse is also prominent in the entire life cycle of a radiating star. Stars form as a result of the gravitational collapse of localised dense regions of molecular clouds in the interstellar medium [2], often referred to as star-forming regions or nurseries. Core contraction continues in a subtle fashion (called gravitational compression) throughout the star's life as it burns off its matter content in the nuclear fusion process. At the end of the luminous phase of the star's life cycle, when all the nuclear fuel has been exhausted and the matter content consists of the mostly very stable atomic nuclei 56 Fe, 56 Si or 60 Zn, it will experience inwardly directed contraction due to gravity. Due to the enormous amount of energy being released in such a process, it remains a paramount case study in general relativity, modified gravity, high energy physics, astrophysics and even, due to the above mentioned stellar nucleosynthesis, nuclear physics. The end result of such a contraction is either a more compact object like a neutron star or, in the event of a more intense electromagnetic field, a pulsar or magnetar. Otherwise, in the case of truly massive stars, with masses above the Tolman-Oppenheimer-Volkoff limit, collapse proceeds in a continual manner without reaching an equilibrium end state. The singularity theorems of Penrose [3] then predict that the end state of gravitationally induced stellar collapse is a spacetime singularity.
The singularity theorems envisage that the end state of gravitational collapse is the formation of spacetime singularities [4]. Their incidence depends on certain circumstances such as gravity's inherent attractive nature, the existence of trapped surfaces and the preservation of causality. These conditions hold within a variety of gravitational theories. There also exist a large family of solutions to Einstein's field equations which are geodesically incomplete. In fact, one of the most important requirements for geodesic incompleteness is the existence of trapped surfaces [3]. These theorems do not take into account the nature of the spacetime singularities, only their formation. Whether or not it is possible for any future directed null geodesics to escape the confines of a black hole to infinity, leaving the singularity naked, is not regarded by the theorems of singularity formation. Penrose [5] proposed the cosmic censorship conjecture to avoid naked singularities. Should any matter distribution that is deemed physically realistic, collapse under its own gravity, the end state of this contraction must be a spacetime singularity which is covered by a trapping horizon for its entire existence. Therefore, the end state of collapse is a black hole with a central curvature singularity covered by the event horizon which acts as an aegis between it and the external universe. The event horizon is defined as the smooth boundary of the causal past of null infinity. As it stands there are no proofs for the conjecture and, in fact, the Vaidya-Papapetrou model [6,7] is one of the earliest models proposed that counters cosmic censorship. Many other counterexamples exist in different contexts for certain matter distributions, for example [8][9][10][11][12][13][14][15][16]. It turns out that an increase in dimension may restore cosmic censorship in general relativity. Under certain physical conditions, it was demonstrated that collapse terminates with the formation of singularities which are no longer naked, in higher dimensions [12,17,18].
The generalisation of the Vaidya spacetime, which contains the additional Type II null string fluid, has been extensively studied in regards to the exteriors of radiating relativistic stars, dynamical black holes and even regular black holes. Various properties of the additional Type II null fluid were examined extensively in [19][20][21][22][23]. Contained within these models were the collapsing monopole, de-Sitter/Anti de-Sitter, the charged Vaidya models and the radiating dyon solution. With regard to regular black holes, several models were analysed in [24][25][26][27][28][29]. It must be mentioned that these models are in a sense contrived since the regularity of the solutions are artificial; coordinate transformations were used to design them. Some examples of similar regularity constructions can be found in [30][31][32].
In modified theories of gravity, Dominguez and Gallo [33] analysed black hole solutions in Einstein-Gauss-Bonnet gravity and Chowdhury et al. [34] found solutions in f (R) gravity using the idea of R-matching, an inescapable feature of f (R) collapse. The reason for modifying general relativity lies in the fact that it is a global theory of gravity and is therefore incomplete. General relativity itself extends Newtonian gravity and so considerations for its own extension are natural. It was shown by Lovelock [35,36] that a polynomial form of the Lagrangian is possible leading to the Lovelock action. Therefore higher order curvature terms appearing in Lovelock gravity are corrections to general relativity. In fact, general relativity can be considered as first order Lovelock gravity, since it is the final limiting case. If the Lovelock polynomial is quadratic in order, we then have the Einstein-Gauss-Bonnet (EGB) action. This is then second order Love-lock gravity, or EGB gravity. Quadratic terms in curvature present as corrections to general relativity. The solution of Boulware and Deser [37] is indeed the higher dimensional EGB analogue of the Schwarzschild solution. As an ansatz to the radiating solution of Vaidya, the Boulware-Deser metric was studied as a radiating solution by Kobayashi [38] and Brassel et al. [39]. Several solutions with equations of state were found as EGB analogues of the Vaidya solutions found by Husain [19] and Wang and Wu [20]. The role of dimensions in any gravitational theory is an important one. It was demonstrated in [40] that naked singularities are inevitable upon the termination of radiation collapse in EGB gravity. For the five dimensional case, a massive timelike naked singularity formed post collapse while in dimensions of six and higher, a massless ingoing naked singularity resulted. The gravitational collapse of the radiating Boulware-Deser spacetime was analysed in detail by [41,42] where it was shown that in five dimensions, the minimum dimension of the EGB theory, the collapse terminates with an extended weak and conical central singularity that is initially naked, before succumbing to a trapping horizon. This is a remarkable feature not present in general relativity. Conical singularities have been discussed in the literature by [43][44][45]. It should be noted that in all dimensions greater than six, collapse terminates with a strong curvature singularity for positive mass, mirroring general relativity, and thus highlighting the importance of the fifth dimension in EGB gravity. Weak and strong curvature singularities were discussed in [46]. Consider a volume element defined by linearly independent spacelike Jacobi fields 1 with zero vorticity propagating along null (timelike) geodesics or is orthogonal to its tangent vector. A gravitational curvature singularity is called strong if such a volume element vanishes at the spacetime singularity. Otherwise it is finite and called weak. The necessary and sufficient conditions for the formation of strong curvature singularities were given by Clarke and Królak [47].
An important feature of gravitational collapse, and one of the bases of this paper, is the notion of the electromagnetic field contribution. Charged collapse is well known in the literature since the discovery of the Reissner-Nordström [48][49][50][51] and Kerr-Newman [52,53] metrics. In general relativity, charged collapse has been analysed in detail, for example see [54][55][56]. In modified theories of gravity, Ghosh [57] studied radiating black holes in five dimensional Einstein-Yang-Mills-Gauss-Bonnet (EYMGB) gravity, discussing the effect of the Yang-Mills gauge charge on the structure and locations of trapping horizons. Hansraj [58] obtained solutions in static Einstein-Gauss-Bonnet-Maxwell (EGBM) theory by utilising a coordinate redefinition which led to several new solutions. The charged Boulware-Deser metric, which is the basis of this treatment, was first studied in detail by Wiltshire [59] where a generalisation of the Birkhoff theorem was proved with the effect of the Gauss-Bonnet term. The conditions for horizon formation were also provided and compared to the general relativistic limit. Abbas and Ahmed [60] studied the contraction of a charged spherically symmetric fluid in f (R, T ) gravity where the trapped surfaces, apparent horizon and singularity structure were addressed. More recently, the effects of the electromagnetic field on conformally flat collapsing stellar configurations with an anisotropic heat flow in f (R) gravity were studied in [61]. It was shown that the effect of the charge contribution and heat flux slowed down the rate of collapse. It should be noted however that in general relativity, the charge compression in collapsing a charged mass is significantly higher (around forty orders of magnitude) than the gravitational attraction. Therefore it is unlikely that a black hole forming in nature will have any significant charge contribution. Whether this is the case in modified theories of gravity is unknown in general.

This paper
We analyse the continual gravitational contraction of a charged spherically symmetric shell of radiation in five dimensional EGBM gravity. We will show that the final fate of such a collapse is a branch-like singularity which separates the physical collapsing spacetime from an unphysical region resulting from the charge contributing to the metric becoming complex. We compare this with the uncharged case where radiation collapse ceases with the formation of a necessarily weak and conical central singularity, under certain conditions for the mass function, which remains initially naked for a time before becoming shrouded by a trapping horizon within a black hole. This paper is organised in the following way: In Sect. 2 we give an overview of units and the role of dimension in gravity. In the following section we highlight some definitions indicative of the EGB theory. The gravitational collapse of a radiating spacetime within an electromagnetic field is considered in Sect. 5, where we highlight the differences in the collapse process with both the uncharged scenario and the general relativistic case. Finally, we briefly look at collapse in arbitrary dimensions, noting that the dimension critically alters the collapse dynamics of radiating matter in EGB gravity. A brief recap of charged radiation collapse in general relativity is given in Appendix B.

Higher dimensional preamble
The role of dimensions in all gravity theories is a germane one. In this research paper, we will utilise units where G = c = 1. In arbitrary dimensions, Einstein's coupling constant is expressed in terms of the factorial function and is given by The surface area of the (N − 2)-sphere is given by the following where we note the presence of the Gamma function Γ (. . .).
Since EGB gravity is only relevant in five dimensions or higher, these terms become important in the analysis, especially with regards to the electromagnetic field. The Lorentzian signature of the spacetime manifold is (−, +, +, . . . , +). These expressions become especially critical when considering a charged radiating stellar mass, or black hole, which is the basis of this paper. The surface gravity of a black hole is not well defined in general. However, it is possible to define it for a black hole whose event horizon is a Killing horizon. 2 Every stationary (non-rotating) black hole has a Killing horizon [62,63]. The surface gravity K of a Killing horizon is the acceleration (as exerted at infinity) required to keep an object at the horizon of a black hole. It can be derived from the Raychauduri equation [64,65] and is given by where k a is a Killing vector. In the above, the semicolon represents covariant differentiation.

Second order Lovelock gravity
The Lovelock (or Lovelock-Zumino) action in arbitrary dimensions is given by where we have R a r b r c r d r , In the context of Lovelock gravity in higher dimensions, some spherically symmetric exact solutions for cosmological backgrounds were found by Bajardi et al. [66]. They also obtained solutions in five dimensional, topological Chern-Simons gravity where it was shown that for particular Lovelock parameters, Anti-de Sitter invariant Chern-Simons gravity is contained in Lovelock gravity in five dimensions. If we consider second order Lovelock gravity, the action reduces to where α 0 = Λ is the cosmological constant term, α 1 is the constant (usually unity) associated with the Einstein-Hilbert action (R = R) and α 2 = α is the constant associated with the Gauss-Bonnet term. It is important to note that in order to avoid any pathologies, α > 0 [37]. In the above, is the second order Lovelock Lagrangian. Varying the second order action (5) with respect to the metric g ab gives the second order Lovelock or EGB field equations Here, G ab is the Einstein tensor, T ab is the energy momentum tensor and H ab is the Lovelock tensor which arises as a consequence of varying the action (5). It is given by The nontrivial nature of this expression requires that the minimum dimension of the spacetime in EGB gravity has to satisfy N ≥ 5. If we have that N < 5, the Lovelock tensor (8) vanishes, and there is no contribution from the higher curvature of the theory. 3 In the limit where α → 0, general relativity will be regained in five dimensions. For the case of a nonvanishing electromagnetic field, the Einstein-Gauss-Bonnet-Maxwell (EGBM) field equations are written as for a relevant energy momentum tensor T ab and the electromagnetic energy tensor E ab . In the above, J a is the current and the fluid N -velocity is u. The electromagnetic energy tensor can be written as and is expressed in terms of the Faraday tensor F ab = Φ b;a − Φ a;b , the gauge potential Φ and the surface area (2). We note that the tensor (10) is trace-free only in four dimensions.

The generalised radiating Boulware-Deser solution
The EGB analogue of the Schwarzschild solution was found by Boulware and Deser [37]. Similar to the Vaidya solution, this metric can be written in Eddington-Finkelstein coordinates and then can be made to radiate [38,39,41,42]. The five dimensional radiating metric is given by where We have that M(v, r ) represents the gravitational mass of the five dimensional hypersurface. In the above expression (12), the negative branch solution has the general relativity limit as α → 0 [37]. The solution (12) comes about as a result of the generalised uncharged Type II matter with a string fluid. The negative branch solution (12) is essentially the EGB analogue of the Vaidya spacetime, to which it reduces in the general relativity limit when α → 0. This motivation is necessary to compare the EGB effects to those in the general relativistic limit. The energy momentum tensor is written as where we have that with the restrictions l c l c = n c n c = 0, l c n c = −1.
The null vector l a is a double null eigenvector of the energy momentum tensor (13). The EGB field equations (7) take the simple form where The Einstein tensor components and Lovelock tensor components will be presented in Appendix A. They are complicated in form but their union in (7) yields the simple expressions above. The energy conditions for the Type II fluid (as shown in [21,67,69]) are -The null energy condition: -The weak energy condition: -The strong energy condition: -The dominant energy condition: These are equivalent in the five-dimensional general relativity limit. Using (14) and (15), the surface gravity K can also be written as

Charged Boulware-Deser collapse
The charged analogue of the Boulware-Deser metric was first analysed in detail by Wiltshire [59]. If we consider the mass function in (12), we make the choice which gives 4 which is analogous to the choice made in the Appendix (B.7). This choice comes about by following the methodology used in Wang and Wu et al. [20]. We then have the metric for the charged Boulware-Deser spacetime where we now have The solution (25) was first found by Wiltshire [59]. In the limit when α → 0, the above reduces to the Vaidya metric 4 The quantities (1) and (2) evaluate to in five dimensions.
(B.8) with nonzero charge q(v) (see Appendix B). The five dimensional gauge potential and five-current are given by The electromagnetic field contribution (10) is built into the above definition (23) in five dimensions. In the above, Q = Q(v) is the charge indicative of the Einstein-Gauss-Bonnet-Maxwell (EGBM) theory. Note that the five dimensional Einstein constant and surface area of the 3-sphere are incorporated into this definition analogously to (B.7).
Using (16), we can write the EGBM field equations (9) with the mass function (23) as We note that in the absence of charge, M(v, r ) = M(v) in the above, and only one field equation (26b) will survive as expected. 5 For the metric (25) the energy conditions (17)-(20) must be satisfied. From the field equations (26), this gives Similarly to the general relativity case (shown in Appendix B), the condition (27) is violated near r = 0 [70,71]. For the case of the metric (25) the reason for this violation is different to the general relativity case, which we explain in the paragraphs below. The collapse dynamics of the charged Boulware-Deser spacetime differ significantly to its uncharged counterpart [41]. We observe from (25) that has a negative sign, there must exist a maximal charge contribution before which this metric becomes complex. Therefore a branch singularity is generic for the charged collapse in EGB gravity [72,73]. Considering the square root term in (25), this only remains real if otherwise it becomes complex. There exists no real spacetime whenever the above inequality is violated. This inequality can be written as which has six solutions for r = r s (v), two of which are real and one of which is positive. There exists a branch singularity at r = r s (v) for any value of the mass function M, if Q = 0, which separates the real spacetime from the complex metric.
The positive real solution to the above inequality is given by where we have set E = 1536M 3 α 3 + 81Q 2 α 2 . This branch singularity forms at r s (0) = 0 = v and extends into the future. The value of the inequality determines the domain of r which is 0 < r s < r < ∞; there is no spacetime at r = 0, which is one reason the weak energy condition (27) is violated. The Kretschmann invariant K = R abcd R abcd can be used to check whether any spacetime manifold has singularities. For the charged metric (25) it is given by where we have set F = 1 + 8α M r 4 − 2α Q 2 r 6 . We note the r 6 + 8α Mr 2 − 2α Q 2 term outside the bracket in (31), which is of the same form as (29). Therefore in the region r 6 s + 8α Mr 2 s − 2α Q 2 = 0, there is a curvature singularity; this branch singularity r = r s (v) is indeed the curvature singularity of the spacetime. It can be observed that in this case the Kretschmann invariant diverges at the branch r = r s (v) (see Figs. 1 and 2) as K ≈ r −24 which is significantly more than the uncharged case (K ≈ r −16 ), and even more so than the general relativistic charged Vaidya analogue (K ≈ r −8 ). This highlights the role of the curvature corrections of EGB gravity in tandem with the charge component Q of the electromagnetic field. There is no physical spacetime in the region 0 < r s . In the limiting case of Q = 0, the scalar (31) reduces to that of the uncharged metric studied in [40,41].  Figure 1 depicts the behaviour of the Kretschmann scalars for the charged Vaidya spacetime and the charged Boulware-Deser spacetime. The charged Vaidya metric diverges strongly near the singularity r = 0 in a steadily decreasing manner with increasing r . However it is apparent that the divergence of the charged Boulware-Deser spacetime is negative nearer to the branch singularity r = r s ≈ 0.258, and then increases as one moves along the radial coordinate. This is an interesting feature since at some finite radial coordinate, the quantity (31) vanishes entirely before increasing steadily with increasing r . Figure 2 shows the behaviour of the Kretschmann invariants for the uncharged and charged Boulware-Deser metrics. We note that for the uncharged Boulware-Deser spacetime, the divergence pattern is similar to the charged Vaidya spacetime. Near the central singularity, the Kretschmann scalar diverges, before steadily decreasing with increasing r . We can clearly see in Figs. 1 and 2 that the charged Boulware-Deser spacetime diverges rapidly near the branch singularity located at r = r s ≈ 0.258.
The apparent horizon forms when which, after simplification, turns out to be a quartic polynomial in r H . It yields four solutions, namely Again, only the positive parts of (33a) and (33b) are applicable so we therefore have Depending on the value of the discriminant three cases arise. 3. IfΔ > 0, then there will exist two real physical apparent horizons for the charged radiating black hole. The solution (34a) will be the outer horizon and (34b) will represent the inner horizon. To an external observer, only the outer horizon will be visible. This is in line with the conclusions drawn in [59].
We now make the following observations: - studied in [41].
Similarly, the mass function from setting (32) to zero can be written as where the horizon radius must satisfy r H > 0. When Q = 0, which further iterates that there is a branch singularity at r = r s for all M and nonzero Q.
In order to quantify this analysis, we assume α = 2, M = 15 and Q = 2 and calculate approximate values of the branch r s , as well as the inner and outer horizons r 1 and r 0 respectively. This gives the three solutions ≈ 0.258196, and it is clear that for these values, the branch singularity is covered by the two horizons within the domain for r which is 0 < r s < r 1 ≤ r 0 < ∞. Table 1 depicts the different scenarios for the formation of singularities and horizons for differing mass values. We have used the values α = Q = 2 throughout; for these values the two critical solutions for   , we also note that for 0 < M < 2 and 2 < M < 6, no horizon will form since the discriminantΔ < 0. When M = 6, we have r 0 = r 1 = 1 and so there is one horizon which forms upon collapse. For all M > 6, two horizons will always form covering the singularity r s . It is interesting to note that as the values of the mass function increase, the value of the branch singularity r s appears to decrease, however this value will never terminate for Q > 0. Figure 3 depicts the possible collapse scenario, from an initially flat spacetime, for charged null matter. We have radiating charged null matter falling into a black hole. The two apparent horizons, of which only the outer is visible to an observer, form at v = r s (0) = 0 unlike in the uncharged scenario, and encompass trapped surfaces in a dense region 0 < v < V 0 . The Gauss-Bonnet constant α does not delay the formation of the two apparent horizons, unlike the single horizon of the uncharged case. The presence of the charge contribution 2α Q(v) 2 r 6 in the metric (25) indicates that for a specific value of r (= r s ) (in tandem with Q), the metric becomes complex and there is no physical spacetime below this value. This contributes to the formation of a curvature singularity (a branch singularity) which separates the unphysical part of the metric from the rest of the collapsing spacetime. This singularity forms at v = r s (0) = 0 and extends into the future, however, covered by the two horizons. The electromagnetic field contribution is significant in the late stages of collapse in EGB gravity. This is unlike the general relativity case where the electromagnetic field contribution appears to be minimal. At a time v = V 0 , the apparent horizons match smoothly to the two event horizons, of which the outer is visible to the external universe. The charged Boulware-Deser exterior is separated from the black hole by this outer horizon. Within the charged Boulware-Deser black hole resides a strong central curvature singularity separating the collapsing matter from the unphysical metric. Incidentally, the collapse scenario remains the same for all N ≥ 6, so gravitational collapse in EGBM gravity is similar for all dimensions greater than five (see Sect. 6).
The collapse scenario differs when Q = 0, as analysed by [40,41]. The uncharged Boulware-Deser spacetime is, first of all, metric regular in five dimensions. The spacetime itself is singular due to the diverging diffeomorphism invariants [41]. Therefore collapse terminates with a sufficiently weak central conical singularity. This conical singularity exists for a period of time where it is naked, and this is due to the Gauss-Bonnet coupling constant α which delays the formation of the Only the outer horizon r = r 0 is visible to an external observer at infinity. In the case of the discriminantΔ = 0, the inner and outer horizons coincide; there will be a single extremal horizon separating the singularity from the external universe, and so this diagram remains the same. It should be noted that this diagram also remains unchanged for all N ≥ 6 apparent horizon. There are also a separate class of curvature singularities in EGB gravity. In the case of uncharged null matter, (25) becomes Upon inspection of the above metric, a zero within the squareroot implies a curvature singularity. Two branches of solutions meet here and so the singularity is also a branch singularity 6 at some finite radius r = r s , however unlike in the charged situation, this can be resolved. Letting the square (37) and solving for M yields and should r s > 0 hold, the domain of r becomes 0 < r s < r < ∞. We note that a branch singularity only forms for the case M < 0 which is unphysical in itself. Therefore for any positive choice of the mass function M(v), collapse terminates with the above mentioned central conical singularity (see Fig. 4). The collapse situation differs in dimensions N ≥ 6, in that collapse will terminate with a strong curva- 6 Lovelock gravity, in general admits more than one class of solutions, therefore a branch singularity at a finite radius r = r s > 0 is common.  37) and (25) ture singularity since the Boulware-Deser spacetime is metric singular for all dimensions greater than five. However, the branch singularity case mentioned above, i.e. M < 0, holds for all dimensions greater than five. These collapse scenarios in EGB gravity are fundamentally different to the higher dimensional Vaidya collapse in general relativity, highlighting the effect of the curvature corrections indicative of the Lovelock theory. In five dimensions, the temperature at the horizon is calculated as where we have f (v, r ) from (25) at the horizon. We note the presence of the surface area term A 3 , due to the horizon surrounding the 3-sphere in the five dimensional collapsing geometry. The surface gravity is equivalent to the definition Note that the temperature at the horizon r = r H is equal to the Hawking temperature 7 so one can write where we have used (40). The surface gravity of the charged metric (25) is given by where we note the charge contribution Q. When the charge vanishes, the above reduces to 7 In arbitrary dimensions, the Hawking temperature is written as where we note the surface area term for the higher dimensional sphere.
In four dimensions this reduces to the well known result T H = K 2π .
The Hawking temperature at the horizon is given by The above quantity gives a measure of the black body radiation emanating from a black hole due to quantum effects near the horizon [75]. In Fig. 5 the Hawking temperatures for all three spacetimes are represented. We again have used the values α = 2, m(v) = M(v) = 15 and q(v) = Q(v) = 2. Now since the Hawking temperature is measured at the horizon, we consider each spacetime separately.
1. Charged Vaidya: Using the numerical values above, the outer horizon (B.13a) evaluates to It can clearly be seen that the Hawking temperature is positive at this value and then steadily decreases with increasing r . 2. Uncharged Boulware-Deser: The value of the apparent horizon r H = √ M − 2α evaluates to and it can be seen that the Hawking temperature corresponding to this value is also positive and slowly decreasing as expected with increasing r . 3. Charged Boulware-Deser: The value corresponding to the outer horizon (34a) is given by Again, we observe that the Hawking temperature at this value is positive and expectedly decreasing as one moves away from the outer horizon to infinity.
We note that below the respective horizons, these temperatures diverge rapidly when nearing the central singularities.
In Table 2 we highlight the differences between the three spacetimes studied in this paper, with regards to the end state of gravitational collapse.
Another point to emphasise is the effect of the electromagnetic field on the gravitational collapse dynamics in EGB gravity compared to conventional Einstein gravity. In general relativity, the effect of the charge contribution appears to be minimal since the collapse dynamics of the Vaidya spacetime are not too different from its charged counterpart. Collapse for both spacetimes ceases with the formation of a strong curvature singularity which may or may not be naked [12,13]. If it is not naked, it is covered within the confines of a black hole. Radiation collapse is markedly different in EGB gravity. Firstly, the gravitational collapse of the uncharged Boulware-Deser spacetime is critically affected by the presence of the Gauss-Bonnet constant α. The metric is regular and well defined for all r , however the spacetime itself is singular due to its diverging diffeomorphism invariants. This implies that collapse terminates with the formation of an extended, weak conical singularity at the centre that is initially naked if M(v) ≥ 0; the Gauss-Bonnet constant α delays its formation. Eventually the apparent horizon forms and matches smoothly with the event horizon forming a single trapping horizon that contains the trapped surfaces within a black hole with a centrally weak conical singularity. Secondly, the charged Boulware-Deser metric is dynamically very different to its uncharged counterpart. The metric is no longer regular or well defined. The presence of the charge contribution changes the type of singularity encountered upon the cessation of collapse; it ceases with the formation of a curvature singularity which is a branch singularity forming at v = r s (0) = 0 separating the spacetime from an unphysical region which is the result of a complex metric. The horizons can also form at v = r s (0) = 0 unlike with the uncharged case. The Gauss-Bonnet constant appears not to affect the collapse dynamics as profoundly as for the uncharged case. Only the outer horizon is visible to an external observer at infinity and covers all trapped surfaces within the confines of a black hole. A further point to highlight is the role of dimensions in this analysis. In EGB gravity, it is well known that N = 5 and N = 6 are the physically important dimensions. We note that the Boulware-Deser spacetime in N ≥ 6 dimensions, whether charged or not, will collapse to a strong, central curvature singularity for all M(v) ≥ 0 (or a branchlike singularity for the charged case). The uncharged case mirrors general relativistic collapse. It is only in five dimensions, where a marked difference can be seen in the collapse dynamics of the uncharged Boulware-Deser spacetime and its counterpart with the electromagnetic field contribution. This highlights that N = 5 is the critical dimension in EGB gravity, and that the electromagnetic field effects are more profound in EGB gravity, than in general relativity.

Higher dimensions
Gravitational contraction in higher dimensions is vital with regards to modified or higher order theories of gravity. Higher dimensional collapse has been extensively studied, for example see [76][77][78] in general relativity and [42,79,80] in modified gravity. As iterated earlier, the collapse process of the Boulware-Deser spacetime in dimensions greater than five differs from the five dimensional case. For all N > 5 the collapse ceases with the formation of a strong curvature singularity, whether the spacetime has a charge contribution or not. If the charge component is there, the result of collapse is a branch singularity separating the physical spacetime from the region where the metric function is complex. The Ndimensional generalised Boulware-Deser metric is given by where we have the (N − 2)-sphere In the above where we have setα = α(N − 3)(N − 4) for neatness. In higher dimensions the choice of the mass function (23) becomes where we note the explicit appearance of Einstein's gravitational constant and the surface area of the (N − 2)-sphere. The electromagnetic stress energy tensor (10) is built into this definition (47) above. The metric (46) then becomes with We note the following: -When N = 5, the above expression (48) reduces to (25) found by Wiltshire [59]. -When N = 5 and Q(v) = 0, the expression (48) reduces to the uncharged metric (37).
We also note that there exists a maximal charge contribution to the higher dimensional metric (48) for which it remains real. This contribution depends on the dimension N , surface area A N −2 , Einstein's gravitational constant κ N and the Gauss-Bonnet constantα = α(N − 3)(N − 4). It is given by for all dimensions N ≥ 5 and α = 0. This can also be written in the form where r = r s (v) > 0 denotes a branch singularity separating the physical spacetime from the complex metric, analogously to the five dimensional case. When N = 5 this reduces to (29). Therefore, the collapse process of null radiating matter within an electromagnetic field is similar for all dimensions N ≥ 5. The dimension appears to focus the effects of the higher order curvature indicative of the theory. If we consider the general N -dimensional metric (46), the formation of the apparent horizon will occur when This simplifies tô It is important to note that whenever the order of the polynomial is odd, there will always be at least one real root. When the order is even, this need not be the case; all roots could possibly be complex. It is indeed possible to incorporate the electromagnetic field into the generalised mass function M(v, r ) in (52). This giveŝ Equivalently, we can solve the above for the mass function to get When N = 5 in (54), this reduces to (35), and so the form of M H (r H ) is qualitatively the same for all dimensions N ≥ 5.
Therefore, as previously mentioned, the gravitational collapse of charged null matter in higher dimensions is not different to the five dimensional case in EGBM gravity. This is markedly different to the uncharged scenario where in five dimensions collapse ceases with a weak conical singularity (for positive mass) and for dimensions six and higher, collapse terminates with a strong curvature singularity. We will present different cases for the polynomials in different dimensions in Table 3.
It is interesting to note that when the electromagnetic field is present, we see in Table 3 that the polynomial equations in r are always of even order. When Q = 0, the order of the polynomial equations alters from even to odd, with odd dimensions corresponding to even ordered polynomials and even dimensions corresponding to odd ordered equations. We deduce from this that for vanishing charge, and for all even dimensions, it is guaranteed that the resulting polynomial equation for the horizon formation will yield at least one real root. If such a root is positive, this means that an apparent horizon will always form to cover the central singularity. We can state the following in a theorem: Theorem 1 Consider an uncharged collapsing, inhomogeneous and radiating Boulware-Deser spacetime in N > 5 dimensions, from a regular epoch with a generalised mass function M(v, r ) that obeys all physically reasonable energy conditions and is at least C 2 in the entire spacetime manifold. Whenever the dimension of spacetime is even (i.e. N = {2k, k ≥ 3, k ∈ Z}), and the real solution to the horizon formation equation is positive, the final fate of collapse will be a strong curvature singularity that will eventually be covered by an event horizon.

Conclusion
In this paper we studied the continual gravitational contraction of a charged null radiation shell in five dimensional EGBM gravity. The final outcome of collapse was the formation of a curvature branch singularity covered by two trapping horizons, of which the outer is visible to the external universe. The charge contribution yielded a situation where the metric is complex in a certain region, separated by this branch singularity. This is very different to uncharged radiating collapse where, for a positive mass, the final end state was an extended weak conical singularity which remains initially naked before succumbing to an apparent horizon at a later time. The electromagnetic field contribution critically alters the dynamics of the collapse process of the radiating Boulware-Deser spacetime. In general relativity, the charge contribution to the collapse process appears to be minimal when compared to EGB gravity. The collapse of the Vaidya spacetime and its charged counterpart both cease with the for-  16 15 π 3 N = 9 κ 9 = 7 18 π 4 r 12 + 30αr 10 − 1 3 Mr 6 + Q 2 36 = 0 mation of a strong curvature singularity which may or may not be naked.
In EGB gravity the electromagnetic field contribution alters the type of singularity formed, post collapse. When Q = 0 and M > 0, the Boulware-Deser spacetime collapses to a weak and conical singularity which is initially kept naked by the presence of the positive Gauss-Bonnet coupling constant α; the curvature corrections directly affect the collapse and post collapse scenarios. While naked, families of trajectories can escape to infinity from the singularity. Eventually an apparent horizon forms which encloses all trapped surfaces and this horizon eventually matches smoothly to an event horizon when the collapse process is complete. The event horizon separates the trapped surfaces from the external five dimensional Boulware-Deser vacuum. When Q = 0, gravitational collapse abates with the formation of a branch curvature singularity covered by two trapping horizons, the outer of which is visible to any observer at infinity. The branch singularity r = r s (v) results from the fact that the metric is complex whenever the inequality is violated. This marked difference only occurs when the dimension of the spacetimes is N = 5. This is due to the fact that the five dimensional uncharged Boulware-Deser metric (12) is regular and everywhere well defined for all values of r . This is not the case for the charged analogue (25). This is fundamentally different to the general relativistic case of radiation contraction. We note also that for the charged Boulware-Deser metric in all dimensions N ≥ 5, there exists a maximal contribution due to the electromagnetic field, for which the metric on the spacetime remains real. Otherwise it becomes complex. This is critically differ-ent to the charged Vaidya analogue in Einstein gravity. This highlights the effects of the higher order corrections as well as those corrections in tandem with the electromagnetic field. For all dimensions N ≥ 6 the collapse of the charged shells of radiation mirror five dimensional collapse. It is only in the uncharged scenario where the collapse dynamics differ from N = 5 and N ≥ 6. This reinforces the notion that the role of dimension in gravitational collapse cannot be underestimated, especially when the gravity theory has modified curvature components and an electromagnetic field affecting the dynamics of the model. What remains to be done is a singularity analysis of the five dimensional charged collapse model. Are naked singularities possible? The conditions for a locally naked singularity; the determining of the existence of outgoing nonspacelike geodesics, need to be analysed. This will allow one to ascertain whether cosmic censorship is violated or not. This will be the subject of future work.
It will also be an interesting endeavour to look for possible observational signatures for compactified higher dimensions, and whether the underlying theory of gravity is indeed EGB gravity. There are numerous qualitative differences between four and higher dimensional gravitational collapse of null or charged null dust in Einstein gravity and these differences escalate in the presence of Gauss-Bonnet term. For example, in five dimensional neutral null dust collapse in EGB gravity, there is a possibility of forming a weak curvature extended naked conical singularity at the centre, that may give rise to strong quantum gravity signatures, that are not shielded from the external observers by the horizon [41]. Furthermore, as we have already seen, after the black hole forms in the radiating EGB collapse, the surface gravity and hence the Hawking temperature of the black holes crucially depend on the EGB coupling constant α. In four dimensional general relativity, we know that the Hawking temperature of solar mass black holes is significantly less than the cosmic microwave background radiation (CMBR) temperature, and hence these are not readily detectable. However in higher dimensional EGB scenarios, there exists a set of non-zero measures of parameter values, that can actually make a solar mass black hole "hotter" than the CMBR, and hence should be, in principle, detectable via a measure of their Hawking radiation.
where subscripts denote differentiation. In the above, μ is the energy density of the null dust, ρ is the energy density of the null string, and P is the null string pressure. The system (B.5) comes about as a result of the generalised two-fluid energy momentum tensor T ab = μl a l b + (ρ + P)(l a n b + l b n a ) + Pg ab , (B.6) where the vectors l a = δ 0 a and n a = 1 2 f (v, r )δ 0 a + δ 1 a are null. For the case of the charged Vaidya (or Vaidya-Bonner) spacetime, the mass function is given by Therefore the charged Vaidya spacetime in five dimensions is then q(v) 2 4r 4 dv 2 + 2dvdr +r 2 (dθ 2 + sin 2 θ dφ 2 + sin 2 θ sin 2 φdψ 2 ). (B.8) We note the presence of the charge contribution q in (B.8).
When q = 0, the above reduces to the Vaidya metric. We can also write (B.5) as the system where the role of the charge contribution q 2 is more clear. When q = 0 we obtain the single field equation of the five dimensional Vaidya solution. This is consistent with [20]. The energy conditions (17)- (20) for (B.8) will be satisfied if we have that and there is a violation of the first condition near r = 0. It was shown in [70] that the energy conditions can still be satisfied in such a scenario if the equation of motion includes a Lorentz force term.