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

Loosely trapped surface and dynamically transversely trapping surface in Einstein-Maxwell system

Kangjae Lee1    Tetsuya Shiromizu1,2    Hirotaka Yoshino3    Keisuke Izumi2,1    Yoshimune Tomikawa4 1Department of Mathematics, Nagoya University, Nagoya 464-8602, Japan 2Kobayashi-Maskawa Institute, Nagoya University, Nagoya 464-8602, Japan 3Advanced Mathematical Institute, Osaka City University, Osaka 558-8585, Japan 4Faculty of Economics, Matsuyama University, Matsuyama 790-8578, Japan
Abstract

We study the properties of the loosely trapped surface (LTS) and the dynamically transversely trapping surface (DTTS) in Einstein-Maxwell systems. These concepts of surfaces were proposed by the four of the present authors in order to characterize strong gravity regions. We prove the Penrose-like inequalities for the area of LTSs/DTTSs. Interestingly, although the naively expected upper bound for the area is that of the photon sphere of a Reissner-Nordström black hole with the same mass and charge, the obtained inequalities include corrections represented by the energy density or pressure/tension of electromagnetic fields. Due to this correction, the Penrose-like inequality for the area of LTSs is tighter than the naively expected one. We also evaluate the correction term numerically in the Majumdar-Papapetrou two-black-hole spacetimes.

\subjectindex

E0, E31, A13

1 Introduction

Since a black hole creates a strong gravitational field, there exists unstable circular orbits for photons. For a spherically symmetric system, the collection of them makes a surface called a photon sphere. In the Schwarzschild spacetime, for example, the photon sphere exists at the surface r=3mr=3m, where mm is the Arnowitt-Deser-Misner (ADM) mass. Furthermore, a generalized concept of a photon sphere, which is called a photon surface, has been proposed Claudel:2000 . However, the definition of the photon surface requires a highly symmetric spacetime (e.g., Yoshino:2016 ). Moreover, the existence of photon surfaces does not necessarily mean that the gravitational field is strong there Claudel:2000 .

A photon sphere is directly related to observational phenomena. The quasinormal modes of black holes are basically determined by the properties of photon spheres cardoso2009 . The black hole shadow, whose direct picture has been taken by the recent radio observations of the Event Horizon Telescope Collaboration akiyama , is also determined by a photon sphere Virbhadra:1999 .

Motivated by the recent observations, the four of the present authors proposed concepts that characterizes strong gravity regions; a loosely trapped surface (LTS) shiromizu2017 , a transversely trapping surface (TTS) Yoshino:2017 , and a dynamically transversely trapping surface (DTTS) Yoshino:2020-1 ; Yoshino:2020-2 (see also Ref. galtsov for an extension of a TTS). For certain cases, we have proved inequalities analogous to the Penrose inequality penrose1973 , that is, their areas are equal to or less than 4π(3m)24\pi(3m)^{2}, where mm is the ADM mass (see also Ref. hod for an earlier work for a photon sphere). The upper bound is realized for a photon sphere in the Schwarzschild black hole, where 3m3m comes from the areal radius of unstable circular photon orbits. For an LTS, the application is restricted to a spacelike hypersurfaces with a positive Ricci scalar. This restriction is natural because it is guaranteed by the positivity of the energy density for maximally sliced initial data. For a DTTS, on the other hand, one requires the non-positivity of the radial pressure on a DTTS in addition to the positivity of the Ricci scalar. This requirement of the non-positive pressure does not significantly restrict the situation because the vacuum cases do work. However, it remains a little mystery why the non-positivity of the radial pressure is required.

Therefore, in this paper, we shall discuss non-vacuum cases. As a first typical example, we will focus on Einstein-Maxwell systems. We adopt Jang’s work jang1979 to show the inequality, that is, we will employ the method of the inverse mean curvature flow geroch ; wald1977 111The resulted inequality for the minimal surface in Ref. jang1979 has a lower bound. However, it is known that the lower bound is violated for multi-black-hole systems weinstein2004 . The proof for multi-black-hole systems in the Einstein-Maxwell theory is given in Refs. khuri2013 ; khuri2017 .. The upper bound is expected to be given by the area of the outermost photon sphere, i.e., the locus of the unstable circular orbits of photons in a spherically symmetric charged black hole spacetime, namely, a Reissner-Nordström spacetime with the same mass and charge. As a consequence, however, we see that the obtained inequalities depend on (the part of) the energy density and the pressure/tension of the electromagnetic fields which give corrections to the naively expected upper bound. This is impressive because the Penrose inequality for apparent horizons does not depend on such quantities.

The rest of this paper is organized as follows. In Sect. 2, we will briefly describe the Maxwell theory in a curved spacetime. Some of the notations will be explained together. In Sect. 3, we will present the definition of the LTS and prove the Penrose-like inequality in the Einstein-Maxwell system. Then, in Sect. 4, we will present the definition of the DTTS and prove the Penrose-like inequality in the Einstein-Maxwell system. In Sect. 5, we will revisit the problem of DTTSs in Majumdar-Papapetrou two-black-hole spacetimes, which was studied in our previous paper Yoshino:2020-2 , from the viewpoint of the current study. We will examine the properties of the correction term of the Penrose-like inequality through numerical calculations. The last section will give a summary and discussions. In Appendix, we will shortly discuss the case of a TTS defined for static/stationary spacetimes. Note that we use following units in which the speed of light c=1c=1, the Newtonian constant of gravitation G=1G=1 and the Coulomb constant 1/(4πε0)=11/(4\pi\varepsilon_{0})=1, where ε0\varepsilon_{0} is the permittivity of vacuum.

2 Einstein-Maxwell theory and setup

In this paper, we consider an asymptotically flat spacelike hypersurface Σ\Sigma in a four-dimensional spacetime with a metric gabg_{ab}. We suppose nan^{a} to be the future-directed unit normal to Σ\Sigma and the induced metric of Σ\Sigma is given by γab=gab+nanb\gamma_{ab}=g_{ab}+n_{a}n_{b}. In Σ\Sigma, we consider a two-dimensional closed surface, an LTS (denoted by S0S_{0}) or a DTTS (denoted by σ0\sigma_{0}), with the induced metric habh_{ab}. The outward unit normal to that surface is rar_{a} (and therefore, γab=hab+rarb\gamma_{ab}=h_{ab}+r_{a}r_{b}).

We assume the presence of electromagnetic fields on Σ\Sigma. The electromagnetic fields are specified by the anti-symmetric tensor, FabF_{ab}, and its Hodge dual, Fab=(1/2)ϵabcdFcd{}^{*}F_{ab}=(1/2)\epsilon_{abcd}F^{cd}, where ϵabcd\epsilon_{abcd} is the Levi-Civita symbol in a four-dimensional spacetime. The tensors FabF_{ab} and Fab{}^{*}F_{ab} follow Maxwell’s equations,

aFab=0,aFab=4πjb,\nabla_{a}{}^{*}F^{ab}=0,\qquad\nabla_{a}F^{ab}=-4\pi j^{b}, (1)

where jbj^{b} is the four-current vector. The electric and magnetic fields, EaE^{a} and BaB^{a}, are defined by

Ea:=Fabnb,Ba:=Fabnb,E_{a}:=F_{ab}n^{b},\qquad B_{a}:=-{}^{*}F_{ab}n^{b}, (2)

respectively. Obviously, these fields are tangent to Σ\Sigma because Eana=Bana=0E_{a}n^{a}=B_{a}n^{a}=0 holds. The electric charge density ρe\rho_{\rm e} and the electric current (Je)a(J_{\rm e})^{a} are defined by

