Wormhole solutions with NUT charge in higher curvature theories

We present wormholes with a Newman-Unti-Tamburino (NUT) charge that arise in certain higher curvature theories, where a scalar field is coupled to a higher curvature invariant. For the invariants we employ i) a Gauss-Bonnet term and ii) a Chern-Simons term, which then act as source terms for the scalar field. We map out the domain of existence of wormhole solutions by varying the coupling parameter and the scalar charge for a set of fixed values of the NUT charge. The domain of existence for a given NUT charge is then delimited by the set of scalarized nutty black holes, a set of wormhole solutions with a degenerate throat and a set of singular solutions.


Introduction
Wormholes are intriguing solutions of numerous theories of gravity. In General Relativity the presence of some form of exotic matter is required in order to construct traversable Lorentzian wormholes, since the energy conditions must be violated (see, e.g., [1][2][3].) The simplest possibility here is to employ a phantom scalar field, i.e., a field whose kinetic term has the opposite sign as compared to an ordinary scalar field. Such phantom fields have been already employed decades ago by Ellis and Bronnikov, when constructing wormhole solution in General Relativity [4][5][6][7].
In alternative theories of gravity, however, classical Lorentzian wormholes may arise without the need for exotic matter. The reason is, that the gravitational interaction itself can give rise to additional terms, that may be interpreted as contributing to an effective stress energy tensor on the right hand side of the Einstein field equations, which then may lead to violation of the energy conditions. Thus it is the modified gravitational interaction itself which provides the necessary violation of the energy conditions (see e.g., [8][9][10][11][12][13][14][15][16]).
In recent years alternative gravity theories have received much attention [17][18][19][20]. To a large extent this interest has been driven by the quest to resolve various cosmological issues. However, alternative theories of gravity may also lead to interesting predictions for compact objects in the strong gravity regime, where the advent of gravitational wave multi-messenger astronomy has opened a new window into the universe with the new observational possibilities and challenges.
Among the numerous alternative gravitational theories we consider higher curvature theories very attractive, where a scalar field is coupled to a higher curvature invariant. On the one hand higher curvature terms arise under various circumstances, when quantum gravity theories are considered, and on the other hand their coupling constants may still be large, thus leading to interesting potentially observable effects in the strong gravity regime, since constraints from the solar system or binary pulsars may easily be satisfied.
Here we consider two such theories: i) Einstein-scalar-Gauss-Bonnet (EsGB) theory and ii) Einstein-scalar-Chern-Simons (EsCS) theory. In both cases we choose the same type of coupling function of the scalar field to the invariant, namely a function quadratic in the scalar field. When such a quadratic coupling is chosen, a vacuum solution of General Relativity supplemented with vanishing scalar field is also a solution of EsGB and EsCS theory. Thus the Schwarzschild black hole and the Kerr black hole remain solutions of these alternative gravity theories. However, these theories may possess additional black hole solutions with a finite scalar field, i.e., spontaneously scalarized black holes.
In the case of EsGB theory spontaneously scalarized black holes were first obtained in the static spherically symmetric case [21][22][23] and later generalized to include rotation [24][25][26][27][28]. Here the Gauss-Bonnet (GB) invariant acts as an effective mass term in the scalar field equation, that may induce a tachyonic instability of the Schwarzschild or Kerr black holes.
In the case of EsCS theory there are no spontaneously scalarized static spherically symmetric black holes. The reason is that the Chern-Simons (CS) invariant vanishes for the Schwarzschild solution. So no tachyonic instability will be induced in this case. For the CS invariant to contribute one would need the presence of a parity-odd source. Such a source would be rotation [29,30], however, the construction of rotating spontaneously scalarized EsCS black holes will be technically rather challenging.
Therefore, as a first step towards obtaining rotating spontaneously scalarized black holes in EsCS theory, a by far simpler case has been studied in [31]. Here a NUT charge has been introduced to obtain a parity odd source term, but retain a system of ordinary differential equations, since the angular dependence factorizes. The CS invariant can then act as an effective mass term in the scalar field equation and induce a tachyonic instability of the Schwarzschild black holes. A NUT charge is, in fact, quite intriguing. It gives, for instance, rise to a Misner string on the polar axis and implies the presence of closed timelike curves (see e.g. [32] and references therein).
Here we consider wormholes with NUT charge in i) EsGB theory and ii) EsCS theory, recalling and extending previous results [33,36,37]. Since one boundary of the domain of existence of wormholes in such theories is typically given by the corresponding set of black holes, we follow closely Brihaye et al. [31], who have obtained the sets of spontaneously scalararized EsGB and EsCS black holes with NUT charge. The presence of the NUT charge does not lead to asymptotically flat solutions in the usual sense. However, the asymptotic fall-off of the functions is of the usual type of an asymptotically flat spacetime.
The paper is organized as follows: We first present the action for both theories and their equations of motion. Then we discuss the conditions for the presence of a throat. We next determine the boundary conditions, the junction conditions and the energy conditions. Subsequently, we present our results for i) EsGB wormholes and ii) EsCS wormholes with NUT charge. We then end with our conclusions.

