A Review of Stable, Traversable Wormholes in f ( R ) Gravity Theories

,


Introduction
The first study of wormhole physics was done by Ludwig Flamm in 1916 [1] during his research into the Schwarzschild solution to the Einstein field equations.The next solution resembling a wormhole, called the "Einstein-Rosen bridge" [2], was an idea that arose during the investigation of blackhole spacetimes by Einstein and Rosen in 1935.They discovered that, at least theoretically, it was possible for a blackhole surface to act as a bridge that connected to a remote patch of spacetime.The putative wormhole in the Einstein-Rosen bridge is colocated with the blackhole's singularity, so it is not traversable.The surface of the blackhole is the event horizon, which cannot be a wormhole.
A good review of the historic development of research into traversable wormholes can be found in [3].A more recent concise summary of wormholes in extended theories of gravity can be found in [4].In 1957, Misner and Wheeler [5] first introduced the term "wormhole" during their analysis of topological issues in General Relativity (GR).They extensively analyzed the Riemannian geometry of manifolds of nontrivial topology.This is where the phrase "physics is geometry" arose.Wheeler also suggested that the geometry of spacetime might be constantly fluctuating, and it may induce a change in topology to form microscopic wormholes.The first traversable wormhole called the "Ellis drainhole" was proposed by both Ellis and independently by Bronnikov in 1973 [6-8].The "Ellis drainhole" spacetime is a static, spherically symmetric solution of the Einstein field equations in a vacuum, and it includes a scalar field ϕ minimally coupled to the spacetime geometry.
The concept of traversable wormholes, which allow inter-and intrauniverse travel by humans, was first introduced by Morris and Thorne in their classic 1988 paper [9].It is well known that the Weak and Null Energy Conditions (WECs and NECs) [10] must be violated by the stress-energy tensor as a minimum for a stable human traversable wormhole.The other two energy conditions are the Strong Energy Condition (SEC) and the Dominant Energy Condition (DEC) [11].The violation of the NEC is needed due to the flare-out condition requirement, i.e., the throat should open up outward as a human travels through the wormhole [9,10].Quantum field theory permits these violations due to the quantum coherent effects found in any quantum field [12].The traversable wormhole solutions have geometries that allow for closed timelike curves and "effective" superluminal travel without surpassing the speed of light locally [3].There is a claim [13] in the published literature about stable traversable wormholes that can be constructed within the framework of Einstein's General Relativity (henceforth called 'Einstein Gravity' in this paper) but without the need for exotic matter.This claim has been refuted by other authors [14].Some of the suggested reasons in [14] are the following: (i) possible error in calculations; (ii) failing to check for violations at the throat, such as divergence of the inverse of the metric; and (iii) failing to check for discontinuities in the exterior curvature in the vicinity of the throat based on the thin-shell formalism.Some simple examples of traversable wormholes such as the polyhedral wormholes that do not require spherical symmetry were given by Matt Visser in 1989 [15,16].These wormholes also require the presence of exotic matter.

Motivation and Equation of State
There are many attempts to use various modified gravity theories to check the existence of stable traversable wormholes in these theories.In many of these studies [17][18][19], the Morris-Thorne (MT) traversable wormhole has been used to see the effects of a modified gravity background on the stability of the wormhole.One of the main motivations for studying wormholes in modified gravity theories is to resolve the problem of the need for exotic matter to stabilize wormholes in GR.Modified gravity theories are also used to construct viable cosmological models of the Universe and to explain singularities encountered in cosmology.In modified gravity theories, the stress energy tensor is replaced by an effective stress energy tensor that contains curvature terms of higher order introduced due to modifications to GR. f (R) gravity [20,21] is one such leading theory, and the review of wormholes in f (R) gravity is one of the main goals of this paper.
While evaluating wormholes in GR and various modified theories, the choice of an equation of state (EOS) becomes critical.The EOS is an equation that relates the radial pressure p r (r) and tangential pressure p t (r) to the energy density ρ(r).One of the most common EOSs is p(r) = ωρ(r).The motivation for using such an equation comes from cosmology.Observation of the accelerated expansion of the Universe can be explained by an EOS with ω < − 1  3 [22].For the case where ω < −1, the EOS is known as the phantom energy EOS.Phantom fluid violates the NEC and is a good candidate to be used in the study of wormholes [23] in various modified gravity theories.According to a study [24], if a source with such phantom fluid dominates the cosmic expansion, the Universe may end up in a Big Rip singularity in which the phantom energy rips apart the galaxies, solar systems, planets, and, eventually, the molecules, atoms, nuclei, and nucleons that we are made of, leading to the death of the Universe.

Astronomical Observational Signatures of Traversable Wormholes
The presence of naturally occurring negative energy regions in space is predicted to produce a unique signature corresponding to lensing, chromaticity, and intensity effects in micro-and macrolensing events on galactic and extragalactic/cosmological scales [25][26][27][28][29][30].It has been shown that these effects provide a specific signature that allows for discrimination between ordinary (positive mass-energy) and negative energy lenses via the spectral analysis of astronomical lensing events.The theoretical modeling of negative energy lensing effects has led to intense astronomical searches for naturally occurring traversable wormholes in the universe.Computer model simulations and comparison of their results with recent satellite observations of gamma ray bursts (GRBs) have shown that putative negative energy (i.e., traversable wormhole) lensing events very closely resemble the main features of some GRBs.
When background light rays strike a negative energy lensing region, they are swept out of the central region, thus creating an umbra region of zero intensity [25].At the edges of the umbra, the rays accumulate and create a rainbow-like caustic with enhanced light intensity.The lensing of a negative energy region is not analogous to a diverging lens because, in certain circumstances, it can produce more light enhancement than does the lensing of an equivalent positive mass-energy region [25].Real background sources in lensing events can have nonuniform brightness distributions on their surfaces and a dependency of their emission on the observed frequency.These complications can result in chromaticity effects, i.e., in spectral changes induced by differential lensing during the event.The quantification of such effects is quite lengthy, somewhat model dependent, and with recent application only to astronomical lensing events.The rest of this paper is organized as follows: In Section 2, we provide an overview of various modified gravity theories with just enough background necessary to analyze traversable wormholes in such theories.In Section 3, we review the basics of Morris-Thorne (MT) wormhole stabilization in GR.We end Section 3 with a summary of the steps used to analyze MT wormholes in modified gravity theories.In Section 4, which forms the core of this paper, we review the analysis of wormhole stabilization in various f(R) theories.In Section 5, we review some closely related topics such as Casimir wormholes, thin-shell wormholes, natural dark matter halo wormholes, and wormholes in nonlocal theories of gravity (NLGs).We end this paper with a summary of key discussion points in this paper.

Modified Gravity Theories
There have been efforts to construct stable traversable wormholes in f (R) theories, including Lovelock gravity.In this section, we give a brief overview of the basics of these theories to provide us with enough background to explore in the next section the properties of traversable wormholes.