ρe:=jana,(Je)a:=γabjb.\rho_{\rm e}:=-j_{a}n^{a},\qquad(J_{\rm e})^{a}:={\gamma^{a}}_{b}j^{b}. (3)

In our paper, we require the charge density to vanish (outside an LTS or a DTTS), i.e. ρe=0\rho_{\rm e}=0, but do not necessarily require the electric current (Je)a(J_{\rm e})^{a} to be zero.

The electric and magnetic fields satisfy Gauss’ laws, DaEa=4πρeD_{a}E^{a}=4\pi\rho_{\rm e} and DaBa=0D_{a}B^{a}=0, where DaD_{a} is the covariant derivative with respect to γab\gamma_{ab}. The total electric and magnetic charges are

4πqe:=SEara𝑑A,4πqm:=SBara𝑑A,4\pi q_{\rm e}:=\int_{S_{\infty}}E^{a}r_{a}dA,\qquad 4\pi q_{\rm m}:=\int_{S_{\infty}}B^{a}r_{a}dA, (4)

where SS_{\infty} is a sphere at spacelike infinity and rar_{a} is the outward unit normal to SS_{\infty}. It would be important to point out that qeq_{\rm e} can have a nonzero value even if the electric charge density is zero throughout the spacetime, as one can understand by imagining a spacelike hypersurface with an Einstein-Rosen bridge and two asymptotically flat regions in the maximally extended Reissner-Nordström spacetime. Similarly, although we assume the absence of magnetic monopoles throughout the paper, the value of qmq_{\rm m} can be nonzero. In a spherically symmetric spacetime, if both electric and magnetic fields are present, the total squared charge defined by

q2:=qe2+qm2q^{2}:=q_{\rm e}^{2}+q_{\rm m}^{2} (5)

appears in the spacetime metric of the Reissner-Nordström solution. We handle the magnetic charge qmq_{\rm m} in the following way. In Sects. 3 and 4, we will assume qm=0q_{\rm m}=0 (and therefore, q2=qe2q^{2}=q_{\rm e}^{2}), and derive the Penrose-like inequalities. In the final section, we will discuss the modifications to those inequalities when qmq_{\rm m} is nonzero.

The energy-momentum tensor for the electromagnetic fields is given by

Tab(em)=14π(FacFbc14gabFcdFcd).T_{ab}^{\rm(em)}=\frac{1}{4\pi}\left({F_{a}}^{c}F_{bc}-\frac{1}{4}g_{ab}F_{cd}F^{cd}\right). (6)

In addition to Tab(em)T_{ab}^{\rm(em)}, we consider the presence of ordinary matters whose energy-momentum tensor is Tab(m)T_{ab}^{(m)}. The total energy-momentum tensor is given by Tab=Tab(em)+Tab(m)T_{ab}=T_{ab}^{\rm(em)}+T_{ab}^{(m)}. The relations of particular importance in this paper are the energy density,

8πρ:=8πTabnanb=EaEa+BaBa+8πρ(m),8\pi\rho:=8\pi T_{ab}n^{a}n^{b}=E^{a}E_{a}+B^{a}B_{a}+8\pi\rho^{(m)}, (7)

where ρ(m):=Tab(m)nanb\rho^{(m)}:=T^{(m)}_{ab}n^{a}n^{b}, and the radial pressure,

8πPr:=8πTabrarb=(EaEb+BaBb)hab(Eara)2(Bara)2+8πPr(m),8\pi P_{r}:=8\pi T_{ab}r^{a}r^{b}=(E_{a}E_{b}+B_{a}B_{b})h^{ab}-(E_{a}r^{a})^{2}-(B_{a}r^{a})^{2}+8\pi P_{r}^{(m)}, (8)

where Pr(m):=Tab(m)rarbP_{r}^{(m)}:=T^{(m)}_{ab}r^{a}r^{b}.

3 Loosely trapped surface in Einstein-Maxwell system

In this section, we review the definition of an LTS following Ref. shiromizu2017 , and show the Penrose-like inequality for it in Einstein-Maxwell systems.

3.1 Definition of an LTS

The definition of an LTS is motivated by the following observation. As an example, we consider a Reissner-Nordström spacetime. The metric is

ds2=fRN(r)dt2+fRN1(r)dr2+r2dΩ22,\displaystyle ds^{2}=-f_{\rm RN}(r)dt^{2}+f^{-1}_{\rm RN}(r)dr^{2}+r^{2}d\Omega_{2}^{2}, (9)

where fRN(r):=12m/r+q2/r2f_{\rm RN}(r):=1-2m/r+q^{2}/r^{2}, and mm and qq are the ADM mass and total charge, respectively. dΩ22d\Omega_{2}^{2} is the two-dimensional metric of the unit round sphere. From the behavior of a null geodesic, one can find unstable circular orbits of photons at

r=rp:=3m+9m28q22,\displaystyle r=r_{p}:=\frac{3m+{\sqrt{9m^{2}-8q^{2}}}}{2}, (10)

where we suppose 9m28q29m^{2}\geq 8q^{2}. Note that a photon sphere exists even if the spacetime possesses a naked singularity at the center for (9/8)m2q2>m2(9/8)m^{2}\geq q^{2}>m^{2}.

Now we define a similar concept to a photon sphere for general setups in terms of geometry. Here, we recall the fact that an apparent horizon is the minimal surface on time-symmetric initial data. Therefore, one possibility to specify a strong gravity region is to employ the mean curvature, that is, the trace of the extrinsic curvature of two-dimensional surfaces. Therefore, we look at the mean curvature for the Reissner-Nordström spacetime. It is easy to see that the mean curvature of an r=r= constant surface on t=t= constant hypersurface is given by

k=2rfRN1/2(r).\displaystyle k=\frac{2}{r}f^{1/2}_{\rm RN}(r). (11)

From the first derivative of kk with respect to rr,

dkdr=2r2(13mr+2q2r2)fRN1/2(r),\displaystyle\frac{dk}{dr}=-\frac{2}{r^{2}}\left(1-\frac{3m}{r}+\frac{2q^{2}}{r^{2}}\right)f^{-1/2}_{\rm RN}(r), (12)

we find that the maximum value of kk exists at r=rpr=r_{p}. This is exactly the same location with that of unstable circular orbits of photons. In the region between the event horizon and the photon sphere at r=rpr=r_{p}, the mean curvature satisfies k0k\geq 0 and dk/dr0dk/dr\geq 0.

From the above argument, one may adopt the following definition of an LTS shiromizu2017 .

Definition 1.

A loosely trapped surface (LTS), S0S_{0}, is defined as a compact two-surface in a spacelike hypersurface Σ\Sigma, and has the mean curvature kk for the outward spacelike normal vector such that k|S0>0k|_{S_{0}}>0 and k|S00k^{\prime}|_{S_{0}}\geq 0, where is the derivative along the outward spacelike normal vector.

3.2 Penrose-like inequality for an LTS

In this section, we present the inequality for the area of an LTS in Einstein-Maxwell systems. Our theorem is as follows:

Theorem 1.

