This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Insights and guidelines on the Cauchy horizon theorems

Xiao Yan Chewa111[email protected] and Dong-han Yeomb,c,d222[email protected] aDepartment of Physics, School of Science, Jiangsu University of Science and Technology, Zhenjiang 212100, Jiangsu Province, People’s Republic of China
bDepartment of Physics Education, Pusan National University, Busan 46241, Republic of Korea
cResearch Center for Dielectric and Advanced Matter Physics, Pusan National University, Busan 46241, Republic of Korea
dLeung Center for Cosmology and Particle Astrophysics, National Taiwan University, Taipei 10617, Taiwan
Abstract

Recently there has been progress to resolve the issue regarding the non-existence of the Cauchy horizon inside the static, charged, and spherically symmetric black holes. However, when we generically extend the black holes’ spacetime, they are not just static but can be dynamical, thus the interior of black holes does not remain the same as the static case when we take into account the dynamical evolution of black holes. Hence, the properties of the Cauchy horizon could behave differently in the dynamical case. Then, our aim in this paper is to provide a few constructive insights and guidelines regarding this issue by revisiting a few examples of the gravitational collapse of spherically symmetric charged black holes using the double-null formalism. Our numerical results demonstrate that the inside of the outer horizon is no longer static even in late time, and the inner apparent horizon exists but is not regular. The inner apparent horizon can be distinguished clearly from the Cauchy horizon. The spherical symmetric property of black holes allows the inner horizon to be defined in two directions, i.e., the differentiation of the areal radius vanishes along either the out-going or the in-going null direction. Moreover, the Cauchy horizon can be generated from a singularity. Still, the notion of the singularity can be subtle where it can have a vanishing or non-vanishing areal radius; the corresponding curvature quantities could be finite or diverge, although the curvatures can be greater than the Planck scale. Finally, we show some examples that the “hair” which is associated with the matter field on the inner horizon is not important to determine the existence of the Cauchy horizon; rather, the hair on the outer horizon might play an important role on the Cauchy horizon. Therefore, the dynamic properties of the interior of charged black holes could shed light for us to understand deeply about the Cauchy horizon for the extensions of no-Cauchy-horizon theorems.

I Introduction

The investigation of the interior of a black hole is a very interesting and important issue in general relativity. According to the singularity theorem Hawking:1970zqf , under very reasonable and natural assumptions, the singularity must be formed as a result of a gravitational collapse. Then, one may wonder about the properties of the singularities, here we have the Reissner-Nordstrom black hole Reissner ; Nordstrom as a typical example which contains a time-like singularity and two horizons which are the outer horizon and the inner horizon. The outer horizon corresponds to the usual black hole horizon and the inner horizon corresponds to the Cauchy horizon. The Cauchy horizon preserves the predictability in general relativity where a system can be evolved with time from its given initial data which is imposed on it. When a time-like singularity exists and an observer can see an event emanating from it, then the so-called cosmic censorship conjecture is violated Penrose:1969pc .

Will the cosmic censorship conjecture be true in generic situations? If the electric charge of the Reissner-Nordstrom black hole is greater than its mass, its time-like singularity can become a naked singularity. However, it is not quite easy to construct such an overcharged black hole using gravitational collapses Hwang:2010im ; Hansen:2013vha . This justifies the weak cosmic censorship where there must be no observer who can see the effects from the time-like singularity; in other words, no observer can cross the Cauchy horizon which can predict the evolution of a system from the imposed initial conditions on it in the strong cosmic censorship. The violation of the strong cosmic censorship also implies the breakdown of predictability in general relativity. Will the strong cosmic censorship conjecture be true, too?

Although we do not have a definite answer yet, we have some evidence to support that the strong version of the cosmic censorship conjecture is still reasonable. In general, there should exist an inner horizon before an observer can reach the time-like singularity. Here, the inner horizon is generally unstable, because the inner horizon corresponds to an infinite blue-shift instability Simpson:1973ua , while the outer horizon corresponds to an infinite red-shift instability. Relating these instabilities, when there is a matter fluctuation and if an observer measures the corresponding fluctuation near the inner horizon, the observed energy density will be exponentially increased. This effect is known as the mass inflation Poisson:1990eh . As long as mass inflation exists, we cannot trust the interior structure of the Reissner-Nordstrom spacetime anymore, and we have to rely on the full dynamical simulations Bonanno:1994ma .

In this context, recently there are interesting papers in the literature that discuss the existence or non-existence of the Cauchy horizons An:2021plu ; Cai:2020wrp ; Hartnoll:2020fhc ; Dias:2021afz ; Devecioglu:2021xug ; Grandi:2021ajl ; Yang:2021civ ; Henneaux:2022ijt ; Luk:2022rgs ; Devecioglu:2023hmn . In the literature, some authors proved the non-existence of the Cauchy horizons by considering some matter fields which support the black holes, i.e., the existence of hairs which refers to the extra global charges (primary hair) associated with the matter fields or only the matter fields themselves (secondary hair). However, to prove a theorem on the non-existence of the Cauchy horizons, one needs to assume some conditions to define the inner horizons or Cauchy horizons. To generalize the corresponding mathematical theorems, the underlying assumptions must be true, but one may wonder whether these assumptions will be still true in fully dynamical situations.

The aim of this paper is not to criticize other literature but to provide useful insight to learn some interesting properties of the Cauchy horizons for the black holes in the dynamical evolution. Indeed, the Cauchy horizon and the inner apparent horizon are distinguishable in dynamic cases, and the situations become more complicated if we take into account the quantum effects. We will report the detailed results that we can learn from numerical computations. We also wish our work can shed light on the development of some advanced theorems about the Cauchy horizons of black holes in the future.

This paper is organized as follows. In Sec. II, we discuss some preliminary topics about the interior structure of the Reissner-Nordstrom black hole, the instability of inner horizons, and the recent progress with several literature that reformulates the no-Cauchy-horizon theorems. In Sec. III, we mainly demonstrate the dynamical properties of the inner horizon by revisiting several models of charged black holes in the gravitational collapse using the double-null formalism which can study the full dynamical process of black holes that begins from their formation and then ends with their evaporation. Finally, in Sec. IV, we provide a few implications and guidelines that are related to the Cauchy horizon based on our numerical results and comment on possible future research directions.

II Preliminaries of inner horizons

In this section, we maximally extend the Reissner-Nordstrom black hole which is a typical example of a static and charged black hole to understand its basic interior structure, particularly the Cauchy horizon inside the outer horizon. We then discuss the instability of the inner horizon. Moreover, we discuss the recent progress that reformulates the non-existence of the Cauchy horizon by providing several references.

Refer to caption
Figure 1: The Penrose diagram of a static charged black hole. There exists a time-like singularity (r=0r=0). The inner and outer apparent horizons correspond to the Cauchy and event horizons. r,v=0r_{,v}=0 horizons (red color) are always parallel to the out-going null direction, while r,u=0r_{,u}=0 (blue color) horizons are always parallel to the in-going null direction.

II.1 Static charged black holes

Here we begin with the Reissner-Nordstrom black hole which is a static, stationary, and charged black hole, with the given metric as shown in the following,