Action and equations of motion
We employ the effective action for Einstein-scalar-higher curvature invariant theories as presented in Brihaye et al. [31] (with geometrized units with G = c = 1) where R is the curvature scalar, φ denotes a massless scalar field without self-interaction, coupled with coupling function F (φ) to an invariant I(g). For the coupling function F (φ) we choose a quadratic φ-dependence with coupling constant α, This is the simplest choice that leads to spontaneous scalarization of black holes. We emphasize, that we follow the conventions of Brihaye et al. [31], since we would like to compare with their results. Obviously, any prefactor of the kinetic term in the action Eq. (1) can be absorbed by re-definitions of the scalar field φ and parameter α of the coupling function F (φ) = αφ 2 . For the invariant I(g) we make the following choices: (i) we choose the Gauss-Bonnet term and (ii) we choose the Chern-Simons term Here we have employed the Hodge dual of the Riemann-tensor * R µ ν ρσ = 1 2 η ρσκλ R µ νκλ which is defined with the 4-dimensional Levi-Civita tensor η ρσκλ = ǫ ργστ / √ −g. We note that both invariants are topological in four dimensions. However, their coupling to the scalar field φ lets them contribute to the equations of motion.
The coupled sets of field equations are then obtained by variation of the action (1) with respect to the scalar field and the metric, Here G µν is the Einstein tensor, as usual, while T (eff) µν denotes the resulting effective stress energy tensor consisting of the standard scalar field contribution and a gravitational contribution from the respective higher curvature invariant I(g). For the two above invariants we obtain (i) and (ii) To obtain wormhole solutions with a NUT charge N we assume the line element to be of the form With this ansatz the angular dependence factorizes, and all three functions, the two metric functions f 0 and f 1 and the scalar field function φ, depend only on the radial coordinate r. Insertion of the above ansatz (11) for the metric and the scalar field into the scalar field equation (5) and the Einstein equations (6) with effective stress-energy tensor (7) leads to five coupled, nonlinear ordinary differential equations (ODEs). However, these are not all independent, and one ODE can be treated as a constraint. We thus retain three coupled ODEs of second order to be solved numerically. We remark that in case (i) we can even obtain one first order and two second order ODEs.
The field equations possess an invariance under the scaling transformation with constant χ > 0.

Throat and boundary conditions
In order to obtain wormhole solutions, the presence of a throat is mandatory. We therefore introduce the circumferential (or spherical) radius of the spacetime. When the circumferential radius has a single finite extremum, a minimum, this corresponds to the single throat of the wormhole, since the throat is a surface of minimal area. While wormholes with more extrema also exist, which then feature an equator surrounded by a double throat, we here focus on single throat wormholes. We now require the presence of an extremum of the circumferential radius at some r = r 0 , in order to obtain the first set of boundary conditions. This leads to Consequently, we will refer to the two-dimensional submanifolds defined by r = r 0 and t =const. as the throat of the wormholes. The presence of the NUT charge endows these wormholes with an interesting feature: the throat metric of the wormholes with NUT charge changes its signature, when the coefficient of dϕ 2 changes sign. This happens at the critical angles θ c and π − θ c . where θ c is obtained by requiring det(g th ) ≥ 0, and thus θ c ≤ θ ≤ π − θ c with The metric has positive signature, as required for a two-dimensional Riemannian surface, only for θ c ≤ θ ≤ π−θ c . The signature change is, however, not a surprise for a spacetime with a NUT charge, since a NUT charge leads to a non-causal structure of a spacetime and allows for closed timelike curves. The second set of boundary conditions is obtained by imposing the usual fall-off as for asymptotically flat spacetimes and thus the usual boundary conditions for r → ∞ [31]. Thus the metric functions and the scalar field satisfy the asymptotic expansions where the constants M and D denote the mass and the scalar charge of the wormholes, respectively, while φ ∞ denotes the asymptotic value of the scalar field. All higher order terms in the expansion can be expressed in terms of the constants M , D and φ ∞ . A solution is therefore uniquely determined by these three constants (and, of course, the parameters of the theory). Since we follow Brihaye et al. [31] in order to compare with their spontaneously scalarized black hole solutions, we need to also employ their chosen asymptotic value φ ∞ = 0.

