Existence of stable wormholes on a non-commutative-geometric background in modified gravity

In this paper, we discuss spherically symmetric wormhole solutions in f(R, T) modified theory of gravity by introducing well-known non-commutative geometry in terms of Gaussian and Lorentzian distributions of string theory. For analytic discussion, we consider an interesting model of f(R, T) gravity defined by f(R,T)=f1(R)+λT\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f(R,T)=f_{1}(R)+\lambda T$$\end{document}. By taking two different choices for the function f1(R)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{1}(R)$$\end{document}, that is, f1(R)=R\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{1}(R)=R$$\end{document} and f1(R)=R+αR2+γRn\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{1}(R)=R+\alpha R^{2}+\gamma R^{n}$$\end{document}, we discuss the possible existence of wormhole solutions. In the presence of non-commutative Gaussian and Lorentzian distributions, we get exact and numerical solutions for both these models. By taking appropriate values of the free parameters, we discuss different properties of these wormhole models analytically and graphically. Further, using an equilibrium condition, it is found that these solutions are stable. Also, we discuss the phenomenon of gravitational lensing for the exact wormhole model and it is found that the deflection angle diverges at the wormhole throat.


Introduction
In modern cosmology, the phenomenon of accelerated cosmic expansion and its possible causes, as confirmed by numerous astronomical probes, have become a focus of interest for the researchers [1][2][3]. In this respect, the first attempt was made by Einstein, by introducing a well-known C DM model; but in spite of its all beauty and success, this model cannot be proved to be problem free [4]. Later on, a bulk of different proposals have been presented by the researchers that can be grouped into two kinds: modia e-mails: mzubairkk@gmail.com;drmzubair@ciitlahore.edu.pk b e-mail: gmustafa3828@gmail.com c e-mail: swaheed@pmu.edu.sa d e-mail: abbasg91@yahoo.com fied matter proposals and modified curvature proposals. The tachyon model, quintessence, the Chaplygin gas and its different versions, phantom, quintom etc., are all obtained by introducing some extra terms in the matter section and hence are members of the modified matter proposal group [5][6][7][8][9][10][11][12][13][14][15]. The other idea is to modify curvature sector of Einstein's general relativity (GR) by including some extra degrees of freedom there. One of the primary alterations was the speculation of the Einstein-Hilbert Lagrangian density with an arbitrary function f (R) instead of the Ricci scalar R. This theory has been widely used in the literature [16] to examine the dark energy (DE) and its resulting speedy cosmic expansion. Moreover, f (R) theory of gravitation provides a unified picture of early stages of cosmos (inflation) as well as the late stages of accelerated cosmos. Some other well-known examples include Brans-Dicke gravity, generalized scalar-tensor theory, f (τ ) gravity, where τ is a torsion, Gauss-Bonnet gravity and its generalized forms like f (G) gravity, f (R, G) gravity, and f (τ, τ G ) theory, etc. [17][18][19][20][21][22][23][24][25][26].
Another significant modification of Einstein gravity namely f (R, T ) gravity was proposed by Harko et al. [27] almost five years ago. In this formulation, a generic function f (R, T ), representing the coupling of the Ricci scalar and energy-momentum tensor trace, replaces the Ricci scalar R for the possible modification of curvature sector. Using metric formalism, they derived the associated field equations for some specific cases. In [28][29][30][31][32], some interesting cosmological f (R, T ) models have been developed by employing various scenarios namely, auxiliary scalar field, dark energy models and anisotropic universe models. In literature [33][34][35][36][37][38][39][40][41][42][43][44][45][46][47], different cosmological applications of f (R, T ) gravity have been discussed like energy conditions, thermodynamics, exact and numerical solutions of field equations with different matter content, phase space perturbation, compact stars and stability of collapsing objects etc.
The existence and construction of the wormhole solutions is one of the most fascinating topics in modern cosmology. Wormholes are topological passage-like structures connecting two distant parts of the same universe or different universes together through a shortcut called tunnel or bridge. Generally, in nature, wormholes are categorized into two sorts namely static wormholes and dynamic wormholes [48]. For the development of wormhole structures, an exotic fluid (hypothetical form of matter) is required which violates the null energy condition (NEC) in GR. This violation of energy condition is regarded as one of the basic requirements for wormhole construction. The existence of the wormhole solutions in GR has always been a great challenge for the researchers. Although GR allows the existence of wormholes but it is necessary to first modify the matter sector by including some extra terms (as the ordinary matter satisfies the energy bounds and hence violates the basic criteria for existence of the wormhole). These extra terms are responsible for energy bound violation and hence permits the existence of wormhole in GR. In 1935, Einstein and Rosen [49] discussed the mathematical criteria of wormholes in GR and they obtained the wormhole solutions known as Lorentzian wormholes or Schwarzchild wormholes. In 1988, it was shown [50] that wormholes could be large enough for humanoid travelers and even permit time travel. In the literature [51][52][53][54][55][56][57][58][59] numerous authors constructed wormholes by including different types of exotic matter like quintom, scalar field models, non-commutative geometry and electromagnetic field etc. and obtained different interesting and physically viable results. Some important and interesting results regarding the stable wormhole solutions without inclusion of any exotic matter are discussed in [60,61]. In Ref. [62], the existence of the wormhole solutions and its different properties in f (R, T ) theory gravity has been discussed.
"On a D-brane, the coordinates may be treated as noncommutative operators", this is one of the most interesting aspect of non-commutative geometry of string theory, providing a mathematical way to explore some important concepts of quantum gravity [63,64]. Basically, noncommutative geometry is an effort to construct a unified platform where one can take the spacetime gravitational forces as a combined form of weak and strong forces with gravity. Non-commutativity has the important feature of replacing point-like structures by smeared objects and hence corresponds to spacetime discretization, which is due to the commutator defined by [x α , x β ] = iθ αβ , where θ αβ is an anti-symmetric second-order matrix. This smearing effect can be modeled by including a Gaussian distribution and a Lorentzian distribution of minimal length √ θ instead of the Dirac delta function. The spherically symmetric, static particle-like gravitational source representing the Gaussian distribution of non-commutative geometry with total mass M has an energy density given by [65][66][67] while, with reference to Lorentzian distribution, we can take the density function of particle-like mass M as follows: Here the total mass M can be considered as wormhole, a type of diffused centralized object and clearly, θ is the noncommutative parameter. The Gaussian distribution source has been utilized by Sushkov to model phantom-energy upheld wormholes [68]. Also, Nicolini and Spalluci [69] used this distribution to demonstrate the physical impacts of shortseparation changes of non-commutative coordinates in the investigation of black holes.
Motivated by this literature, in this manuscript, we will construct spherically symmetric static wormholes in the presence of curvature-matter coupling with non-commutative geometry. In the next section, we will describe the basic mathematical formulation of f (R, T ) gravity and the corresponding field equations for static spherically symmetric spacetime. In Sect. 3, we shall discuss the wormhole solutions for both Gaussian and Lorentzian distributions of non-commutative geometry by taking the linear model of f (R, T ) gravity, i.e., f (R, T ) = R + λT . Section 4 provides wormhole solutions for both these distributions of non-commutative geometry where the model f (R, T ) = R + α R 2 + γ R n + λT will be taken into account. In Sect. 5, the stability of these obtained wormhole solutions will be discussed through graphs. Section 6 will be devoted to an investigation of the gravitational lensing phenomenon for the exact model of Sect. 3 by exploring deflection angle at the wormhole throat. The last section will summarize the whole discussion by highlighting the major achievements.