An Overview of Modified Gravity Theories
The various modified gravity theories that we survey here are f (R) theories, f (T ) theories, f (R, T ) theories, f (G) theory, f (R, L m ) theory, f (Q) theory, f (Q, T ) theory, Lovelock gravity (a special case of f (R) theory), Einstein-Gauss-Bonnet (EGB) gravity (another special case of f (R) theory), Brans-Dicke theory, and Kaluza-Klein (KK) theory.
Lovelock gravity and EGB gravity are f (R) theories that include higher-order curvature terms, and they are also applicable to higher dimensional spacetimes.For example, in EGB gravity [10], the authors replace the 2-sphere in the MT wormhole with an (n-2)-sphere, and thus the MT wormhole line element is modified as follows: The standard Einstein-Hilbert action of GR is In f (R) theories [31], the Ricci scalar R is replaced by a function of the Ricci Scalar f (R) as follows: where f (R) is a function of the Ricci scalar, and g is the determinant of the metric.Examples of commonly used functions f (R) are f (R) = R, f (R) = R + αR n , and f (R) = R + αe (βR) .
In teleparallel gravity [32,33], the Ricci scalar R is replaced by the torsion scalar T as follows: In GR, we assume spinless particles to follow the geodesic of the underlying spacetime, and hence, we have only R in the action and no T .In teleparallel gravity, T replaces R.This interpretation holds true in the case where we see teleparallel gravity as a gauge theory [34] for the translation group.
The f (T ) gravity [35] is an additional modification to teleparallel gravity, in which T is replaced by the torsion function-which is a function of T , namely f (T ), as follows: and we can obtain the modified field equations by varying the action S with respect to the metric in the same way it is done in Section 2.2 of this paper for f (R) theory.
T is a fundamental geometric quantity in the context of theories of gravity that involve spacetime torsion.Unlike GR, where spacetime curvature plays a central role in describing gravity, theories that incorporate torsion consider the twisting or nonmetricity of spacetime, and the torsion scalar is a measure of this twisting.The torsion scalar is given by the equation where S αβµ is the contorsion tensor, which is the difference between the affine (Levi-Civita) connection's components and the Weitzenböck connection's components.The contorsion tensor quantifies nonmetricity, and it describes how spacetime is twisted.The indices α, β, and µ refer to spacetime coordinates.In teleparallel gravity theory, the teleparallel connection is used instead of the metric affine connection that is used in GR.
In the f (R, T ) theory, T is the trace of the energy-momentum tensor T µν .Similar to f (R) theories [17], the presence of f (R, T ) in the action leads to changes in gravitational dynamics as compared to Einstein gravity.These modifications can have implications on the behavior of gravitational fields in various contexts.They also have consequences in cosmology and gravitational lensing.An example of an f (R, T ) function used to stabilize a wormhole is where α and λ are constants.In f (R, T ) theory, wormhole solutions with normal matter are feasible when appropriate shape functions are used.The coupling parameters α and λ in the action of f (R, T ) theory play an important role in accommodating the composition of matter.According to [17], when α < 0, wormholes exist in the presence of exotic matter, and when α > 0, wormholes exist even in the absence of exotic matter.
In Brans-Dicke theory [36], a scalar field ϕ is introduced to modify the standard Einstein-Hilbert action as follows [37]: where ϕ is the Brans-Dicke scalar field, ϕ ;µ is the covariant derivative of ϕ, and ω is a coupling constant that couples the Brans-Dicke scalar field ϕ with the gravitational field.L m is the Lagrangian density for the matter field(s).
The field equations in Brans-Dicke theory can be obtained as always by varying the action.These field equations are modified as compared to Einstein's field equations in GR, and they can be expressed in terms of the scalar field ϕ and its derivatives, as well as the metric tensor, the curvature tensors, and the scalars.The value of the coupling parameter ω affects the behavior of the theory, and it can influence the stability of wormholes.Kaluza-Klein (KK) theory [38] extends the usual four-dimensional spacetime of GR to include one or more extra dimensions.These extra dimensions are compactified or 'curled up' so small that they are not perceptible on macroscopic scales.The total spacetime is a product of the usual four-dimensional spacetime and the compactified extra dimensions.
In KK theory, the metric tensor describes the geometry of the higher-dimensional spacetime.It has components corresponding to both the usual four dimensions denoted by µ and ν and the extra dimensions denoted by a and b.The metric tensor can be decomposed into a four-dimensional part and an extra-dimensional part.The extra-dimensional part manifests as an additional vector or scalar fields, and these fields are typically associated with electromagnetic interactions.This process of decomposing the higher dimensional metric tensor field is called dimensional reduction.
The action of KK theory after dimensional reduction is where d 4 x is the four-dimensional spacetime volume element, R (4) is the four-dimensional scalar curvature, F ab is the electromagnetic field strength tensor, and L m is the usual Lagrangian density for the matter field(s).
There are currently three sets of geometric theories of gravity.The first one is the general theory of relativity based on curvature.The second is teleparallel gravity based on torsion (T), as we have seen in f (T ) theory.The third set of geometric theories is based on nonmetricity (Q), as described here.The origin of f (Q) theory is from 'symmetric teleparallel gravity', which is based on the nonmetricity scalar Q. Nonmetricity is defined as the covariant derivative of the metric tensor g µν , i.e., Q αµν ≡ ∇ α g µν .It vanishes in the case of Riemannian geometry, and it can be used to study non-Riemannian spacetimes.f (Q) gravity has inspired research in blackholes, wormholes, and cosmology.In cosmology, f (Q) models can be used to explain phenomena related to both early and late time cosmology [39], without dark energy, dark matter, or the inflation field.The action of f (Q) gravity is given by [40] Here, as usual, g is the determinant of the metric g µν , and L m is the matter Lagrangian density.f (Q) is an arbitrary function of the nonmetricity scalar Q [41] given by where Q α ≡ Q µ αµ and Qα ≡ Q µα µ are two independent traces of the nonmetricity tensor Q αµν , and these are obtained by contracting the nonvanishing tensor Q αµν ≡ ∇ α g µν .
In f (Q, T ) theories [42], the nonmetricity is coupled minimally to the trace of the matter energy-momentum tensor.The coupling between Q and T leads to nonconservation of the energy-momentum tensor, which has important physical implications such as changes to the thermodynamics of the Universe [43], the nongeodesic motion of test particles, and the appearance of an additional force.
In f (R, L m ) theory [44], f is a function of both the Ricci scalar and the matter Lagrangian.It is possible to have both an additive function such as or an exponential function such as where Λ > 0 is an arbitrary constant.This function becomes in the limit The observed late time acceleration of the Universe can be described by f (R, L m ) gravity [45,46].In [47], the MT wormhole solution was studied in f (R, L m ) gravity assuming three different types of EOSs, namely linear barotropic, anisotropic, and isotropic EOSs.The wormhole solutions obey the flare-out condition for both the barotropic and the anisotropic cases.For the isotropic case, the shape function does not follow the flatness condition.The NEC is violated in the vicinity of the throat.The amount of exotic matter required near the wormhole throat is minimized in f (R, L m ) gravity as compared to GR.
In the next two subsections, we take a deeper dive into f (R) and Lovelock theories of gravity, which will set the required background to study wormholes in these theories.In Section 3, we provide a synopsis of the study of wormholes in each of these f (R) theories of gravity.
In f (R) theories, there are two formalisms [20,49] to derive the Einstein field equations from the action.They are (i) the metric formalism, in which a matter term S m (g µν , ψ) is added to the action, where ψ represents the matter field(s).The action is then varied with respect to the metric by not treating the connections Γ µ αβ independently to obtain the field equations.Then, there is (ii) the Palatini formalism, in which an independent variation is done with respect to the metric and the connection.The action is the same, but the curvature tensors and Ricci scalar are constructed with this independent connection.
Here, we will follow [21,48] and show the derivation of the Einstein field equations in much more detail using the metric formalism.The connection coefficient Γ α βγ and the components of the Riemann curvature tensor R α βγδ are calculated using the standard equations and where the comma denotes partial derivatives.Before we vary the action, we first vary each of the quantities in the action.The variation of the determinant is The Ricci Scalar is R = g µν R µν , and the variation of R with respect to g µν is where Γ ρ µν is the Christoffel symbol representing the Levi-Civita connection, and ∇ µ is a covariant derivative.Now, by substituting (21) in (20), we obtain where By varying the action (15), we obtain Substituting δf (R) and δ √ −g from ( 24) and ( 18) into (25), we obtain Now, we substitute δR from ( 22) into ( 26) to obtain After integrating by parts and factoring out δg µν , we obtain Finally, by requiring that the action remain invariant with the variation of the metric, we obtain the field equations for f (R) modified gravity theory: where and L m is the Lagrangian for matter.

