{'fraction_non_alphanumeric': 0.03639498722172615, 'fraction_numerical': 0.012193244449918404, 'mean_word_length': 4.510349507974211, 'pattern_counts': {'":': 0, '<': 42, '<?xml version=': 0, '>': 17, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 5, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Traversable wormholes in Einsteinian cubic gravity

In the present work we investigate wormhole configurations described by a constant redshift function in Einstein-Cubic gravity ({{\sf ECG}}). We derive analytical wormhole geometries by assuming a particular equation of state ({{\sf EoS}}) and investigate the possibility that these solutions satisfy the standard energy conditions. We introduce exact asymptotically flat and anti-de Sitter (AdS) spacetimes that admit traversable wormholes. These solutions are obtained by imposing suitable values for the parameters of the theory so that the resulted geometries satisfy the weak energy condition ({\sf WEC}) in the vicinity of the throat, due to the presence of higher-order curvature terms. Moreover, we find that AdS solutions satisfy the {\sf WEC} throughout the spacetime. A description of the geodesic motion of timelike and null particles is presented for the obtained wormhole solutions. Also, using gravitational lensing effects, observational features of the wormhole structure are discussed.


Introduction
A wormhole is a hypothetical topological feature of the spacetime which creates shortcuts between two distinct spacetimes. The concept of "wormhole" for these structures was invented for the first time in 1957 within the seminal papers of Misner and Wheeler [1] and Wheeler [2], in order to provide a mechanism for having "charge without charge". The study of Lorentzian wormholes in the context of General Relativity GR dates back to the fundamental works of Morris and Thorne in 1988 [3] where, introducing a static spherically symmetric line element, they showed that, exact solutions representing wormhole geometries could be found by solving the Einstein field equation. According to GR, the fundamental faring-out condition of throat causes to the violation of null energy condition (NEC) and thus, traversable wormholes are only possible if the so called "exotic matter" [4] exists at their throat, which involves an energy-momentum tensor (EMT) violating the NEC. This condition is in turn a part of the WEC that the physical meaning of which is that the energy density be non-negative in any reference frame. In this respect, traversable wormhole geometries have been obtained using the exotic matter distribution, e.g., with the help of phantom energy distribution [5]. This type of matter, though exotic in the laboratory context, is of observational interest in cosmological scenarios [6]. Phantom energy possesses peculiar attributes, namely, a divergent cosmic energy density in a finite time [7], prediction of existence of a new long range force [8], and the appearance of a negative entropy and negative temperature [9].
One of the most significant issues in wormhole geometries is the fulfillment of standard energy conditions. In this regard, many attempts have been devoted in the literature in order to find some realistic material sources that support the wormhole configuration or minimize the employment of exotic matter. Work along this line have been done to investigate the construction of thin-shell, dynamical, and rotating wormholes [10]. However, we concentrate, in the herein model, on modified theories of gravity where the effective scenario could provide static spherically symmetric solutions representing traversable wormholes without resorting to exotic matter source. In this context, it has been shown that the presence of higher order terms in curvature would allow for constructing thin-shell wormholes supported by ordinary matter [11]. Moreover, wormholes in modified gravity involving higher-order curvature invariants can satisfy the energy conditions at least at the throat [12] and throughout the spacetime [13]. The study of wormhole solutions in the framework of modified gravity has become a main focus of interest in modern cosmology. This topic has attracted much attention in recent years and a large amount of work has been devoted to this issue among which we quote: spacetimes admitting wormhole solutions in Brans-Dicke theory [14], f (R) gravity [15], Born-Infeld theory [16], Einstein-Gauss-Bonnet theory [17], Kaluza-Klein gravity [18], Rastall gravity [19], scalar-tensor gravity [20] wormhole solutions in the presence of a cosmological constant [21], non-commutative geometry [22] and other modified gravity theories [23].
Over recent years there has been a remarkable interest in the subject of higher curvature gravity, that the main focus of much has been motivated by the wish to account for quantum phenomena within gravitational field. In [24], it was shown that including higher curvature terms within the Einstein-Hilbert action can provide a setting towards a renormalizable theory of gravity. It is generally expected that higher order terms show up themselves, e.g., within the renormalization of quantum field theory in curved spacetime [25], or in the construction of low-energy effective action of string theory [26] such as Lovelock theory [27]. Moreover, in order that the gravitational action be free of ghost terms, the quadratic curvature corrections to the Einstein-Hilbert action should be proportional to the Gauss-Bonnet term [28]. Also, gravity theories with higher curvature corrections have attracted a great deal of interest in holography [29] and considered in the context of cosmology [30].
A particular proper feature of a specified higher-order gravity is its linearized spectrum, that is, the set of physical degrees of freedom propagating through metric perturbations in vacuum. For instance, in the framework of holography, useful information about the corresponding holographic CFT stress-energy tensors can be extracted by the virtue of linearized equations of a given higherorder gravity as these equations are dual to the metric perturbation, see e.g., [31]. There are some higher-order gravity theories such as quasi-topological gravity [32] and f (Lovelock) [33] which are equivalent to GR at the linearized level in the vacuum, i.e., the only physical mode propagated by the metric perturbations is a transverse and massless graviton. However, the inconvenient characteristic of most of higher-order gravity theories is the dependency of the couplings of the different curvature invariants on spacetime dimension and in fact, these theories are different gravity theories in different dimensions. Nevertheless, the only known theories that their couplings are independent of spacetime dimension and share spectrum with GR are Lovelock theories [27]. In this regard, it is shown that, up to cubic order corrections in curvature, there is only one additional theory which satisfies this criterion. Furthermore, this theory is non-trivial in four dimensions, unlike the quadratic and cubic Lovelock theories. Recently, a new model of higher order gravity has been presented in [34] to which, when quadratic and cubic Lovelock terms are added, is the unique cubic model for gravitational interaction that shares its graviton spectrum with GR. Moreover it is shown that this theory has dimensionindependent coupling constants and is free of massive gravitons in general dimensions [35]. In the context of this theory, which is known as Einsteinian cubic gravity (ECG), black hole solutions in fourdimensions [36] and in higher dimensional spacetimes have been studied [37]. The asymptotically-AdS black brane solutions of ECG have also been considered in [38]. Motivated by the above considerations we search for static spherically symmetric solutions representing wormhole configurations in ECG and study the effects of higher curvature terms within the wormhole structure. This paper is organized as follows: In section 2, we give a brief review on gravitational field equations in ECG, and proceed with analyzing the the energy conditions using a general form for spacetime metric of a worm hole. In section 3, we take an EoS for the radial and tangential pressures and search for exact wormhole solutions along with checking the satisfaction of energy conditions. Our conclusion is drawn in section 5.

Action and Field equations
The action of pure Einsteinian Cubic gravity in four dimensions is given by [34] S =ˆd 4 x √ −g (L g + L m ) , where L m and L g are the material and gravitational parts of the action, respectively, κ = 8πG is the gravitational coupling constant and λ is a dimensionless coupling constant. Also, R is the Ricci scalar and the new Einsteinian cubic gravity term P contributes to the action as cubic curvature terms in four dimensions, defined as where R abcd is the Riemann curvature tensor constructed out of Christoffel symbols. The field equations of ECG can then be written as [34,35] where P abcd = ∂L ∂R abcd , and T ab is the EMT of matter fields. In this work, we are interested in wormhole solutions in fourdimensions. Therefore, we proceed with employing a general static spherically symmetric line element which represents a wormhole and is given by where ψ(r) being the redshift function and b(r) is the wormhole's shape function. The shape function must satisfy the flare-out condition at the throat, i.e., we must have b ′ (r 0 ) < 1 and b(r) < r for r > r 0 in the whole spacetime, where r 0 is the throat radius. In the present work, we consider ψ(r) = 0 in order to ensure the absence of horizons and singularities throughout the spacetime. The EMT of matter is given by the following diagonal form where ρ(r) is the energy density and p r (r) and p t (r) are the radial and tangential pressures, respectively. Using field equation (3) along with taking into account the metric (5) we get where a prime denotes derivative with respect to the radial coordinate r. It is notable that the radial and tangential pressure profiles are exactly as in the GR case and the higher order curvature terms within the action and field equations (3) of the background theory contribute only to the energy density. This is due to the assumption for vanishing redshift function and for ψ(r) = 0, one gets more complicated expressions for pressure profiles. In view of the local energy conditions, we examine the WEC, which asserts that T ab U a U b ≥ 0 where U a is a timelike vector field. For the diagonal EMT (6), the WEC implies ρ ≥ 0, ρ + p r ≥ 0 and ρ + p t ≥ 0. Note that the WEC implies the NEC. Using equations (7)- (9), we arrive at the following relations One can easily show that for λ = 0 the NEC, and consequently the WEC are violated at the throat (ρ + p r < 0), due to the flaring-out condition [1,2]. Note that at the throat, one verifies which shows that for λ = 0 the NEC, and consequently the WEC, are violated at the throat. In order to impose ρ + p r > 0 in ECG, it is now possible to find an adequate range for the model parameters such that the NEC is satisfied at the throat. In the following section, we search for exact wormhole solutions in ECG and investigate in detail the conditions under which the WEC is respected.

Wormhole Geometries
In this section we present exact wormhole solutions in the context of ECG. Here, we have a system of differential equations (7)- (9) with four unknown functions ρ(r), p r (r), p t (r) and b(r). In order to determine spacetimes admitting wormhole structures, we must adopt a strategy for specifying the shape function b(r). We may also consider a pre-determined form for the functionality of b(r) and consequently obtain the EMT components. Also, it is common to consider a specific form for the equation of state (EoS), which provides a relation between the EMT components, namely, a linear EoS [39], the traceless [40] and EoS in general form [41]. However a linear EoS involving energy density leads to a complicated differential equation. This is due to the presence of third order derivative for b(r). Thus, we proceed with a less complicated differential equation for the shape function and consider an EoS relating the tangential and radial pressures, so that the energy density is determined through Eq. (7). Here, we consider an EoS of the form p t = αp r [42], whereby using Eqs. (8)- (9) we obtain the following differential equation for the shape function for which the solution reads where the integration constant C 0 can be determined from the condition b(r 0 ) = r 0 at the throat and is given by The flare-out condition at the throat results in the following inequality and the energy density for this solution is found as where We also get the following expressions for radial and tangential pressure profiles, as From the expressions (17), (22) and (23) we can find the energy conditions at the throat as where In what follows, we discuss the properties of the solution (14) in detail.

Specific case: Λ = 0
Setting Λ = 0 in Eq. (14), we get the metric function as We note that condition b ′ (r 0 ) < 1 in (16) implies that α < 0. It is then clear that these solutions are asymptotically flat, i.e., b(r)/r tends to zero as r → ∞. For Λ = 0 and in the limit of large values of r coordinate, the ξ 0 coefficient in (17) can be ignored in comparison to the last term and hence the energy density and pressure profiles for this case read This is effectively the GR limit of ECG which can be achieved in the limit of λ → 0. We therefore observe that for large values of r, the quantities ρ(r) and ρ + p t are positive for α > − 1 2 and α > −1, respectively. Hence, the tangential component of WEC is satisfied for −1/2 < α < 0. However the value of ρ + p r is always negative in the limit of large r and thus, the radial component of WEC is violated in this limit. Substituting Λ = 0 into Eqs. (24)- (26) we obtain, at the throat, (4α 2 − 5α + 1) > 0 for α < 0; thus the condition ρ + p r > 0 at the throat is respected for λ > 0. Therefore, in order to satisfy the WEC at the throat, we can choose suitable values of positive λ such that ρ, ρ + p r and ρ + p t get positive values. It can also be seen that ρ, ρ + p r and ρ + p t have real roots (r i ), where their values are positive in the interval r 0 < r < r i . The values of r i correspond to positive real roots of the following equations where In this case, we obtain r 1 = r 2 = r 3 = r 0 4α−3 5α−3

1/2α
when λ is very large. The quantities ρ, ρ + p r and ρ + p t are sketched in Fig. (1). Note that the components of EMT tend to zero as r tends to infinity. The left panel shows that for λ > 0 the WEC (and also NEC) is violated in the vicinity of the wormhole throat, but for λ < 0, it can be satisfied near the wormhole throat as it is shown in the right panel. It is worth mentioning that for black hole solutions in ECG, the coupling parameter λ can affect the mass, entropy, Hawking temperature and horizon radius of the black hole, e.g., the horizon radius with λ > 0 is greater than the Schwarzschild value [36]. In [43], possible observational implications of ECG have been explored and it is shown that, for spherically symmetric black holes, the corrections due to ECG are remarkable on the scales near the horizon. However, the λ parameter here affects only on the behavior of energy density (and consequently WEC) and faraway from the wormhole configuration the ECG corrections in energy density disappear. It is therefore reasonable to choose those values of λ parameter for which the WEC is satisfied, at least at the wormhole throat.

General case: Λ = 0
Let us proceed with the general case in which Λ = 0. Firstly, let us rewrite the metric function (14) as We note that this solution does not correspond to an asymptotically flat spacetime, however, one may match this solution to an exterior vacuum geometry [44]. The flare-out condition (16) implies that we must have Λ > α (α−1)r 2 0 for α > 1 and Λ < α (α−1)r 2 0 for α < 1. In what follows, we analyze the physical properties and characteristics of these wormholes along with checking the conditions under which the WEC can be satisfied.
For α > 1, it is seen that for a suitable value of λ < 0, the WEC is satisfied throughout the spacetime. For this case, the right panel of Fig. (2) displays the behavior of ρ , ρ + p r and ρ + p t that, as is seen, these quantities are positive throughout the spacetime, implying that the WEC is satisfied for all values of r. In case we take α < 0 the value α (α−1)r 2 0 is positive and therefore for −∞ < Λ < Λ c = α (α−1)r 2 0 the condition b ′ (r 0 ) < 1 is satisfied. To be a solution of a wormhole, the condition 0 < r − b(r) is also imposed. For the range of values 0 < Λ < Λ c , condition b(r) = r leads to two real and positive roots given by r − = r 0 and r + = r 0 1 + 1 1) , and thus the spatial extension of this type of wormhole solution cannot be arbitrarily large. We then have a finite wormhole within the range r 0 < r < r 0 1 + 1 The left panel in Fig. (3) shows that an increase in the value | α | enlarges the wormhole spatial extension. Note that for −∞ < Λ < 0, we have Ads wormhole solution. In order to study energy conditions for these class of solutions we proceed with obtaining the behavior of quantities ρ, ρ + p r and ρ + p t at infinity. All of the quantities tend to the value Λ as r tends to infinity. Therefore, in the large limit of r, for 0 < Λ < Λ c , the WEC is satisfied. The right panel of Fig. (3) shows that for λ < 0 the WEC (and also NEC) is violated in the vicinity of wormhole throat, but for λ > 0, these conditions can be satisfied near the throat. In order to have normal matter at infinity, one can choose the positive value for Λ. We also note that the radius of throat can be determined by imposing energy conditions at the throat. To do so, one can use expressions (24)-(26) to find three inequalities. Depending on the model parameters, these inequalities then decide the minimum and maximum values for r 0 . Finally, we can see that for 0 < α < 1 the term α (α−1) is negative and therefor we have wormholes for negative value of Λ. We plot the quantities ρ , ρ + p r and ρ + p t in the left panel of Fig. (2) . All the quantities tend to Λ c = −0.2 as r tends to infinity. We have considered a suitable value for λ < 0 in order to have normal matter distribution in the vicinity of the throat.

Particle Trajectories Around the Wormhole
In this section we analyze geodesic equations in ECG wormhole spcetime described by the metric Eq. (14), using the Lagrangian formalism [45]. Due to the spherical symmetry, it suffices to consider the equatorial plane θ = π/2. The corresponding Lagrangian for metric (14) is then found as where a dot denotes derivative with respect to the affine parameter η. As the Lagrangian is constant along a geodesic we can consider L(x µ ,ẋ µ ) = ǫ so that time-like and null geodesics correspond to ǫ = −1 and ǫ = 0, respectively. Using the Euler-Lagrange equation, one can readily identify the following constants of motion − 2ṫ = −2E and 2r 2φ = 2L, where E is the energy and L the angular momentum of the test particle. Inserting these constants of motion into (35) we getṙ It is convenient to rewrite Eq. (38) in terms of the proper radial distance which is finite for all finite values of r. Note that the spacetime is extended in such a way that l monotonically increases from −∞ to +∞. l < 0 or l > 0 correspond to two parallel universes joined together via a throat at l = 0. Using the proper radial distance, Eq. (38) takes the simple forṁ where the effective potential is defined as In what follows, we discuss the trajectory of particles around the wormhole, using the above form for the effective potential. In fact, geodesic equation (40) can be interpreted as a classical scattering problem with a potential barrier V eff (L, l). Moreover, using Eq. (37) we can rewrite Eq. (40) as an ordinary differential equation for orbital motion We note that, in traversable wormhole spacetimes, particles can travel through the throat of the wormhole from one asymptotically flat part of the manifold to other one. Then, a geodesic can pass through the throat into the other universe if Similarly, for a geodesic reflected back on the same universe by the potential barrier, we have E 2 < V eff (L, 0). In this case, there is a turning point at l = l tu which is obtained by solving the following equation A generic feature of this effective potential is that it possesses a global maximum at the throat The flaring out condition leads to d 2 V eff dℓ 2 < 0 at the throat. This clearly has an unstable orbit since it occurs at the maximum of the potential for E 2 = V eff (L, l 0 ). We note that these conditions are independent of whether the geodesic is null or timelike. We now consider the wormhole solutions presented in section (3.1) and restrict ourselves to the class of wormholes with Λ = 0. Setting the shape function b(r) = r (r/r 0 ) 2α in Eq. (39), we find There is not an explicit representation for the effective potential V eff (L, l) as a function of l, since Eq. (46) cannot in general be solved for the radial coordinate as a function of proper distance. However, this issue can be numerically treated using standard techniques. Nevertheless, for α = −1 we can obtain the behavior of effective potential as a function of the proper distance for timelike and null geodesics, see Fig. (4). We are interested in two kinds of trajectories for the obtained wormhole solutions. In the the first kind, test particles are not allowed to pass through the throat and remain within the original universe, whereas for the geodesics of the second kind the test particle passes through the throat from lower universe to the upper one. In order to visualize the trajectory of particles, we must solve numerically Eq. (42) with specified initial conditions. We note that the constants of motion E and L can be expressed by appropriate initial conditions. In Figs. (5) and (6) we have plotted trajectories of particles in the embedding diagram of the wormhole for timelike and null geodesics. In these plots, radius of the throat has been set as r 0 = 2 and the initial value for proper distance is fixed at l i = 4. In each graph, two kinds of trajectories for null Fig. (5) and timelike Fig. (6) geodesics are plotted with different initial conditions. We would also like to mention that although, for the present wormhole solutions, the effects of coupling parameter may not appear directly in the behavior of particle trajectories around the wormhole, these effects can be detected using the Raychaudhuri equation [46], where the λ parameter, by the virtue of energy density profile, would play its own role in the expansion rate of a congruence of timelike or lightlike geodesics. Not straightforwardly, but rather presenting additional information on the role of ECG coupling parameter, more general wormhole solutions can be obtained considering nonzero redshift function. In this case, the effects of λ parameter on particle trajectories can be better detected and the corresponding results can be compared with those of GR. We hope to address these issues and other questions in future work. A wormhole configuration has its own influences on its environment, an important aspect through which observational signatures of the wormhole can be surveyed. It is well known that the gravitational lensing is now one of useful tools to seek not only for dark and massive objects, but also wormholes. In this regard, we end this section by studying lensing effects of the obtained wormhole solutions. Let us consider a beam of light incoming from infinity, reaching the closest approach distance r c from the center of the gravitating body and then emerges in another direction. Therefore, a light ray approaching the wormhole with the shape function b(r) = r 0 (r/r 0 ) 2α+1 is scattered and bends with the deflection angle [47] Θ(r c , α) = −π + 2ˆ∞ rc r −(1+α) r c dr   where use has been made of Eq. (38) for null rays (ǫ = 0) along with the second part of Eq. (37). The quantity r c is defined as the distance at which dr/dφ = 0 and for the present model, it is related to the impact parameter β = L/E as r c = β. Figure 7 shows how the light deflection angle depends upon the ratio of tangential to radial pressures, i.e., the EoS parameter α, and the distance of closest approach (or correspondingly the impact parameter). We then observe that for r 0 < r c < ∞ the quantity Θ(r c , α) is positive and finite. The deflection angle tends to zero as r c → ∞, i.e., light ray is unaffected by the gravitating object. Moreover, the deflection angle increases as r c → r 0 and diverges at the wormhole throat where an unstable photon sphere is present 1 [49]. Consequently, the wormhole structure may produce infinitely many relativistic images of an appropriately placed light source. This infinite sequence corresponds to infinitely many light rays whose limit curve asymptotically spirals towards the unstable photon sphere [50]. Since the photon sphere is located at the throat, such a sphere may be detectable providing thus a setting to search for observational evidences of the wormhole. The integral (47) cannot be exactly solved for arbitrary values of EoS parameter, however we can evaluate it in the weak field limit, i.e., far away enough from the gravitating body where r 0 ≪ r c . Introducing new variable y = 1 − r c /r the integral (47) can be re-expressed in the following form The throat radius has been set as r 0 = 1.
In the weak field regime where ζ ≪ 1, we expand the numerator and perform the integration to get From the above solution we find that when α → −1, the deflection angle Θ(α) → πζ 2 /4 + 9πζ 4 /64. Hence, in the limit where radial and tangential pressures cancel each other, the lensing characteristics of the wormhole structure mimic those of Ellis wormhole [51].

Concluding Remarks
In the present study we constructed models of static wormholes with constant redshift function within the frame work of ECG by considering ordinary matter distribution. Taking the EoS of the matter fluid as a linear relation between radial and tangential pressures, we obtained exact traversable wormhole solutions which can be asymptotically flat, de-Sitter (dS) or AdS. For asymptotically flat solutions, the energy conditions in the vicinity of the throat can be satisfied by choosing suitable positive values of λ parameter. Moreover, we showed that for λ > 0, the WEC can be satisfied near the throat for r 0 < r < r 3 and a growth in the value of λ parameter leads to increasing this range for radial coordinate. For dS solutions, the energy conditions are fulfilled throughout the spacetime for α > 1 and λ < 0. For AdS solutions, the energy conditions are satisfied for 0 < α < 1 and negative values of λ parameter. Also, for λ > 0 the energy conditions can be satisfied at the wormhole throat for both dS and AdS solutions provided that α < 0. It is noteworthy that asymptotically AdS wormhole solutions can be constructed in the context of holography. In [52], wormhole solutions with two asymptotic AdS boundaries are introduced and the effects of multi-trace deformations of the boundary quantum field theory on the wormhole state have been investigated. As discussed in [3], the imposition of causality and energy conditions prevents the Einstein-Rosen bridge from being traversable. However, as a result of the relationship between quantum information and quantum gravity in holography [53] and particularly in the framework of Ads/CFT correspondence [54], traversable wormhole configurations can be built in Ads/CFT by inserting a double-trace deformation on the boundary CFTs [55]. In this scenario quantum matter fields can provide the necessary negative energy (generated by explicit couplings between the two dual boundary CFTs) to retain the throat of the wormhole open and thus achieve its traversability. Nevertheless, traversable wormhole solutions presented here are obtained for classical matter fields that respect energy conditions for specific values of the model parameters.
We also obtained embedding diagrams for wormhole solutions and the behavior of timelike and null geodesics in wormhole configuration were discussed. The study of wormhole geodesics is of great importance as it could help us to figure out possible observational effects that arise as a result of the scattering of particles on wormholes [56]. Likewise, such an investigation provides a promising way to test modified theories of gravitation as well as sufficient incentives for scientists towards probing wormhole structures in the universe. Work along this line has attracted many researchers to investigate observational aspects of a wormhole such as gravitational lensing [51,49] and microlensing [57] effects, accretion disks around wormholes [58], wormhole shadows [59] and several observables such as rotation curves [60]. Traversable wormhole solutions in GR can be obtained by considering some form of exotic fluid as the supporting matter for wormhole geometry. However, the situation could be completely different in the context of modified gravity theories where a modified EMT could provide a suitable setting for traversable wormhole solutions, without the need of introducing any form of exotic matter. This objective was pursued within the framework of ECG gravity and as we saw, depending on model parameters, a domain of existence of traversable wormhole solutions which respect energy conditions can be obtained. Beside the present model some other modified gravity theories have been investigated, whose geometrical attributes (not present in GR) may provide a setting for investigating wormhole solutions without resorting to exotic matter. In this regard, much attempt has been made during the last decades to investigate the structure of wormholes in alternative theories of gravity such as, Einstein-Cartan theory [61], modified gravity theories with curvature-matter coupling [62] and multimetric gravity theories [63]. It is also worth mentioning that the solutions obtained in the present work can be generalized to the case of charged wormhole configuration. In this case, the Lagrangian for Maxwell field has to be added to the material part of the action (1). Therefore, adding the extra source of matter (electromagnetic field) could enhance the degrees of freedom of the model (e.g., charge distribution may support the model against violation of WEC and NEC at the throat), providing then, a new class of traversable wormholes that respect energy conditions. Moreover, adding rotation to the present wormhole configurations may provide a different scenario in which, depending on the model parameters, wormhole structures with different properties (as compared to those of GR) including their mass, angular momentum, quadrupole moment, and ergosphere can be constructed [64]. However, dealing with the resultant field equations for the rotating case may not be an easy task and a more detailed study on this issue would be desirable.
Finally, we would like to point out that, the stability or instability of celestial bodies is an important issue in gravitational physics. A traversable wormhole structure is of physical importance if it is stable against perturbations. In this regard, the stability of static wormholes has been surveyed utilizing specific equation of state or by considering a linearized radial perturbations around a static solution. Work along this line has been carried out for a thin-shell wormhole constructed from the Schwrazschild spacetime [65]. Stability of thin-shell wormhole in Einstein-Maxwell theory with a Gauss-Bonnet term has been studied in [66] and stable and unstable static solutions of spherically symmetric thin shell wormholes supported by generalized Chaplygin gas has been investigated in [67], see also [64,68] and references therein. An alternative technique for exploring the stability issue of wormhole structure is the thermodynamic stability. In this method, if the thin shell wormhole fulfills the first law of thermodynamics such that it respects the thermodynamic stability conditions in a successful way, and admits a meaningful heat capacity, then it can be considered as a stable wormhole configuration. Such an investigation has been performed for thin shell structures [69] and thin shell wormholes [70]. In this regard, thermodynamic properties of the herein wormhole configuration could provide a testbed for thermodynamic stability of the wormhole solutions, however, this study is beyond the scope of the present article and the results of future investigations will be reported as an independent work.