Field equations of f (R, T ) gravity and spherically symmetric wormhole geometry
In this section, we shall discuss the basic formulation of f (R, T ) gravity and its corresponding field equations for spherically symmetric spacetime in the presence of ordinary matter. For this purpose, we take the following action of this modified gravity [27]: where f (R, T ) is an arbitrary function of the Ricci scalar R and the trace of the energy-momentum tensor T = g μν T μν .
Here L m represents the Lagrangian density of ordinary matter. By taking the variation of the above action, we have By contracting the above equation, we have a relation between the Ricci scalar R and the trace T of the energymomentum tensor as follows: 8π These two equations involves covariant derivative and d'Alembert operator denoted by ∇ and , respectively. Furthermore, f R (R, T ) and f T (R, T ) correspond to the function derivatives with respect to R and T , respectively. Also, the term μν is defined by The energy-momentum tensor for an anisotropic fluid is given by where V μ is the 4-velocity vector of the fluid given by V μ = e −a δ μ 0 and χ μ = e −b δ μ 1 , which satisfy the relations V μ V μ = −χ μ χ μ = 1. Here we choose L m = ρ, which leads to the following expression for μν : We relate the trace equation (3) with Eq. (2), then the Einstein field equations take the form given by The spherically symmetric wormhole geometry is defined by the spacetime: where (r ) and b(r ) both are functions of radial coordinate r and represent redshift and shape functions, respectively [50,70]. In the subsequent discussion, we shall assume the red shift function to be constant, i.e., (r ) = 0. Here the radial coordinate r is non-monotonic as it decreases from infinity to a minimum value r 0 , representing the location of wormhole throat, i.e., b(r 0 ) = r 0 , then it increases back from r 0 to infinity. The most important condition for existence of the wormhole is the flaring out property where the shape function satisfies the inequality: (b−b r )/b 2 > 0, while at the wormhole throat, it satisfies b(r 0 ) = r 0 . Further, the property b (r 0 ) < 1, is also a necessary condition to be satisfied for the wormhole solutions. Basically these conditions lead to NEC violation in classical GR. Furthermore, another condition that needs to be satisfied for wormhole solutions is 1 − b(r )/r > 0. All these conditions collectively provide basic criteria for the existence of a physically realistic wormhole model. In order to find the relations for ρ, p r and p t , we substitute the corresponding quantities for the metric (5) in Eq. (4) and then by rearranging the resulting equations, we have b where The curvature scalar R is given by and f R (R, T ) obeys the following expression: Since the above system, involving higher-order derivatives with many unknowns, is very complicated to solve for the quantities ρ, p r and p t therefore, for the sake of simplicity in calculations, we assume a particular form of the function f (R, T ) given by the relation f (R, where λ is a coupling parameter. After inserting this form of f (R, T ) and then by simplifying the corresponding Eqs. (6)-(8), we get