Lovelock Gravity Theory
Introduced in 1971, Lovelock's theory of gravity [50] is considered to be the most generalized extension to the theory of gravitation in D dimensions, because it satisfies the requirements of GR that the field equations be covariant and not include more than the second-order derivatives of the metric.
We consider a generally covariant theory of gravity in D dimensions.The Lagrangian is a functional of independent variables (g µν , R µ νρσ ), and the action is where S m is the action of the matter.Note that, depending on the purpose, different pairs of independent variables might be chosen to describe the system, such as (g µν , R µνρσ ) or (g µν , R µν ρσ ).However, the pair of (g µν , R µ νρσ ) is the most appropriate one when deriving the field equations.For later convenience, we define the following: where, by definition, P µν is symmetric, and the entropy tensor P µνρσ follows the symmetries We further construct the generalized form of the Ricci tensor as where it can be shown that [51,52] To obtain equations of motion, one needs to vary the action under the variation of the metric g µν → g µν + δg µν that leads to the following variations By varying the action (31), plugging in (36), and doing the partial integration, we obtain where we have used the symmetry properties of P µνρσ as well.Finally, the field equations of motion read as E µν = (1/2)T µν , where and T µν is the energy-momentum tensor of matter.Moreover, P µν is replaced by the generalized Ricci tensor R µν via (35).
Given the fact that P µνρσ already contains second derivatives of the metric by definition, one realizes that the last term in the field equations (38) will involve up to fourth derivatives of g µν .To ensure that we do not include derivatives higher than the second order, we impose the following condition: Even though one might think of imposing ∇ λ ∇ γ P µνρσ = 0 as a more general condition, it turns out that this condition does not make any difference.Imposing (39) reduces the field equations to which are nonlinear in second derivatives of the metric.In fact, imposing the linearity condition leads to the Einstein-Hilbert action.Here, one must be careful, since the form of the field equations is analogous to Einstein's field equations.However, R µν here is the generalized Ricci tensor defined in (34).
The next step is finding the form of the generalized Lagrangian that leads to field equations with no derivatives higher than the second order.This task reduces to finding those scalar functions of the metric and the Riemann tensor that satisfy (39).It can be shown that, at each order, such functions are uniquely constructed as which determines the order m Lagrangian in D dimensions.It contains m factors of the Riemann tensor, and is a completely antisymmetric determinant tensor defined as The entire Lanczos-Lovelock Lagrangian is then given by a sum over various orders of The expansion coefficients c m are initially arbitrary.The detailed explanation of the construction procedure of Lanczos-Lovelock (LL) Lagrangian and the proof of its uniqueness are out of the scope for this review.So, we refer the interested reader to the original papers by Lanczos [53,54] and Lovelock [50].
Note that the order m Lagrangian in (41) has been written in terms of R µν ρσ .This change of variable enables us to use the following identity that simplifies our guesses for a generalized form of the Lagrangian, where This identity implies that the Lagrangian must be independent of the metric tensor, and hence, indices of the Riemann tensor R µν ρσ must be contracted only by the Kronecker deltas in the form that we have in (41).
The zeroth order Lagrangian has m = 0, which simply leads to a constant term (e.g., cosmological constant).In the first order, we obtain So, the first order of the LL Lagrangian gives back the Einstein-Hilbert Lagrangian.As expected, the generalized equations of motion (40) simply reduce to the Einstein fields equations for L (0) + L (1) , i.e., where Λ is the cosmological constant, and G µν = R µν − 1 2 Rg µν is the well-known Einstein tensor.
To study the wormhole solutions of generalized gravity discussed in [18,55], we go up to the third order in the LL Lagrangian: where the next two orders are determined as and α m are the Lovelock coefficients.The second order term ( 48) is known as the Gauss-Bonnet Lagrangian.
The generalized field equation (40) corresponding to various orders of the Lagrangian can be written separately as well, and T µν is the energy-momentum tensor.
The field equations corresponding to the second and third order Lagrangian are respectively determined by R µν which are R (2)  µν ≡ and 3 Fundamentals of Morris-Thorne (MT) Wormhole In Section 3.1, we give a brief background of how the stability of a wormhole (usually an MT wormhole) is studied in GR.In Section 3.2, we summarize a general methodology to study wormholes in f (R) gravity theories based on the various papers we reviewed on this subject.This background will be useful to follow Section 4, where we review these calculations in greater detail for general f (R) theories and Lovelock gravity theory.