Let Σ\Sigma be an asymptotically flat spacelike hypersurface with the Ricci scalar R(3)2(EaEa+BaBa)+16πρ(m){}^{(3)}R\geq 2\left(E^{a}E_{a}+B^{a}B_{a}\right)+16\pi\rho^{(m)}, where ρ(m)\rho^{(m)} is the non-negative energy density for other matters.222From Eq. (7) and the Hamiltonian constraint of the Einstein equations, this condition is equivalent to KabKabK2K_{ab}K^{ab}\geq K^{2}, where KabK_{ab} is the extrinsic curvature of Σ\Sigma. This condition is obviously satisfied by maximally sliced hypersurfaces, on which K=0K=0 holds. We assume that Σ\Sigma is foliated by the inverse mean curvature flow, and a slice of the foliation parameterized by yy, SyS_{y}, has topology S2S^{2}. We also suppose the electric charge density ρe\rho_{\rm e} to vanish outside the LTS, S0S_{0}. Then, the areal radius r0=(A0/4π)1/2r_{0}=(A_{0}/4\pi)^{1/2} of the LTS, S0S_{0}, in Σ\Sigma satisfies the inequality

m13(1+Φ0+)r0+2q23r0,\displaystyle m\geq\frac{1}{3}\Bigl{(}1+\Phi_{0}^{+}\Bigr{)}r_{0}+\frac{2q^{2}}{3r_{0}}, (13)

where mm is the ADM mass and qq is the total charge. Φ0+\Phi_{0}^{+} is defined by

Φ0+:=18πS0[(Bara)2+(EaEb+BaBb)hab]𝑑A,\displaystyle\Phi^{+}_{0}:=\frac{1}{8\pi}\int_{S_{0}}[(B_{a}r^{a})^{2}+(E_{a}E_{b}+B_{a}B_{b})h^{ab}]dA, (14)

where habh_{ab} is the induced metric of S0S_{0} and rar^{a} is the outward-directed unit normal vector to S0S_{0} in Σ\Sigma.

Proof.

On Σ\Sigma, the derivative of the mean curvature kk along rar^{a} is given by

raDak=φ1𝒟2φ12R(3)+12R(2)12(k2+kabkab),r^{a}D_{a}k\ =\ -\varphi^{-1}{\cal D}^{2}\varphi-\frac{1}{2}{}^{(3)}R+\frac{1}{2}{}^{(2)}R-\frac{1}{2}(k^{2}+k_{ab}k^{ab}), (15)

where DaD_{a} is the covariant derivative of Σ\Sigma, 𝒟a{\cal D}_{a} is the covariant derivative of S0S_{0}, R(2){}^{(2)}R is the Ricci scalar of S0S_{0}, R(3){}^{(3)}R is the Ricci scalar of Σ\Sigma, kabk_{ab} is the extrinsic curvature of S0S_{0} and φ\varphi is the lapse function for yy, that is, ra=φDayr_{a}=\varphi D_{a}y. Then, the integration of Eq. (15) over S0S_{0} gives us

12S0R(2)𝑑A\displaystyle\frac{1}{2}\int_{S_{0}}{}^{(2)}RdA =\displaystyle= S0[raDak+φ2(𝒟φ)2+12R(3)+12k~abk~ab+34k2]\displaystyle\int_{S_{0}}\left[r^{a}D_{a}k+\varphi^{-2}\left(\cal{D}\varphi\right)^{2}+\frac{1}{2}{}^{(3)}R+\frac{1}{2}\tilde{k}_{ab}\tilde{k}^{ab}+\frac{3}{4}k^{2}\right] (16)
\displaystyle\geq 34S0[k2+43(EaEa+BaBa)]𝑑A\displaystyle\frac{3}{4}\int_{S_{0}}\left[k^{2}+\frac{4}{3}(E^{a}E_{a}+B^{a}B_{a})\right]dA
=\displaystyle= 34S0[k2+43(Eara)2]𝑑A+8πΦ0+,\displaystyle\frac{3}{4}\int_{S_{0}}\left[k^{2}+\frac{4}{3}(E_{a}r^{a})^{2}\right]dA+8\pi\Phi^{+}_{0},

where k~ab:=kab(1/2)habk\tilde{k}_{ab}:=k_{ab}-(1/2)h_{ab}k. Note that Φ0+0\Phi^{+}_{0}\geq 0. Using the Gauss-Bonnet theorem and Cauchy-Schwarz inequality, we can derive the following inequality for the mean curvature

S0k2𝑑A16π343(4πq)2A0323πΦ0+,\displaystyle\int_{S_{0}}k^{2}dA\leq\frac{16\pi}{3}-\frac{4}{3}\frac{(4\pi q)^{2}}{A_{0}}-\frac{32}{3}\pi\Phi^{+}_{0}, (17)

where we used Gauss’ law for the electric field, S0Eara𝑑A=SEara𝑑A=4πq\int_{S_{0}}E_{a}r^{a}dA=\int_{S_{\infty}}E_{a}r^{a}dA=4\pi q. Here, SS_{\infty} denotes the two-sphere at spacelike infinity.

Let us consider Geroch’s quasilocal energy geroch ; wald1977 ; jang1979

E(y):=A1/2(y)64π3/2Sy(2R(2)k2)𝑑A,\displaystyle E(y):=\frac{A^{1/2}(y)}{64\pi^{3/2}}\int_{S_{y}}(2{}^{(2)}R-k^{2})dA, (18)

where A(y)A(y) is the area of SyS_{y}. Here, we suppose that the surfaces y=0y=0 and y=y=\infty correspond to the LTS and a sphere at spacelike infinity, respectively. Under the inverse mean curvature flow generated by the condition kφ=1k\varphi=1, the first derivative of E(y)E(y) is computed as

dE(y)dy=A1/2(y)64π3/2Sy[2φ2(𝒟φ)2+R(3)+k~abk~ab]𝑑A.\frac{dE(y)}{dy}\ =\ \frac{A^{1/2}(y)}{64\pi^{3/2}}\int_{S_{y}}\Bigl{[}2\varphi^{-2}({\cal D}\varphi)^{2}+{}^{(3)}R+\tilde{k}_{ab}\tilde{k}^{ab}\Bigr{]}dA. (19)

Using R(3)2(EaEa+BaBa){}^{(3)}R\geq 2(E^{a}E_{a}+B^{a}B_{a}), we can derive

dE(y)dyA1/2(y)32π3/2[(4πq)2A(y)+8πΦy+],\displaystyle\frac{dE(y)}{dy}\geq\frac{A^{1/2}(y)}{32\pi^{3/2}}\left[\frac{(4\pi q)^{2}}{A(y)}+8\pi\Phi^{+}_{y}\right], (20)

with the same procedure as the derivation of the inequality of Eq. (17), where 8πΦy+=Sy[(Bara)2+(EaEb+BaBb)hab]𝑑A8\pi\Phi_{y}^{+}=\int_{S_{y}}[(B_{a}r^{a})^{2}+(E_{a}E_{b}+B_{a}B_{b})h^{ab}]dA.

The integration of the inequality of Eq. (20) over yy in the range 0y<0\leq y<\infty implies us

m\displaystyle m \displaystyle\geq A01/24π1/2(1116πS0k2𝑑A)+π1/2q2A01/2+14π1/20Φy+A1/2(y)𝑑y\displaystyle\frac{A_{0}^{1/2}}{4\pi^{1/2}}\left(1-\frac{1}{16\pi}\int_{S_{0}}k^{2}dA\right)+\frac{\pi^{1/2}q^{2}}{A_{0}^{1/2}}+\frac{1}{4\pi^{1/2}}\int^{\infty}_{0}\Phi^{+}_{y}A^{1/2}(y)dy (21)
\displaystyle\geq A01/26π1/2+43π1/2q2A01/2+A01/26π1/2Φ0+.\displaystyle\frac{A_{0}^{1/2}}{6\pi^{1/2}}+\frac{4}{3}\frac{\pi^{1/2}q^{2}}{A_{0}^{1/2}}+\frac{A_{0}^{1/2}}{6\pi^{1/2}}\Phi^{+}_{0}.

