Non-Singular Model of Magnetized Black Hole Based on Nonlinear Electrodynamics

universe

Article Non-Singular Model of Magnetized Black Hole Based on Nonlinear Electrodynamics

Sergey I. Kruglov 1,2

1 Department of Physics, University of Toronto, 60 St. Georges St., Toronto, ON M5S 1A7, Canada; [email protected] 2 Department of Chemical and Physical Sciences, University of Toronto Mississauga, 3359 Mississauga Road North, Mississauga, ON L5L 1C6, Canada

Received: 15 October 2019; Accepted: 10 December 2019; Published: 13 December 2019

Abstract: A new modified Hayward metric of magnetically charged non-singular black hole spacetime in the framework of nonlinear electrodynamics is constructed. When the fundamental length introduced, characterising quantum gravity effects, vanishes, one comes to the general relativity coupled with the Bronnikov model of nonlinear electrodynamics. The metric can have one (an extreme) horizon, two horizons of black holes, or no horizons corresponding to the particle-like solution. Corrections to the Reissner–Nordström solution are found as the radius approaches infinity. As r → 0 the metric has a de Sitter core showing the absence of singularities, the asymptotic of the Ricci and Kretschmann scalars are obtained and they are finite everywhere. The thermodynamics of black holes, by calculating the Hawking temperature and the heat capacity, is studied. It is demonstrated that phase transitions take place when the Hawking temperature possesses the maximum. Black holes are thermodynamically stable at some range of parameters.

Keywords: modiﬁed Hayward metric; magnetically charged black hole; nonlinear electrodynamics; thermodynamics

1. Introduction It is well-known that General Relativity (GR) is ultraviolet (UV) incomplete. In addition, there is a problem of singularities in the classical Einstein theory of gravity. Thus, solutions of the Einstein equations for charged (the Reisner–Nordström metric) black holes (BHs) have curvature singularities in the center (r = 0). Therefore, GR should be modified when the curvature is large. There are some attempts to overcome problems in the classical Einstein theory of gravity. So, if one adds curvature terms of the higher order or terms with higher derivatives, the UV behaviour of the Einstein gravity will be improved [1,2]. But the price for this is the existence of ghosts (non-physical degrees of freedom). A ghost free modification of the GR, which is UV-complete, was considered in References [3–5] but such a theory is non-local and has an infinite number of derivatives. Because the fundamental quantum gravity theory (UV-complete) is absent, some phenomenological models can be useful for solving problems of singularities. Following References [6–8] (see also Reference [9]), we assume that there is a critical energy µ and the corresponding length l = µ−1 in such a way that the metric is modified when the spacetime curvature is in the order of l−2. The length scale, characterising quantum gravity effects, is smaller than l and one may use the classical metric gµν. In addition, we suppose that the limiting curvature condition R ≤ cl−2 (R is one of the curvature invariants, c is dimensionless constant depending on the curvature invariant) is satisfied [6–8]. A simple metric satisfying the above conditions was proposed by Hayward [10] for a neutral BH. This is the phenomenological approach that we explore here.

Universe 2019, 5, 225; doi:10.3390/universe5120225 www.mdpi.com/journal/universe Universe 2019, 5, 225 2 of 9