MT Wormhole Stabilization in GR
In [9], they explain in great detail why blackholes and Schwarzschild wormholes are not traversable.However, the MT wormhole is designed to be made traversable if it has the following properties that ensure wormhole stability for traversability.The MT metric allows for the realization of faster-than-light interstellar space travel that does not violate the special relativistic light speed limit.The metric should be spherically symmetric and static (time independent).The following metric has these properties [9]: The solution must obey the Einstein field equations as does (53).The solution must have a throat that connects two asymptotically flat regions of spacetime.The spatial geometry must have a wormhole shape consistent with the well-known Flamm diagram for a spherically symmetric throat.This puts the following constraints on the shape function b(r) and redshift function Φ(r): • The throat is at minimum of r, which is specified as r 0 .
• Proper radial distance is defined by and should be real and finite for r > r 0 .• As l − → ±∞ (asymptotically flat regions of the Universe), b/r − → 0, and so r ∼ |l|.• There should be no horizons, since it will prevent two-way travel through the wormhole.
There are no singularities.This implies that Φ is finite, continuous, and differentiable everywhere, as well as the fact that τ measuring proper time in asymptotically flat regions implies Φ − → 0 as l → ± ∞.
and so rb ′ < b.That is, the throat of the wormhole must expand outward from the central point.The throat of the wormhole must open up as one travels through it.
The tidal gravitational forces experienced by a traveler must be ≤g, where g is the acceleration due to Earth's gravity.This condition is not a rigid requirement.
The procedure to check the stability of this traversable wormhole in GR involves the following steps: Compute Curvature tensors: Here, we give as briefly as possible the method to calculate the curvature tensors.First, using the MT metric written in the form with x 0 = t, x 1 = r, x 2 = θ, and x 3 = ϕ, the connection coefficient Γ α βγ and the components of the Riemann curvature tensor R α βγδ are calculated using the standard Equations ( 16) and (17).By applying these equations to the metric, we can obtain the 24 nonzero components of the Riemann tensor, as shown in Equation ( 5) of [9].These were obtained using the basis vectors (e t , e r , e θ , e ϕ ).
To further simplify the calculations, we can switch to the following orthonormal basis vectors: In this basis, the metric coefficients become the same as those of flat (Minkowski) spacetime, and the 24 nonzero components of the Riemann tensor take a much simpler form, as seen in Equation ( 8) of [9].
Contraction: Next, by contracting the Riemann tensor, we obtain the Ricci tensor, and again, by contracting the Ricci tensor, we obtain the Ricci scalar R.

Compute Einstein tensor:
We finally obtain the Einstein tensor G µν from the metric, Ricci tensor, and the Ricci scalar.This forms the left-hand side of the Einstein field equations, as given in (46).This yields the nonzero components of G µν in terms of the shape function and redshift function, namely and Compute stress-energy tensor: The next step involves computing the stress-energy tensor (the right-hand side of the Einstein field equations).Based on Birkhoff's theorem, which states that "any spherically symmetric solution of the vacuum field equations must be static and asymptotically flat", the exterior solution must be given by the Schwarzschild metric (with certain modifications), which is a spherical wormhole.Therefore, we cannot have a vacuum solution for a traversable wormhole, which implies that our wormhole must be threaded by matter with a nonzero stress-energy tensor.Based on Einstein's field equations, T µν must have the same algebraic structure as G µν in the orthonormal basis that we have chosen.Similar to G µν , only the four components T tt , T rr , T θθ , and T ϕϕ are nonzero.Based on a remote static observer's measurement, each of the components have a simple physical interpretation as follows: is the total rest energy density that the static observer measures, is the radial tension measured per unit area, and is the tangential pressure measured per unit area in a direction orthogonal to the radial tension τ (r).This give us the final stress-energy tensor: Engineer the traversable wormhole: In the next step, we need to 'engineer' the traversable wormhole to obtain the properties enumerated earlier in this section.This can be done by controlling the shape function b(r) and the redshift function Φ(r) based on a suitable T µν that we require.We substitute the G µν and T µν found in the previous steps into (60) to solve for ρ(r), τ (r), and p(r) in terms of b(r) and Φ(r) to obtain Our strategy for the stabilization of the MT wormhole in Einstein gravity will involve tailoring b(r) and Φ(r) to build a wormhole with the required properties.Our choice of b(r) gives us ρ(r).Our choice of b(r) and Φ(r) will give us the tangential pressure τ (r).We next find p(r) using ρ(r), τ (r), and Φ(r).We then ensure that the averaged NEC and WEC are violated, which means that the source of matter for traversable wormholes must be exotic.
With this process for analyzing the stability of the MT wormhole in GR, many authors in our review have used similar techniques to analyze the stability of the Morris-Thorne wormhole in f (R) gravity theories.We will hone into the nuances of these analysis in the following sections.

A General Methodology for MT Wormhole Stabilization in f(R) Gravity Theories
Based on Section 4, we briefly summarize the method that can be used to stabilize the MT wormhole in f (R) gravity theories.First, we require that matter threading the wormhole satisfy the NEC and WEC.We include the required flare-out condition for traversable wormholes.We delegate the required energy condition violations and sustenance of the nonstandard wormhole geometry to the higher-order curvature terms, including the derivative terms (T (c) µν ).Consider a redshift function (Φ = c, Φ ′ = 0), where c is a constant, which simplifies calculations and still provides physically relevant solutions.This condition defines an ultrastatic traversable wormhole and is a very narrow and specific subclass of solutions called zero tidal force (ZTF) solutions.We specify b(r).Some examples used in [48] √ r 0 r, and b(r) = r 0 + γ 2 r 0 (1 − r 0 r ), with 0 < r < 1.For each shape function b(r), we assume an equation of state such as p r = p r (ρ) or p t = p t (ρ).We find F (r) (see Section 4.1) from the modified gravitational field equations, with the Ricci curvature scalar R(r) obtained from the MT metric.Finally, we obtain the exact f (R) that we need from the trace Equation (68).
4 Wormholes in f(R) Gravity Theories

Wormholes in f(R) Gravity Theory
In this section, we mainly follow the references [21,22,31,48] to analyze the stability of an MT wormhole in f (R) gravity theories.In Section 2.2, we obtained (29) for f (R) gravity theory.Contracting this field equation, we obtain the trace of T µν , where F = f ′ (R) is used for convenience, R is the Ricci scalar, and T = T µ µ is the trace of the stress-energy tensor T µν .By substituting ( 68) into ( 29) and the rearranging terms, we obtain an updated field equation: We will call (69) the effective field equation for f (R) gravity theory, where µν is the curvature stress-energy tensor for a higher-order curvature given by and Here, we write the energy-momentum tensor in terms of an anisotropic distribution of matter as follows: such that u µ u µ = −1 and x µ x µ = −1, where u µ is a four-velocity vector, r is a unit spacelike vector in the radial direction, ρ(r) is the rest energy density, p r (r) is the radial tension, and p t (r) is the tangential pressure orthogonal to x µ .This gives us and By substituting this value of T in the trace form of the field equations for f (R) gravity (68), we obtain F R − 2f + 3□F = (−ρ + 2p r + p t ).( 75) We now use the MT wormhole metric (53) just as we did in GR, where Φ(r) is the redshift function, and b(r) is the shape function as before.The radial coordinate r decreases from ∞ to a minimum value r 0 at the wormhole throat.At the throat, b(r 0 ) = r 0 , and it then increases from r 0 back to ∞.We recall the flare-out condition, which is important for traversability: b At the throat, b(r 0 ) = r = r 0 b(r 0 ) > b ′ (r 0 )r or r > b ′ (r 0 )r, b ′ (r 0 ) < 1 is the flareout condition that we need for a traversable wormhole solution.For the wormhole to be traversable, no horizon should be present.Horizons are surfaces with e 2Φ − → 0, so we want Φ(r) to be finite everywhere.This implies that we can use the constant redshift functions Φ(r), Φ = c, and Φ ′ = 0.This will help to simplify calculations and avoid fourthorder differential equations.The effect of using a variable redshift function is discussed in Section 4.3.The following steps can be used to calculate ρ, p r , and p t in terms of the shape and redshift functions.Substitute T µν in the effective field Equation (69).Use the MT metric (53) Note that F ′ = d dR ( df (R) dR ), and dR .Thus, we obtain the following simplified equations from the effective field Equation ( 69): and These are the required simplified equations for matter threading the wormhole as a function of b(r) and F (r).We will use them later in this paper.We can now determine the matter content by choosing an appropriate shape function and a specific form of F (r).At this point, let us recap the strategy for stabilizing a wormhole in f (R) gravity.We choose a shape function b(r).Then, we specify an equation of state such as p r = p r (ρ) or p t = p t (ρ).This will let us compute F (r) from the effective field equation of Equation (69).We can also find the Ricci scalar R(r) from the MT wormhole metric (53).Then, we obtain T = T µ µ as a function of r.Finally, we compute the function f (R) from the trace of the field Equation (68).