In the above, we have used the well-known relation A(y)=A0exp(y)A(y)=A_{0}\exp(y) that holds in the inverse mean curvature flow at the first step, and the non-negativity of Φy+\Phi_{y}^{+} and the inequality of Eq. (17) in the second step. Then, we find the inequality of Eq. (13).

There are four remarks. First, the minimum value of the right-hand side of the inequality of Eq. (13) implies

m2231+Φ0+|q|.\displaystyle m\geq\frac{2{\sqrt{2}}}{3}{\sqrt{1+\Phi_{0}^{+}}}\ |q|. (22)

Setting Φ0+=0\Phi_{0}^{+}=0, this inequality is reduced to m223|q|m\geq\frac{2{\sqrt{2}}}{3}|q| which corresponds to the condition for the existence of a photon sphere in the Reissner-Nordström solution [see Eq. (10)].

Next, under the condition given by Eq. (22), the rearrangement of the inequality of Eq. (13) gives us

4πrLTS2A04πrLTS+2,\displaystyle 4\pi r_{{\rm LTS}-}^{2}\leq A_{0}\leq 4\pi r_{{\rm LTS}+}^{2}, (23)

where

rLTS±:=3m±9m28(1+Φ0+)q22(1+Φ0+).\displaystyle r_{{\rm LTS}\pm}:=\frac{3m\pm\sqrt{9m^{2}-8\left(1+\Phi^{+}_{0}\right)q^{2}}}{2\left(1+\Phi^{+}_{0}\right)}. (24)

This inequality must be interpreted carefully in the sense that the lower bound would not hold in general. In the case of the ordinary Riemannian Penrose inequality weinstein2004 , it has been pointed out that the lower bound is expected to be incorrect for multi-black-hole systems. We consider this also may be the case for an LTS with multiple components. This is not a contradiction to our proof since an LTS with multiple components is out of the application of our theorem. Therefore, in general, we would have only the upper bound for the area,

A04πrLTS+2.\displaystyle A_{0}\leq 4\pi r_{{\rm LTS}+}^{2}. (25)

Note that if we restrict out attention on an LTS with a single component, the lower bound must hold true. The physical reason is as follows. Let us consider a Reissner-Nordström spacetime in the parameter region (9/8)m2>q2>m2(9/8)m^{2}>q^{2}>m^{2} which possesses both a naked singularity and two photon spheres. The singularity of a Reissner-Nordström spacetime is known to be repulsive. This is because the energy of electromagnetic fields outside a small sphere near the singularity (e.g., in the sense of the Komar integral) exceeds the ADM mass mm, and hence, the gravitational field is generated by negative energy in that region. The repulsive gravitational field is, of course, not strong. This is the reason why an LTS with a single component cannot exist in the vicinity of a naked singularity with an electric charge and a lower bound exists for its area in the Einstein-Maxwell theory.

The third remark is that there is no contribution from Φ0+\Phi_{0}^{+} in the Riemannian Penrose inequality for the Einstein-Maxwell system jang1979 , whereas in our theorem it appears. This is because the Riemannian Penrose inequality discusses the minimal surface with k=0k=0, for which the inequality of Eq. (17) is unnecessary.

Finally, the presence of Φ0+\Phi_{0}^{+} makes the inequality tighter than the case of Φ0+=0\Phi^{+}_{0}=0. From Eq. (15), one can see that the quantity Φ0+\Phi^{+}_{0} appears from (the part of) the energy density in the Hamiltonian constraint. The electromagnetic energy density increases if Φ0+\Phi_{0}^{+} is turned on, and through the relation given by Eq. (15), the positive energy density tends to make the formation of an LTS more difficult. This means that the area of S0S_{0} will become smaller.

4 Dynamically transversely trapping surface in Einstein-Maxwell system

In this section, we first explain the observation that motivates the definition of a DTTS, and introduce the definition of a DTTS. Then we prove an inequality for its area in Einstein-Maxwell systems.

4.1 Definition of a DTTS

The concept of a DTTS is inspired by the induced geometry of a photon surface in a spherically symmetric spacetime. The photon surface is defined as a timelike hypersurface SS such that any photon emitted tangentially to SS at an arbitrary point of SS remains in SS forever Claudel:2000 . Let us consider spherically symmetric spacetimes with the metric,

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

Solving the null geodesic equations, we can find that SS satisfies dr/dt=±f(1b2f/r2)1/2dr/dt=\pm f(1-b^{2}f/r^{2})^{1/2}, where bb is the impact parameter. The induced metric of the photon surface SS is obtained as

ds2|S=α2(r)dt2+r2dΩ22,\displaystyle ds^{2}|_{S}=-\alpha^{2}(r)dt^{2}+r^{2}d\Omega_{2}^{2}, (27)

where α(r):=bf(r)/r\alpha(r):=bf(r)/r. The mean curvature of t=t= constant surface σt\sigma_{t} in SS is given by

k¯=2bf(r)drdt,\displaystyle\bar{k}=\frac{2}{bf(r)}\frac{dr}{dt}, (28)

and the Lie derivative along n¯a\bar{n}^{a} for k¯\bar{k} is computed as

£¯n¯(3)k¯=2r2(f12rf),\displaystyle{}^{(3)}\bar{\mbox{\pounds}}_{\bar{n}}\bar{k}=\frac{2}{r^{2}}\left(f-\dfrac{1}{2}rf^{\prime}\right), (29)

where n¯a\bar{n}^{a} is the future-directed unit normal to σt\sigma_{t} in SS. For the Reissner-Nordstro¨\rm{\ddot{o}}m black hole, it becomes

£¯n¯(3)k¯=2r2(13mr+2q2r2).\displaystyle{}^{(3)}\bar{\mbox{\pounds}}_{\bar{n}}\bar{k}=\frac{2}{r^{2}}\Bigl{(}1-\frac{3m}{r}+\frac{2q^{2}}{r^{2}}\Bigr{)}. (30)

Thus, £¯n¯(3)k¯{}^{(3)}\bar{\mbox{\pounds}}_{\bar{n}}\bar{k} is negative and positive in the inside and outside regions of the photon sphere, respectively. Hence, the non-positivity of £¯n¯(3)k¯{}^{(3)}\bar{\mbox{\pounds}}_{\bar{n}}\bar{k} is expected to indicate the strong gravity.

We now give the definition of a DTTS Yoshino:2020-1 .

Definition 2.

A closed orientable two-dimensional surface σ0\sigma_{0} in a smooth spacelike hypersurface Σ\Sigma is a dynamically transversely trapping surface (DTTS) if and only if there exists a timelike hypersurface SS that intersects Σ\Sigma at σ0\sigma_{0} and satisfies k¯=0\bar{k}=0, max(K¯abkakb)=0\mathrm{max}(\bar{K}_{ab}k^{a}k^{b})=0, and £¯n¯(3)k¯0{}^{(3)}\bar{\mbox{\pounds}}_{\bar{n}}\bar{k}\leq 0 at every point in σ0\sigma_{0}, where k¯\bar{k} is the mean curvature of σ0\sigma_{0} in SS, K¯ab\bar{K}_{ab} is the extrinsic curvature of SS, n¯a\bar{n}^{a} is the unit normal vector of σ0\sigma_{0} in SS, £¯n¯(3){}^{(3)}\bar{\mbox{\pounds}}_{\bar{n}} is the Lie derivative associated with SS and kak^{a} is arbitrary future-directed null vectors tangent to SS.