ds2=f(r)dt2+1f(r)dr2+r2dΩ22,\displaystyle ds^{2}=-f(r)dt^{2}+\frac{1}{f(r)}dr^{2}+r^{2}d\Omega_{2}^{2}, (1)

where

f(r)=12Mr+Q2r2,\displaystyle f(r)=1-\frac{2M}{r}+\frac{Q^{2}}{r^{2}}, (2)

where MM is the ADM mass and QQ is the asymptotically defined electric charge Reissner ; Nordstrom . Since f(r)=0f(r)=0 can yield two solutions, one can notice that, if M>QM>Q, this metric has two apparent horizons which can be denoted as r+r_{+} and rr_{-} with r+>rr_{+}>r_{-}, where r+r_{+} is the outer horizon and rr_{-} is the inner horizon.

Let us define that uu denotes the in-going null direction and vv denotes the out-going null direction. After we maximally extend the spacetime in Kruskal-Szekeres coordinates dInverno:1992gxs , i.e.,

ds2=4r+4(rr)1+r2/r+2r2(r+r)2err+r+2rdudv+r2dΩ2\displaystyle ds^{2}=-\frac{4r_{+}^{4}\left(r-r_{-}\right)^{1+r_{-}^{2}/r_{+}^{2}}}{r^{2}\left(r_{+}-r_{-}\right)^{2}}e^{\frac{r_{-}-r_{+}}{r_{+}^{2}}r}dudv+r^{2}d\Omega^{2} (3)

with r=r(u,v)r=r(u,v), one can notice that there are two directions of the horizons, where r,v=0r_{,v}=0 horizons (red color) are always parallel to the out-going null direction, while r,u=0r_{,u}=0 horizons (blue color) are always parallel to the in-going null direction as depicted in Fig. 1. In addition, there exists a time-like singularity at r=0r=0; the inner horizon is the boundary where we cannot decide the geometry from the past data, and hence, this horizon is known as the Cauchy horizon Hawking:1973uf . The outer horizon is the boundary between the interior and the asymptotic infinity, and hence, this horizon is known as the event horizon. So, in the static and charged black hole case, the inner apparent horizon is the Cauchy horizon and the outer apparent horizon is the event horizon.

Here, we list down several common senses of the charged black holes:

  • r,v=0r_{,v}=0 apparent horizon is parallel to the out-going null direction; r,u=0r_{,u}=0 apparent horizon is parallel to the in-going null direction.

  • The r=0r=0 singularity is time-like.

  • The inner apparent horizon is the Cauchy horizon; the outer apparent horizon is the event horizon.

However, these common senses are valid only if the metric is static. In the dynamic cases, all three assertions are not valid in general. On the other hand, if one applies such common sense to some proofs about the non-existence of inner horizons, the corresponding results might not be rigorous or even wrong because, in a realistic situation, the metric should be dynamic.

II.2 Instability of inner horizons

One of the most important issues about the inner horizon is the instability of inner horizons. When there exists an energy flux along the in-going null direction, we can define an observer who is moving along the out-going null direction and measures the energy flux. Near the inner horizon, the observed energy flux is approximately

ρe+κv,\displaystyle\rho\propto e^{+\kappa v}, (4)

where κ\kappa is the surface gravity of the inner horizon Poisson:1990eh . Therefore, the observed energy must be exponentially amplified near the inner horizon. This effect is known as the mass inflation.

Theoretically, we can have a static, stationary, and charged black hole solution that is independent of time. However, in reality, there is no such an eternal object in the universe; every black hole must be generated from gravitational collapses. Therefore, energy fluctuations necessarily exist inside a charged black hole. Thus, the consequence of mass inflation is that, in realistic black holes that are generated from gravitational collapses, there must be a strong back-reaction due to the exponentially increased energy flux Bonanno:1994ma .

Hence, it is unlikely to precisely study the interior structures of a charged black hole from the perturbative analysis on the background of static solution. The only possible way to precisely understand the interior structure is to numerically evolve the corresponding equations of motion. In this paper, we apply the double-null formalism to study the formation of an inner horizon in the gravitational collapse of charged black holes.

II.3 No-Cauchy-horizon theorem

In this context, there have been some mathematical investigations and proofs regarding the inner horizon of black holes. In the literature, these investigations are mainly to formulate the no-Cauchy-horizon theorems. For example, we specify several mathematical assertions from the literature as follows:

  • Ref. An:2021plu has shown numerically the Cauchy horizon does not exist in the interior of the scalar-charged hairy black holes with spherical and planar symmetries in the Born-Infield theory.

  • Ref. Cai:2020wrp has shown numerically the Cauchy horizon does not exist in the interior of the scalar charged hairy black holes with spherical and planar symmetries in Einstein-Maxwell-scalar theory of a complex scalar field with an arbitrary form of the scalar potential.

  • Ref. Devecioglu:2021xug has studied the non-existence of the Cauchy horizon for some specific types of black holes analytically in the Einstein-Maxwell-Horndeski theories.

  • Ref. Grandi:2021ajl has studied numerically the non-existence of Cauchy horizon for scalar charged hairy black holes in the 5-dimensions of Einstein-Gauss-Bonnet theory with a massive complex scalar field.

  • Ref. Yang:2021civ has studied the possibility of putting a constraint on the number of horizons for the static black holes which can only possess at most one non-degenerated inner Killing horizon inside the event horizon when they satisfy the strong or null energy condition.

  • Ref. Devecioglu:2023hmn has studied numerically the non-existence of Cauchy horizon for charged black holes with planar configuration in the Einstein-Maxwell-Gauss-Bonnet-scalar theory.

Here we list down some common features in the literature as follows.

  • 1. Some authors assumed that the black hole geometry will follow a regular and static solution. If a gravitational collapse occurs and then after a long time has sufficiently passed, a black hole solution will approach a static solution, but this is in the perspective of the outside observer. One can ask whether this assumption is still true for an inside observer.

  • 2. Some authors assumed hairs to prove the theorems for some specific cases, i.e., the dynamics of matter fields near the inner horizon. If any inconsistency can arise from the corresponding dynamics, then they can prove the non-existence of the Cauchy horizon. However, the no-hair theorem has not been established near the inner apparent horizon, and thus there is no systematic way to evade the no-hair theorem inside the outer horizon; hence, it would be risky to assume the existence of such a conserved quantity associated with the symmetry of the theory for the inner horizon.

In this paper, we will show some examples that satisfy the following conditions. (1) The inner horizon is formed and approaches a (weak) curvature singularity. Hence, the interior geometry is no longer static or regular. (2) Near the inner horizon, there can be non-trivial field dynamics. This is not surprising because the interior geometry is no longer static. In the next section, we will discuss more issues that go beyond common sense.

III Lessons from numerical investigations

In this section, we revisit several models of charged black holes in the gravitational collapse using the double-null simulation by demonstrating the properties of the inner horizon in the gravitational collapse. For the sake of readers, we briefly introduce the double-null formalism as the first step. Then we also briefly define the properties on the inner horizon and singularity in the double-null coordinate since they can behave differently in the gravitational collapse. Lastly, we will present and discuss the properties of the inner horizon in the gravitational collapse.