Violation of Energy Conditions
Just as the motion of a single particle is governed by the geodesic equation, the equation of motion of a family of particles, also known as a congruence, is governed by the Raychaudhuri equation [56].From this equation, the following focusing condition in terms of the Ricci tensor arises: where k µ is a null vector.If this condition is satisfied, then the geodesic congruences focus into a finite value of the parameter for labeling points on the geodesics.In GR, this is written as the NEC, in terms of the stress-energy tensor, as follows: In modified gravity, in particular in f (R) theories, we could first assume that T (m) µν satisfies this energy condition, and violation of the energy condition can be assumed to come from the higher-order curvature terms T µν .Note that this condition applied to f (R) gravity theory does not mean geodesics are focused, as required by the Raychaudhuri equation.In terms of the radial null vector, violation of the NEC requires and takes the form where ρ eff + p eff r < 0.
Using the gravitational field equation of Equation ( 69), we obtain By applying the flare-out condition, At the throat, with r = r 0 , At the throat, this gives us We have been highlighting that matter threading the wormhole should obey the WEC as well, i.e., ρ ≥ 0, and ρ + p r ≥ 0. Applying these WECs to the simplified Equations ( 79)-( 81) for ρ, p r , and p t , we obtain The four inequalities above ( 86)-( 89) must be obeyed by the function f (R) for traversable wormholes in f (R) theories, thus taking into account that matter threading the wormhole satisfies the NEC and WEC, and the flare-out conditions are satisfied at the wormhole throat.The task of maintaining the wormhole geometry (violating the NEC) is delegated to the higher-order curvature terms T (c) µν .

Effect of Using a Variable Redshift Function
The effect of using a variable redshift function for wormholes has been studied by [57] in κ(R, T ) gravity with Φ(r) = α r , where α is a constant.The authors of [57] obtained solutions that require exotic matter.Wormholes in f (R) gravity were studied by [58][59][60][61][62] with a variable redshift function Φ(r) = ln r 0 r + 1 , and the authors in [58][59][60][61][62] found solutions that require nonexotic matter.The authors in [63] studied the effect of a variable redshift function Φ(r) = − α r with α ≥ 0, where α is a constant, and [64] used Φ(r) = α r .[21] in particular studied wormholes in f(R) gravity with a variable redshift function Φ(r) = 1 r in great detail.Here, we discuss the solution obtained by them.[21] used the function f (R) = R − µR c ( R Rc ) p , where µ, R c , and p are constants, with µ > 0, R c > 0, and 0 < p < 1. [21] also used the shape function b(r) = r e r−r 0 .Their motivation for selecting the above f (R) function is related to obtaining viable dark energy models, and it is also related to finding explanations for the accelerated expansion of the universe.Based on past studies, this f (R) function seems to be a good candidate for exploring wormholes in f (R) gravity with both a constant and a variable redshift function.
With the redshift function Φ = 1 r , they found ρ > 0 for r ≥ 1.2, ρ + p r > 0 for r ≥ 1.2, ρ + p t > 0 for r > 1.8, ρ + p r + 2p t < 0 for r > 1.2, ρ − |p r | > 0 for r ≥ 1.2, and ρ − |p t | > 0 for r > 1.8, and the NEC, WEC, and DEC were satisfied for r > 1.8 with the variable redshift function.So in conclusion, for wormholes with r 0 > 1.8-where r 0 is the radius of the throat-and with variable redshift function 1 r , they obtained wormhole geometries free of exotic matter.They also did similar calculations in GR with the variable redshift function 1 r and found that there was no solution in GR without exotic matter for any value of r.

The Stability Condition and the Speed of Sound
In [65], the authors analyzed the stability of a wormhole in addition to traversability in f (R) gravity.The authors in [65] used a power law f (R) = f 0 R 1+ϵ , which can be written in the form f (R) = R + ϵRln(R) + O(ϵ 2 ), for ϵ ≪ 1, to consider small deviations from Einstein gravity.Such an approach was used in the study of compact objects such as neutron stars and blackholes [66,67], where departure from GR can be useful to explain observations.For the analysis, the authors used the MT wormhole metric with a redshift function 2Φ(r) = r 0 r and two different shape functions b(r) r = r 0 r 1+β , where β is a real number, β+1 > 0, and b(r) r = r 0 1+αr -with α ≡ r 0 −1 r 0 -to ensure that the wormhole is not singular at r 0 .They obtained the lengthy Equations ( 13)-( 16) in [65] for ρ(r), p r (r), p t (r), and the average pressure p(r) = 1  3 [p r (r)+2p t (r)], respectively.Next, they used ideas from fluid dynamics [68][69][70] and defined the adiabatic speed of sound: To obtain a stability condition, we need a vanishing speed of sound at the throat: By inserting Equations ( 13) and ( 14) in [65] for ρ(r) and p r (r) into Equation ( 91) above, they obtained the required stability conditions, namely Equations ( 24) and ( 25) in [65], for b(r) r = ( r 0 r ) β+1 and for b(r) r = r 0 1+αr , respectively.Their primary finding was that for a specific range of values of ϵ in the assumed f (R), and using the condition for vanishing speed of sound, stable wormhole solutions exist in the presence of nonexotic matter (a perfect fluid).