Since the emitted photons do not form a photon sphere in general without spherical symmetry, here we emit photons in the transverse direction (to satisfy the condition k¯=0\bar{k}=0), and adopt the location of the outermost photons as the surface SS [the condition max(K¯abkakb)=0\mathrm{max}(\bar{K}_{ab}k^{a}k^{b})=0]. Then, we judge that the surface σ0\sigma_{0} exists in a strong gravity region if the condition £¯n¯(3)k¯0{}^{(3)}\bar{\mbox{\pounds}}_{\bar{n}}\bar{k}\leq 0 is satisfied. See our previous papers Yoshino:2020-1 ; Yoshino:2020-2 for more details.

4.2 Penrose-like inequality for a DTTS

We present the following theorem on a Penrose-like inequality for a DTTS:

Theorem 2.

Consider a spacetime which satisfies the Einstein equations with the energy-momentum tensor TabT_{ab} composed of the Maxwell field part, Tab(em)T_{ab}^{(\rm em)} given by Eq. (6), and the ordinal matter part, Tab(m)T_{ab}^{(m)}. We suppose that an asymptotically flat spacelike hypersurface Σ\Sigma is time-symmetric and foliated by the inverse mean curvature flow. We also assume that a slice of the foliation parameterized by yy, where each of the y=constanty=\mathrm{constant} surfaces, σy\sigma_{y}, has topology S2S^{2}, and σ0\sigma_{0} is a convex DTTS. We further assume the electric charge density ρe\rho_{\rm e} to vanish outside σ0\sigma_{0}. Then, if ρ(m):=Tab(m)nanb0\rho^{(m)}:=T_{ab}^{(m)}n^{a}n^{b}\geq 0 in the outside region of σ0\sigma_{0} and Pr(m):=Tab(m)rarb0P_{r}^{(m)}:=T_{ab}^{(m)}r^{a}r^{b}\leq 0 on σ0\sigma_{0}, where nan^{a} is the future-directed timelike unit normal to Σ\Sigma and rar^{a} is the outward spacelike unit normal to σy\sigma_{y} in Σ\Sigma, the areal radius r0=(A0/4π)1/2r_{0}=(A_{0}/4\pi)^{1/2} of the convex DTTS σ0\sigma_{0} satisfies the inequality

m13(1+Φ0)r0+2q23r0,\displaystyle m\geq\frac{1}{3}\Bigl{(}1+\Phi_{0}^{-}\Bigr{)}r_{0}+\frac{2q^{2}}{3r_{0}}, (31)

where mm is the ADM mass and qq is the total charge. Here, Φ0\Phi^{-}_{0} is defined by

Φ0:=18πσ0[(Bara)2(EaEb+BaBb)hab]𝑑A,\displaystyle\Phi^{-}_{0}:=\frac{1}{8\pi}\int_{\sigma_{0}}[(B_{a}r^{a})^{2}-(E_{a}E_{b}+B_{a}B_{b})h^{ab}]dA, (32)

where habh_{ab} is the induced metric of σ0\sigma_{0}.

Proof.

On σ0\sigma_{0}, the Lie derivative of the mean curvature k¯\bar{k} along nan^{a} is given by Yoshino:2020-1

£¯n(3)k¯=12R(2)8πPr+12(2kkL+k2kabkab),\displaystyle{}^{(3)}\bar{\mbox{\pounds}}_{n}\bar{k}=-\frac{1}{2}{}^{(2)}R-8\pi P_{r}+\frac{1}{2}(2kk_{\rm L}+k^{2}-k_{ab}k^{ab}), (33)

where R(2){}^{(2)}R is the Ricci scalar of σ0\sigma_{0}, kabk_{ab} is the extrinsic curvature of σ0\sigma_{0} in Σ\Sigma, that is, kab=(1/2)£r(3)habk_{ab}=(1/2){}^{(3)}\mbox{\pounds}_{r}h_{ab}, kk is its trace, kLk_{\rm L} is the largest value of the eigenvalues of kabk_{ab} and PrP_{r} is the radial pressure defined in Eq. (8). Then, the condition £¯n(3)k¯|σ00{}^{(3)}\bar{\mbox{\pounds}}_{n}\bar{k}|_{\sigma_{0}}\leq 0 gives us

34k212R(2)+(EaEb+BaBb)hab(Eara)2(Bara)2+8πPr(m)\frac{3}{4}k^{2}\ \leq\ \frac{1}{2}{}^{(2)}R+(E_{a}E_{b}+B_{a}B_{b})h^{ab}-(E_{a}r^{a})^{2}-(B_{a}r^{a})^{2}+8\pi P_{r}^{(m)} (34)

on σ0\sigma_{0}, where we used Eq. (8) and the inequality for a convex DTTS Yoshino:2020-1 , 2kkL+k2kabkab(3/2)k22kk_{\rm L}+k^{2}-k_{ab}k^{ab}\geq(3/2)k^{2}. With the condition Pr(m)0P_{r}^{(m)}\leq 0, the integration of the above over σ0\sigma_{0} gives us 333Note that the origin for Φ0\Phi^{-}_{0} in the case of a DTTS is different from that for Φ0+\Phi^{+}_{0} in the case of an LTS. As we can see from the derivation, Φ0±\Phi^{\pm}_{0} come from the energy density and the radial pressure of electromagnetic fields in the cases of the LTS and the DTTS, respectively. This is essential reason why we have different results for the LTS and the DTTS.

34σ0k2𝑑A\displaystyle\frac{3}{4}\int_{\sigma_{0}}k^{2}dA \displaystyle\leq 12σ0R(2)𝑑A8πΦ0σ0(Eara)2𝑑A\displaystyle\frac{1}{2}\int_{\sigma_{0}}{}^{(2)}RdA-8\pi\Phi_{0}^{-}-\int_{\sigma_{0}}(E_{a}r^{a})^{2}dA (35)
\displaystyle\leq 4π(4πq)2A08πΦ0,\displaystyle 4\pi-\frac{(4\pi q)^{2}}{A_{0}}-8\pi\Phi^{-}_{0},

where we used the Gauss-Bonnet theorem, the Cauchy-Schwarz inequality and Gauss’ law as in the proof for Theorem 1.

We now consider the inverse mean curvature flow in which the foliation is given by y=constanty=\mathrm{constant} surfaces. Similarly to the proof for Theorem 1, the surfaces y=0y=0 and y=y=\infty are set to be the DTTS and a sphere at spacelike infinity, respectively. With the same procedure, one can derive the inequality of Eq. (20) again. Then, the integration of the inequality of Eq. (20) over the range 0y<0\leq y<\infty shows us

m\displaystyle m \displaystyle\geq A01/24π1/2(1116πσ0k2𝑑A)+π1/2q2A01/2\displaystyle\frac{A_{0}^{1/2}}{4\pi^{1/2}}\left(1-\frac{1}{16\pi}\int_{\sigma_{0}}k^{2}dA\right)+\frac{\pi^{1/2}q^{2}}{A_{0}^{1/2}} (36)
\displaystyle\geq A01/26π1/2+43π1/2q2A01/2+A01/26π1/2Φ0,\displaystyle\frac{A_{0}^{1/2}}{6\pi^{1/2}}+\frac{4}{3}\frac{\pi^{1/2}q^{2}}{A_{0}^{1/2}}+\frac{A_{0}^{1/2}}{6\pi^{1/2}}{\Phi}^{-}_{0},

where we used the inequality of Eq. (35) at the last step. Then, we arrive at the inequality given by Eq. (31).