Wormhole solutions: Gaussian and Lorentzian distributions for f 1 (R) = R model
In this section, we shall consider a specific and interesting f (R) model [71,72] that is given by the linear function of the Ricci scalar: Using this relation in Eqs. (12)- (14) and after doing some simplifications, we get the following set of field equations: Here we include the smearing effect mathematically by substituting Gaussian distribution of insignificant width √ θ in place of the Dirac delta function, where θ is a noncommutative parameter of Gaussian distribution. Here we consider the mass density of a static, spherically symmetric, smeared, particle-like gravitational source given by The particle mass M, rather than being splendidly restricted at the point, diffused on a region of the direct estimate √ θ . This is because of the fact that the uncertainty is encoded in the coordinate commutator.
Comparing Eqs. (16) and (19), and then solving the resulting differential equation, we get the shape function b(r ) in terms of error function as follows: where Here C 1 is a constant of integration. Also, λ = −8π which clearly leads to b(r ) = 0. Using Eq. (20) in (16)-(18), we get the following relations for the ordinary energy density, tangential and radial pressures, which will be helpful to discuss the energy bounds: In the case of non-commutative geometry with reference to the Lorentzian distribution, we take the density function as follows: where M is a mass which is diffused centralized object such as a wormhole and θ is a non-commutative parameter. Comparing (16) and (24) and then solving the resulting differential equation, we get the following form of the shape function: where C 2 is an integration constant. Again using Eq. (25) in (16) and (18), we get a new set of equations which help us to discuss the energy conditions for the existence of wormhole structure. In this case, the expressions for energy density, radial and tangential pressures are given by Now we will present the graphical illustration of the obtained shape functions as well as the conditions that need to be fulfilled for existence of the wormhole. For this purpose, we take different suitable choices for the involved free parameters. Firstly, we check the behavior of the shape function b(r ) for the Gaussian distribution where the red shift function has been taken as a constant. The left graph of Fig. 1 indicates  the positive increasing behavior of the shape function and its right graph corresponds to the behavior of the shape function ratio to radial coordinate, i.e., b(r ) r which shows that as the radial coordinate gets larger values, the ratio b(r ) r approaches zero, and hence confirms the asymptotic behavior of the shape function. The left part of Fig. 2 indicates the behavior of b(r ) − r , which shows that the wormhole throat for this model is located at r 0 = 0.2 where b(r 0 ) = r 0 . In the right part of this figure, we check the flaring out condition for this model by plotting b (r ). It shows that at the wormhole throat r 0 = 0.2, clearly the condition b (r 0 ) < 1 is satisfied. The graphical behavior of density function as well as the null energy conditions ρ + p r and ρ + p t are shown in Figs. 3 and 4, respectively. It is clear from these graphs that the energy density function and the function ρ + p t indicate the positive but decreasing behavior versus radial coordinate while ρ + p r shows negative and increasing behavior and hence violates the NEC. Thus it can be concluded that the obtained wormhole solutions are acceptable in this modified gravity.  In the case of the wormhole solution with a Lorentzian distribution, the graphical behavior of the shape function as well as its corresponding properties is given in Figs. 5, 6 and 7. The left curve of Fig. 5 corresponds to the behavior of the shape function while the right graph shows the behavior of b(r ) r . It is clear from the curves that the shape function is positive and increasing satisfying the asymptotic flatness condition as r → 0. Figure 6 indicates the location of wormhole throat and the flaring out condition. It is seen that the wormhole throat is located at r 0 = 0.1 where the function b(r ) − r crosses the radial coordinate axis. Also, at this wormhole throat, the flaring out condition b (r 0 ) < 1 is satisfied for this case as shown in the right part of Fig. 6. The behavior of the energy density profile and the functions ρ + p r and ρ + p t is presented in Figs. 7 and 8. These show that the energy density remains positive and increasing with increasing values of r . Similarly, the function ρ + p t indicates the Here, we fix θ = 0.9, M = 0.0001, C 2 = 0.1 and λ = 2 positive but decreasing behavior whereas the function ρ + p r shows the negative increasing behavior versus r . This confirms the violation of NEC in this case and hence allows the existence of the wormhole. Thus in both cases, all necessary and important characteristics of the shape function for the existence of the wormhole are satisfied and thus it can be concluded that the obtained solutions are physically viable.
A comparison of Eqs. (19) and (30) (Gaussian distribution) yields the following non-linear differential equation: which is complicated and hence we solve it numerically for the shape function b(r ).
In a similar way, by comparing Eqs. (24) and (30) (Lorentzian distribution), we get the following non-linear differential equation: whose analytic solution is also not possible, thus we evaluate the possible form of the shape function by solving this equation numerically. Now we will discuss the behavior of the shape functions that are obtained by numerical approach as well as their corresponding important and necessary properties for the existence of wormhole structure for both Gaussian and Lorentzian distributions. For this purpose, we utilize a fixed value n = 3 for the modified model (29) which results in the cubic form given by f (R) = R + α R 2 + γ R 3 . For the other higher values, i.e., n > 3, it is observed that the resulting form of the shape function is not physically viable. For graphical illustration of the shape functions and their other properties, we will take different feasible values of the free parameters. The left part of Fig. 9 indicates that the shape function remains positive and increasing for the Gaussian distribution (obtained numerically), while its right part shows behavior of the function b(r ) r versus radial coordinate. Clearly, it indicates that as the radial coordinate increases, the function tends to zero and hence leads to the asymptotic behavior of the shape function. In Fig. 10, the left curve corresponds to   the function b(r ) − r , which provides the location of wormhole throat at r 0 = 0.001 where it cuts the r -axis. Its right curve provides information as regards the flaring out condition, i.e., b (r 0 ) < 1 which is clearly compatible at the obtained wormhole throat. Furthermore, the graphical illustration of energy density, tangential and radial pressures is given in Figs. 11 and 12. The left part of Fig. 11 corresponds to the energy density which shows a positive but decreasing behavior, while the right curve shows the graph of ρ + p t , which is also positive decreasing. Figure 12 indicates the behavior of ρ + p r , which is clearly negative and increasing versus r and hence violates the NEC. Thus all the conditions are satisfied allowing the existence of physically viable wormhole solution.
Similarly, for the Lorentzian distribution, the graphical behavior of the shape function and its properties like asymptotic behavior, wormhole throat and the flaring out condition are shown in Figs. 13 and 14. It can easily be observed that the obtained shape function is positive increasing and is compatible with all conditions. Further, the graphs for the resulting    16 This shows the development of ρ + p r versus r for the Lorentzian distribution density profile, ρ + p t and ρ + p r , are given in Figs. 15 and 16, respectively, which confirm the violation of NEC for this wormhole model. Thus it can be concluded that the obtained wormhole solutions for this cubic polynomial f 1 (R) model are physically interesting for both non-commutative distributions.