III.1 Brief summary of the double-null formalism

Let us consider the most generic double-null coordinates with spherical symmetry:

ds2=α2(u,v)dudv+r2(u,v)dΩ22.\displaystyle ds^{2}=-\alpha^{2}(u,v)dudv+r^{2}(u,v)d\Omega_{2}^{2}. (5)

In the the double-null formalism, one can obtain the equations of motion and solve them numerically based on this metric with appropriate boundary conditions (see Appendix C) Nakonieczna:2018tih . So far, plenty of numerical simulations have been carried out to investigate the dynamics of the formation of inner horizons during the gravitational collapse of the charged black holes using the double-null formalism Hod:1998gy ; Hwang:2011mn ; Sorkin:2001hf ; Hong:2008mw ; Sorkin:2001hf ; Hansen:2014rua ; Hansen:2015dxa ; Nakonieczna:2015umf ; Nakonieczna:2016iof ; Chen:2015nma . In this section, we summarize some important results that are relevant to the issue of the inner horizon by revisiting some of our previous works.

  • The authors investigated the properties of the Cauchy horizon inside the formation of a charged black hole during a gravitational collapse caused by an input pulse consisting of a complex scalar field in the Einstein-Maxwell theory Hod:1998gy ; Hwang:2011mn . The numerical results show that the inner horizon is formed after the gravitational collapse. If we turn off Hawking radiation, the inner horizon becomes a Cauchy horizon.

  • The authors extended the above investigation to study the semi-classical description for the emission of Hawking radiation from a charged black hole Sorkin:2001hf ; Hong:2008mw . The inclusion of Hawking radiation allows that the inner horizon and the Cauchy horizon can be distinguished.

  • The authors generalized the above two investigations from the Einstein gravity to the Brans-Dicke theory that can cover a diverse range of theories, including string-inspired models Hansen:2014rua ; Hansen:2015dxa ; Nakonieczna:2015umf ; Nakonieczna:2016iof , f(R)f(R) gravity Hwang:2011kg ; Chen:2015nma , etc. These results show that the inner and outer horizons can be varied dynamically due to the Brans-Dicke field. Also, we can see the existence or non-existence of inner horizons relating couplings between the Maxwell field and the Brans-Dicke field.

Before we proceed to present the above results, we need to define clearly the inner horizon and singularity in the next subsection, since their appearance is highly sophisticated in the double-null simulation.

III.2 Subtleties to define inner horizon and singularity

The Reissner-Nordstrom black hole in the double-null coordinate contains two kinds of inner horizons: r,v=0r_{,v}=0 and r,u=0r_{,u}=0, where the former is parallel to a constant uu line and another is parallel to a constant vv line. The two horizons appear with the same areal radius, and they are well coincide with the causally defined Cauchy horizon.

However, we need to distinguish three kinds of horizons for the dynamical cases:

  • H1. r,v=0r_{,v}=0 inner apparent horizon.

  • H2. r,u=0r_{,u}=0 inner apparent horizon.

  • H3. Cauchy horizon (boundary where one can causally determine from the past data).

To define the third type of horizon (H3) precisely, one needs to take into account the subtleties to define the singularity and its properties. Thus, let us grasp several basic concepts about singularity as follows.

  • Type 1: it lies at r=0r=0 and the corresponding curvature quantities diverge; this kind of singularity is in the classical sense.

  • Type 2: it lies in r>0r>0 and the corresponding curvature quantities diverge; this kind of singularity is also in the classical sense which is so-called a weak singularity Ori:1991zz .

  • Type 3: the radius of singularity is non-vanishing but at the Planck scale rPr\simeq\ell_{\mathrm{P}} and the corresponding curvature quantities are also the Planck scale.

  • Type 4: the radius of singularity is greater than the Planck scale r>Pr>\ell_{\mathrm{P}} but the curvature quantities are at the Planck scale. This kind of singularity is known as the quantum version of the weak singularity.

Therefore, we can classify that Type 1 and 2 are classical singularities, while Type 3 and 4 are defined when there exists the Planck constant \hbar for the quantum regime. Type 1, 2, and 3 will dynamically generate the Cauchy horizon because one cannot extend computations inside the singularity; however, the property of Type 4 is less known, because one may extend general relativistic computations inside the singularity. The only problem is that such an extension is inconsistent in terms of quantum gravity except for some extreme cases (e.g., introducing the large NN-rescaling Yeom:2009zp ; Kim:2013fv ).

Refer to caption
Refer to caption
Figure 2: Left: Dynamical formation of a charged black hole, where P=0P=0 Hong:2008mw . There exists a Type 1 singularity r=0r=0 (thick black curve) and an outer apparent horizon r,v=0r_{,v}=0 (red curve). As vv goes to infinity, there appears an H1 horizon (r,v=0r_{,v}=0) which is parallel to the ingoing null direction; this will approach a Type 2 curvature singularity (r>0r>0). Here, the grid spacing is 11 for black contours and 0.10.1 for green contours. Right: The Penrose diagram for the dynamical formation of a charged black hole.
Refer to caption
Refer to caption
Figure 3: Left: The evaporating charged black hole, where P=0.1P=0.1 Hong:2008mw . There exists a Type 3 singularity with rPr\simeq\ell_{P} and an outer apparent horizon at r,v=0r_{,v}=0. At the inner apparent horizon, an H1 horizon (r,v=0r_{,v}=0) and an H2 horizon (r,u=0r_{,u}=0) are overlapped. Between the Type 3 singularity and H1/H2 horizon, there appears an H3 Cauchy horizon; we cannot decide beyond the null surface. Beyond the H1/H2 inner apparent horizon, the curvature rapidly increases, and hence, the gray-colored region is the super-Plankian region which goes beyond the scope of the semi-classical description. Therefore, the gray-colored region corresponds to Type 4 singularity (r>Pr>\ell_{P}). Right: The Penrose diagram for an evaporating charged black hole.

III.3 Formation of the inner horizon in Einstein-Maxwell theory

In this subsection, we discuss the formation and evaporation process of charged black holes based on Ref. Hong:2008mw . Hence, we begin with the typical model where a complex scalar field is coupled with a U(1)U(1) gauge field that is given as follows:

S=g𝑑x4[16π(ϕ;μ+ieAμϕ)(ϕ¯;μieAμϕ¯)18πFμνFμν],\displaystyle S=\int\sqrt{-g}dx^{4}\left[\frac{\mathcal{R}}{16\pi}-\left(\phi_{;\mu}+ieA_{\mu}\phi\right)\left(\bar{\phi}^{;\mu}-ieA^{\mu}\bar{\phi}\right)-\frac{1}{8\pi}F_{\mu\nu}F^{\mu\nu}\right], (6)

where \mathcal{R} is the Ricci scalar, ϕ\phi is a complex scalar field, AμA_{\mu} is the U(1)U(1) gauge field, ee is the gauge coupling, and Fμν=Aν;μAμ;νF_{\mu\nu}=A_{\nu;\mu}-A_{\mu;\nu}. In the double-null formalism, one can observe a charged black hole that is formed in this model due to the appearance of the input pulse ϕ(u,v)\phi(u,v) at u=0u=0 to cause a gravitational collapse to occur in the spacetime,