Junction conditions
Our aim is to obtain wormholes, which are symmetric with respect to the throat. Inspection of the functions at the throat shows, that if the solutions would be simply continued to the other side of the throat, the solutions would not be symmetric. Worse, however, is that further integration beyond the throat will lead to a curvature singularity, and asymptotic infinity on that side cannot be reached. Therefore we cut the solutions at the throat and paste a symmetric copy on the other side of the throat. In order to have continuous and differentiable wormholes, we therefore need to impose junction conditions at the throat.
To that end we introduce the radial coordinate η, where r 0 is a constant, and define the constant η 0 via η 0 = 2r 0 . We then reexpress the metric in terms of the new radial coordinate η and the new metric function F 1 , Thus the throat is located at η = 0, and the solution for η ≤ 0 is obtained from the solution for η ≥ 0 by imposing the conditions Since these conditions generically introduce jumps in the derivatives of the functions f 0 and φ at the throat, we impose the presence of a thin shell of matter at the throat. The proper procedure here is the make use of an appropriate set of junction conditions [38,39], namely we evaluate the jumps in the Einstein and scalar field equations that arise for η → −η, Here the stress energy tensor of the matter at the throat has been denoted by s µ ν , while the source term for the scalar field has been named s scal , and dF (φ)/dφ =Ḟ denotes the derivative of the coupling function. Since the thin shells should be composed of ordinary (non-exotic) matter, we assume a perfect fluid at the throat with pressure p and energy density ǫ c . Moreover, we assume a scalar charge density ρ scal together with a gravitational source employed before for EsGB wormholes without NUT charge [13,14,33]. Here λ 1 and λ 0 are constants,h ab is the induced metric at the throat, andR the corresponding Ricci scalar. The junction conditions are obtained by substituting the metric into the set of equations (22). Now we evaluate the junction conditions for both invariants successively. Note, that all functions and derivatives are evaluated at the throat. For the Gauss-Bonnet invariant we obtain where the prime denotes the derivative with respect to η. These conditions follow from the ( t t ), t ϕ , and we obtain For all of our solutions f ′ 0 > 0, therefore the energy density ǫ c is always positive for this choice of the constants λ 0 and λ 1 .
For the Chern-Simons invariant we obtain again from the ( t t ), t ϕ , ϕ ϕ components of the Einstein equations and the scalar field equation, respectively. Considering again pressureless matter, p = 0, with we now obtain

Energy conditions
In order to obtain wormhole solutions the null energy condition (NEC) T µν n µ n ν ≥ 0 (36) must be violated. Here n µ may be any null vector (n µ n µ = 0). Consequently, all one has to show is that null vectors exist, such that T µν n µ n ν < 0 somewhere in the spacetime. Choosing the null vector n µ n µ = 1, −g tt /g ηη , 0, 0 , and thus n µ = g tt , √ −g tt g ηη , 0, 0 , the NEC leads to This shows that the NEC is violated when Alternatively, choosing the null vector n µ = 1, 0, −g tt /g θθ , 0 , leads to NEC violation when These conditions have been investigated before for various wormhole solutions (see e.g., [1,13,14,33]).