Wormholes in Lovelock Gravity Theory
We will now discuss how the modifications in Lovelock gravity help with the stabilization of a wormhole.In this section, we mainly follow [18,55].The second-and third-order curvature terms in Lovelock gravity are important to be considered, especially in the case of wormholes with smaller throat diameters, where the curvature is very high.First, the MT wormhole metric is modified for n dimensions as follows: where Φ(r) is the redshift function, and b(r) is the shape function.The WEC requires that matter threading the wormhole have positive energy density ρ(r) and positive ρ(r) − τ (r) , where τ (r) is the radial tension, as well as positive ρ(r) + p(r), where p(r) is the tangential pressure orthogonal to the radius.If the WEC ρ = T µν u µ u ν ≥ 0 is satisfied everywhere, then the wormhole can be constructed with normal matter without the need for exotic matter.Here, u µ is a timelike velocity of the observer.Note that the stress tensor T µν used in the WEC is calculated using (50), which includes all the Lovelock gravity modifications to Einstein's field equations.
The NEC is The authors in [18] used the orthonormal basis set and to calculate the components of the energy-momentum tensor, which gives us ρ, τ , and p.
That is, and Then, by setting the Φ(r) = c = 0, where c is a constant, (ρ − τ ) and (ρ + p) are calculated as follows: and where α 2 and α 3 are the second-order and third-order Lovelock coefficients.In order to obtain positive ρ and (ρ + p) values, the rest of the analysis was done by studying three types of shape functions b(r): the power law, the logarithmic, and the hyperbolic shape functions.The power law shape function is given by b with positive m.The functions ρ and (ρ + p) for the power law above are positive for r > r 0 , provided that r 0 > r c , and r c is the largest positive real root of certain equations given in [18].For the logarithmic shape function, we have where ρ and (ρ + p) are positive for r > r 0 , with the condition r 0 ≥ r c , and r c is the largest real root of a second set of equations in [18].For the hyperbolic shape function, we have where ρ and (ρ + p) are positive if r 0 > r c , with the condition that r c is the largest root of a third set of equations in [18].
To ensure that (ρ − τ ) is positive, we need where α 2 and α 3 are the Lovelock coefficients.For certain combinations of the Lovelock coefficients contributing to the throat radius, with a negative value for either α 2 or α 3 , the above condition is satisfied only in the vicinity of the throat for all three types of shape functions discussed above.The throat radius must be in the range r − < r 0 < r + with and The following points summarize why the higher-order curvature terms are needed to stabilize a wormhole, as well as how the GB (second-order) and Lovelock (third-order) curvature terms actually help with using normal matter as opposed to exotic matter in the vicinity of the throat.Higher-order curvature corrections are used in wormholes with smaller throat radii, since the curvature near the throat is very large for such wormholes.The matter near the throat can be normal for the region r 0 ≤ r ≤ r max , where r max depends on the Lovelock coefficients and the shape function.We obtain a larger radius (more region) with normal matter in the case of third-order Lovelock gravity with z negative coupling constant α 3 as compared to the second-order Gauss-Bonnet wormholes.There is a lower limit on the throat radius r 0 imposed by the positivity of ρ and (ρ + p) in the case of Lovelock gravity but not in the case of Einstein gravity.The lower limit depends on the Lovelock coefficients, the dimensionality of the spacetime, and the shape function.