ϕ(0,v)=A4πsin2(vvf)exp(2πivvf),\displaystyle\phi(0,v)=\frac{A}{\sqrt{4\pi}}\sin^{2}\left(\frac{v}{v_{\mathrm{f}}}\right)\exp\left(2\pi i\frac{v}{v_{\mathrm{f}}}\right), (7)

where the pulse is defined in 0vvf0\leq v\leq v_{\mathrm{f}}, A=0.25A=0.25, and e=0.1e=0.1.

To turn on the Hawking radiation, we introduce the renormalized energy-momentum tensor T^μν\langle\hat{T}_{\mu\nu}\rangle in an approximate form Davies:1976ei , where the tensor is proportional to PNP2P\propto N\ell_{P}^{2} (the quantity PP is precisely defined in Appendix A), where NN is the number of scalar field that contribute to the Hawking radiation; this is the natural cutoff if NN fields contribute to Hawking radiation Dvali:2012uq . If P=0P=0, we turn off Hawking radiation; we choose P=0.1P=0.1 for the evaporating case.

Solving the equations of motion derived from Eq. (6) numerically using the double-null formalism (see Appendices A and C for more details), our numerical results demonstrate that a charged black hole is formed and evaporated during the process of gravitational collapse. Regarding this, we address several important remarks from the numerical results in Hong:2008mw :

  • 1. There exist two kinds of singularities, which are Type 1 and Type 2. On the left of Fig. 3, one can see a Type 1 singularity at r=0r=0 (thick black curve) and an outer apparent horizon r,v=0r_{,v}=0 (red curve).

    In the late time limit, i.e., as vv goes to infinity, there appears the H1 horizon (r,v=0r_{,v}=0) which is parallel to the ingoing null direction. This will eventually approach a Type 2 curvature singularity (r>0r>0) due to the mass inflation effect.

    Therefore, the Type 2 singularity coincides with the horizon H1, although it is parallel with constant vv lines. This is somehow counter-intuitive with the common sense which has been mentioned in Sec. II. However, this will help to preserve the strong cosmic censorship conjecture, because no observer can penetrate the Cauchy horizon due to the weak singularity.

  • 2. If we turn on the Hawking radiation as well as the incoming negative energy flux, the situation becomes more complicated. In the left of Fig. 3, there exists the Type 3 singularity rPr\simeq\ell_{P} and an outer apparent horizon r,v=0r_{,v}=0. At the inner apparent horizon, the H1 horizon (r,v=0r_{,v}=0) and the H2 horizon (r,u=0r_{,u}=0) are overlapped. Hence, as an observer penetrates inside the inner apparent horizon, the observer will experience that the areal radius increases, which is analogous to a wormhole structure (see more details in Hwang:2010im ). Due to the existence of the Type 3 singularity, there appears the H3 Cauchy horizon which is the future affected region from the singularity. Beyond the H1/H2 inner apparent horizon, the curvature rapidly increases, approximately me+M2m\sim e^{+M^{2}}, where mm is the Misner-Sharp mass Hong:2008mw ; hence, the gray-colored region represents the super-Plankian region which goes beyond the scope of the semi-classical description. Therefore, the gray-colored region (inside the H1 and H2 inner horizons) corresponds to Type 4 singularity (r>Pr>\ell_{P}).

    To summarize, there appears both Type 3 singularity and Type 4 singularity. Type 3 singularity generates the Cauchy horizon, while it is hidden by the Type 4 singularity. On the other hand, if one can trust the computation beyond the Type 4 singularity, e.g., assuming a large number of matter fields (eM2\sim e^{M^{2}}, see Yeom:2009zp ; Kim:2013fv ), the strong cosmic censorship might be violated even with the semi-classical effects Poisson:1997my . The Type 4 singularity very well coincides with an inner apparent horizon, although there is an ambiguity in defining the exact location of the singularity because it depends on the choice of the cutoff.

Therefore, if we only define the apparent horizon using H1, one may neglect the important and interesting behaviors relating to the inner horizon and the Cauchy horizon. We can conclude that in dynamical situations, H1, H2, and H3 horizons are distinguishable; also, Type 1, Type 2, Type 3, and Type 4 singularities all appear in different situations.

As we have mentioned, Fig. 3 confirms that the interior of a black hole does not remain as the interior of the static solution. The inner horizon becomes the r,v=0r_{,v}=0 horizon and is parallel to the in-going null direction, which is inconsistent with the static solution. Moreover, in the late time (vv\rightarrow\infty) limit, the inner horizon becomes a curvature singularity, and hence, the regularity condition is also not satisfied.

III.4 Do hairs remove the Cauchy horizon?: A case of the Brans-Dicke theory

Some literature heavily relied on the assumption that there may exist hairs associated with the matter fields to establish some theorems about the non-existence of Cauchy horizons. In this subsection, we focus on the cases in which some hairs may help to remove the Cauchy horizon, but surprisingly there exists a counter-example that there is the Cauchy horizon with a scalar hair. Here, a crucial factor to support the counter-example is the existence of scalar hair at the outer horizon, although the non-trivial scalar hair can exist on the inner horizon.

To investigate this issue, we consider the following model with a non-minimally coupled scalar field Hansen:2014rua ; Hansen:2015dxa ; Nakonieczna:2015umf ; Nakonieczna:2016iof :

S=g𝑑x4[116π(ΦωΦΦ;μΦ;μV(Φ))+Φβ(12(ϕ;μ+ieAμϕ)(ϕ¯;μieAμϕ¯)116πFμνFμν)],S=\int\sqrt{-g}dx^{4}\left[\frac{1}{16\pi}\left(\Phi\mathcal{R}-\frac{\omega}{\Phi}\Phi_{;\mu}\Phi^{;\mu}-V(\Phi)\right)+\Phi^{\beta}\left(-\frac{1}{2}\left(\phi_{;\mu}+ieA_{\mu}\phi\right)\left(\bar{\phi}^{;\mu}-ieA^{\mu}\bar{\phi}\right)-\frac{1}{16\pi}F_{\mu\nu}F^{\mu\nu}\right)\right], (8)

where Φ\Phi is the Brans-Dicke field, ω\omega is the Brans-Dicke parameter, V(Φ)V(\Phi) is the potential of the Brans-Dicke field, and β\beta is a constant. Similarly, to observe the formation of the inner horizon for charged black hole during the gravitational collapse using the double-null formalism we send an input pulse of a charged scalar field (0vvf0\leq v\leq v_{\mathrm{f}}) at u=0u=0, where it is given by

ϕ(0,v)=A4πsin4(vvf)exp(2πivvf).\displaystyle\phi(0,v)=\frac{A}{\sqrt{4\pi}}\sin^{4}\left(\frac{v}{v_{\mathrm{f}}}\right)\exp\left(2\pi i\frac{v}{v_{\mathrm{f}}}\right). (9)

For simplicity, we fix e=0.3e=0.3 and β=0\beta=0 or β=1\beta=1.

