On a New Approach for Constructing Wormholes in Einstein-Born-Infeld Gravity

We study a new approach for the wormhole construction in Einstein-Born-Infeld gravity, which does not require exotic matters in the Einstein equation. The Born-Infeld field equation is not modified from"coordinate independent"conditions of continuous metric tensor and its derivatives, even though the Born-Infeld fields have discontinuities in their derivatives at the throat in general. We study the relation of the newly introduced conditions with the usual continuity equation for the energy-momentum tensor and the gravitational Bianchi identity. We find that there is no violation of energy conditions for the Born-Infeld fields contrary to the usual approaches. The exoticity of energy-momentum tensor is not essential for sustaining wormholes. Some open problems are discussed.


I. INTRODUCTION
Wormholes are non-singular space-time structures connecting two (or more) different universes or parts of the same universe [1]. However, in the conventional approaches for (traversable) wormholes, there are problems in the "naturalness". First, the hypothetical exotic matters, which violate the energy conditions in the standard general relativity (GR) but are essential for sustaining the throats of wormholes, have never been discovered [2]. Second, the usual cuts and pastes of the throats and subsequent putting of the hypothetical matters to the throats by hand to satisfy the Einstein equation [3] are too artificial to be considered as a natural process. Moreover, we don't know much about the formation mechanism of wormholes due to the gravitationally repulsive nature of the exotic matters.
To circumvent the obstacle by exotic matters, one may consider other extended theories of gravity so that the exotic matters are naturally generated. For example, one may consider non-minimally coupled scalar-tensor gravity theories [4], higher curvature gravities [5], higher dimensional gravities [6], and brane-worlds theories [7]. (For a more complete and modern review, see Ref. [8] and references therein.) On the other hand, to cure the problem of artificiality of the conventional cuts and pastes construction of throats of wormholes, a new approach for the wormhole construction has been proposed recently [9,10]. In the new approach, the throat is defined as the place where the solutions are smoothly joined. There, the metric and its derivatives are continuous so that the exotic matters are not introduced at the throats.
From the new definition, throats can not be constructed arbitrarily contrary to the conventional cuts and pastes approach. For example, we consider a four dimensional spherically symmetric wormhole connecting two remote parts of the same universe (or mathematically equivalently, the reflection symmetric two universes). The metric is described by ds 2 = −N ± (r) 2 c 2 dt 2 + dr 2 f ± (r) + r 2 dθ 2 + sin 2 θdφ 2 (1) with two coordinate patches, each one covering the range [r 0 , +∞). Then, the radius of throat r 0 is defined as with the usual matching condition by demanding that the metric and the Schwarzschild coordinates be continuous, though not generally smooth, across the throat. If there exists a coordinate patch M + in which the singularities are absent for all values of r ≥ r 0 , one can construct a smooth regular wormhole-like geometry, which may or may not be traversable depending on the location of r 0 , by joining M + and its mirror patch M − at the throat r 0 . Here, it is important to note that, in the new approach, f ± (r 0 ) needs not to be vanished in contrast to the Morris-Thorne approach [2], while the quantities dN ± (r 0 )/dr, df ± (r 0 )/dr in (2) need not to be vanished in both the Morris-Thorne approach [2] and the cuts and pastes approach of Visser [3]. So, in the new approach, we have two options for exact solutions of the same Einstein equation when the throat point r 0 exists. We may consider either the usual solutions for compact objects, like the black holes which may have singularities (but usually shielded by the horizons) generally, or wormhole-like objects which do not have singularities in the whole space-time domain. We do not have a priori reason to choose only one of the two options so that we can consider both options equally. This is in contrast to the Morris-Thorne-type traversable wormholes, which can not co-exist with the solutions of black holes with the same conserved charges. In this context, we call the new type of wormholes as "natural wormholes" and the their throats as the natural throats to distinguish these from the previously known wormholes.
Originally the new wormhole solutions [9] have been studied in Hořava (or Hořava-Lifshitz) gravity, which has been proposed as a renormalizable quantum gravity based on different scaling dimensions between the space and time coordinates [11]. But the new approach can be applied to other gravity theories as well if the natural throat exists. However, the existence of the throat for black holes with rotation would be very difficult unless some coincidences occur. For the Kerr black hole in GR, there are more metric functions and there is no solution for the throat where all the metric functions join smoothly simultaneously. This implies that the wormhole throat in the spherically symmetric configurations could be easily destroyed by adding other conserved charges. Conversely, the wormhole throat could be formed only after losing all charges except the mass [10]. It has been also argued that the situation with electromagnetic charges would be similar because there are additional gauge fields to be joined smoothly at the same throat. One can easily show that there is no solution for the throat where the electromagnetic fields F ± µν join smoothly, dF ± µν (r 0 )/dr = 0, F + µν (r 0 ) = F − µν (r 0 ), in addition to the conditions for the metric (3). The purpose of this paper is to see whether the concept of natural wormholes could be also generalized to the charged cases or not. In particular, in order to see the generic role of charges for non-linear electromagnetic fields at short distances, we consider the Born-Infeld (BI) action, instead of the usual Maxwell's, without modifying the Einstein-Hilbert action at short distances.
The organization of this paper is as follow. In Sec. II, we set up the basic equations of Einstein-Born-Infeld gravity with a cosmological constant and study the exact solutions for spherically symmetric black holes and their physical properties, including the black hole thermodynamics and phase structures. In Sec. III, we study the natural wormhole geometry as the solution of the Einstein equation without introducing additional matters at the throat based on the construction sketched in Sec. I. We show that the BI field equation without additional terms at the throat is still valid even with the throat from some reasonable conditions of "coordinate independent" continuity of metric tensors, the preservation of the usual gravitational Bianchi identity, and the continuity equation of energy-momentum tensors for the BI fields. In Sec. IV, we conclude with several discussions. Throughout this paper, we use the conventional units for the speed of light c and the Boltzman's constant k B , c = k B = 1, but keep the Newton's constant G and the Planck's constanth unless stated otherwise.

