Construction of a traversable wormhole from a suitable embedding function

In this work, we construct a traversable wormhole by providing a suitable embedding function ensuring the fulfilling of the flaring--out condition. The solution contains free parameters that are reduced through the study of the acceptable conditions of a traversable wormhole. We compute both the quantifier of exotic matter and the quasi--normal modes through the $13^{th}$ order WKB as a function of the remaining free parameters. We obtain that the wormhole geometry can be sustained by a finite amount of exotic matter and seems to be stable under scalar perturbations.

From a technical point of view, any traversable wormhole could be constructed by providing its geometry in terms of free parameters which should be constrained based on the acceptability conditions it must satisfy: i) the existence of a throat connecting two asymptotically flat regions, ii) small tidal forces bearable by a human being iii) finite proper time to traverse the throat, among others. However, the fulfilling of all the requirements is not always possible. For example, the solution should not be asymptotically flat or should require an infinite amount of exotic matter supporting it. It should be emphasized that the above mentioned requirements are not universal but sufficient for a wormhole to be a physically viable and suitable for interstellar travel by human beings. For example, we can construct asymptotically AdS wormholes which are suitable for interstellar travels. [25]. Nevertheless, in this work, it is our main goal to construct an asymptotically flat traversable wormhole supported by a finite amount of exotic matter by assuming a general embedding function. Besides, we explore its stability thorough its response to scalar perturbation.
The response of a wormhole to perturbations is dominated by damped oscillations called quasi-normal modes. The computation of the QNM modes can be performed through a variety of methods (for an incomplete list see [26][27][28][29][30][31][32][33][34][35][36][37][38][39][40][41][42] and references therein, for example). However, * econtreras@usfq.edu.ec in this work, we shall use the recently developed WKB approximation to the 13 th order which has brought the attention of the community [43]. It should be emphasized that, for the application of the method in the context of traversable wormholes, a bell-shape potential must be ensured. In this work, we study the QNM for the model after providing the suitable sets of parameters that ensure a bell-shaped potential.
This work is organized as follows. In section II we review the mains aspects related to traversable wormholes. Next, in section III we propose the embedding function, obtain the shape of the wormhole and analyse the quantifier of the exotic matter. In section IV we implement the 13 th order WKB approximation to compute and interpret the quasinormal modes associated to the scalar perturbations of the wormhole. Finally, in the last section we conclude the work.

II. TRAVERSABLE WORMHOLES
Let us consider the spherically symmetric line element with φ = φ(r) and b = b(r) the redshift and shape functions respectively. Assuming that (1) is a solution of Einstein's equations arXiv:2208.09048v1 [gr-qc] 18 Aug 2022 with κ = 8πG/c 4 1 , sourced by T µ ν = diag(−ρ, p r , p t , p t ) we arrive at In what follows we shall describe the main aspects of a traversable wormhole by its embedding in the three dimensional Euclidean space. First, note that as our solution is spherically symmetric, we can consider θ = π/2 without loss of generality. Now, considering a fixed time, t = constant, the line element reads The surface described by (6) can be embedded in R 3 where the metric in cylindrical coordinates (r, φ, z) reads Next, as z is a function of the radial coordinate we have from where ds 2 = 1 + dz dr 2 dr 2 + r 2 dφ 2 .
Finally, from (6) and (9) we obtain where is clear that b > 0 for r ∈ [r 0 , ∞). At this point some comments are in order. First, the wormhole geometry must be endowed with minimum radius which leads to dz/dr → ∞ as r → b 0 (that occurs when b = r). Accordingly, the existence of a minimum radius requires that at r = b 0 the shape function must be b = b 0 . Second, we demand that the solution is asymptotically flat which implies both, b/r → 0 (from where dz/dr → 0) and φ → 0 as r → ∞. Third, as as the conditions lim r→∞ dz dr = 0, 1 In this work we shall assume c = G = 1. must be satisfied, the smoothness of the geometry is ensured whenever the embedding surface flares out at or near the throat, namely which corresponds to the flaring-out condition. It is worth mentioning that, the flaring out condition (14) leads to the violation of the null energy condition (NEC) as we shall see in what follows. Let us define the quantity which can be written as Now, as (r − b) → 0 at the throat, we have so that Note that if ρ > 0 the above condition implies p r < 0 which entails that T 1 1 should be interpreted as a tension. Furthermore, if we define τ = −p r the flaring out condition leads to which implies that, for this exotic matter, the throat tension must be greater than the total energy density which violates the NEC, as we stated before. Although there is not evidence of exotic matter in the universe, we can minimize the amount required to construct a traversable wormhole by demanding that the quantifier [44] is finite.
From a technical point of view, the construction of traversable wormholes require solving the system (3)-(5), namely three equations with five unknowns. The strategy should be either supplying the metric functions that satisfy the geometric constraints listed above or giving one of the metrics and an auxiliary condition, namely an equation of state or a metric constraint. In this work we follow an alternative route which consists of proposing a suitable embedding function.