Equilibrium condition
In this segment, we explore the stability of obtained solutions using equilibrium conditions in the presence of Gaussian and Lorentzian distributions of non-commutative geometry. For this purpose, we take the Tolman-Oppenheimer-Volkov equation [59], which is given by where σ (r ) = 2 (r ). This equation determines the equilibrium state of configuration by taking the gravitational, hydrostatic as well as the anisotropic forces (arising due to anisotropy of matter) into account. These forces are defined by the following relations: and thus Eq. (33) takes the form given by Since we assumed the red shift function to be a constant so that (r ) = 0, we have F g f = 0 and hence the equilibrium condition reduces to the following form: We shall discuss the stability condition for both exact and numerical solutions in the presence of both distributions of non-commutative geometry. Firstly, we calculate F a f and F h f for the Gaussian distribution as follows: while, for a Lorentzian distribution, these are given by F a f = 6π 2 C 2 (θ + r 2 ) 2 − √ θ(8π + λ)Mr(3θ + 5r 2 ) + 3M(8π + λ)(θ + r 2 ) 2 tan −1 r √ θ 2π 2 (8π + λ)r 4 (θ + r 2 ) 2 , The graphical behavior of these forces is given in Figs. 17 and 18. The left graph indicates the behavior of these forces for the Gaussian distribution while the right graph corresponds to Lorentzian distribution for simple f 1 (R) model. It is clear from the graph that both these forces show the same but opposite behavior and hence cancel each other's effect and thus leaving a stable wormhole configuration. Similarly, we investigate the stability of numerical solutions for modified cubic f (R) model using both the Gaussian and the Lorentzian distributions. The graphical behavior of resulting forces is given by Fig. 18. Its left part corresponds to behavior of these forces for the Gaussian distribution whereas the right graph provides the behavior for the Lorentzian distribution which clearly indicates that these forces are also balancing each other's effect and thus leading to a stable wormhole structure.