The first pioneering work representing a regular BH in GR is in Reference [11]. It was shown in Reference [12] that the Bardeen model can describe the gravitational field of a nonlinear magnetic monopole. In References [13,14] the regular electrically charged BH solution in GR was presented, where the source is a nonlinear electrodynamics (NED) field satisfying the weak energy condition. It worth noting that in accordance with Reference [15], regular electric solutions with the Maxwell weak-field limit can be described only by different NED theories in different parts of spacetime. Thus, there is a significant shortcoming in the models of References [13,14]. In this paper, we consider the spherically symmetric non-singular model of the magnetically charged BH based on NED. In some NED, the electric field in the center of point-like charges is finite [16–20] and the self-energy of charges is finite unlike classical electrodynamics. It is worth mentioning that quantum corrections to Maxwell’s electrodynamics, within QED , lead to NED [21]. The universe inflation also can be explained in the framework of the GR coupled with NED [22–29]. Here, with the help of the modified Hayward metric, we study regular magnetically charged BH solutions within NED considered in Reference [15]. The BH thermodynamics and phase transitions are investigated. In References [30,31] the authors also considered BH solutions with the modified Hayward metric based on NED proposed in References [32,33], respectively. The thermodynamics for a magnetically charged regular BH, which comes from the action of GR and NED, was investigated in Reference [34]. These authors also used NED proposed in Reference [15]. A similar study was performed in References [35,36]. The work in Reference [37] analyzes the minimal model proposed by Hayward for an uncharged BH within GR. The authors introduced an anisotropic fluid and postulated the expressions for the energy density and pressure but the Lagrangian corresponding to the stress tensor was not obtained. In the present study we use the NED of Reference [15], explore a phenomenological extension of GR by introducing a fundamental length l using the modified Hayward metric, and investigate the magnetically charged BH. This paper is organized as follows. In Section2 the modified Hayward metric is studied and we obtain the asymptotic of the metric and mass functions as r → 0 and r → ∞. Corrections to the Reissner–Nordström (RN) solution are found. The asymptotic of the Ricci and Kretschmann scalars are calculated and we show that curvature singularities are absent. In Section3 we calculate the Hawking temperature and the heat capacity of BHs. We demonstrate that the second-order phase transitions occur. It is shown that in some range of parameters BHs are stable. Section4 is a conclusion.

2. A Regular Magnetized BH Solution To describe the magnetically charged BH solution we consider the Lagrangian density of NED [15]:

F L = − p , (1) cosh2 4 |βF|

µν 2 2 where F = (1/4)FµνF = (B − E )/2 and the ﬁeld tensor is Fµν = ∂µ Aν − ∂ν Aµ. The parameter β in Equation (1) is positive and it possesses the dimension of (length)4. At the weak ﬁeld limit the Lagrangian density (1) becomes: L → −F βF 1, (2) that is, the correspondence principle holds. We will derive the metric function representing the static magnetic regular BH. Let us consider the spherically symmetric line element, which is given by:

1 ds2 = − f (r)dt2 + dr2 + r2(dϑ2 + sin2 ϑdφ2). (3) f (r)

The Hayward metric function [10] is given by:

2GMr2 f (r) = 1 − , (4) r3 + 2GMl2 Universe 2019, 5, 225 3 of 9 where G is the Newton constant, M = constant and l is the fundamental length. We interpret this metric in the framework of an extension of GR for an uncharged source and replace the Schwarzschild metric. It should be noted that in GR the metric function (4) may be obtained as a solution within NED with a nonzero magnetic charge [38]. But in this case the NED Lagrangian is ill-deﬁned. At the weak-ﬁeld limit the NED Lagrangian does not approach to the Maxwell Lagrangian. One can consider and investigate other geometries of the horizon in Equation (3). At l = 0 we come to the Schwarzschild metric of a BH which is a solution to Einstein’s equation without sources. Now we suppose that the BH is magnetically charged. Then the mass function of a BH varies with r and is:

Z r Z ∞ Z ∞ 2 2 2 M(r) = m0 + ρ(r)r dr = m0 + ρ(r)r dr − ρ(r)r dr, (5) 0 0 r

R ∞ 2 where m0 is the Schwarzschild mass, ρ(r) is the magnetic energy density and mM = 0 ρ(r)r dr is the magnetic mass of the BH. In Reference [15] the mass m0, which can be considered as a constant of integration, was not introduced. But the case m0 6= 0 allows us to consider the uncharged BH when the charge q = 0. Indeed, if q = 0 (ρ(r) = 0) in Equation (5), the mass function M becomes constant and we come to the Hayward metric function (4). At E = 0 the magnetic energy density, corresponding to Equation (1), is given by: F ρ(r) = −L = p , (6) cosh2 4 |βF| where F = B2/2 = q2/(2r4), and q is a magnetic charge. Then the mass function (5) becomes: √ ! q3/2 β1/4 q M(r) = m0 + mM − tanh , (7) 23/4β1/4 21/4r where the BH magnetic mass is given by [36]:

Z ∞ 3/2 2 q mM = ρ(r)r dr = . (8) 0 23/4β1/4

The total BH mass is M(∞) ≡ m = m0 + mM. Then the metric function corresponding to a charged BH is: 2GM(r)r2 f (r) = 1 − , (9) r3 + 2GM(r)l2 where M(r) is given by (7). √ For a convenience we introduce the dimensionless parameter x = 21/4r/(β1/4 q). Then from Equations (7)–(9) one obtains the metric function:

Ax2g(x) f (x) = 1 − , (10) x3 + Bg(x) where: √ 2Gq 2Gl2 23/4β1/4m 1 A = , B = , C = 0 , g(x) = C + 1 − tanh . (11) pβ β q3/2 x From Equations (10) and (11) we ﬁnd the asymptotic of the metric function as r → ∞ and r → 0:

p ! 2Gm Gq2 G βq3 f (r) = 1 − + − √ − 4Gl2m2 + O(r−5) r → ∞, (12) r r2 r4 3 2

r2 r5 ( ) = − + + O( 6) → f r 1 2 4 r r 0. (13) l 2Gm0l Universe 2019, 5, 225 4 of 9

Equation (12) shows the corrections to the RN solution that are in the order of O(r−4). At l = 0 and m0 = 0 (when the total BH mass is the magnetic mass) Equation (12) is converted into the equation obtained in Reference [36]. As r → ∞ we have f (∞) = 1, and the spacetime becomes ﬂat. According to Equation (13) limr→0 f (r) = 1. Thus, the spacetime has a smooth de Sitter core and the BH is regular. If β = 0, l = 0 one has the RN solution. The plot of the function f (x) is depicted in Figure1. This plot is typical for any regular solution described by the metric (3) (see, e.g., Reference [13]).

0.8

0.6

0.4 A=2 A=3.93 0.2 f(x) A=6

−0.2

−0.4

−0.6 0 5 10 15 x

Figure 1. The plot of the function f (x) for B = 1 and C = 0 (m0 = 0). The dashed-dotted line corresponds to A = 6, the solid line corresponds to A = 2 and the dashed line corresponds to A = 3.93.

In accordance with Figure1 at A < 3.93 (B = 1, C = 0) horizons are absent and we have particle-like solution. At A ≈ 3.93 the horizons shrink into one horizon (the extreme solution). If > 3.93, we have two horizons of a BH. The horizon radii xh are roots of the equation f (xh) = 0. 1/4 1/4√ From Equation (10), at B = 1, C = 0 one ﬁnds the inner x− (x− = 2 r−/(β q)) and outer x+ 1/4 1/4√ (x+ = 2 r+/(β q)) horizon radii of the BH that are given in Table1.

Table 1. The BH inner and outer horizon radii (B = 1, C = 0).

A 4 5 6 7 8 9 10 15

x− 1.56 1.04 0.88 0.78 0.72 0.68 0.64 0.54 x+ 2.11 3.58 4.71 5.78 6.82 7.84 8.87 13.92

The asymptotic of the Ricci and Kretschmann scalars can be obtained from the relations:

4 f (r) − 1 R = − f 00(r) − f 0(r) − 2 , (14) r r2

2 f 0(r) 2 4 ( f (r) − 1)2 K = f 002(r) + + . (15) r r4 From Equations (14) and (15) we ﬁnd:

12 21r3 = − + O( 4) → R 2 4 r r 0, (16) l Gm0l Universe 2019, 5, 225 5 of 9

p ! 86Gm 12Gq2 34G βq3 R = − + − √ − 4Gl2m2 + O(r−7) r → ∞, (17) r3 r4 r6 3 2

24 84r3 = − + O( 4) → K 4 6 r r 0, (18) l Gm0l 48G2m2 96G2mq2 K = − + O(r−8) r → ∞. (19) r6 r7 As r → ∞ the Ricci and Kretschmann scalars vanish and the spacetime becomes ﬂat. Equations (16)–(19) indicate that solutions obtained are regular.

3. Thermodynamics and Phase Transitions Let us study the thermal stability of magnetized BHs and the possible phase transitions. The Hawking temperature is given by [39]:

κ f 0(r ) T = = h , (20) H 2π 4π where κ is the surface gravity and rh is the horizon radius. Making use of Equations (10) and (20) we obtain the Hawking temperature: √ qG T = −2x g(x ) H 5/4 3/4 3 h h 2 πβ (xh + Bg(xh))

4 2 1 g(xh)(B + 3x cosh (1/xh)) − + h . (21) 2 3 2 cosh (1/xh) (xh + Bg(xh)) cosh (1/xh)