II. BLACK HOLE SOLUTIONS IN EINSTEIN-BORN-INFELD GRAVITY
Born-Infeld action of non-linear electrodynamics was originally introduced to solve the infinite self-energy of a point charge in Maxwell's linear electrodynamics [12]. When this action is combined with Einstein action to study the gravity of charged objects, the shortdistance behavior of the metric is modified [13]. With the recent development of D-branes, the Born-Infeld-type action has attracted renewed interests as an effective action for lowenergy superstring theory [14]. There have been extensive studies about black hole solutions in the Einstein gravity coupled to Born-Infeld electrodynamics (Einstein-Born-Infeld (EBI) gravity) [15][16][17][18][19][20][21][22][23][24][25]. In this section, we briefly describe EBI gravity and its black hole solutions.
Since our results do not depend much on the space-time dimensions D ≥ 4, we consider the EBI gravity action with a cosmological constant Λ in D = 3+1 dimensions for simplicity, where L(F ) is the BI Lagrangian density given by Here, the parameter β is a coupling constant with dimensions [length] −2 which needs to flow to infinity to recover the usual Maxwell electrodynamics with L(F ) = −F µν F µν ≡ −F 2 at low energies. The correction terms from a finite β represent the effect of the nonlinear BI fields. This situation is similar to the Lagrangian for a relativistic free particle mc 2 (1 − 1 − v 2 /c 2 ) which reduces to the Newtonian Lagrangian 1 2 mv 2 for c → ∞ limit. As the speed of light c gives the upper limit for the speed v of a particle, the parameter β gives the upper limit for the strength of fields − 1 2 F 2 if F 2 < 0. Taking 16πG = 1 for simplicity, the equations of motion are obtained as where the energy momentum tensor for BI fields is given by Let us now consider a static and spherically symmetric solution with the metric ansatz For the static 1 electrically charged case where the only non-vanishing component of the field strength tensor is F rt ≡ E, one can find Then, the general solution for the BI electric field is obtained from the equation (6) as where Q is the integration constant that represents the electric charge localized at the origin. Note that the electric field is finite in the limit r → 0, E(r) ≈ ±|β| for |β| < ∞ due to the non-linear effect of BI fields while it reduces to the usual Coulomb field E ≈ Q/r 2 in the limit |β| → ∞ (Fig. 1). Now, from (10) and (11), the Einstein equation (7) reduces to and a simple integration over r gives where C is an integration constant. With a straightforward integration, one can express the solution in a compact form in terms of the incomplete elliptic integral of the first kind, For small r, this may be expanded as which shows a mass-like term −2C/r. The parameter C is called an intrinsic mass [16] and represents the mass-like behavior near the origin. However, it is important to note that this is not the usual ADM mass defined at large r. In order to see this, it is convenient to reorganize the integral in (13) as r 0 = ∞ 0 + r ∞ so that the integral r ∞ is convergent in the limit r → ∞. Then, we have another expression for the solution in terms of the hypergeometric function [17][18][19][20]23]. Here, M is another mass parameter, defined by where M 0 comes from the ∞ 0 part of the integral. One can find that the total mass M, which is the sum of the intrinsic mass C and the (finite) self-energy of a point charge M 0 , is the usual ADM mass defined by the asymptotic behavior of the metric at large r This shows that the physical parameters, like the mass and the cosmological constant, are shifted at the short distances due to non-linear effects of BI fields: The fourth term on the right hand side of (15) represents the space with an effective cosmological constant Λ eff = Λ − 2β 2 , which can behave like an anti de-Sitter space (Λ eff < 0) near the origin even though it behaves like a flat (Λ = 0) or a de-Sitter (Λ > 0) space at the asymptotic infinity. The solution (14) or (16) has a curvature singularity at the origin, r = 0, Note that the leading singularity near r = 0 is of the order O(r −6 ) in the Riemann tensor square, which is the same as that of Schwarzschild (Sch) black hole. But this singularity is absent for the marginal case of C = 0 (M = M 0 ) due to the regular behavior of the solution near the origin (15), which is the same as that of Reissner-Nodstrom (RN) black hole 2 . On the other hand, it can be shown that there is no curvature singularity at the horizon, which is defined as the solution for f (r) = 0 in our case, as expected [17]. The solution can have two horizons generally and the Hawking temperature for the outer horizon r + is given by One can not directly obtain the singular behavior of RN black hole near r = 0 from (20) as one can not obtain RN black hole solution from (15) by taking |β| → ∞ limit because |β| < ∞ has been implicitly assumed in the expansions of (15) and (20). One can obtain RN results only from (19) which is valid There exists an extremal black hole limit of vanishing temperature, where the outer horizon r + meets the inner horizon r − at for a non-vanishing cosmological constant Λ and for vanishing Λ. At the extremal point, the ADM mass M gets the minimum The first law of black hole thermodynamics is found as with the usual Bekenstein-Hawking entropy formula by and the scalar potential [18] A 0 (r) = 1 Now, from the result (22), one can classify the black holes in terms of the values of βQ and the cosmological constant Λ .
(i) βQ > 1/2: In this case, the black hole may have one non-degenerate horizon (type I), or two non-degenerate (non-extremal) horizons or one degenerate (extremal) horizons (type II), depending on the mass for the asymptotically AdS (Λ < 0) or flat (Λ = 0), i.e., the type I for M ≥ M 0 and the type II for M < M 0 (Fig. 2a, 2b). There can be a phase transition from Schwarzschild (Sch)-like type I black hole to Reissner-Nodstrom(RN)-like type II black hole when the mass becomes smaller than the marginal mass M 0 , which is the mass value at r + = 0 in Fig. 3. When the mass is smaller than the extremal mass M * , which is the mass value at the extremal point in Fig. 3 (top curve), the singularity at The top two curves represent M 0 > M * , M 0 = M * for βQ > 1/2, βQ = 1/2, respectively, with the extremal mass M * , whereas the bottom two curves represent the cases where M * is absent for βQ < 1/2. We consider βQ = 2/3, 1/2, 2/4.5, 1/3 (top to bottom) with β = 2. The effect of cosmological constant is important for large r + (right): (Λ < 0 (thin curve), Λ = 0 (medium curve), Λ > 0 (thick curve)). We consider Λ = ±1/5 for the (A)dS cases. r = 0 becomes naked as usual. In the limit β → ∞, only the RN-like type II black holes are possible as in GR. On the other hand, for asymptotically dS (Λ > 0) case, the situation is more complicated due to the cosmological horizon. In this case, there can exist maximally three horizons (two smaller ones for black holes and the largest one for the cosmological horizon), depending on the mass and the cosmological constant (Fig. 2c).
(ii) βQ = 1/2: In this case, only the Sch-like type I black holes are possible (Fig. 4). It is peculiar that the horizon shrinks to zero size for the marginal case M = M 0 , even though its mass has still a non-vanishing value (Fig. 3). This does not mean that the metric is a flat Minkowskian. This rather means that all the gravitational energy is stored in the BI electric field as a self-energy. When the mass is smaller than the marginal mass M 0 , the singularity at r = 0 becomes naked. Note that this is the case where the GR limit β → ∞ does not exist.  (iii) βQ < 1/2: This case is similar to the case (ii), except that the marginal case has no (even a point) horizon so that its singularity is naked always (Fig. 5).