Gravitational lensing effect of wormhole for simple f (R) = R model
In this section, we will explore the possible detection of traversable wormhole through gravitational lensing phenomena. For this purpose, we consider the static spherical symmetric metric involving x = r 2M representing the radius in Schwarzschild units and given by ds 2 = −A(x)dt 2 + B(x)dx 2 + C(x)(dθ 2 + sin 2 θ dφ 2 ). (34) Here the closest path taken by the light ray isx =r 2M . Here we will consider the obtained exact form of the shape function in the case of a simple linear f 1 (R) = R model (Sect. 3). The integration of this shape function from the wormhole throat r 0 to r is given as follows: Here clearly the coupling constant satisfies λ = −8π . Basically, we consider the form of static spherically symmetric wormhole metric (5) where e 2 (r ) = ( where v φ indicates the rotational velocity. In [75][76][77][78], it is pointed out that m = 0.000001, which is a very small value (nearly zero) and hence leaving the red shift function as a constant (as we assumed in previous sections). A comparison of these metrics leads to the following relations: The deflection angle for the light ray is given by Here d represents the mouth of wormhole because of exterior Schwarzschild line element, while the internal metric contribution is provided by I (x), which implies that the closest path taken by the ray of light is bigger than the mouth of the wormhole. This is defined by the relation In our case, this integral take the following form: representing the closest approach for the light ray to be inside the mouth of the wormhole. Here the function G(x) is given by In order to investigate the convergence/divergence of this integral, we can redefine the variable as y = x x for the sake of simplicity in calculations. Thus the integral takes the following form: In the integrand of the above integral, we can assume that Taylor's series can be used to expand the function H (y) around y = 1 as follows: Here we truncate the Taylor expansion up to second order where O(y − 1) 3 indicates the cubic and higher-order terms of factor (y − 1). It can easily be observed that the integral I (x) converges or diverges because of the leading term in the above expression. Integral can be convergent if the first (y − 1) 1/2 leads the expression where g(1) = 0. If g(1) = 0, then second term will lead the expression and whose integration will be ln(y − 1). Since y = 1, it turns out be undefined there and hence the integral diverges. If we choose the nearest approach of light ray as the wormhole throat, i.e.,r = r 0 , then consequently, we have y 0 = x 0 x and thus y 0 = 1. Using these values in f (y), it can easily be verified that f (1) = 0. Hence a photon sphere with radiusr (closest path taken by the light ray) equal to radius of the throat r 0 , can be found.

Conclusions
The existence and construction of the wormhole solutions in GR with some exotic matter has always been of great interest for the researchers. The presence of exotic matter is one of the most important requirement for wormhole construction as it leads to NEC violation and hence permits the existence of the wormhole. In the case of modified theories, construction of wormholes has become more fascinating topic as these include the effective energy-momentum tensor that violates NEC without inclusion of any exotic matter separately. In the present paper, we have constructed spherically symmetric wormhole solutions in the presence of two interesting Gaussian and Lorentzian distributions of non-commutative geometry in f (R, T ) modified gravity. For this purpose, in order to make system of equations closed, we assumed the function f (R, T ) = f 1 (R) + λT with two different forms of f 1 (R), i.e., the linear form f 1 (R) = R and f 1 (R) = R + α R 2 + γ R n , n ≥ 3.
Firstly, we talked about the possible wormhole construction for the linear f 1 (R) model with both Gaussian and Lorentzian distributions. For Lorentzian and Gaussian distribution, we found the exact solution. In order to examine the physical behavior of these obtained solutions, we plotted b(r ) versus radial coordinate. It is observed that shape functions show positive increasing behaviors for both these noncommutative distributions. Further we found the location of wormhole throats and analyzed some important characteristics of the shape functions namely asymptotic behavior, the flaring out condition and the violation of NEC using graphs. This discussion has been given in Figs. 1, 2, 3, 4, 5, 6, 7 and 8. It is concluded from these graphs that the obtained shape functions show asymptotic behavior, i.e., b(r ) r → 0 as r → ∞. Also, for both cases, wormhole throats are located at r 0 = 0.2 and r 0 = 0.1. Furthermore, the obtained shape functions are compatible with the flaring out condition and NEC as the function ρ + p r indicated negative behavior for both distributions. Thus the obtained solutions are viable permitting wormhole to exist in non-commutative f (R, T ) gravity.
Secondly, we checked the existence of the wormhole for model f 1 (R) = R + α R 2 + γ R n , n ≥ 3 by taking both non-commutative distributions into account. In this case, we obtained very complicated non-linear differential equations for b(r ) whose analytic solutions are not possible, therefore we solved them numerically. It is worthwhile to mention here that we fixed n = 3 for numerical solutions and their graphical behaviors as it is found that, for n > 3, the obtained solutions are not physically interesting (not meeting the necessary criteria for existence of the wormhole). In the left parts of Figs. 9 and 13, it is shown that the numerical solutions for b(r ) indicate increasing positive behavior. Other necessary conditions like asymptotic behavior of the shape function, flaring out condition as well as NEC have been given in Figs. 9, 10, 11, 12, 13, 14, 15 and 16. The wormhole throat for solutions in both distributions are located at r 0 = 0.001. Also, ρ + p r shows negative behavior and hence NEC is incompatible for this solution. Thus it is concluded that all the conditions are satisfied for the chosen specific values of free parameters and hence the obtained wormhole solutions are viable. It is also interesting to mention here that, for a different selection of free parameters θ, M, λ etc. (other than the used values in the present paper), all the functions show a similar graphical behavior as presented in the figures. Thus all the necessary conditions for existence of the wormhole will also be satisfied in these cases and hence the wormhole solution still exist.
Further, we examined the stability of obtained solutions using equilibrium condition given by Tolman-Oppenheimer-Volkov equation. Here we explored the stability for both models of f 1 (R) in the presence of Gaussian and Lorentzian distributions. After evaluating the possible expressions of anisotropic and hydrostatic forces for these cases, we examined them graphically as shown in Figs. 17 and 18. It can easily be observed from the graphs that these forces are almost equal in magnitude but opposite in behavior, therefore canceling each other's effect and hence leaving a balanced final wormhole configuration. Furthermore, we explored the possible detection of photon sphere at the wormhole throat. For this purpose, we followed the procedure given in Refs. [75][76][77][78] and explored the convergence of deflection angle. It is observed that, for the obtained exact solution for f 1 (R) = R, the resulting integral diverges at wormhole throat and hence it is concluded that a photon sphere, with radius r 0 (closest path taken by the light ray) equal to radius of the throat, can be detected.