The horizon radii rh (and xh) are deﬁned as roots of the equation f (rh) = 0 (and f (xh) = 0). From Equation (10) we obtain: 3 Gq x + Bg(xh) = √h . (22) p 2 β 2xhg(xh)

According to Equation (22) the horizon radius rh (and xh) depends on the magnetic charge q and p the model parameter β. Substituting Gq/ β from Equation (22) into Equation (21) we obtain the ﬁnal equation for the Hawking temperature:

1 T = −2x g(x ) H 7/4 √ 1/4 2 h h 2 π qβ g(xh)xh

4 2 1 g(xh)(B + 3x cosh (1/xh)) − + h . (23) 2 3 2 cosh (1/xh) (xh + Bg(xh)) cosh (1/xh) √ 1/4 The plot of the reduced Hawking temperature TH qβ is depicted in Figure2 for different values of the parameter B (or l). Universe 2019, 5, 225 6 of 9

0.02

0.01

−0.01 B=1 B=5

) B=10 h −0.02 (x H T −0.03 1/4 β 1/2

q −0.04

−0.05

−0.06

−0.07

−0.08 1 2 3 4 5 6 7 8 9 10 x h

√ 1/4 Figure 2. The plot of the function TH qβ vs. horizons xh for C = 0 (m0 = 0). The dashed-dotted line corresponds to B = 10, the solid line corresponds to B = 1 and the dashed line corresponds to B = 5.

It follows from Figure2 that at the bigger value of l (or B) the maximum of the Hawking temperature shifts to the bigger value of the horizon radius. The temperature curve has one extremum (maximum) resulting in one phase transition during the evaporation. Similar form of the temperature curve for a BH takes place in the models studied in [36]. The heat capacity at the constant charge is deﬁned by the relation [40]: ∂S TH∂S/∂rh 2πrhTH Cq = TH = = . (24) ∂TH q ∂TH/∂rh G∂TH/∂rh

2 The entropy obeys the Hawking area low S = A/(4G) = πrh/G. When the Hawking temperature has the extremum (∂TH/∂rh = 0) the heat capacity is singular and the second-order phase transition p takes place. In Figure3 the function GCq/( βq) vs. the horizon radius xh for different values of B for C = 0 (m0 = 0) is presented. Figure2 shows the similarity in the considered thermodynamics of our model and the thermodynamics of a neutral regular BH. It is seen from Figure3 that second-order phase transitions at the discontinuity points occur between negative and positive heat capacities. Figure2 shows that the maximum of the temperature (the Davies point), where the phase transitions take place, separates areas with increasing and decreasing BH temperatures. The unstable point between the positive and negative heat capacities has a discontinuity. The positive heat capacity corresponds to the late stage and the negative capacity to the early stage of the thermodynamics process. Thus, there is an interval of the horizon radius where the heat capacity is positive and the BH is stable. In accordance with Figures2 and3 , the heat capacity possesses a discontinuity at the horizon where the Hawking temperature possesses a maximum. When the parameter B is bigger, the second-order phase transition of the BH occurs at the larger value of the horizon radius rh (xh). For the large value of xh the BH is unstable (Cq < 0). Universe 2019, 5, 225 7 of 9

4 x 10 5

4 B=1 B=5 3 B=10

1 q) 1/2 β 0 /( q

GC −1

−2

−3

−4

−5 3 3.5 4 4.5 5 5.5 6 6.5 x h

p Figure 3. The plot of the function CqG/( βq) vs. xh for C = 0 (m0 = 0). The dashed-dotted line corresponds to B = 10, the solid line corresponds to B = 1 and the dashed line corresponds to B = 5.

4. Conclusions Solutions of a magnetically charged regular BH in the new model were obtained. This model is of interest because of its simplicity. We found the mass and metric functions possessing simple analytical structures. The BH can have one (an extreme horizon), two horizons (trapping horizons), or no horizons (untrapped surface, see Figure1). These plots are typical for any regular solution described by the metric (3). One can find the same behavior of the metric function in Reference [13] for another BH model. Corrections to the RN solution that are in the order of O(r−4) were obtained as the radius approaches to infinity. As r → ∞ the spacetime becomes flat. The model of a electrically charged BH [13] was formulated in so-called P-frame (the Hamiltonian framework). But the Lagrangian dynamics is specified in F-framework. It was shown in Reference [15] that the regular electric solution in P-frame corresponds to different Lagrangians in different parts of the space if the function P(F) is not monotonic. But this problem is absent for magnetic solutions. Thus, in the model [13] the problem of singularities was not solved completely [41]. We calculated the asymptotic of the Ricci and Kretschmann scalars as r → ∞ and r → 0 showing the absence of singularities. It was shown that the spacetime as r → 0 has a de Sitter core (the flatness at the center). Thus, the singularity at r = 0 has been smoothed out. Our solution describes nonsingular BH with the finite curvature everywhere including r = 0. The regular behavior of the Ricci and Kretschmann scalars also was observed in Reference [13]. The Hawking temperature and the heat capacity of the BH were found demonstrating that second-order phase transitions take place. It was shown that second-order phase transitions separate areas between negative and positive heat capacities and areas with increasing and decreasing BH temperatures. For small values of the horizon radius, depending on the parameters of the model, the Hawking temperature is negative (see Figure2). The thermodynamic stabilities of black holes were studied and it was shown that in some range of horizon radii the BHs are stable (the heat capacity is positive) (see Figure3). The long standing problem of singularities is solved in this model, at r = 0 curvature invariants are finite and the BH is regular. In addition, at the weak field limit NED (1) becomes the Maxwell electrodynamics, that is, the correspondence principle holds. It is worth noting that, in the Bardeen model, the correspondence principle breaks out [12].

Funding: This research received no external funding. Conﬂicts of Interest: The authors declare no conﬂict of interest. Universe 2019, 5, 225 8 of 9

References

1. Stelle, K.S. Classical gravity with higher derivatives. Gen. Rel. Grav. 1978, 9, 353–371. [CrossRef] 2. Stelle, K.S. Renormalization of Higher Derivative Quantum Gravity. Phys. Rev. D 1977, 16, 953–969. [CrossRef] 3. Biswas, T.; Gerwick, E.; Koivisto, T.; Mazumdar, A. Towards singularity and ghost free theories of gravity. Phys. Rev. Lett. 2012, 108, 031101. [CrossRef][PubMed] 4. Biswas, T.; Conroy, A.; Koshelev, A.S.; Mazumdar, A. Generalized ghost-free quadratic curvature gravity. Class. Quant. Grav. 2014, 31, 015022. [CrossRef] 5. Modesto, L. Super-renormalizable Quantum Gravity. Phys. Rev. D 2012, 86, 044005. [CrossRef] 6. Markov, M.A. Limiting density of matter as a universal law of nature. JETP Lett. 1982, 36, 266. 7. Markov, M.A. Problems of a perpetually oscillating universe. Ann. Phys. 1984, 155, 333. [CrossRef] 8. Polchinski, J. Decoupling Versus Excluded Volume or Return of the Giant Wormholes. Nucl. Phys. B 1989, 325, 619–630. [CrossRef] 9. Frolov, V.P. Notes on nonsingular models of black holes. Phys. Rev. D 2016, 94, 104056. [CrossRef] 10. Hayward, S.A. Formation and evaporation of nonsingular black holes. Phys. Rev. Lett. 2006, 96, 031103. [CrossRef] 11. Bardeen, J. Non-singular general-relativistic gravitational collapse. In Proceedings of the GR5, Tbilisi, U.S.S.R., 9–13 September 1968. 12. Ayon-Beato, E.; Garcia, A. The Bardeen model as a nonlinear magnetic monopole. Phys. Lett. B 2000, 493, 149. [CrossRef] 13. Ayon-Beato, E.; Garcia, A. Regular black hole in general relativity coupled to nonlinear electrodynamics. Phys. Rev. Lett. 1998, 80, 5056. [CrossRef] 14. Ayon-Beato, E.; Garcia, A. New regular black hole solution from nonlinear electrodynamics. Phys. Lett. B 1999, 464, 25. [CrossRef] 15. Bronnikov, K.A. Regular magnetic black holes and monopoles from nonlinear electrodynamics. Phys. Rev. D 2001, 63, 044005. [CrossRef] 16. Born, M.; Infeld, L. Foundations of the new field theory. Proc. Royal Soc. (Lond.) A 1934, 144, 425. [CrossRef] 17. Gitman, D.M.; Shabad, A.E. A note on “Electron self-energy in logarithmic electrodynamics” by P. Gaete and J. Helayël-Neto. Eur. Phys. J. C 2014, 74, 3186. [CrossRef] 18. Kruglov, S.I. A model of nonlinear electrodynamics. Ann. Phys. 2015, 353, 299. [CrossRef] 19. Kruglov, S.I. Nonlinear arcsin-electrodynamics. Ann. Phys. (Berlin) 2015, 527, 397. [CrossRef] 20. Kruglov, S.I. Modified nonlinear model of arcsin-electrodynamics. Commun. Theor. Phys. 2016, 66, 59. [CrossRef] 21. Heisenberg, W.; Euler, H. Folgerungen aus der diracschen theorie des positrons. Z. Physik 1936, 98, 714. [CrossRef] 22. García-Salcedo, R.; Breton, N. Born–Infeld cosmologies. Int. J. Mod. Phys. A 2000, 15, 4341. [CrossRef] 23. Camara, C.S.; de Garcia Maia, M.R.; Carvalho, J.C.; Lima, J.A.S. Nonsingular FRW cosmology and nonlinear electrodynamics. Phys. Rev. D 2004, 69, 123504. [CrossRef] 24. Novello, M.; Bergliaffa, S.E.P.; Salim, J.M. Nonlinear electrodynamics and the acceleration of the universe. Phys. Rev. D 2004, 69, 127301. [CrossRef] 25. Novello, M.; Goulart, E.; Salim, J.M.; Bergliaffa, S.E.P. Cosmological effects of nonlinear electrodynamics. Class. Quant. Grav. 2007, 24, 3021. [CrossRef] 26. Vollick, D.N. Homogeneous and isotropic cosmologies with nonlinear electromagnetic radiation. Phys. Rev. D 2008, 78, 063524. [CrossRef] 27. Kruglov, S.I. Universe acceleration and nonlinear electrodynamicsv. Phys. Rev. D 2015, 92, 123523. [CrossRef] 28. Kruglov, S.I. Nonlinear electromagnetic fields as a source of universe acceleration. Int. J. Mod. Phys. A 2016, 31, 1650058. [CrossRef] 29. Kruglov, S.I. Acceleration of universe by nonlinear electromagnetic fields. Int. J. Mod. Phys. D 2016, 25, 1640002. [CrossRef] 30. Ali, A.; Saifullah, K. Asymptotic magnetically charged non-singular black hole and its thermodynamics. Phys. Lett. B 2019, 792, 276. [CrossRef] 31. Mazharimousavi, S.H.; Halilsoy, M. Note on regular magnetic black hole. Phys. Lett. B 2019, 796, 123. [CrossRef] Universe 2019, 5, 225 9 of 9

32. Kruglov, S.I. Black hole as a magnetic monopole within exponential nonlinear electrodynamics. Ann. Phys. 2017, 378, 59. [CrossRef] 33. Kruglov, S.I. Nonlinear Electrodynamics and Magnetic Black Holes. Ann. Phys. (Berlin) 2017, 529, 1700073. [CrossRef] 34. Myung, Y.S.; Kim, Y.-W.; Park, Y.-J. Thermodynamics of regular black hole. Gen. Rel. Grav. 2009, 41, 1051. [CrossRef] 35. Tharanath, R.; Suresh, J.; Kuriakose, V.C. Phase transitions and geometrothermodynamics of regular black holes. Gen. Rel. Grav. 2015, 47, 46. [CrossRef] 36. Kruglov, S.I. On a model of magnetically charged black hole with nonlinear electrodynamics. Universe 2018, 4, 66. [CrossRef] 37. Myung, Y.S.; Kim, Y.-W.; Park, Y.-J. Quantum cooling evaporation process in regular black holes. Phys. Lett. B 2007, 656, 221. [CrossRef] 38. Fan, Z.Y.; Wang, X.B. Construction of Regular Black Holes in General Relativity. Phys. Rev. D 2016, 94, 24027. [CrossRef] 39. Hawking, S.W. Particle creation by black holes. Commun. Math. Phys. 1975, 43, 199. [CrossRef] 40. Novikov, I.D.; Frolov, V.P. Physics of Black Holes; Kluver Academic Publishers: Dordrecht, The Netherlands, 1989. 41. Bronnikov, K.A. Nonlinear electrodynamics, regular black holes and wormholes. Int. J. Mod. Phys. D 2018, 27, 1841005. [CrossRef]

c 2019 by the author. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).