Wormholes in Einstein Dilaton Gauss-Bonnet (EdGB) Gravity
The authors in [71] reported stable wormhole solutions in EdGB gravity.Other authors reported in [72] that the same wormhole solution in EdGB gravity was indeed unstable.Since these two papers are excellent examples of stability analysis, in this section, we will review their key findings and methods.
In order to avoid using the required exotic matter necessary to support traversable wormholes, modified gravity theories were used.One such modified gravity theory is EdGB gravity.In this theory, the low-energy heterotic-string-theory-based effective action in four dimensions is given by [73,74] where, in addition to the scalar curvature term R, we have a quadratic GB curvature term L (2) given by ( 48) and a scalar field term ϕ known as the dilaton field, with a coupling constant γ. α is a positive numerical constant given in terms of the Regge slope parameter.
In [71], the authors analyzed the stability of a spherically symmetric wormhole solution in EdGB gravity.With a coordinate transformation r 2 = l 2 + r 2 0 , where r 0 is the radius at the throat, the spherically symmetric wormhole solution becomes pathology-free and is given by In terms of the new coordinates, the expansion at the throat (l = 0) gives where f i , ν i , and ϕ i are constant coefficients, and all curvature invariants, including the GB term, remain finite for l → 0. Next, we analyze the stability of the wormhole in terms of radial perturbations.The metric and dilaton functions depend on both l and t.We then decompose the functions into unperturbed and perturbed parts.The perturbations are σ is defined in [71] as an unspecified eigenvalue that appears in their Ordinary Differential Equations (ODEs).Substituting these perturbations into the Einstein and Dilaton equations, as well as linearizing in ϵ, we obtain a system of linear ODEs for the functions δν(l), δf (l), and δϕ(l).Rearranging the ODEs, the authors arrived at a coupled second-order equation in δϕ as follows: (δϕ where q 0 , q 1 , and q σ depend on the unperturbed solution.To be able to normalize, δϕ → 0 as l → ∞.At l = 0, q σ is bounded, and q 0 and q 1 diverge as 1 l .To avoid this singularity, a transformation is made as follows: where F (l) satisfies the equation This implies that where 4 + q 0 is bounded at l = 0.Then, ψ → 0 as l → ∞ and for l = 0. Solving these ODEs for several values of α r 2 0 and f 0 , a solution exists only for certain values of the eigenvalues σ 2 .The negative modes obtained for families of wormhole solutions indicate that the wormhole is unstable.In addition, there are positive modes in a region of parameter space where the EdGB wormhole solutions are linearly stable under radial perturbations.
In [72], the same wormhole was proved to be unstable for the following reason.In [71], the perturbation function is where the factor A(l) diverges at the throat as 1 l .Vanishing boundary conditions were imposed at the throat to obtain finite perturbation at the throat.This disconnects the regions of space on both sides of the throat for purely radial modes of perturbations.δϕ is not a gauge invariant quantity.For general spherically symmetric perturbations of the wormhole, the gauge invariant quantity χ(t, l) ∝ γδϕ(t, l) − r 2 l ϕ ′ (l)δr(t, l) satisfies the wavelike equation where V ef f is a finite effective potential.When χ is finite at the throat, δϕ does not diverge unless the throat size is fixed by chosing δr = 0. Therefore, the wormhole is unstable with perturbation at the throat radius.The behavior of the dilaton field ϕ at the throat is a pure artifact of the chosen gauge and therefore can be safely ignored.
5 Other Wormhole Studies

Casimir Wormholes
Traversable wormhole solutions to Einstein's field equations require the existence of exotic matter (matter violating the weak, null, and dominant energy conditions).The most common type of exotic matter proposed is that which has negative energy density ρ < 0. Classically, negative energy density is thought to be impossible.However, in quantum field theory, negative energy densities can and do occur.The most famous example of this is the Casimir vacuum energy.The usual derivation of it follows quantizing the electromagnetic field between two neutral parallel plates and renormalizing the energy.This results in an energy density that is negative or less than the energy of the undisturbed vacuum: where a is the separation between the plates.As a result, there has been much speculation and study undertaken into the possibility of using the Casimir vacuum energy to construct traversable wormholes.This section will serve as an introduction to the basic concepts and results of attempting to use the Casimir vacuum energy to construct a traversable wormhole.Morris, Thorne, and Yurtsever [9,75] were the first to publish the idea of utilizing the Casimir effect to support traversable wormholes.They noted that there were two distinct possibilities for utilizing the Casimir effect.The first is to consider the Casimir vacuum energy as is, i.e., with a representing the fixed plate separation.The second is to promote the plate separation to the radial coordinate r.
The case of fixed plate separation was analyzed by Garattini [76], in which he studied the connection between the Casimir energy and absurdly benign traversable wormholes.He utilized the semiclassical Einstein field equations: where ⟨T µν ⟩ is the renormalized quantum expectation value of the stress-energy tensor, and G µν is the classical Einstein tensor.The justification for this equation was first given by Hawking [77] in his derivation for particle creation by black holes, in which he argued that quantum gravity contributions to the Einstein field equations can be assumed negligible above the Planck length.The same assumption is made here.Upon solving the Einstein field equations using the metric for a static spherically symmetric traversable wormhole, defined by Morris and Thorne, the shape function is found to be where r 0 is the throat radius.It is clear from this expression that this spacetime is not asymptotically flat, and in fact, it is asymptotically de Sitter.As stated in [76], the Casimir vacuum energy can then be viewed as a cosmological constant.However, in close proximity to the throat, the shape function can be approximated as where it has been assumed that r is very close to r 0 , and l p is the Planck length.This shape function has the same form as the shape function for "absurdly benign traversable wormholes" defined by Morris and Thorne [9], namely where d is given by (128), and exotic matter is confined in the region (r 0 ≤ r ≤ r 0 + d).
In order to make the connection between the Casimir shape function and this shape function, the zero tidal force (ZTF) condition must be imposed, where Φ(r) = 0.The leading order for this shape function close to the throat is which gives us the required connection between these two shape functions: The absurdly benign traversable wormhole is defined so that the exotic matter is confined close to the throat (d is small).In order to fit this requirement, there are two possibilities: the first is r 0 > 10 34 and the second is a ∼ 10 −15 m for a wormhole with r 0 ∼ 10 10 m.Neither of these conditions are physical.It is noted that the identification made for the relationship between d and a is a physically meaningless assumption.Therefore, Garattini attempted to refine and reduce the throat radius by examining various examples of inhomogeneous equations of state, and he managed to establish a throat radius of r 0 ∼ 10 17 m for a plate separation of 1 nm.However, these results still present a nonphysical result, and they point to the likelihood that the semiclassical Einstein field equations may be ill suited for the description of Casimir traversable wormholes.
Promoting the plate separation to radial coordinates was discussed in [78] as follows.Once again, semiclassical quantum gravity was used but now with the Casimir vacuum energy density given by where r is the radial coordinate.This gives the following shape function Here, . First, observe that this cannot be transformed into an absurdly benign traversable wormhole.In fact, attempting to impose the ZTF condition produces disastrous results.Instead, we use the Einstein field equation and the shape function to solve the redshift function.We also use the EOS p r = ωρ.In solving the redshift function, it is necessary to impose the condition ωr 2  1 − r 2 0 = 0.This condition is required for the existence of the wormhole, since if ωr 2  1 − r 2 0 > 0, we obtain a black hole.With this condition, we obtain the following shape function and redshift function, respectively, and These result in a wormhole of Planckian size.
We now return to the question of using the semiclassical Einstein field equations.As stated above, Hawking originally justified the use of this approximation in the derivation of Hawking radiation.In this scenario, it makes sense that quantum gravitational effects can be ignored in the Einstein tensor, since Hawking never considered Planckian-sized black holes, as the black holes would be expected to have evaporated before reaching this size.So, any such effects could be assumed to be hidden by the formation of the event horizon.
When considering the use of the Casimir effect in constructing traversable wormholes, we find that we might not be able to ignore quantum gravity.First, as we have just seen above, promoting the plate separation to a radial coordinate results in a Planckian-sized wormhole negating Hawking's assumption of not having to deal with any system of this size.For the fixed plate separation, we encountered solutions that produce wormholes with throat radii on the order of the Sun's radius for plate separation on the order of nanometers and smaller.This makes no physical sense, as the Casimir vacuum energy is confined within the plates and would have no way of interacting with a body of such size.Lastly, we do not have any event horizon by definition and therefore cannot assume that quantum gravity effects are hidden away.So far, no one has made an attempt to incorporate any such quantum gravity effects into the analysis of Casimir wormholes, and it is therefore left as the subject for future publications.However, other studies have been done on the subject of Casimir traversable wormholes.For example, Garattini [79] studied the effects of an electric charge on a Casimir wormhole.He found that the throat radius becomes directly dependent on the strength of the electric charge, meaning that there exists the possibility of creating a throat larger than the Planck length.However, it should be noted that, in his analysis, it was found that there exists a radius at which the energy density becomes positive.It is not discussed whether this negates the existence of the wormhole.
Casimir wormholes have also been studied in general D dimensions [80].It has been shown that, even in these scenarios, the radius is found to be on the order of the Planck length or a nonphysical size.It is worth reiterating that these studies have been constructed primarily using the semiclassical Einstein field equations and therefore do not include any possible effects of quantum gravity, which is a fact that is most likely the source of many of the nonphysical results.

Thin-Shell Wormholes in Modified Gravity
There is a set of wormhole solutions in modified gravity theories that use special techniques to satisfy the energy conditions and thus remove the need for exotic matter.We summarize the analysis of these wormholes here.In [81], the authors derived asymptotically flat traversable wormhole solutions that satisfy the NEC in a quadratic form of the gravitational hybrid metric-Palatani gravity theory.Their solution has an interior part with a nonexotic perfect fluid near the throat, an exterior Schwarzschild solution, and a double gravitational layer thin shell at the junction hypersurface Σ between the interior and exterior solution.
The generalized hybrid metric-Palatani gravity theory has the action Here, κ 2 ≡ 8πG c 4 , Ω is the spacetime manifold in which a set of coordinates x a are defined, g is the determinant of the metric g ab , and f (R, R) is a well-behaved function of R and R. R ≡ g ab R ab is the Ricci scalar, R ≡ g ab R ab is the Palatini scalar curvature, and R ab is the Palatini Ricci tensor written in terms of an independent connection Γ c ab .L m is the matter Lagrangian density minimally coupled to the metric g ab .The junction conditions at the separation hypersurface Σ were defined to be [R] = 0, [R] = 0, (137) and where h αβ = e a α e b β g ab is the induced metric at the hypersurface Σ, K αβ = e a α e b β ∇ a n b is the extrinsic curvature, K = K α α is the trace of the extrinsic curvature, and S β α is the stressenergy tensor of the thin shell arising at the hypersurface Σ.It was shown in [81] that the general set of junction conditions shown above can be simplified for particular forms of f (R, R).In particular, they selected an f (R, R) that is quadratic in R and linear in R. For this f (R, R), the junction conditions [R] = 0 and [R] = 0 are not mandatory conditions anymore.This gives rise to additional terms in the stress-energy tensor S ab of the thin shell, which then forms a double layer thin-shell distribution at the junction hypersurface Σ.
Unlike the general case [82], in which the solutions obtained using scalar-tensor representation are scarce, the simplified set of junction conditions in [81] implies that it is possible to obtain numerous solutions for a wide variety of metrics and actions.There is also the advantage that asymptotic flatness can be preserved in this case, unlike the general case where the NEC can be guaranteed only for the asymptotically flat AdS spacetime.
In [83], traversable wormholes with double layer thin shells in quadratic gravity were analyzed, where f (R) ≡ R + αR 2 , in which R is the Ricci scalar.The NEC is satisfied in this case at the throat, as well as the whole wormhole interior.However, the NEC is not satisfied for the double layer stress-energy distribution component at the thin shell.
In [84], traversable wormhole solutions were analyzed in f (R, T ) gravity theory for a linear model, i.e., f (R, T ) = R + γT , where T = T ab T ab , and T ab is the energy-momentum tensor.It was shown that there are a large set of wormhole solutions in which the matter field satisfies all the energy conditions (NEC, WEC, SEC, and DEC).Since the field equations are quadratic in the matter quantities ρ, p r , and p t , as well as complex, they used an analytical recursive algorithm to extract the nonexotic wormhole solutions.The solutions obtained were not naturally localized, so the junction conditions were derived for this theory.It was proven that a matching between two spacetimes must always be smooth and does not allow thin shells at the boundary.Traversable localized static and spherically symmetric wormhole solutions satisfying all energy conditions were obtained by matching the interior wormhole spacetime to an exterior vacuum Schwarzschild solution.

Dark Matter Halo Wormholes
It is well known that an exotic type of matter with negative energy density can stabilize wormholes.Dark matter halos are vast, invisible regions of space that surround galaxies.They are composed of dark matter [85], a substance that does not emit, absorb, or reflect light, so they cannot be directly detected.Their presence is inferred through their gravitational effect on visible matter such as stars and gas clouds.In recent studies of galaxy formation, it was found that every galaxy forms within a dark matter halo [86].The formation and growth of galaxies over time is related to the growth of the halos in which they form.
In recent years, a few research groups have been investigating traversable wormholes supported by the dark matter halo [87][88][89][90]-both in GR and modified theories of gravity such as f(Q,T).Dark matter halo wormholes are usually studied using the Milky Way galaxy (MWG) halo profiles, pseudothermal, Navarro-Frenk-White (NFW) models I and II, and Universal Rotation Curves (URCs).A sample investigation of dark matter profiles in the Milky Way galaxy can be found in [91].In [90], they used the "Einasto dark matter density profile" to produce suitable redshift and shape functions.The NEC is violated for redshift functions Φ = C and Φ = α r , and the shape function satisfies the flare-out condition.The anisotropic dark matter content within the wormhole creates the appropriate environment to stabilize the wormhole structure by violating the NEC.The global monopole charge η [92][93][94] plays an important role in the violation of the NEC.The probability of violation of the NEC decreases for an increasing value of η, and so it is important to minimize the value of η.
In [95], wormholes supported by galactic halos have been investigated in 4D EGB gravity.This analysis was done for three different dark matter profiles, namely URC, NFW, and Scalar Field Dark Matter (SFDM).The NEC was violated at the neighborhood of the wormhole throat.The Gauss-Bonnet coefficient α has an influence on the NEC.α > 0 gives negative energy.The contribution to the violation increases with α.The SEC is also violated for each of the DM profiles.

Wormholes in Nonlocal Theories of Gravity (NLGs)
There is a class of nonlocal integral kernel theories of gravity where the inverse of the d'Alembert operator in the gravitational action is taken into account.These are called nonlocal theories of gravity (NLGs).A spherically symmetric MT wormhole solution satisfying NLG field equations was studied in [96].Linear and exponential NLG correction terms were selected due to the existence of Noether symmetry [97].It was found that the nonlocal gravity contributions allow for the stability and traversability of a wormhole without considering exotic matter.

Summary and Discussion
In this review paper, our goal is to review stable, traversable wormholes in f (R)-like gravity theories.We started with a historic development of ideas about wormholes in the Introduction section and then provideed some motivation for this paper and the need for an equation of state (EOS) in Section 1.1.We then concluded the Introduction with a discussion of the astronomical observational signatures of natural wormholes in Section 1.2.In Section 2.1, we started with a brief discussion of each modified gravity theory such as f (R), f (R, T ), Lovelock, EGB, Brans-Dicke, and KK theory, and we followed with a brief discussion of nonmetricity theories such as f (Q) and f (Q, T ).We ended Section 2.1 with a review of f (R, L m ) theory and its applications.In Sections 2.2 and 2.3, we gave a more detailed treatment of f (R) gravity theory and Lovelock gravity theory by deriving the EOM from the action for the respective theory.In Section 3, we discussed the MT wormhole in greater detail, since it is the standard traversable wormhole geometry used in most of the studies of wormhole in modified gravity theories, as seen in the literature.In Section 4, we reviewed the study of wormholes in f(R) gravity and Lovelock gravity based on the flareout condition and energy condition violation requirements.We concluded Section 4 with a detailed review of stability analysis under a perturbation using EdGB gravity as an example.In Section 5, we took the reader on a tour of other wormholes such as Casimir wormholes, thin-shell wormholes in modified gravity, dark matter halo wormholes, and wormholes in nonlocal theories of gravity.
In the examples that we discuss, it is possible to have stable wormholes without the use of exotic matter in many of the modified gravity theories.In some cases, it is shown that the amount of exotic matter that is needed can be minimized.In their analyses, the authors varied the redshift function, shape function, the type of fluid used, the equation of state, and the primary function used in the corresponding modified gravity theory.
The key takeaways for wormholes in f (R) gravity theory are as follows: In GR, violation of the NEC is required for static traversable wormholes.In the f (R) theory of gravity, modified field equations are obtained by varying the modified action with respect to the metric.Using these modified field equations, we require that matter (stress-energy tensor) threading the wormhole satisfy the NEC, and the required violation of the NEC can be enabled by the total stress-energy tensor, which includes higher-order curvature terms.The higher-order curvature terms, interpreted as a gravitational fluid, support the nonstandard wormhole geometries.A constant redshift function is assumed in many of these analyses to reduce the complexity of the calculations.It is possible to use a variable redshift function as well.Based on the review of papers related to wormholes in Lovelock gravity, we have the following key takeaways: Higher-order curvature terms become useful for the analysis of wormholes with smaller throat radius.There exists a lower limit for the throat radius in Lovelock gravity imposed by the requirement that ρ and (ρ + p) be positive.There is no such limit in Einstein's gravity.The radius of the region with normal matter is higher for wormholes in third-order Lovelock gravity, with a negative coupling constant (α 3 ), compared to wormholes in second-order Gauss-Bonnet gravity.
In Section 5, we discussed "Other Wormholes" such as Casimir, thin-shell, dark matter halo, and nonlocal gravity wormholes.In Casimir wormholes, the required negative energy density for a stable traversable wormhole is provided by the Casimir effect.Casimir wormholes are currently being studied in several modified gravity theories.Thin-shell wormholes use special techniques to satisfy the energy conditions and thus remove the need for exotic matter.For example, there are asymptotically flat traversable wormhole solutions that satisfy the NEC in the quadratic form of the gravitational hybrid metric-Palatani gravity.These solution have an interior portion with a nonexotic perfect fluid near the throat, an exterior Schwarzschild solution, and a double gravitational layer thin shell at the junction hypersurface Σ between the interior and exterior solutions.In dark matter halo wormholes, the dark matter content within the wormhole creates the appropriate environment to stablize the wormhole structure by violating the NEC.In NLG wormholes, it was found that the nonlocal gravity contributions allow for the stability and traversability of a wormhole, without considering exotic matter.