III. A WORMHOLE MODEL
In this section, we construct a traversable wormhole geometry by providing a general embedding function with the aim to integrate Eq. (10) and obtain the shape function b. Note that, although Eq. (10) can always be inverted numerically, in this work we look for analytical solutions so we propose where a, c and d are free parameters. The free parameters can be constrained by imposing both the existence of a throat (a minimum radius r 0 ) and the flaring out condition given by Eqs. (11) and (12) which lead to In figure 1 we show the embedding function for different values of the parameters involved. Note that the profiles flare out slowly as both c and d increases.
The shape function is obtained by replacing (21) in (10). As a result we obtain where Note that the conditions z(r 0 ) = 0 and b(r 0 ) = r 0 hold, as expected. In figure 2 we show the embedding function and the shape function for different values of the parameters involved. Note that, the solution is asymptotically flat as required.
In order to specify the wormhole metric completely, we must propose a suitable redshift function. In this work we shall take the simplest choice, namely φ = 0, so the solution has a vanishing radial tidal force. At this point, the matter sector can be completely specified but the expressions of the density and pressure are too long to be shown here. Nevertheless, a more interesting issue is the analysis of the quantifier of the exotic matter which is finite for all the parameters under consideration as shown in Fig. (3). Besides, for c = 0.4, the amount of exotic matter required to sustain the wormhole increases as d grows. In contrast, for c = {0.8, 1.2, 1.6} the decreasing of the amount of exotic matter is associated with an increasing of the free parameter d.