There are four remarks. Similarly to Theorem 1, in general, the minimum value of the right-hand side of the inequality of Eq. (31) implies us the lower bound for mm. However, unlike Φ0+\Phi_{0}^{+} in the case of an LTS, the quantity Φ0\Phi^{-}_{0} does not have a definite signature. For Φ0<1\Phi_{0}^{-}<-1, there is no such restriction for mm, whereas, for Φ01\Phi_{0}^{-}\geq-1, mm has a lower bound,

m2231+Φ0|q|.\displaystyle m\geq\frac{2{\sqrt{2}}}{3}{\sqrt{1+\Phi_{0}^{-}}}\ |q|. (37)

Next, under the condition of the inequality of Eq. (37), a short rearrangement of the inequality of Eq. (36) gives

4πrDTTS2A04πrDTTS+2,\displaystyle 4\pi r_{{\rm DTTS}-}^{2}\leq A_{0}\leq 4\pi r_{{\rm DTTS}+}^{2}, (38)

where

rDTTS±:=3m±9m28(1+Φ0)q22(1+Φ0).\displaystyle r_{{\rm DTTS}\pm}:=\frac{3m\pm\sqrt{9m^{2}-8\left(1+\Phi^{-}_{0}\right)q^{2}}}{2\left(1+\Phi^{-}_{0}\right)}. (39)

However, from the same reason to the remark in Theorem 1, we expect that the lower bound is not correct for a DTTS with multiple components, and just the inequality

A04πrDTTS+2\displaystyle A_{0}\leq 4\pi r_{{\rm DTTS}+}^{2} (40)

would hold in a general context. For a DTTS with a single component, the lower bound must hold true with the same physical reason as the one given in Sect. 3.2.

As a third remark, in a similar way to the case of an LTS, the obtained inequality depends on the electromagnetic field. Interestingly, if Φ0\Phi_{0}^{-} is negative, the contribution from Φ0\Phi_{0}^{-} makes the inequality weaker than the cases of Φ0=0\Phi^{-}_{0}=0. Furthermore, for the case of Φ01\Phi^{-}_{0}\leq-1, the upper bound disappears. Let us discuss the effect of Φ0\Phi_{0}^{-} physically. It is known that there are two kinds of pressure for magnetic fields. One is the negative pressure in the direction of magnetic field lines, called the magnetic tension. The other is repulsive interaction (i.e., positive pressure) between two neighboring magnetic field lines, called the magnetic pressure. A similar thing happens also to electric field lines (say, the electric tension and the electric pressure). We recall the formula for 8πPr8\pi P_{r}, Eq. (8). In that formula, (Eara)2(Bara)2-(E_{a}r^{a})^{2}-(B_{a}r^{a})^{2} is the contribution of the electric/magnetic tension, while (EaEb+BaBb)hab(E_{a}E_{b}+B_{a}B_{b})h^{ab} is the contribution of the electric/magnetic pressure. Then, Eq. (33) tells that the electric/magnetic tension makes the formation of a DTTS difficult, while the electric/magnetic pressure helps the formation of a DTTS. Therefore, in the presence of the electric/magnetic pressure, the area of a DTTS tends to be larger. This is the reason why the upper bound of the area of a DTTS becomes larger when Φ0\Phi_{0}^{-} is negative. Nevertheless, the negativity of Φ0\Phi_{0}^{-} would not change the situation so much in the following reason. If Φ0\Phi_{0}^{-} is negative, the upper bound for the DTTS becomes weaker and the DTTS can exist at farther outside. However, Φ0\Phi_{0}^{-} depends on the position of the DTTS and we naively expect that it is sharply decreasing according to the distance from the center, if the electromagnetic field is intrinsic to the compact object; namely, monopole or multi-pole fields. Therefore, when we take a farther surface, Φ0\Phi_{0}^{-} becomes immediately negligible. Then, the area of the DTTS cannot be large. On the other hand, Φ0\Phi_{0}^{-} could be large at some point by extrinsic effects, such as, external fields and/or dynamical generation of fields.

The final remark is on the relation between an LTS and a convex DTTS σ0\sigma_{0} on time-symmetric initial data. Recall Proposition 1 in Ref. Yoshino:2020-1 , that is, a convex DTTS with k>0k>0 in time-symmetric initial data is an LTS as well if ρ+Pr=0\rho+P_{r}=0 is satisfied on σ0\sigma_{0}. Since

8π(ρ(em)+Pr(em))=2(EaEb+BaBb)hab,\displaystyle 8\pi\left(\rho^{({\rm em})}+P_{r}^{({\rm em})}\right)=2(E_{a}E_{b}+B_{a}B_{b})h^{ab}, (41)

the presence of Φ0±\Phi^{\pm}_{0} disturbs the equivalence between an LTS and a convex DTTS in general. This feature is reflected in the two inequalities obtained in this paper.

5 Numerical examination of the Majumdar-Papapetrou spacetime

In Sect. 4.2, we have obtained the Penrose-like inequality for a DTTS. There, the quantity Φ0\Phi^{-}_{0} appears, and this quantity depends on the configuration of electromagnetic fields. The purpose of this section is to examine the values of Φ0\Phi^{-}_{0} in an explicit example. Specifically, in our previous paper Yoshino:2020-1 , we numerically solved for marginally DTTSs in systems of two equal-mass black holes adopting the Majumdar-Papapetrou solution. We revisit this problem from the viewpoint of our current work.

A Majumdar-Papapetrou spacetime is a static electrovacuum spacetime. The metric is

ds2=U2dt2+U2[dr~2+r~2(dθ~2+sin2θ~dϕ~2)],ds^{2}=-U^{-2}dt^{2}+U^{2}[d\tilde{r}^{2}+\tilde{r}^{2}(d\tilde{\theta}^{2}+\sin^{2}\tilde{\theta}d\tilde{\phi}^{2})], (42)

where the spatial structure is conformally flat, and we span the spherical-polar coordinates here. The electromagnetic four-potential is

Aa=U1(dt)a=na.A_{a}=U^{-1}(dt)_{a}=-n_{a}. (43)

Any solution to the Laplace equation ¯2U=0\bar{\nabla}^{2}U=0 gives an exact solution, where ¯2\bar{\nabla}^{2} is the flat space Laplacian. In this situation, EaE^{a} and BaB^{a} are calculated as

Ea=DaUU,Ba=0.E^{a}=-\frac{D^{a}U}{U},\qquad B^{a}=0. (44)

Setting E2:=habEaEbE_{\parallel}^{2}:=h_{ab}E^{a}E^{b}, we have 8πΦ0=σ0E2𝑑A8\pi\Phi^{-}_{0}=-\int_{\sigma_{0}}E_{\parallel}^{2}dA. Therefore, the value of Φ0\Phi^{-}_{0} is non-positive.

In our previous paper Yoshino:2020-1 , we chose the solution

U=1+m/2r~2+z022r~z0cosθ~+m/2r~2+z02+2r~z0cosθ~U=1+\frac{m/2}{\sqrt{\tilde{r}^{2}+z_{0}^{2}-2\tilde{r}z_{0}\cos\tilde{\theta}}}+\frac{m/2}{\sqrt{\tilde{r}^{2}+z_{0}^{2}+2\tilde{r}z_{0}\cos\tilde{\theta}}} (45)

that represents the system in which two extremal black holes with the same charge are located with the coordinate distance 2z02z_{0}. Then, assuming the functional form r~=h(θ~)\tilde{r}=h(\tilde{\theta}), we numerically solved for a marginally DTTS that surrounds both black holes for each value of z0z_{0}. The solution was found in the range 0z0/m0.793530\leq z_{0}/m\leq 0.79353. We refer readers to our previous paper Yoshino:2020-1 for explicit shapes of the obtained solutions.