III. CONSTRUCTION OF NATURAL WORMHOLE SOLUTIONS
The strategy to construct a natural wormhole is to find the throat radius r 0 defined by (2) with the matching condition (3) [9,10]. Based on this definition, cases (i) and (ii) in Fig.  2 and 4, respectively, and the case (iii) in Fig. 5a (AdS case) show the throats depending on the mass M for given values of βQ and Λ. More details are as follows.
For the case (i) of βQ > 1/2, there is no throat for the Sch-like type I black holes with M > M 0 down to the marginal case of M = M 0 , where a point-like throat is generated at the origin. The zero-size (r 0 = 0) throat can grow as the black hole loses its mass. But, down to the extremal case of M = M * (< M 0 ), where the inner and outer horizons meet at r * + and Hawking temperature vanishes (Fig. 6), the throat is inside the outer horizon and non-traversable. 3 As the black hole loses its mass further, there emerges and grows a 3 In the dS case (Fig. 2c), regardless of the existence of the inside throat, there is another outside radius which satisfies the throat condition (2). But this can not be considered as the wormhole throat because it does not meet the requirement that the throat occupies the minimum radial coordinate, which has been assumed implicitly in the definition of Sec. I. In principle, it could also be possible to have a throat at the maximum radius so that one universe is smoothly connected to another by just traveling toward the traversable throat from the extremal black hole horizon. 4 For the AdS case (Fig. 2a), the throat grows and persists until its mass M vanishes. However, for the cases of flat (Fig. 2b) and dS (Fig. 2c), the throat does not grow to an indefinitely large size. There is a maximum size at a certain mass M before M vanishes, i.e., M * > M > 0. If we further reduce the mass below M , only the space-time with a naked singularity is a possible solution. So, the usual naked singularity solution for M < M * may be replaced by the non-singular wormhole geometry until a certain (non-negative) mass M is reached in flat and dS cases. This may suggest the generalized cosmic censorship [9,10] by the existence of a wormhole-like structure around the origin, which has been thought to be a singular point of black hole solutions where all the usual physical laws could break down and might be regarded as an indication of the incompleteness of general relativity.
For the case (ii) of βQ = 1/2, a point-like throat can be generated at the origin contrary to the generation of a finite-size throat for the case (i). But its growing process is the same as the case (i), showing the existence of a maximum size of the throat with the mass M > 0 for flat and dS cases, and M = 0 for AdS case (Fig. 4).
For the case (iii) of βQ < 1/2, the natural wormhole geometry exists only for AdS case (Fig. 5a) and its growing behavior is similar to the case (ii). On the other hand, for flat and dS cases (Fig. 5b, 5c), only the Sch-like type I black hole solution for M > M 0 or the naked singularity solution for 0 < M < M 0 can exist.
So far, we have considered geometric aspects in constructing the natural wormhole. Here, it is important to note that the right-hand side of the Einstein equation depends only on the values of the BI field strength themselves, not on their derivatives. This means that Eq. (7) is still valid even in the presence of the throat r 0 , where the BI fields do not join smoothly, i.e., their (spatial) derivatives are discontinuous (Fig. 1). This would be valid for any kind of matter field including gauge field if the coupling with the matter and gravity does not contain derivatives. In other words, the same solution of the Einstein equation (7) without the throat can also be the exact solution for each patch of the wormhole geometry cosmological horizon. But, though interesting, this is not our main concern and we will not consider this possibility in this paper. 4 For the discussion about a dynamical generation of this unusual process, which is beyond the usual linear perturbation analysis, see Ref. [10] and references therein. And this process is different from the dynamical wormhole picture suggested by Hayward where wormholes emerge from bifurcating, i.e. non-extremal black hole horizons [10,26]. even when considering matter coupled system. However, the above argument would not hold for the BI field equation (6), which depends on the derivatives of BI fields. Equation (6) for the wormhole geometry with a throat can be modified due to the discontinuities of derivatives of BI fields. To find the possible modification term, let us introduce another coordinate η which spans the whole wormhole space-times with a geodesically complete ranges (−∞, +∞) by joining two coordinate patches, each one covering [r 0 , +∞) with the throat at η = 0 for r = r 0 . For example, we may consider the Ellis-type metric ansatz for a traversable wormhole, i.e., a throat-type geometry that physical objects can propagate through, where l is the proper length from the throat given by Then, in the region close to the throat η = 0, one may expand the original coordinate r ± for each patch as Here, we note that (dr/dη)| + r 0 ≥ 0, (dr/dη)| − r 0 ≤ 0, and (dr 2 /dη 2 )| ± r 0 > 0 in order that the minimum radius r 0 exists.
In general, there are two options for the choice of η coordinate.
(b) Second, if we consider the case of (dr/dη)| ± r 0 = 0 (bottom curve in Fig. 7), we have a singularity -delta function singularity due to discontinuity of (dr/dη) ± 0 -in the third term of (31) so that the coordinate transformation r ± (η) becomes singular at the throat.
The second option was not considered before [2,3] but it seems that both options are possible on the general ground. For example, if we consider the derivative of the metric tensor g ± µν (r) with respect to the proper length l in the Schwarzschild coordinate (1) one can join the metric smoothly, dg + µν (r)/dl| r 0 = dg − µν (r)/dl| r 0 = 0, by considering either f ± (r 0 ) = 0 (Morris-Thorne approach [2]) or dg + µν /dr| r 0 = dg − µν /dr| r 0 = 0 (new approach in this paper [9,10]).
But, independently of the choice of the two options, we may generally obtain the radial component of the BI electric field and its derivative, where we have used the continuity of the field E + (r 0 ) = −E − (r 0 ) and dǫ(±η)/dη = ±δ(η) with the step function ǫ(η) (Fig. 8).
In general, singularities may enter into the BI field equation when it involves the radial derivative d dr for the wormhole geometry with a throat, regardless of the choice of η coordinate unless we fine-tune the coordinate η [2,3] so that dη dr Note that the condition (34) has been implicitly assumed in the Morris-Thorne construction for the first option, though (dη/dr) ± r 0 becomes infinity [2]. Now, in order to understand the condition (34) in our construction for the second option, let us consider the second derivatives of the metric tensor g ± µν (r) with respect to η d 2 g ± µν (r) dη 2 = dr ± dη of the throat. However the right hand side of (40) depends on dE(r)/dr, as well as E(r), which introduces the additional delta-function term at the throat. This means that, in order to have a "consistent" equation (40), we need a counter term on the right hand side to cancel the singular term at the throat implying the violation of the continuity equation Because the Einstein equation (7) is not affected by the existence of the natural wormhole, this violation of the continuity equation means that the Bianchi identity for the curvature tensor fails also at the throat, with the Einstein tensor G µν = R µν − 1 2 Rg µν . In electrodynamics, the Bianchi identity ∂ µF µν = 0, with the dual field strengthF µν = 1 2 ǫ µνσρ F σρ , is violated at the location of the magnetic monopole. In our case, this electromagnetic Bianchi identity works as one can easily check from the solution (11) sinceF µν = 0. But we may have the violation of the gravitational Bianchi identity at the throat if we do not require the condition (34). 5 This can be considered as another support of the condition (34) in our new approach.