Refer to caption
Figure 4: The real part (left) and the imaginary part (right) of the complex scalar field when ω=1000\omega=1000 and β=0\beta=0. Hence, the dynamics of the Brans-Dicke field can be negligible and this corresponds to the Einstein gravity Nakonieczna:2015umf ; Nakonieczna:2016iof . Here, we used A=0.25/2A=0.25/\sqrt{2}. The thin white contours are the contours of the fields and the red curves are apparent horizons (r,v=0r_{,v}=0). The blue-colored regions are where the contours are space-like; the red-colored regions are where the contours are time-like. This directly shows that the scalar field has non-trivial dynamics in the Cauchy horizon limit.
Refer to caption
Figure 5: The Brans-Dicke field Φ\Phi and horizons of charged black holes with Φ\Phi when A=0.15A=0.15 and ω=1.4\omega=-1.4 Hansen:2014rua ; Hansen:2015dxa . Left: β=0\beta=0 and hence there is no Brans-Dicke hair at the outer horizon. Inside the outer apparent horizon, r,v=0r_{,v}=0 (red curve) and r,u=0r_{,u}=0 (light blue curve) horizons appear in a highly sophisticated manner because of the Jordan frame. In the vv\rightarrow\infty limit, there exists a Type 2 singularity, but there exists a gradient of the Brans-Dicke field along the uu direction. Right: β=1\beta=1 and there exists Brans-Dicke hair along the outer apparent horizon. The corresponding hair causes neither the inner apparent horizon nor Type 2 singularity to exist inside the outer horizon.

In this model, we can consider the Brans-Dicke field which is non-minimally coupled with the Maxwell field. Solving the equations of motion derived from Eq. (8) numerically in the double-null formalism (see Appendices B and C for more details), we demonstrate our numerical results from Hansen:2014rua ; Hansen:2015dxa ; Nakonieczna:2015umf ; Nakonieczna:2016iof in the following.

  • 1. In Fig. 5, we considered β=0\beta=0 and ω=1000\omega=1000. The large value of ω\omega implies that the dynamics of the Brans-Dicke field can be switched off effectively and this result is consistent with Einstein gravity.

    Interestingly, the dynamics of the real and imaginary parts of the complex scalar field do not disappear in the large vv limit. We find that the derivative of a complex scalar field with respect to uu might vanish but might not vanish with respect to vv. Therefore, the dynamics of a non-trivial complex scalar field can exist in the inner horizon.

    If we turn on the Hawking radiation, one can still see the dynamics of the scalar fields that penetrate beyond the inner apparent horizon: see Fig. 15 of Hong:2008mw .

  • 2. We considered β=0\beta=0 in the left of Fig. 5. The minimal coupling between the Maxwell field and the Brans-Dicke field in the Einstein frame does not form Brans-Dicke hair on the outer horizon.

    Inside the outer apparent horizon, the appearance of r,v=0r_{,v}=0 (red curve) and r,u=0r_{,u}=0 (light blue curve) horizons are highly sophisticated because we perform the calculations in the Jordan frame. Hence, both H1 and H2 can appear in various places in the Jordan frame.

    In the vv\rightarrow\infty limit, there exists a Type 2 singularity, but there exists a gradient of the Brans-Dicke field along the uu direction. Therefore, one can interpret that there exists the Cauchy horizon in the vv\rightarrow\infty limit, and at the same time, there exists a non-trivial scalar field dynamics along the Cauchy horizon.

  • 3. In the right of Fig. 5, we considered the β=1\beta=1 case. The non-minimal coupling between the Brans-Dicke field and the Maxwell field gives rise to the existence of the Brans-Dicke hair along the outer apparent horizon.

    One interesting observation is that there is no inner apparent horizon or Type 2 singularity inside the outer horizon. So, this implies that if there exists a hair at the outer horizon, the inner Cauchy horizon structure disappears, which is very different from pure Einstein gravity. Of course, at this moment, this is just a phenomenological observation and we need further mathematical investigations.

Therefore, the scalar hair at the outer apparent horizon plays a more important role than the gradient of the (scalar or vector) field along the inner apparent horizon to determine the existence of the Cauchy horizon.

IV Conclusion

In this paper, we investigated the properties of inner horizons. First, we distinguished two notions of horizons, where one is the quasi-local notion (e.g., apparent horizons) and another is the global notion (Cauchy horizon). For the quasi-local notion, we can distinguish two horizons, i.e., r,v=0r_{,v}=0 or r,u=0r_{,u}=0. For the global notion, we can distinguish the horizons whether the singularity was generated from the classical singularity or the quantum (Planckian) singularity. In fully dynamic cases, one can construct several models so that all the notions can be distinguished.

Our results provide a few important implications for the cosmic censorship conjecture. How can we define the singularity inside a black hole, for example, in the classical or quantum sense? What are the conditions to rescue the cosmic censorship conjecture, for example, in the classical or quantum sense? If there is no Cauchy horizon, is it related to the quasi-local or the global notion of the horizons? To preserve the cosmic censorship conjecture, should the horizon be disappeared at all or the existence of the inner horizon itself is still fine? These are the diverse questions that we need to ask before we proceed with more proof of the non-existence of the Cauchy horizon. We will leave these fascinating topics for future investigation and hope that our paper can shed light on the future advances for this issue.

In addition, we can provide some guidelines for the Cauchy horizon theorems. Our simulations demonstrate that when we consider black holes in the dynamical cases, their interior structure is no longer to be described by the static and regular metric. Also, the scalar hair on the outer horizon is more important than the scalar hair on the inner horizon to determine the existence of the Cauchy horizon. We need to investigate these issues using not only numerical but also analytic approaches. Will this behavior generally be true in most cases? We leave this interesting question to be investigated in the future.

Acknowledgment

We would like to thank Mu-In Park for his valuable comments. DY is supported by the National Research Foundation of Korea (Grant no.: 2021R1C1C1008622, 2021R1A4A5031460). XYC acknowledges the support from the starting grant of Jiangsu University of Science and Technology (JUST).

Appendix A: Charged black holes in Einstein gravity

We start from the model with a complex massless scalar field ϕ\phi and a gauge field AμA_{\mu} Hong:2008mw :

S=g𝑑x4[16π(ϕ;μ+ieAμϕ)(ϕ¯;μieAμϕ¯)18πFμνFμν],\displaystyle S=\int\sqrt{-g}dx^{4}\left[\frac{\mathcal{R}}{16\pi}-\left(\phi_{;\mu}+ieA_{\mu}\phi\right)\left(\bar{\phi}^{;\mu}-ieA^{\mu}\bar{\phi}\right)-\frac{1}{8\pi}F_{\mu\nu}F^{\mu\nu}\right], (10)

where Fab=Ab;aAa;bF_{ab}=A_{b;a}-A_{a;b} and ee is the gauge coupling. We introduce the double-null coordinates

ds2=α2(u,v)dudv+r2(u,v)dΩ22,\displaystyle ds^{2}=-\alpha^{2}(u,v)dudv+r^{2}(u,v)d\Omega_{2}^{2}, (11)

and define variables:

hα,uα,dα,vα,fr,u,gr,v,w4πϕ,u,z4πϕ,v.\displaystyle h\equiv\frac{\alpha_{,u}}{\alpha},\quad d\equiv\frac{\alpha_{,v}}{\alpha},\quad f\equiv r_{,u},\quad g\equiv r_{,v},\quad w\equiv\sqrt{4\pi}\phi_{,u},\quad z\equiv\sqrt{4\pi}\phi_{,v}. (12)

In addition, we assume the Coulomb gauge Aμ=(a(u,v),0,0,0)A_{\mu}=(a(u,v),0,0,0) which allows us to define the charge function q(u,v)2r2a,v/α2q(u,v)\equiv 2r^{2}a_{,v}/\alpha^{2}.

The semi-classical Einstein equation is defined as

Gμν=8π(TμνC+T^μνH),\displaystyle G_{\mu\nu}=8\pi\left(T^{\mathrm{C}}_{\mu\nu}+\langle\hat{T}^{\mathrm{H}}_{\mu\nu}\rangle\right), (13)

where TμνCT^{\mathrm{C}}_{\mu\nu} is the classical energy-momentum tensor and T^μνH\langle\hat{T}^{\mathrm{H}}_{\mu\nu}\rangle is the renormalized energy-momentum tensor. We introduce the SS-wave approximation of the renormalized energy-momentum tensor:

T^uuH\displaystyle\langle\hat{T}^{\mathrm{H}}_{uu}\rangle =\displaystyle= P4πr2(h,uh2),\displaystyle\frac{P}{4\pi r^{2}}\left(h_{,u}-h^{2}\right), (14)
T^uvH=T^vuH\displaystyle\langle\hat{T}^{\mathrm{H}}_{uv}\rangle=\langle\hat{T}^{\mathrm{H}}_{vu}\rangle =\displaystyle= P4πr2d,u,\displaystyle-\frac{P}{4\pi r^{2}}d_{,u}, (15)
T^vvH\displaystyle\langle\hat{T}^{\mathrm{H}}_{vv}\rangle =\displaystyle= P4πr2(d,vd2),\displaystyle\frac{P}{4\pi r^{2}}\left(d_{,v}-d^{2}\right), (16)

with PNP2/12πP\equiv N\ell_{\mathrm{P}}^{2}/12\pi, where NN is the number of scalar field and P\ell_{\mathrm{P}} is the Planck length.

The semi-classical Einstein equations, as well as the scalar field equation (s4πϕs\equiv\sqrt{4\pi}\phi) and the Maxwell equation, are summarized as follows:

d,u=h,v\displaystyle d_{,u}=h_{,v} =\displaystyle= 11P/r2[fgr2+α24r2α2q22r412(wz¯+w¯z)iea2(sz¯s¯z)],\displaystyle\frac{1}{1-P/r^{2}}\left[\frac{fg}{r^{2}}+\frac{\alpha^{2}}{4r^{2}}-\frac{\alpha^{2}q^{2}}{2r^{4}}-\frac{1}{2}(w\overline{z}+\overline{w}z)-\frac{iea}{2}(s\overline{z}-\overline{s}z)\right], (17)
g,v\displaystyle g_{,v} =\displaystyle= 2dgrzz¯Pr(d,vd2),\displaystyle 2dg-rz\overline{z}-\frac{P}{r}\left(d_{,v}-d^{2}\right), (18)
f,u\displaystyle f_{,u} =\displaystyle= 2fhrww¯iear(w¯sws¯)e2a2rss¯Pr(h,uh2),\displaystyle 2fh-rw\overline{w}-iear\left(\overline{w}s-w\overline{s}\right)-e^{2}a^{2}rs\overline{s}-\frac{P}{r}\left(h_{,u}-h^{2}\right), (19)
f,v=g,u\displaystyle f_{,v}=g_{,u} =\displaystyle= fgrα24r+α2q24r3Prd,u,\displaystyle-\frac{fg}{r}-\frac{\alpha^{2}}{4r}+\frac{\alpha^{2}q^{2}}{4r^{3}}-\frac{P}{r}d_{,u}, (20)
a,v\displaystyle a_{,v} =\displaystyle= α2q2r2,\displaystyle\frac{\alpha^{2}q}{2r^{2}}, (21)
q,v\displaystyle q_{,v} =\displaystyle= ier22(s¯zsz¯),\displaystyle-\frac{ier^{2}}{2}\left(\overline{s}z-s\overline{z}\right), (22)
z,u=w,v\displaystyle z_{,u}=w_{,v} =\displaystyle= fzrgwriearzrieagsrie4r2α2qs.\displaystyle-\frac{fz}{r}-\frac{gw}{r}-\frac{iearz}{r}-\frac{ieags}{r}-\frac{ie}{4r^{2}}\alpha^{2}qs. (23)

Appendix B: Charged black holes in Brans-Dicke gravity

We consider the model for the charged black holes in Brans-Dicke gravity Hansen:2014rua ; Hansen:2015dxa

S=g𝑑x4[116π(ΦωΦΦ;μΦ;μV(Φ))+Φβ(12(ϕ;μ+ieAμϕ)(ϕ¯;μieAμϕ¯)116πFμνFμν)],S=\int\sqrt{-g}dx^{4}\left[\frac{1}{16\pi}\left(\Phi\mathcal{R}-\frac{\omega}{\Phi}\Phi_{;\mu}\Phi^{;\mu}-V(\Phi)\right)+\Phi^{\beta}\left(-\frac{1}{2}\left(\phi_{;\mu}+ieA_{\mu}\phi\right)\left(\bar{\phi}^{;\mu}-ieA^{\mu}\bar{\phi}\right)-\frac{1}{16\pi}F_{\mu\nu}F^{\mu\nu}\right)\right], (24)

where Φ\Phi is the Brans-Dicke field, ω\omega is the Brans-Dicke parameter, V(Φ)V(\Phi) is the potential of the Brans-Dicke field, and β\beta is a constant. We define the additional variables:

WΦ,u,ZΦ,v.\displaystyle W\equiv\Phi_{,u},\quad Z\equiv\Phi_{,v}. (25)

The Einstein equations, as well as the Brans-Dicke field equation, the complex scalar field equation, and the Maxwell equation, are given by the following:

d,u=h,v\displaystyle d_{,u}=h_{,v} =\displaystyle= 𝔄𝔅r2rΦ,\displaystyle\mathfrak{A}-\frac{\mathfrak{B}}{r}-\frac{\mathfrak{C}}{2r\Phi}, (26)
f,v=g,u\displaystyle f_{,v}=g_{,u} =\displaystyle= 𝔅2Φ,\displaystyle\mathfrak{B}-\frac{\mathfrak{C}}{2\Phi}, (27)
Z,u=W,v\displaystyle Z_{,u}=W_{,v} =\displaystyle= r,\displaystyle\frac{\mathfrak{C}}{r}, (28)
r,uu\displaystyle r_{,uu} =\displaystyle= 2fhr2Φ(W,u2hW)rω2Φ2W24πrΦβ1TuuC,\displaystyle 2fh-\frac{r}{2\Phi}\left(W_{,u}-2hW\right)-\frac{r\omega}{2\Phi^{2}}W^{2}-4\pi r\Phi^{\beta-1}{T}^{\mathrm{C}}_{uu}, (29)
r,vv\displaystyle r_{,vv} =\displaystyle= 2gdr2Φ(Z,v2dZ)rω2Φ2Z24πrΦβ1TvvC,\displaystyle 2gd-\frac{r}{2\Phi}\left(Z_{,v}-2dZ\right)-\frac{r\omega}{2\Phi^{2}}Z^{2}-4\pi r\Phi^{\beta-1}{T}^{\mathrm{C}}_{vv}, (30)
a,v\displaystyle a_{,v} =\displaystyle= α2q2r2,\displaystyle\frac{\alpha^{2}q}{2r^{2}}, (31)
q,v\displaystyle q_{,v} =\displaystyle= ier22(s¯zsz¯)βqZΦ,\displaystyle-\frac{ier^{2}}{2}\left(\overline{s}z-s\overline{z}\right)-\beta q\frac{Z}{\Phi}, (32)
a,vv\displaystyle a_{,vv} =\displaystyle= α2r2(dgr)qieα24(zs¯sz¯)βqα2Z2r2Φ,\displaystyle\frac{\alpha^{2}}{r^{2}}\left(d-\frac{g}{r}\right)q-\frac{ie\alpha^{2}}{4}\left(z\overline{s}-s\overline{z}\right)-\beta q\frac{\alpha^{2}Z}{2r^{2}\Phi}, (33)
q,u\displaystyle q_{,u} =\displaystyle= ier22(s¯wsw¯)r2e2ass¯βqWΦ,\displaystyle\frac{ier^{2}}{2}\left(\overline{s}w-s\overline{w}\right)-r^{2}e^{2}as\overline{s}-\beta q\frac{W}{\Phi}, (34)
a,uv\displaystyle a_{,uv} =\displaystyle= α2r2(hfr)q+ieα24(ws¯sw¯)α22e2ass¯βqα2W2r2Φ,\displaystyle\frac{\alpha^{2}}{r^{2}}\left(h-\frac{f}{r}\right)q+\frac{ie\alpha^{2}}{4}\left(w\overline{s}-s\overline{w}\right)-\frac{\alpha^{2}}{2}e^{2}as\overline{s}-\beta q\frac{\alpha^{2}W}{2r^{2}\Phi}, (35)
s,uv\displaystyle s_{,uv} =\displaystyle= fzrgwriearzrieagsrie4r2α2qsβ2Φ(Wz+Zw+iesaZ),\displaystyle-\frac{fz}{r}-\frac{gw}{r}-\frac{iearz}{r}-\frac{ieags}{r}-\frac{ie}{4r^{2}}\alpha^{2}qs-\frac{\beta}{2\Phi}\left(Wz+Zw+iesaZ\right), (36)

where

𝔄\displaystyle\mathfrak{A} \displaystyle\equiv 2πα2r2Φβ1TθθC12rΦ(gW+fZ)ω2Φ2WZ,\displaystyle-\frac{2\pi\alpha^{2}}{r^{2}}\Phi^{\beta-1}T^{\mathrm{C}}_{\theta\theta}-\frac{1}{2r\Phi}\left(gW+fZ\right)-\frac{\omega}{2\Phi^{2}}WZ, (37)
𝔅\displaystyle\mathfrak{B} \displaystyle\equiv α24rfgr+4πrΦβ1TuvC1Φ(gW+fZ),\displaystyle-\frac{\alpha^{2}}{4r}-\frac{fg}{r}+4\pi r\Phi^{\beta-1}T^{\mathrm{C}}_{uv}-\frac{1}{\Phi}(gW+fZ), (38)
\displaystyle\mathfrak{C} \displaystyle\equiv fZgW2πrα2Φβ3+2ω(TC2βEM),\displaystyle-fZ-gW-\frac{2\pi r\alpha^{2}\Phi^{\beta}}{3+2\omega}\left(T^{\mathrm{C}}-2\beta{\mathcal{L}}^{\mathrm{EM}}\right), (39)
TuuC\displaystyle T^{\mathrm{C}}_{uu} \displaystyle\equiv 14π[ww¯+iea(w¯sws¯)+e2a2ss¯],\displaystyle\frac{1}{4\pi}\left[w\overline{w}+iea(\overline{w}s-w\overline{s})+e^{2}a^{2}s\overline{s}\right], (40)
TuvC\displaystyle T^{\mathrm{C}}_{uv} \displaystyle\equiv (a,v)24πα2,\displaystyle\frac{{(a_{,v})}^{2}}{4\pi\alpha^{2}}, (41)
TvvC\displaystyle T^{\mathrm{C}}_{vv} \displaystyle\equiv 14πzz¯,\displaystyle\frac{1}{4\pi}z\overline{z}, (42)
TθθC\displaystyle T^{\mathrm{C}}_{\theta\theta} \displaystyle\equiv r24πα2[(wz¯+zw¯)+iea(z¯szs¯)+2(a,v)2α2],\displaystyle\frac{r^{2}}{4\pi\alpha^{2}}\left[(w\overline{z}+z\overline{w})+iea(\overline{z}s-z\overline{s})+\frac{2{(a_{,v})}^{2}}{\alpha^{2}}\right], (43)
TC\displaystyle T^{\mathrm{C}} \displaystyle\equiv 4α2TuvC+2r2TθθC,\displaystyle-\frac{4}{\alpha^{2}}T^{\mathrm{C}}_{uv}+\frac{2}{r^{2}}T^{\mathrm{C}}_{\theta\theta}, (44)
EM\displaystyle\mathcal{L}^{\mathrm{EM}} \displaystyle\equiv 14πα2(wz¯+zw¯)+iea4πα2(z¯szs¯)+a,v22πα4.\displaystyle\frac{1}{4\pi\alpha^{2}}\left(w\bar{z}+z\bar{w}\right)+\frac{iea}{4\pi\alpha^{2}}\left(\bar{z}s-z\bar{s}\right)+\frac{{a_{,v}}^{2}}{2\pi\alpha^{4}}. (45)

Appendix C: Boundary conditions, integration algorithm, and numerical convergence test

Continuing from Appendix B, we need to provide the initial conditions for all functions (α,r,g,f,h,d,s,w,z,a,q,Φ,W,Z\alpha,r,g,f,h,d,s,w,z,a,q,\Phi,W,Z) on the initial surfaces u=0u=0 and v=0v=0 to solve all equations.

Using the gauge freedom to choose the initial rr function, we define r(0,v)=r0+0.5vr(0,v)=r_{0}+0.5v and r(u,0)=r00.5ur(u,0)=r_{0}-0.5u, where r0r_{0} is an arbitrary constant. We used r0=10r_{0}=10 in the Einstein gravity case and r0=20r_{0}=20 in the Brans-Dicke gravity case. This provides g(0,v)=0.5g(0,v)=0.5 and f(u,0)=0.5f(u,0)=-0.5. We have specified the form of the scalar field profile s(0,v)s(0,v) along the u=0u=0 surface, while it vanishes along the v=0v=0 surface. In addition, we fix Φ(u,0)=Φ(0,v)=1\Phi(u,0)=\Phi(0,v)=1 in the Brans-Dicke model. This automatically fixes z(0,v)=ϕ,v(0,v)z(0,v)=\phi_{,v}(0,v), w(u,0)=ϕ,u(u,0)w(u,0)=\phi_{,u}(u,0), Z(0,v)=0Z(0,v)=0, and W(u,0)=0W(u,0)=0. Hence, q(u,0)=a(u,0)=0q(u,0)=a(u,0)=0 and α(u,0)=1\alpha(u,0)=1. Other functions will be determined by the equations with uuuu- or vvvv- derivation equations (Eqs. (29), (30), (33)). Interested readers can refer to Nakonieczna:2018tih on the generic procedure to assign the boundary conditions.