Refer to caption
Refer to caption
Figure 1: The quantities related to Φ0\Phi^{-}_{0} of a marginally DTTS in a Majumdar-Papapetrou two-black-hole spacetime. Left panel: The value of E2E_{\parallel}^{2} as a function of θ~\tilde{\theta} on a marginally DTTS for z0/m=0.25z_{0}/m=0.25, 0.3750.375, 0.50.5, 0.6250.625, 0.750.75, and 0.793530.79353. Right panel: The value of Φ0\Phi^{-}_{0} for a marginally DTTS as a function of z0/mz_{0}/m.

We examine the value of Φ0\Phi^{-}_{0}. The left panel of Fig. 1 presents the behavior of E2E_{\parallel}^{2} as a function of θ~\tilde{\theta} on a marginally DTTS for z0/m=0.25z_{0}/m=0.25, 0.3750.375, 0.50.5, 0.6250.625, 0.750.75, and 0.793530.79353. The value of E2E_{\parallel}^{2} is generally nonzero, but is less than 10310^{-3}. The right panel of Fig. 1 plots the value of Φ0\Phi^{-}_{0} as a function of z0z_{0}. It is negative and its absolute value is less than 5×1045\times 10^{-4}.

Refer to caption
Figure 2: The area AA of a common marginally DTTS in a Majumdar-Papapetrou spacetime as a function of z0/mz_{0}/m. The cases of two kinds of normalizations are shown. One is A/4πrDTTS+2A/4\pi r_{\rm DTTS+}^{2} and the other is A/4πrp2A/4\pi r_{p}^{2}. See text for details.

Figure 2 shows the relation between the area AA of the marginally DTTS and z0z_{0}. We normalize the value of AA in two ways: One is A/4πrDTTS+2A/4\pi r_{\rm DTTS+}^{2} (a red solid curve), where rDTTS+r_{\rm DTTS+} is defined in Eq. (39), and the other is A/4πrp2A/4\pi r_{p}^{2} (a green dotted curve), where 4πrp24\pi r_{p}^{2} is the area of a photon sphere with the same mass and charge [see Eq. (10)]. Because the value of Φ0\Phi^{-}_{0} is small, the difference is scarcely visible. Both of these values are in agreement with the Penrose-like inequalities.

Refer to caption
Figure 3: The marginally DTTS for z0/m=0.79353z_{0}/m=0.79353 (the blue curve) and the contour surfaces of UU in the Majumdar-Papapetrou two-black-hole spacetime (gray and red curves). The red curve depicts the contour surface of U=2.14U=2.14.

The reason why the value of Φ0\Phi^{-}_{0} is so small is that the electric field is approximately perpendicular to the marginally DTTS. From the formula for EaE^{a} given in Eq. (44), this means that the marginally DTTS approximately coincides with a contour surface of UU. Figure 3 confirms this feature for z0/m=0.79353z_{0}/m=0.79353. Here, the blue curve depicts the marginally DTTS, and the red curve shows the contour surface of U=2.14U=2.14. They agree well.

The lesson from this numerical experiment is that if a spacetime is static, the quantity Φ0\Phi^{-}_{0} is small and does not play an important role in the Penrose-like inequality for a DTTS. Of course, it is expected that the absolute value of Φ0\Phi^{-}_{0} may become large if dynamical situations are considered. For example, if two black holes have opposite charges, the contribution of the electric pressure would become important, although such a situation is more difficult to study. Exploring such issues is left as a remaining problem.

6 Summary and discussion

In this paper, we have examined the properties of LTSs and DTTSs for Einstein-Maxwell systems, particularly focusing on the derivation of Penrose-like inequalities on their area. Similarly to the Riemannian Penrose inequality for charged cases, the electric charge comes into the inequalities, but there are additional contributions from the density or the pressure/tension of electromagnetic fields in general. This is a rather interesting result because one naively expects that the upper bound for the area of an LTS and a DTTS is that of the photon sphere 4πrp24\pi r_{p}^{2}, where rpr_{p} is the radius of an unstable circular orbit of a photon in the Reissner-Nordström spacetime given in Eq. (10). For an LTS, we have a tighter inequality than the naive one. For a DTTS, the obtained inequality can become both stronger and weaker depending on the configuration of electromagnetic fields. We have numerically examined the value of the correction term, represented by Φ0\Phi^{-}_{0} in Eq. (32), for a Majumdar-Papapetrou two-black-hole spacetime. Although the correction term makes the inequality weaker, we have checked that the value of Φ0\Phi^{-}_{0} is very small in that situation.

Up to now, we have assumed that the magnetic charge qmq_{\rm m} vanishes. Here, we consider what happens to our inequalities when nonzero qmq_{\rm m} is present. We first consider the case of an LTS. Since the total squared charge q2:=qe2+qm2q^{2}:=q_{\rm e}^{2}+q_{\rm m}^{2} appears in the metric of the Reissner-Nordström solution when both the electric and magnetic fields are present, we would like to present the Penrose-like inequality in terms of q2q^{2}. For this reason, we have to use the inequality

S0[(Eara)2+(Bara)2]𝑑AA01[(4πqe)2+(4πqm)2]=(4πq)2A0\int_{S_{0}}\left[(E_{a}r^{a})^{2}+(B_{a}r^{a})^{2}\right]dA\geq A_{0}^{-1}\left[(4\pi q_{\rm e})^{2}+(4\pi q_{\rm m})^{2}\right]=\frac{(4\pi q)^{2}}{A_{0}} (46)

in the calculations of Eqs. (16) and (17). As a result, we must introduce the quantity

Φ0:=18πS0[(EaEb+BaBb)hab]𝑑A\Phi_{0}:=\frac{1}{8\pi}\int_{S_{0}}\left[(E_{a}E_{b}+B_{a}B_{b})h^{ab}\right]dA (47)

instead of Φ0+\Phi^{+}_{0} of Eq. (14). The resultant inequality is the one of Eq. (13) but Φ0+\Phi^{+}_{0} being replaced by Φ0\Phi_{0}. Next, we consider the case of a DTTS. Similarly to the case of an LTS, we must use the inequality of Eq. (46) (but S0S_{0} being replaced by σ0\sigma_{0}) in the calculations of Eqs. (34) and (35). As a result, instead of Φ0\Phi_{0}^{-}, we must introduce Φ0=(1/8π)σ0[(EaEb+BaBb)hab]𝑑A-\Phi_{0}=-(1/8\pi)\int_{\sigma_{0}}\left[(E_{a}E_{b}+B_{a}B_{b})h^{ab}\right]dA. The resultant inequality is the same as the one of Eq. (31), but Φ0\Phi^{-}_{0} being replaced by Φ0-\Phi_{0}. These results are summarized as follows:

Corollary 1.

In the presence of nonzero qmq_{\rm m}, Theorems 1 and 2 hold by changing from Φ0±\Phi^{\pm}_{0} to ±Φ0\pm\Phi_{0}, where Φ0\Phi_{0} is defined by Eq. (47).

In the main article of this paper, we have not considered a TTS for the static and stationary spacetimes defined in our previous paper Yoshino:2017 . We note that the concepts of a TTS and a DTTS are related but independent of each other in the sense that no inclusion relationship can be found Yoshino:2020-1 . In Appendix A, we present a theorem on the Penrose-like inequality for a TTS in a static spacetime, which is very similar to Theorem 2.

Throughout this study, we have not used the property of Maxwell’s equations except for Gauss’ law. The information from Maxwell’s equation may further restrict the properties of LTSs, DTTSs, and TTSs, especially for static/stationary spacetimes with static/stationary electromagnetic fields.