IV. ENERGY CONDITIONS
In the usual set-up of wormholes, exotic matters violating the energy conditions are essential to sustain the throats. 6 Even when we look for wormhole solutions without the explicit introduction of exotic matters, we need some effective matter terms in the Einstein equation which play the role of exotic matters. Wormhole solutions in higher-derivative gravities is one example. In that case, higher-derivative terms act as the effective energymomentum tensor which can violate the energy conditions generally. However, this is not the case in our construction of natural wormholes and all the energy conditions can be normal.
To see this, let us consider where we have used ρ + p r = 0 from (39). Then, it is easy to see that the quantities in (43)∼(45), as well as ρ in (39), are all non-negative for the whole range of E(r) ≤ β. In 9, we have plotted the case for two horizons in the RN-like type II black hole but the result is the same as the Sch-like type I case. This shows that the strong energy condition (SEC: ρ + p i ≥ 0, ρ + i p i ≥ 0), which includes the null energy condition (NEC: ρ + p i ≥ 0), is satisfied as well as the weak energy condition (WEC: ρ ≥ 0, ρ + p i ≥ 0) and the dominant energy condition (DEC: ρ ≥ |p i | or equivalently ρ ± p i ≥ 0).
In order to understand the physical origin of this result, which is unusual in the conventional wormhole physics, let us consider the usual Morris-Thorne ansatz instead of (1) [2,3], In the usual approach of the first option, the quantity in the bracket vanishes at the throat r = r 0 since, when we consider the proper length l, for example, (dr/dl)| r 0 = ± 1 − b(r 0 )/r 0 = 0, i.e., b(r 0 ) = r 0 . But, in order that b(r) < r for r > r 0 so that r is the space-like coordinate for an observer outside the throat, the quantity in the bracket needs to be positive, i.e., it has "positive" derivative 7 d dr which reduces to so that all the energy conditions are violated. However, this is not the case in our approach of the second option, where the quantity in the bracket does not vanish at the throat for the traversable wormhole. The radial derivative of the quantity in the bracket, which is a positive quantity for r > r 0 , can have any values. Actually, in our example, we have ρ + p r = 0. One can also obtain similar results for other combinations of the energy and pressures as in Fig. 9. This shows that the exoticity of the energy-momentum tensor for matters or effective matters from the modified gravities is not essential in the new approach.