Results
We solve the coupled Einstein and scalar field equations numerically. To this end we introduce the inverse radial coordinate x = 1/r. Consequently, the asymptotic region r → ∞ corresponds to the region x → 0. Here the expansion of the metric functions and the scalar field then reads (Eqs. (17)-(19)) To solve the system of ODEs, we treat it as an initial value problem, and employ the fourth order Runge Kutta method. The above expansion then yields the initial values, where the prime denotes the derivative with respect to x. We start the computation at spatial infinity, x = 0, and end at the throat at some finite value x = x 0 , where we reach the condition (14). By following the outlined numerical procedure we have constructed numerous sets of wormhole solutions with NUT charge for both invariants I(g). In the following we demonstrate our results for both invariants successively. The metric functions e f0 and e F1 , and the scalar field function φ are shown in Fig. 1(a) versus the scaled radial coordinate η/M . Also shown is the scaled circumferential radius R c /η 0 versus η/M , where η 0 /M = 0.517. The minimum of R c at the coordinate origin η = 0 corresponds to the location of the throat.
As required the wormhole solutions violate the energy conditions. We demonstrate the violation of the NEC conditions −T t t + T η η ≥ 0 and −T t t + T θ θ ≥ 0 for the same solution in Fig. 1(b). Clearly, both expressions are negative in a region close to the throat. Note, that we have omitted in the notation of the energy conditions, that we refer to the effective stress energy tensor. We emphasize again that the violation is caused by the gravitational part of the effective stress energy tensor.
We next turn to the domain of existence of these EsGB wormholes with NUT charge. The domain of existence is shown in terms of the scaled coupling constant 8α/M 2 and the scaled scalar charge D/M in Fig. 1(c) for several values of the scaled NUT charge, n = N/M = 1, 2, and 3, and compared with the case of vanishing NUT charge, n = 0 [33]. Note, that we only show the right hand side D ≥ 0 of the domain, since the domain is symmetric with respect to D → −D, which is a consequence of the quadratic coupling function and the chosen boundary condition φ ∞ = 0.
For any given value of the NUT charge the domain of existence is delimited by three sets of solutions. The first set corresponds to the scalarized black holes with NUT charge of Brihaye et at. [31], which we fully reproduce. These black holes are shown by a solid green curve. The second set is shown by a dotted red curve and corresponds to a set of singular solutions. Here a cusp singularity is encountered at some value η ⋆ of the radial coordinate. Cusp singularities are a recurring feature of scalarized wormholes (see [33][34][35][36]). They form, when the determinant of the coefficients of the second order terms of the ODEs vanishes.
The third set of boundary solutions is shown by a dashed blue curve and referred to as degenerate set of solutions. To clarify the physical meaning of this set, we recall that we identified the throat by looking for a minimum of the circumferential radius. However, it may happen that this extremum becomes a degenerate extremum. Here, in principle, a new set of wormholes emerges, which possess an equator and a double throat, i.e., the degenerate extremum will split into a maximum (equator) and two minima (throats). We remark, that for D = 0 the scalar field vanishes, thus vacuum solutions of General Relativity are recovered. Overall we note, that with increasing NUT charge the domain of existence of wormhole solutions increases rapidly.
We now turn to the causal structure of the EsGB wormhole spacetimes. The relevant quantity to consider here is the polar angle, where the g ϕϕ component changes sign, and thus closed timelike curves will be allowed. In particular, we are interested in the critical polar angle θ c , where this change of sign occurs at the wormhole throat. We show the boundary of the domain with changed causal structure at the throat in Fig. 1(d), where we exhibit the critical polar angle θ c versus the scaled coupling constant 8α/M 2 for several values of the scaled NUT charge, n = N/M = 1, 2, and 3. For a fixed value of the NUT charge the dotted blue curve shows the singular boundary set, while the dashed green curve shows the degenerate boundary set. Only for the limiting case of black holes the critical polar angle θ c goes to zero. Thus all the wormholes with NUT charge possess a region with closed timelike curves that extends to their throat.