IV. QNM BY THE WKB APPROXIMATION
The perturbations of the TW can be carried out by adding test fields (scalar or vectorial) to the background or through perturbations of the space-time itself. However, independently of how it is performed, the equation governing the evolution of the perturbation can be reduced to a like-Schrödinger equation given by where r * is the tortoise radial coordinates defines as and V (r) is an effective potential. Notices that the tortoise coordinate r * is defined in the interval (−∞, ∞) in such way that the spatial infinity at both sides of the wormhole corresponds to r * = ±∞ and the wormhole throat is located at r * = 0. In this work, we will focus on the study of scalar perturbations, so the effective potential takes the form where L is called the multipole number or fundamental tone. The solution of (25) with the boundary conditions corresponding to purely out-going waves at infinity, are the QNM with frequency ω = Re(ω) + iIm(ω). The real part of of the QNM frequencies corresponds to the frequency of oscillation, while the imaginary part Im(ω) relates with the damping factor due to the loss of energy produced by the gravitational radiation. It is worth noticing that when Im(ω) > 0, the perturbation grows exponentially meaning an instability in the system. For the system to be stable, is required that Im(ω) < 0. Also consider that a complete perturbation will be a superposition of different tones L, which means that we need that every possible frequency satisfy Im(ω) < 0. The QNM frequencies are usually obtained by numerical methods but, in this work, we shall implement the WKB approach taking advantage of the similarity of Eq. (25) with the one-dimensional Schrödinger equation with a potential barrier. This method was first used by Schutz and Will in Ref. [45] to study scattering around black holes and has been extended to higher orders around the top of the bell-shaped potential [40]. Specifically, the 13 th order formula reads being V 0 the maximum height of the potential and V 0 its second derivative with respect to the tortoise coordinate evaluated at the radius where V 0 reaches a maximum which for a TW with a bell-shaped potential occurs at the throat (r = r 0 , r * = 0). The higher order corrections are encoded in Λ j which depend on the value of the potential its derivatives evaluated at the maximum.
The exact expressions for this corrections can be found in [43].
In what follows, we shall show numerical results for the QNM associated to scalar perturbations of the model. As we stated previously, the implementation of the WKB method requires a bell-shaped potential as a function of the tortoise coordinate as shown in Fig. (4) for different values of the parameters involved. We note that the peak of the potential decreases as l decreases. Besides, the potential spread out and its peak decreases as c grows.
In Fig. 5 we show the imaginary part of the frequency as a function of the parameter d. For c = 0.4 (first row), we note that the profile is monotonously increasing and approaches asymptotically to zero. In contrast, for n = 2 and n = 3, Im(ω) decreases with d, reach a minimum and grows again approaching asymptotically to zero. Particularly, for l = 6, the frequency has positives values for d ∈ (0.2, 0.3) which means that the wormhole is unstable for this interval. For c = 0.8 (second row), c = 1.2 (third row) and c = 1.6 (fourth row), Im(ω) is an increasing function and is always negative for the values of d under consideration so the solution can be considered as stable under scalar perturbation for this parameters. We also note that the value of Im(ω) decreases as the overtone n increases.
In Fig. 6 we show the Re(ω) as a function of d for different values of the parameter c. For c = 0.4 (first row), the profile reaches a minimum located at different values of d depending of the overtone. More precisely, as n increases, the location of the minimum shift to larger values of d. Besides, except for n = 3, the Re(ω) converge to the same value as d grows which means that the oscillatory behaviour of the signal in indistinguishable for each overtone for large d. The asymptotic line shift to bigger values of Re(ω) as l increases. For c = 0.8 (second row) except for n = 0 which remains constant, the Re(ω) reaches a maximum in contrast to the previous case. Moreover, the signal approaches asymptotically to a constant value which is different for each overtone. For c = 1.20 (third row) and l = 6, the signal is constant for n = 0 and increases monotonously for n = 1 and n = 2 and reach a minimum for n = 3. In contrast to the previous cases, the Re(ω) do not approach asymptotically to any value in the interval under consideration. For l = 7 and l = 8, except for the lowest overtone, Re(ω) increases monotonously. For c = 1.6, the signal is constant for n = 0, and increases monotonously for n = 1 and n = 2. For n = 3, the signal increases with but undergoes some oscillatory behaviour for large d. This behaviour can be associated to numerical instabilities.
Based on the previous results, at this point some comments are in order. First, for c = 0.4, the behaviour of signal for large d is almost independent of the overtone. In particular, the damping factor, given by e Im(ω) , is almost the same for each n under consideration. Similarly, the oscillatory behaviour is almost monochromatic: Re(ω) converge to the same value as l grows. In this regard, if the main goal is to differentiate the behaviour from overtones for c = 0.4, our model must be tested in the interval d ∈ (0.2, 1). It is worth mentioning that the model seems unstable in this interval for n = 3. However, it is a well known fact that the model works well for lower overtones so that such instability could be associated to inaccuracy of the method. Second, for c = 0.8 the behaviour clearly depends on the value of the overtone in the whole interval of d under consideration. In particular, the damping is both stronger for large n and weaker a d grows. Regarding the Re(ω), the behaviour of the oscillatory part depends on the value of d. Indeed, Re(ω) decreases as n grows in d ∈ (0.2, 0.4) but increases for large n in d ∈ (0.4, 2). Interestingly, all the frequencies of the oscillatory part coincide for d = 0.4. Finally, for c = 1.2 and c = 1.6, the pattern is clear: the damping is stronger as n increases and d decreases and the frequency of the oscillations grow as n grows.

V. CONCLUSIONS
In this work we obtained a traversable wormhole with vanishing radial tidal force by proposing a general embedding function with some free parameters. The parameters were reduced by imposing the basic requirements that must be satisfied by a wormhole geometry. In particular, we demanded the existence of a minimum radius (which defines the throat of the hole) and the flaring out condition. In order to explore how the geometry behaves in terms of the remaining parameter we analyzed both the quantifier of the exotic matter and the quasi normal modes of the solution. As a results, we observed that the solution requires a finite amount of exotic matter that decreases for certain values of the parameters involved. Besides, we obtained that in general the solution seems stable after scalar perturbations. Indeed, the imaginary part of the quasinormal frequencies remains negative which leads to a suitable damping factor for the signal.   5. Imaginary part of the frequency as a function of d for n = 0 (black), n = 1 (blue), n = 2 (green), n = 3 (red). We have set c = 0.4 (first row), c = 0.8 (second row), c = 1.2 (third row) and c = 1.6 (fourth row). For each row we have l = 6 (left panel), l = 7 (center panel) and l = 8 (right panel).  6. Real part of the frequency as a function of d for n = 0 (black), n = 1 (blue), n = 2 (green), n = 3 (red). We have set c = 0.4 (first row), c = 0.8 (second row), c = 1.2 (third row) and c = 1.6 (fourth row). For each row we have l = 6 (left panel), l = 7 (center panel) and l = 8 (right panel).