V. DISCUSSION
We have studied a new approach to construct a natural wormhole geometry in EBI gravity without introducing additional exotic matters at the throat. BI field equations are not modified at the throat for coordinate independent conditions of continuous metric tensor and its derivatives, even though BI fields have discontinuities in their derivatives generally. If we do not require the newly introduced conditions, the modification term exists and it produces the violation of the gravitational Bianchi identity. This supports the necessity of the new conditions in our construction. Remarkably, there is no violation of energy conditions and the energy-momentum tensors are normal. It is shown that this is the crucial effect of the new approach and the exoticity of the energy-momentum tensors for the matters is not essential for natural wormholes. It would be a challenging problem to see whether the non-exoticity in the energy-momentum tensors is a sign for the stability of wormholes or not. Several further remarks are in order.
First, the new approach for wormholes has been first studied in Hořava gravity [9], which has been proposed as a renormalizable quantum gravity model, where the effective energymomentum tensors from the higher spatial derivative contributions violate the usual energy conditions in general. There, it has been argued that the Hořava gravity wormhole is the result of microscopic wormholes created by negative energy quanta which have entered the black hole horizon in Hawking radiation process 8 . In other words, the existence of wormholes may reflect the quantum gravity effect of black holes. However, we have found that the 7 At the throat r = r 0 , this is equivalent to (d 2 r/dl 2 )| r0 = 1 2 d dr (1 − b/r)| r0 > 0, which is known as the "flaring-out" condition [2]. 8 Hawking radiation process involves virtual pairs of particles near the event horizon, one of the pair enters into the black hole while the other escapes. The outgoing particle can be observed as a real particle with natural wormhole construction is still valid in the classically charged black holes in GR and the usual energy conditions can be satisfied too. This indicates that the formation mechanism for the natural charged wormholes has the classical origin in contrast to vacuum wormholes in Hořava gravity and other higher curvature gravities. So, there are two possible origins for wormholes, one for the classical and the other for the quantum mechanical ones, and both effects need to be considered generally. Actually, in our wormholes for EBI gravity, there is the upper bound for the electric field which may reflect the quantum mechanical pair creation of photons 9 at short distances so that both the classical and quantum effects are involved in the natural wormholes.
Second, we note that the new approach may have some difficulties when extending to higher derivative gravity theories due to possible singularities at the throat. For example, when third-order spatial derivatives appear in the Einstein equation, we need to consider In the usual Morris-Thorne approach, only the third term survives at the horizon due to (dr ± /dl)| r 0 = 0 and it remains finite from the continuity condition (d 2 r + /dl 2 )| r 0 = (d 2 r − /dl 2 )| r 0 . Similarly, one may generalize this argument to arbitrarily higher-order derivatives in the Einstein equation. On the other hand, in the new approach the third term vanishes but the first two terms survive at the throat due to (dg ± µν /dr ± )| r 0 = 0. However, the second term can cause a singularity due to the discontinuity of (dr ± /dl)| r 0 unless we consider the case of (d 2 g ± µν /dr ± 2 )| r 0 = 0, which does not seem to be quite generic. This might imply that the new approach is not complete to describe all the generic higher-derivative gravity theories. Or this might imply that the natural wormholes are unstable for the case where the higher derivatives are involved like rotating black holes in Hořava gravity or other higher curvature gravities, as has been pointed out by one of us recently [10].
Note added: After finishing this paper, a related paper [31] appeared whose analysis of black hole in four dimensions is partly overlapping with ours.