Case ii): EsCS wormholes
We now discuss the properties of EsCS wormholes with NUT charge. We start again with a typical wormhole solution, that has been obtained for the parameter choice α/M 2 = 4, D/M = 1 and n = N/M = 3. We show its metric functions e f0 and e F1 , and its scalar field function φ versus the scaled radial coordinate η/M in Fig. 2(a), along with its scaled circumferential radius R c /η 0 , where η 0 /M = 0.495. We note, that the functions are rather similar to those of the EsGB wormholes shown in Fig. 1(a).
We demonstrate the violation of the NEC conditions −T t t + T η η ≥ 0 and −T t t + T θ θ ≥ 0 in Fig. 2(b) for this solution. Again we find a region close to the throat, where both conditions are violated.
The domain of existence of the EsCS wormholes is shown in terms of the scaled coupling constant 8α/M 2 and the scaled scalar charge D/M in Fig. 2(c) for several values of the scaled NUT charge, n = N/M = 0.5, 1, 2, and 3. The boundaries of the domain of existence are completely analogous to the previously discussed case of EsGB wormholes with NUT charge.
The first set of boundary solutions corresponds to the scalarized black holes with NUT charge of Brihaye et at. [31] and is shown by a solid green curve. The second set of boundary solutions corresponds to the set of solutions with a cusp singularity and is shown by a dotted red curve. Finally the third set of boundary solutions corresponds to wormhole solutions with a degenerate throat, beyond which wormholes with an equator and a double throat arise. They are shown by a dashed blue curve. The domain of existence is again symmetric with respect to D → −D, and for D = 0 solutions of General Relativity are recovered.
However, we also notice a new feature of the EsCS wormholes as compared to the EsGB wormholes. Whereas for EsGB wormholes the domain of existence moves continuously to smaller values of the coupling constant α, as the NUT charge is decreased towards zero, this is not the case for the EsCS wormholes. Here the domain of existence exhibits the same dependence as the scalarized EsCS black holes, namely after reaching a minimum of α with decreasing NUT charge, a further decease of the NUT charge leads to increasing values of α. This is seen in Fig. 2(c), where the domain of existence for n = 0.5 has moved to larger values of α as compared to the domain for n = 1. As the NUT charge tends towards zero, the coupling constant α tends towards infinity. However, in the limit a finite source term for the scalar field results, since the NUT charge and the coupling constant conspire appropriately.
We finally turn to the causal structure of the EsCS wormhole spacetimes. We show the critical polar angle θ c , where the metric component g ϕϕ changes sign at the wormhole throat, versus the scaled coupling constant 8α/M 2 in Fig. 2(d) for the same set of values of the scaled NUT charge, n = N/M = 0.5, 1, 2, and 3. Again, for a fixed NUT charge, the boundaries are given by the dotted blue curve showing the singular boundary set, and the dashed green curve showing the degenerate boundary set. Only for the limiting case of black holes the critical polar angle θ c goes to zero.

Conclusions
We have considered scalarized wormholes with NUT charge in two alternative gravity theories, EsGB theory and EsCS theory. In both cases a scalar field has been coupled with a quadratic coupling function to the respective invariant, since this choice of coupling function allows for spontaneously scalarized black holes with NUT charge, i.e., solutions obtained before by Brihaye et al. [31].
In these alternative gravity theories wormholes arise without the need of introducing exotic matter. The reason is the presence of gravitational terms in the stress energy tensor that result, respectively, from the higher curvature GB term or CS term. These terms allow for violations of the energy conditions and therefore the presence of wormholes.
We have mapped out the domain of existence of these wormholes for various values of the NUT charge. In all cases, for a given NUT charge, the domain has one boundary consisting of scalarized black holes with NUT charge, one boundary formed by solutions with a cusp singularity, and one boundary composed of wormholes with a degenerate throat, where a new type of wormhole arises, that features an equator and a double throat, i.e., wormholes to be addressed in the future.
In order to have symmetric wormholes without singularities (apart from the Misner string caused by the NUT charge), we have cut the wormhole spacetimes at the throat, retaining only the part extending to infinity, and pasted the symmetrically reflected copy at the throat. We have shown that the resulting junction conditions at the throat can be fulfilled by a thin shell of ordinary matter such as dust, for example.
Since these wormholes carry NUT charge, they possess closed timelike curves. Therefore we have investigated, in particular, the causal structure at the throat. Interestingly, all these wormholes with NUT charge possess closed timelike curves that extend all the way to the throat. However, for the limiting set of scalarized black holes with NUT charge the closed timelike curves do not extend to the horizon, but can get arbitrarily close to it.
We have considered wormholes with a NUT charge as a toy model for rotating wormholes in EsGB theory and EsCS theory. In this toy model the angular dependence factorizes, leaving a set of nonlinear coupled ODEs to be solved numerically. In the case of rotating wormholes, we will have to solve sets of nonlinear coupled partial differential equations, which will be rather challenging for these alternative gravity theories. Moreover, the proper conditions for the throat and the proper set of junction conditions will have to be formulated and solved in this case, as well.