\ack

T. S. and K. I. are supported by Grant-Aid for Scientific Research from Ministry of Education, Science, Sports and Culture of Japan (No. 17H01091). H.Y. is supported by the Grant-in-Aid for Scientific Research (C) (No. JP18K03654) from Japan Society for the Promotion of Science (JSPS). The work of H.Y. is partly supported by Osaka City University Advanced Mathematical Institute (MEXT Joint Usage/Research Center on Mathematics and Theoretical Physics).

Appendix A Transversely trapping surface

The four of the present authors also proposed the concept of a TTS Yoshino:2017 . This concept is applicable only to static or stationary spacetimes. The definition is as follows :

Definition 3.

A static/stationary timelike hypersurface SS is a transversely trapping surface (TTS) if and only if arbitrary light rays emitted in arbitrary tangential directions of SS from arbitrary points of SS propagate on SS or toward the inside region of SS.

The necessary and sufficient condition for a surface SS to be a TTS (the TTS condition hereafter) is expressed as K¯abkakb0\bar{K}_{ab}k^{a}k^{b}\leq 0, where K¯ab\bar{K}_{ab} is the extrinsic curvature of SS and kak^{a} are arbitrary null vectors tangent to SS. For a static spacetime, there is the Killing time coordinate tt whose basis vector tat^{a} is orthogonal to the t=constantt=\mathrm{constant} hypersurface, denoted by Σ\Sigma. The lapse function α\alpha is defined by ta=αnat^{a}=\alpha n^{a}, where nan^{a} is a future-directed unit normal to Σ\Sigma. We denote the two-dimensional section of SS and Σ\Sigma by σ0\sigma_{0}. As shown in Eq. (15) of Ref. Yoshino:2017 , the TTS condition is reexpressed in terms of α\alpha as

kLraDaαα,k_{\rm L}\leq\frac{r^{a}D_{a}\alpha}{\alpha}, (48)

where kLk_{\rm L} is the largest value among the two eigenvalues of the extrinsic curvature kabk_{ab} of σ0\sigma_{0} in Σ\Sigma, rar_{a} is a spacelike unit normal to SS, and DaD_{a} is the covariant derivative with respect to Σ\Sigma. In this situation, it is possible to derive the relation

R(2)=16πPr+2α𝒟2α+2kraDaαα+k2kabkab.{}^{(2)}R=-16\pi P_{r}+\frac{2}{\alpha}\mathcal{D}^{2}\alpha+2k\frac{r^{a}D_{a}\alpha}{\alpha}+k^{2}-k_{ab}k^{ab}. (49)

For a convex TTS, we can derive the inequality

2kraDaαα+k2kabkab32k2.2k\frac{r^{a}D_{a}\alpha}{\alpha}+k^{2}-k_{ab}k^{ab}\geq\frac{3}{2}k^{2}. (50)

These relations are presented as Eqs. (24) and (25) in Ref. Yoshino:2017 . We now consider the Einstein-Maxwell system. Using Eq. (8), we find

34k212R(2)+(EaEb+BaBb)hab(Eara)2(Bara)2+8πPr(m)1α𝒟2α.\frac{3}{4}k^{2}\ \leq\ \frac{1}{2}{}^{(2)}R+(E_{a}E_{b}+B_{a}B_{b})h^{ab}-(E_{a}r^{a})^{2}-(B_{a}r^{a})^{2}+8\pi P_{r}^{(m)}-\frac{1}{\alpha}\mathcal{D}^{2}\alpha. (51)

Compare this inequality with the one of Eq. (34). Integrating over σ0\sigma_{0}, we obtain exactly the same inequality as the one of Eq. (35). Therefore, a TTS in a static spacetime satisfies the same inequality as the Penrose-like inequality for a DTTS in time-symmetric initial data. This result is summarized as the following theorem:

Theorem 3.

The static time cross section of a convex TTS, σ0\sigma_{0}, in an asymptotically flat static spacetime has topology S2S^{2} and its areal radius r0=A0/4πr_{0}=\sqrt{A_{0}/4\pi} satisfies the inequality of Eq. (31) [with Φ0\Phi^{-}_{0} defined in Eq. (32)] if Pr(m)<0P_{r}^{(m)}<0 holds on σ0\sigma_{0}, k>0k>0 at least at one point on σ0\sigma_{0}, and ρ(m)0\rho^{(m)}\geq 0 in the outside region.

As remarked in the final section, this theorem applies to the case that the magnetic charge qmq_{\rm m} is zero. When qmq_{\rm m} is nonzero, Φ0\Phi^{-}_{0} must be replaced by Φ0-\Phi_{0}, where Φ0\Phi_{0} is given in Eq. (47).

References

  • (1) C. M. Claudel, K. S. Virbhadra and G. F. R. Ellis, J. Math. Phys.  42, 818 (2001). [gr-qc/0005050].
  • (2) H. Yoshino, Phys. Rev. D 95, 044047 (2017) [arXiv:1607.07133 [gr-qc]].
  • (3) V. Cardoso, A. S. Miranda, E. Berti, H. Witek and V. T. Zanchin, Phys. Rev. D 79, 064016 (2009) [arXiv:0812.1806 [hep-th]].
  • (4) K. Akiyama et al. [Event Horizon Telescope Collaboration], Astrophys. J.  875, L1 (2019) [arXiv:1906.11238 [astro-ph.GA]].
  • (5) K. S. Virbhadra and G. F. R. Ellis, Phys. Rev. D 62, 084003 (2000). [astro-ph/9904193].
  • (6) T. Shiromizu, Y. Tomikawa, K. Izumi and H. Yoshino, Prog. Theor. Exp. Phys. 2017, 033E01 (2017) [arXiv:1701.00564 [gr-qc]].
  • (7) H. Yoshino, K. Izumi, T. Shiromizu and Y. Tomikawa, Prog. Theor. Exp. Phys.  2017, 063E01 (2017) [arXiv:1704.04637 [gr-qc]].
  • (8) H. Yoshino, K. Izumi, T. Shiromizu and Y. Tomikawa, Prog. Theor. Exp. Phys. 2020 023E02 (2020) [arXiv:1909.08420 [gr-qc]].
  • (9) H. Yoshino, K. Izumi, T. Shiromizu and Y. Tomikawa, Prog. Theor. Exp. Phys. 2020 053E01 (2020) [arXiv:1911.09893 [gr-qc]].
  • (10) D. V. Gal’tsov and K. V. Kobialko, Phys. Rev. D 99, 084043 (2019) [arXiv:1901.02785 [gr-qc]].
  • (11) R. Penrose, Annals N. Y. Acad. Sci.  224, 125 (1973).
  • (12) S. Hod, Phys. Lett. B 727, 345-348 (2013). [arXiv:1701.06587 [gr-qc]].
  • (13) P. S. Jang, Phys. Rev. D 20, 834 (1979).
  • (14) R. Geroch, Ann. N. Y. Acad. Sci. 224, 108 (1973).
  • (15) P. S. Jang and R. M. Wald, J. Math. Phys. 18, 41 (1977).
  • (16) G. Weinstein and S. Yamada, Commun. Math. Phys. 257, 703-723 (2005) [arXiv:math/0405602 [math.DG]].
  • (17) M. Khuri, G. Weinstein and S. Yamada, Contemp. Math. 653, 219-226 (2015) [arXiv:1308.3771 [gr-qc]].
  • (18) M. Khuri, G. Weinstein and S. Yamada, J. Diff. Geom. 106, no.3, 451-498 (2017) [arXiv:1409.3271 [gr-qc]].