Regarding the convergence test of numerical calculations, interested readers can also refer to the appendices in the following references for detailed description Hong:2008mw ; Hansen:2014rua ; Hansen:2015dxa . Here Fig. 6 shows the error for the case of charged black holes in Brans-Dicke gravity, the error inside the horizon (u=10,15u=10,15) is still sufficiently small, while the convergence is the second order (this depends on the detailed algorithm).

Refer to caption
Figure 6: An example of the convergence test of Fig. 3 Hong:2008mw . Here, plots of errors with different step sizes are denoted. We compared |r(1)r(2)|/r(2)|r_{(1)}-r_{(2)}|/r_{(2)} (solid curves) and 4|r(2)r(4)|/r(4)4|r_{(2)}-r_{(4)}|/r_{(4)} (dashed curves) along a few constant uu lines, where r(n)r_{(n)} is calculated in an nn-times finer simulation than r(1)r_{(1)}.

References

  • (1) S. W. Hawking and R. Penrose, Proc. Roy. Soc. Lond. A 314, 529-548 (1970).
  • (2) H. Reissner, Ann. Physik 50, 106 (1916).
  • (3) G. Nordstrom, Proc. Kon. Ned. Akad. Wet. 20, 1238 (1918).
  • (4) R. Penrose, Riv. Nuovo Cim. 1, 252-276 (1969).
  • (5) D. Hwang and D. Yeom, Phys. Rev. D 84, 064020 (2011) [arXiv:1010.2585 [gr-qc]].
  • (6) J. Hansen, B. H. Lee, C. Park and D. Yeom, Class. Quant. Grav. 30, 235022 (2013) [arXiv:1307.0266 [hep-th]].
  • (7) M. Simpson and R. Penrose, Int. J. Theor. Phys. 7, 183-197 (1973).
  • (8) E. Poisson and W. Israel, Phys. Rev. D 41, 1796-1809 (1990).
  • (9) A. Bonanno, S. Droz, W. Israel and S. M. Morsink, Phys. Rev. D 50, 7372-7375 (1994) [arXiv:gr-qc/9403019 [gr-qc]].
  • (10) Y. S. An, L. Li and F. G. Yang, Phys. Rev. D 104 (2021) no.2, 024040 [arXiv:2106.01069 [gr-qc]].
  • (11) R. G. Cai, L. Li and R. Q. Yang, JHEP 03 (2021), 263 [arXiv:2009.05520 [gr-qc]].
  • (12) S. A. Hartnoll, G. T. Horowitz, J. Kruthoff and J. E. Santos, SciPost Phys. 10 (2021) no.1, 009 [arXiv:2008.12786 [hep-th]].
  • (13) O. J. C. Dias, G. T. Horowitz and J. E. Santos, JHEP 12 (2021), 179 [arXiv:2110.06225 [hep-th]].
  • (14) D. O. Devecioglu and M. I. Park, Phys. Lett. B 829 (2022), 137107 [arXiv:2101.10116 [hep-th]].
  • (15) N. Grandi and I. Salazar Landea, JHEP 05 (2021), 152 [arXiv:2102.02707 [gr-qc]].
  • (16) R. Q. Yang, R. G. Cai and L. Li, Class. Quant. Grav. 39 (2022) no.3, 035005 [arXiv:2104.03012 [gr-qc]].
  • (17) M. Henneaux, JHEP 03 (2022), 062 [arXiv:2202.04155 [hep-th]].
  • (18) J. Luk, S. J. Oh and Y. Shlapentokh-Rothman, Annales Henri Poincare 24 (2023) no.2, 363-411 [arXiv:2201.12294 [gr-qc]].
  • (19) D. O. Devecioglu and M. I. Park, [arXiv:2307.10532 [hep-th]].
  • (20) R. d’Inverno, “Introducing Einstein’s relativity,” Oxford University Press (1992).
  • (21) S. W. Hawking and G. F. R. Ellis, “The Large Scale Structure of Space-Time,” Cambridge University Press (2023).
  • (22) A. Nakonieczna, Ł. Nakonieczny and D. Yeom, Int. J. Mod. Phys. D 28, no.03, 1930006 (2018) [arXiv:1805.12362 [gr-qc]].
  • (23) S. Hod and T. Piran, Phys. Rev. Lett. 81, 1554-1557 (1998) [arXiv:gr-qc/9803004 [gr-qc]].
  • (24) D. Hwang, H. B. Kim and D. Yeom, Class. Quant. Grav. 29, 055003 (2012) [arXiv:1105.1371 [gr-qc]].
  • (25) E. Sorkin and T. Piran, Phys. Rev. D 63, 124024 (2001) [arXiv:gr-qc/0103090 [gr-qc]].
  • (26) S. E. Hong, D. Hwang, E. D. Stewart and D. Yeom, Class. Quant. Grav. 27, 045014 (2010) [arXiv:0808.1709 [gr-qc]].
  • (27) J. Hansen and D. Yeom, JHEP 10, 040 (2014) [arXiv:1406.0976 [hep-th]].
  • (28) J. Hansen and D. Yeom, JCAP 09, 019 (2015) [arXiv:1506.05689 [hep-th]].
  • (29) A. Nakonieczna and D. Yeom, JHEP 02, 049 (2016) [arXiv:1512.06712 [gr-qc]].
  • (30) A. Nakonieczna and D. Yeom, JHEP 05, 155 (2016) [arXiv:1604.04718 [gr-qc]].
  • (31) D. Hwang, B. H. Lee and D. Yeom, JCAP 12, 006 (2011) [arXiv:1110.0928 [gr-qc]].
  • (32) P. Chen and D. Yeom, JCAP 10, 022 (2015) [arXiv:1506.06713 [gr-qc]].
  • (33) A. Ori, Phys. Rev. Lett. 67, 789-792 (1991).
  • (34) D. Yeom and H. Zoe, Int. J. Mod. Phys. A 26, 3287-3314 (2011) [arXiv:0907.0677 [hep-th]].
  • (35) W. Kim, B. H. Lee and D. Yeom, JHEP 05, 060 (2013) [arXiv:1301.5138 [gr-qc]].
  • (36) P. C. W. Davies, S. A. Fulling and W. G. Unruh, Phys. Rev. D 13, 2720-2723 (1976).
  • (37) G. Dvali, C. Gomez and D. Lust, Fortsch. Phys. 61, 768-778 (2013) [arXiv:1206.2365 [hep-th]].
  • (38) E. Poisson, Annals Israel Phys. Soc. 13, 85 (1997) [arXiv:gr-qc/9709022 [gr-qc]].