Quasi-cosmological Traversable Wormholes in $f(R)$ Gravity

TL;DR

This paper explores traversable wormholes within $f(R)$ gravity, finding exact solutions that can satisfy energy conditions unlike in Einstein gravity, and matching various FRW metrics.

Contribution

It provides new exact wormhole solutions in $f(R)$ gravity that respect energy conditions, expanding possibilities beyond Einstein gravity.

Findings

Wormhole solutions can satisfy NEC at the throat in $f(R)$ gravity.

Some models respect WEC throughout the space.

Solutions match spherical, flat, and hyperbolic FRW metrics.

Abstract

In this paper we study traversable wormholes in the context of $f(R)$ gravity. Exact solutions of traversable wormholes are found by imposing the nonconstant Ricci scalar. These solutions asymptotically match spherical, flat and hyperbolic FRW metric. By choosing some static $f(R)$ gravity models, we verify the standard energy conditions for the asymptotically spherical, flat and hyperbolic wormhole solutions. Unlike the Einstein gravity, we find that in the context of $f(R)$ modified gravity, the asymptotically spherical, flat and hyperbolic wormhole solutions can respect the null energy condition (NEC) at the wormhole throat and near that. We find that in some static $f(R)$ models, asymptotically flat and hyperbolic wormholes respect the weak energy condition (WEC) through the whole space.

Equations

$$S=\frac{1}{2\kappa }\int d^{4}x\sqrt{-g}f(R)+\int d^{4}x\sqrt{-g}{\mathcal{L}}_m(g_{\mu \nu },\psi ),$$

$$F{R}_{\mu \nu }-\frac{1}{2}f{g}_{\mu \nu }-{\nabla }_{\mu }{\nabla }_{\nu }F+g_{\mu \nu }\square F=T_{\mu \nu }^m,$$

$$FR-2f+3\square F=T,$$

$$G_{\mu \nu }\equiv R_{\mu \nu }-\frac{1}{2}R{g}_{\mu \nu }=T_{\mu \nu }^c+{\tilde{T}}_{\mu \nu }^m.$$

$$T_{\mu \nu }^c$$

$${\tilde{T}}_{\mu \nu }^m$$

$$d{s}^2=-e^{2\Phi (r)}d{t}^2+\frac{d{r}^2}{1-b(r)/r}+r^2(d{\theta }^2+{\sin }^2\theta d{\varphi }^2).$$

$$d{s}^2=-d{t}^2+\frac{d{r}^2}{1-b(r)/r}+r^2(d{\theta }^2+{\sin }^2\theta d{\varphi }^2).$$

$$R=\frac{2b^{\prime }(r)}{r^2},$$

$$T_{\mu \nu }=(\rho +p_t)U_{\mu }{U}_{\nu }+p_t{g}_{\mu \nu }+(p_r-p_t){\chi }_{\mu }{\chi }_{\nu },$$

$$\frac{b^{\prime }}{r^2}$$

$$-\frac{b}{r^3}$$

$$\frac{b-b^{\prime }r}{2r^3}$$

$$\rho$$

$$p_r$$

$$p_t$$

$$R=6c_1+\frac{6c_2}{r^n},$$

$$b(r)=r^3\left(c_1+\frac{3c_2}{3-n}{r}^{-n}\right)+c_3,$$

$$i)b(r_0)=r_0,$$

$$ii)b^{\prime }(r_0)<1,$$

$$iii)1-\frac{b(r)}{r}>0.$$

$$b(r)=\left[\left(-r_0^n{c}_1+r_0^{n-2}\right)r^{3-n}+c_1{r}^3\right].$$

$$\frac{b(r)}{r}=c_1{r}^2+(r_0^{n-2}-c_1{r}_0^n)r^{2-n},$$

$$WEC1:\left[3c_1+(n-3)r^{-n}(c_1{r}_0^n-r_0^{n-2})\right]F>0,$$

$$\displaystyle\!\!\!\!\!\!\!{\rm WEC2\!:}\,\frac{1}{2r}\big{[}F\big{(}24\,c_{1}r\!+\!(4\!-\!n)(r_{0}^{n-2}\!\!-6c_{1}r_{0}^{n})r^{1-n}\big{)}$$

$$\displaystyle\hskip 19.91692pt+2F^{\prime}\!\left(c_{1}r^{2}\!\!-\!1\!+\!(r_{0}^{n-2}\!-6c_{1}r_{0}^{n})r^{2-n}\right)\big{]}\!>0,$$

$$WEC3:F\left[12c_1-(n-2)\left(r_0^{n-2}-6c_1{r}_0^n\right)r^{-n}\right]$$

$$+\frac{1}{2}{F}^{\prime }\left[12c_{1}r+(6c_1{r}_0^n-r_0^{n-2})\left(n-2\right)r^{1-n}\right]$$

$$+F^{\prime \prime }\left[6c_1{r}^2-1+(r_0^{n-2}-6r_0^n{c}_1)\right]r^{2-n}>0,$$

$$L=\frac{1}{2}\left[f(\varphi )R+h(\varphi )g^{\mu \nu }{\varphi }_{,\mu }{\varphi }_{,\nu }-2U(\varphi )\right]+L_m,$$

$$F=\frac{df(R)}{dR}>0,\frac{d^{2}f(R)}{d{R}^2}\ne 0.$$

$$f(R)=R-\mu R_{\ast }\tanh (\frac{R}{R_{\ast }}),$$

$$f(R)=R-\mu R_{\ast }\frac{(R/R_{\ast }{)}^{2m}}{(R/R_{\ast }{)}^{2m}+1},$$

$$f(R)=R+\lambda R_{*}\left[\big{(}1+\frac{R^{2}}{R_{*}^{2}}\big{)}^{-q}-1\right],$$

$$f(R)=R+\frac{R^{m}\big{(}aR^{m}-b\big{)}}{1+cR^{m}}\,,$$

$$f(R)=R-\mu R_{\ast }(R/R_{\ast }{)}^p.$$

Full text

1 Introduction

Traversable wormholes are throatlike geometrical structures which connect two separate and distinct regions of spacetimes and have no horizon or singularity. There are some hints in the paper of L. Flamm fl , to such spatial geometry. In 1935, Einstein and Rosen eros firstly obtained the wormhole solutions as bridge model of a particle . The study of Lorentzian wormholes in the context of general relativity (GR) was stimulated by the significant paper of Morris and Thorne in 1988 mtor . A fundamental ingredient in wormhole physics is the flaring-out condition of the throat, which in GR entails the violation of the NEC. Matter that violates the NEC is denoted exotic matter.

A very important challenge in wormhole scenarios is the establishment of standard energy conditions. In this regard, various methods have been proposed in the literature that deal with the issue of energy conditions within wormhole settings. Work along this line has been done in dynamical wormholes kar2 and thin-shell wormholes Mehdizadeh:2015dta , where the supporting matter is concentrated on the wormholes throat. In the context of modified theories of gravity, the presence of higher-order terms in curvature would allow for building thin-shell wormholes supported by ordinary matter thin1 . Recently, many people try to build and study wormhole solutions within the framework of modified gravity for instance, wormhole solutions in BransDicke theory bran1 , Born-Infeld theory brn1 , Einstein-Gauss-Bonnet theory gaus1 , Kaluza-Klein gravity kul1 and scalar-tensor gravity scal1 and Einstein-Cartan theory Bronnikov:2015pha ; Bronnikov:2016xvj ; ecart1 .

In recent years, theories of modified gravity have been studied involving cosmological objects like black holes, gravastars, strange stars and wormholes. Such theories are able to explain accelerating expansion of the universe and also solve the dark matter problemaccel1 . A well- known theory of modified gravity is f(R) gravity which modifies GR by replacing the gravitational action R by an arbitary function f(R), where R is the Ricci scalar. These models can be demonstrated to cause acceleration Starobinsky:1980te ; mod1 . It was found that higher-order curvature invariants can satisfy the energy conditions in f(R) gravity lobener . The power-law $R^{m}$ gravity is also studied in deb , the new static wormholes in f(R) theory using the non-commutative geometry are built in rahm1 and the cosmological evolution of wormhole solutions is investigated in char . Recently, the junction conditions in f(R) applied to the construction of thin-shell wormholes ero2 and pure double layer bubbles ero1 . The static wormhole geometries have also been theorized within curvature-matter coupling theories such as f(R, T) aziz1 .

Lorentzian wormhole solutions were also investigated in viable f(R) modified theories of gravity pmarko , which are consistent with observations of the solar system and cosmological evolution in solar . These wormholes were found in specific f(R) models which are considered to reproduce realistic scenarios of cosmological evolution based on the accelerated expansion of the universe. it was shown that the WEC holds at the vicinity of the throat for certain ranges of the free parameters of the theory. These asymptotically flat solutions analysed with the simple choices of shape function which WEC is satisfied at a specific point in space.

A large class of static black hole solutions with non constant Ricci scalar and traversable wormholes with constant Ricci scalar have also investigated in the background of $f(R)$ gravity saba , Aclas . The existence of wormhole solutions in scalar-tensor and $f(R)$ gravity is studied in Bronnikov:2006pt ; Bronnikov:2010tt and it is shown that static wormhole solutions in these theories obtain when the effective gravitational constant is negative (which means the anti-gravitational property of the solutions) in some region of space.

In this work, our aim is to obtain the wormhole solutions with non constant Ricci scalar in the context of viable f(R) gravity. This paper is organized as follows: After this introduction, we will present the basic field equations of f(R) gravity . In Sec. III, we will derive the properties of wormholes solutions within the framework of f(R) gravity and discuss about the energy conditions. The final section is devoted to our concluding remarks.

2 f(R) gravity wormholes

The action of f(R) modified theories of gravity is given by

[TABLE]

where $\kappa=8\pi G$ and ${\cal L}_{m}$ is the matter Lagrangian density, in which the matter field $\psi$ is minimally coupled to the metric $g_{\mu\nu}$. Throughout this work for notational simplicity we consider $\kappa=1$. One can find the field equation by varying the action (1) with respect to metric $g^{\mu\nu}$ as

[TABLE]

and the trace of the field equation

[TABLE]

where $F=df/dR$. Substituting the trace equation in (2) and re-organizing the terms, one can write the field equation as Lobo:2009ip

[TABLE]

In the above, the two terms in the right hand side of the field equation are in the form

[TABLE]

In this background, we are going to study the static and spherically symmetric traversable wormholes which is given by the following metric111The static and spherically symmetric wormhole geometry is given by

$ds^{2}=-e^{2\Phi(r)}dt^{2}+\frac{dr^{2}}{1-b(r)/r}+r^{2}\,(d\theta^{2}+\sin^{2}{\theta}\,d\phi^{2})\,.$

The solution is traversable when $\Phi(r)$ is finite. In the following we set $\Phi(r)=0$.

[TABLE]

The Ricci scalar for the above wormhole geometry is

[TABLE]

where “ ′ ” denotes derivative with respect to the radial coordinate $r$. As in Lobo:2009ip we choose the energy-momentum tensor of anisotropic distribution of matter

[TABLE]

where $U^{\mu}$ is the four-velocity and $\chi^{\mu}=\sqrt{1-b(r)/r}\,\delta^{\mu}_{r}$ is the unit spacelike vector in the radial direction. $\rho(r)$ is the energy density and $p_{r}(r),\,p_{t}(r)$ are radial and transverse pressure respectively, so the energy-momentum tensor takes to the form $T^{\mu}_{\,\,\,\,\,\nu}={\rm diag}[-\rho,p_{r},p_{t},p_{t}]\,$. Thus the field equation (4) leads to the following relationships

[TABLE]

where $H(r)=\frac{1}{4}\,(FR+\Box F+T)$. By solving the above system, one finds the following expressions Lobo:2009ip

[TABLE]

We are interested in studying inhomogeneous spacetimes which merge smoothly to the cosmological background. Therefore, we may consider the Ricci scalar of the wormhole geometry as

[TABLE]

where $c_{1},\,c_{2}$ and $n$ are free parameters. In the following section we consider some choices for these free parameters and study the corresponding wormhole solutions in some different $f(R)$ backgrounds.

3 Traversable wormhole solutions in $f(R)$ models

We start this section by finding wormhole solutions with Ricci scalar in the form (13). Combining (7) and (13) one can find the shape function as

[TABLE]

where $c_{3}$ is the integration constant. This wormhole solution by choosing $c_{2}=0$ is obtained in Aclas . In the following we set the integration constant to zero. The shape function $b(r)$ for a wormhole solution, should satisfy the following conditions:

[TABLE]

Inserting condition i), one can find $c_{2}$ in terms of $r_{0},\,c_{1}$ and $n$; so it is possible to rewrite the shape function as

[TABLE]

We consider the wormhole solutions with three values $c_{1}=0,\pm 1$. The shape function for these values are depicted in Fig 1. It is easy to check that in the case of $c_{1}=0,-1$ the condition ii) in (3) is satisfied for $n>2$ and in the case of $c_{1}=1$ condition ii) is satisfied when $n>2.4$ .

Setting the shape function into the form (16), one can write

[TABLE]

note that due to the condition ii) we should insert $n>2$ in the above, so the second term in (17) falls down at large $r$. Remembering the wormhole metric (6), one deduces that the wormhole solutions with $c_{1}=-1,\,c_{1}=0$ and $c_{1}=1$ at large $r$ match the hyperbolic, flat and spherical FRW universe respectively, so we call them asymptotically flat ($c_{1}=0$), asymptotically hyperbolic ($c_{1}=-1$) and asymptotically spherical ($c_{1}=1$) wormhole solutions.

In the framework of Einstein gravity, the traversable wormhole geometry (6), violates the weak energy condition WEC1: $\rho>0$ , WEC2: $\rho+p_{t}>0$ , WEC3: $\rho+p_{r}>0$  and the null energy condition WEC2: $\rho+p_{t}>0$ , WEC3: $\rho+p_{r}>0$ . Substituting $\rho,p_{r}$ and $p_{t}$ given in (10)-(12) and the shape function (16), one can rewrite the requirements of WEC as

[TABLE]

where $F=df(R)/dR$ and “ ′ ” denotes derivative with respect to the radial coordinate $r$.

The conditions of stability and absence of ghosts are discussed in Bronnikov:2006pt ; Bronnikov:2010tt . It is shown that in a scalar-tensor theory of gravity formulated in the Jordan frame with Lagrangian

[TABLE]

where $f$, $h$ and $U$ are arbitrary functions and $L_{m}$ is the matter Lagrangian, there is no static wormhole solution that satisfy the null energy conditions when $f(\phi)>0$ and f(\phi)\,h(\phi)+\frac{3}{2}\big{(}\frac{df}{d\phi}\big{)}^{2}>0. In the case of $f(R)$ theories of modified gravity $h(\phi)=0$ and $f(\phi)=F=df(R)/dR$ Bronnikov:2010tt , so the above nonexistence theorem holds for

[TABLE]

In the following we consider the asymptotically flat ($c_{1}=0$), asymptotically hyperbolic ($c_{1}=-1$) and asymptotically spherical ($c_{1}=1$) wormhole solutions in the background of $f(R)$ modified gravity. We are interested to the wormholes that respect the requirements of WEC, in the other word, we look for the solutions which satisfy (3), the conditions ii), iii) in (3) and violate the nonexistence theorem (20) simultaneously222Note that by $d^{2}f(R)/dR^{2}=0$ in (20), the nonexistence theorem violates just in one or some points instead, the $F<0$ violates this theorem in a region at space. Hence in the following we search for $F<0$ to check the violation of this nonexistence theorem.. Note that the regions where $df(R)/dR<0$ are anti-gravitational, in the sense that the effective gravitational constant is negative Bronnikov:2006pt ; Bronnikov:2010tt .

3.1 The case of $c_{1}=0$

For a wormhole solution, as we mentioned above, the second term in (17) falls down at large $r$, therefore in the case of $c_{1}=0$ the traversable wormhole solution (6) matches the flat FRW metric at large $r$. At this stage we consider such asymptotically flat wormholes in the background of some static $f(R)$ gravity models. As we will see, the $f(R)$ modified gravity wormhole solutions in many cases satisfy the null or weak energy condition.

1) The Tsujikawa model Tsujikawa:2007xu ; DeFelice:2010aj ; luca

One of $f(R)$ models which is consistent with cosmological and local gravity conditions, proposed by Tsujikawa. In this viable model, the $f(R)$ considered as

[TABLE]

where $0.905<\mu<1$ and $R_{*}$ is a free small positive parameter Tsujikawa:2007xu . For this model one can easily finds $F=1-\mu+\mu\tanh^{2}(R/R_{*})$. It is obvious that $F$ takes positive values when $\mu<1$, this means that by choosing $\mu$ in the consistent range $0.905<\mu<1$, due to (20), there is no static wormhole solution in this model that satisfy the NEC (and so WEC). However by taking $\mu$ at the neighborhood of this range, one can find wormholes that respect WEC. In figure 2.a we have depicted $F$ in (20) and the requirements of WEC (3) for a solution with $\mu=1.02$ and $n=5$. It is obvious in this figure that $F<0$ , therefore the condition of nonexistence theorem (20) is violated. The figure also shows that WEC is respected at the wormhole throat $r_{0}$.

It is worth to mention that in this model, there is a 4-dimensional parameter space for the wormhole solutions which is made of $\mu,R_{*},r_{0}$ and $n$. Among these parameters, $\mu,R_{*}$ are limited by the model Tsujikawa:2007xu , we also measure the radial distance $r$, in the unit of wormhole throat $r_{0}$ (the WEC1,2,3 in figure 2.a are depicted in term of $r/r_{0}$). Now we are going to find the parameter $n$ in a manner that wormhole solutions respect WEC through the space. In the other word we search for regions in $n-r/r_{0}$ diagram that (3) and conditions ii), iii) in (3) are satisfied and $F$ in (20) takes a negative value simultaneously.

Figure 2.b shows $n-r/r_{0}$ diagram. The blue region in this figure corresponds to asymptotically flat traversable wormhole solutions in the background of (21), which respect the WEC. It is obvious that by choosing $4<n<5.2$, the WEC is respected at the throat $r_{0}$ and through the whole space outside that. Note that unlike the Einstein gravity, considering $f(R)$ gravity model in the form (21), it is possible to find traversable wormhole solutions that respect the WEC at the throat and beyond that.

2) Hu-Sawicki model

The second $f(R)$ model which satisfies both cosmological and local gravity constraints proposed by Hu and Sawicki Hu:2007nk ; luca . This model is in the form

[TABLE]

where $m,R_{*}$ and $\mu$ are positive parameters. It is shown in Tsujikawa:2007xu that $R_{*}$ takes a small positive value and for $m=1$ one should insert $\mu\geq 8\sqrt{3}/9$. Similar to the previous case, we are looking for the wormhole solutions of this $f(R)$ model which respect the NEC or WEC. In figure 3.a by setting $\mu=5$ and $n=2.8$, we have depicted $F$ in (20) and WEC2,3 (3). It is clear from the figure that around the throat $r_{0}$, $F<0$ and WEC2,3 are positive, which means that this wormhole solution respect the null energy conditions.

Figure 3.b shows the $n-r/r_{0}$ diagram for this model. The blue zone in this figure is the region that $F<0$, ${\rm WEC2>0}$, ${\rm WEC3>0}$ and ii), iii) in (3) are satisfied simultaneously. This means that traversable asymptotically flat wormhole solutions in this model respect the NEC around the throat $r_{0}$ when $2<n<3$ and by setting $n\sim 2.8$ the NEC is respected at the throat.

3) The Starobinsky model Starobinsky:2007hu ; Amendola:2006kh ; luca

Another viable modified $f(R)$ gravity model that we consider, has three free parameters $\lambda,R_{*}$ and $q$ as:

[TABLE]

where $0.944<\lambda<0.966$ for $q=2$ according to Tsujikawa:2007xu and $R_{*}$ takes small positive values. Similar to the previous case, we are looking for the wormhole solutions of this $f(R)$ model which respect the WEC (3). One can check that by choosing $0.944<\lambda<0.966$, $F$ in (20) takes positive values and there is no static wormhole solution in this model. However by taking $\lambda$ at the neighborhood of this range, it is possible to obtain negative values for $F$. In figure 4.a we draw $F$ and (3) by choosing $\lambda=0.98,n=2.5$. It is clear that by this choice, $F$ is negative at the throat however the null (and so weak) energy condition is violated at the wormhole throat.

We also tried -by varying the parameter $n$- to find wormhole solutions that respect NEC, but there is no region in $n-r/r_{0}$ diagram that $F<0$ and the NEC is satisfied simultaneously (the $n-r/r_{0}$ diagram is blank). In the other words, our results show that there is no asymptotically flat wormhole in this model that respect the NEC.

4) Nojiri-Odintsov model n-o

Another viable $f(R)$ model suggests that

[TABLE]

where $m\geq 2$ and $a,b$ and $c$ are free positive parameters which satisfy the condition $bc>>a$. Setting $m=2,\,a=2,\,b=20,\,c=20$ and $n=2.8$, we checked that $F<0$ and WEC2, WEC3 are respected around the wormhole throat $r_{0}$ (figure 5.a). Figure 5.b shows that the $f(R)$ model (24) contains asymptoticly flat wormhole solutions with different values for the parameter $n$, that respect the null energy condition.

5) Amendola-Gannouji-Polarski-Tsujikawa model

In the fifth model that we consider, the $f(R)$ is in the form luca ; agp

[TABLE]

This viable model contains three parameters $0<p<1$, $\mu>0$ and $R_{*}$ which takes a small positive value. The same analysis by setting $p=1/2,\,\mu=5$ and $n=2.75$, shows that $F$ in (20) becomes negative and the obtained wormhole solution satisfies NEC at the throat (figure 6.a). It is also obvious in Fig(6.b) that there are asymptoticly flat traversable wormhole solutions in this model which respect NEC outside the throat $r_{0}$, if one choose adequate values for $n$.

6) Exponential gravity model Cognola:2007zu ; Elizalde:2010ts

For the last model that we consider in this section, the $f(R)$ is in the form

[TABLE]

where $\lambda,R_{*}$ are free positive parameters of the model. Similar to the previous cases, one can find a wormhole in this model which respect NEC at the throat (figure 7.a). In figure 7.b the $n-r/r_{0}$ diagram is depicted. It is clear that there is a little possibility for traversable asymptotically flat wormhole solutions in this model to satisfy the null energy conditions around the throat $r_{0}$.

3.2 The cases of $c_{1}\neq 0$

We considered traversable wormhole solution as (6), where $b(r)/r$ given by (17). Remember that the wormhole shape function $b(r)$ satisfies the conditions in (3) which leads to $n>2$ in the case of $c_{1}=-1$ and $n>2.4$ in the case of $c_{1}=1$, therefore the second term in (17) falls down at large $r$. In this subsection we choose $c_{1}=\pm 1$, for which the wormhole metric (6) at large $r$ matches hyperbolic ($c_{1}=-1$) and spherical ($c_{1}=1$) FRW metric, so we call them asymptotically hyperbolic and asymptotically spherical wormhole solutions respectively. In the following we verify the null and weak energy condition for these asymptotically hyperbolic and spherical wormhole solutions in the background of static $f(R)$ modified gravity models that we considered in the previous subsection.

1) The Tsujikawa model

For this model the $f(R)$ is given in (21). We are interested in solutions for which $F$ in (20) takes negative values, conditions ii), iii in (3) are satisfied and (3) are respected simultaneously. Figures 8.a,c show the location of these asymptotically closed ($c_{1}=1$) and open ($c_{1}=-1$) wormhole solutions in the $n-r/r_{0}$ diagram respectively. In figure 8.b (8.d) we have explicitly depicted $F$ and requirements of the WEC for an asymptotically spherical ( hyperbolic) wormhole solution with $n=4.5$ ($n=2.02$). It is obvious that by choosing adequate values for the parameters, one can find traversable asymptotically closed and open wormhole solutions in this $f(R)$ gravity model that respect NEC around the wormhole throat.

2) Hu-Sawicki model

The second model that we studied its wormhole solutions at the previous subsection was the Hu-Sawicki model with the $f(R)$ in the form (22). We verify the NEC for asymptotically spherical and hyperbolic wormhole solutions of this model in Fig(9a,c) respectively. In the blue regions in these figures $F=df(R)/dR$ becomes negative, ii), iii) of (3) are satisfied and (3) are respected simultaneously. We have plotted the NEC for an asymptotically spherical and hyperbolic wormhole solution in figure 9.b,d respectively.

Similar to the previous model, it is clear that by setting the parameters in adequate values, traversable asymptotically spherical and hyperbolic wormholes in this model respect the null energy condition around the throat $r_{0}$.

3) The Starobinsky model

In the case of Starobinsky model which $f(R)$ is given in (23), the $df(R)/dR<0$ could be satisfied for both asymptotically spherical and hyperbolic solutions by choosing $R_{*}$ and $r_{0}$ appropriately, however the null energy condition is violated around the throat by these wormhole solutions (Fig10). In this model the $n-r/r_{0}$ diagram is blank which means that there is no wormhole solution in this case which respect the NEC.

4) Nojiri-Odintsov model

Now we return to the Nojiri-Odintsov modified gravity model where the form of $f(R)$ is given in (24). The $n-r/r_{0}$ diagram for asymptotically hyperbolic ($c_{1}=-1$) solutions in this model by setting $m=2$, is depicted in Fig(11.a). The blue zone in this diagram corresponds to the region that $F<0$, ii), iii) of (3) and (3) are satisfied simultaneously. It is clear that the asymptoticly hyperbolic wormholes with $n>2$ in this model, respect the WEC almost through the whole space outside the throat $r_{0}$. In figure (11.b) we have drawn $df(R)/dR$ and WEC1,2,3 explicitly for a wormhole solution with $n=3.2$ .

In the case of $c_{1}=1$, it is possible to set the parameters in a manner that $F<0$ is satisfied, however these asymptotically spherical wormhole solutions do not respect the NEC (see Fig(11.c) for instance), in the other words, the $n-r/r_{0}$ diagram in this case is blank.

5) Amendola-Gannouji-Polarski-Tsujikawa model

The $f(R)$ in this model is given by (25). In the case of $c_{1}=-1$ we check the NEC for wormhole solutions of this model in figure (12.a) which shows that it is possible to find asymptotically hyperbolic wormholes that respect NEC around the throat, by choosing appropriate parameters in this static model. In figure (12.b) we have explicitly drawn the NEC for a wormhole by setting $n=2.3$ .

In the case of $c_{1}=1$ our results show that, like the previous model, there is no asymptoticly closed wormhole solution in this model which respect NEC (the $n-r/r_{0}$ diagram is blank). As an example we have depicted $F=df(R)/dR$ and WEC2,3 for a wormhole by setting $n=5$ in figure (12.c).

6) Exponential gravity model

For the last model that we considered in the previous subsection, the $f(R)$ is given by (26). The same analysis shows that in the case of $c_{1}=-1$, there are asymptotically hyperbolic wormhole solutions which respect the WEC almost through the whole space outside the wormhole throat (figure 13.a,b). It is also possible in the case of $c_{1}=1$, to find asymptoticly spherical wormhole solutions by choosing parameters in a manner that $F<0$. However the NEC are violated by these wormholes as we can see in the Fig(13.c).

4 Conclusion

In this paper we considered wormhole solutions in the context of some $f(R)$ modified gravity models. In the first step, by solving the field equations similar to Lobo:2009ip , we found the energy density, radial and transverse pressure (10)-(12) of the solutions in terms of $F=df(R)/dR$ and the shape function $b(r)$. In the second step we considered the Ricci scalar of the wormhole solutions in the form (13) and reorganized the shape function into the form $b(r)=\left[\left(-r_{0}^{n}\,c_{1}+r_{0}^{n-2}\right){r}^{3-n}+{c_{1}\,r^{3}}\right]$.

As the next step we inserted the physical conditions on $b(r)$ and found that $n$ must be greater than 2. This convey us to choose the values $c_{1}=-1,0,+1$ for which the wormhole solution (6) at large radial coordinate $r$, matches the hyperbolic, flat and spherical FRW universe respectively (asymptotically hyperbolic, flat and spherical solutions). In the final step we chose six static $f(R)$ modified gravity models and checked whether the null and weak energy conditions (3) hold for the wormhole solutions in the background of these $f(R)$ models.

Our results show that by choosing adequate parameters, it is possible to find traversable asymptotically hyperbolic and flat wormhole solutions that respect NEC around the throat $r_{0}$, in all these $f(R)$ models except that the Starobinsky model. In the case of the Tsujikawa model, asymptoticly flat wormholes can also respect the weak energy condition through the whole space outside the throat. In the case of exponential gravity and Nojiri-Odintsov models, asymptoticly hyperbolic wormholes can respect WEC in the whole space outside the $r_{0}$ too. We also found that asymptotically spherical wormhole solutions respect the NEC around $r_{0}$ just in the Tsujikawa and Hu-Sawicki models. It is worth to mention that due to $F=df(R)/dR<0$ for these wormhole solutions, the “remote mouths” would be located in anti-gravitational regions because of the effective gravitational constant is negative Bronnikov:2006pt ; Bronnikov:2010tt .

As a final remark we emphasis that unlike the Einstein gravity, our results show that it is possible to find traversable wormhole solutions which respect the null and weak energy conditions in the background of some $f(R)$ modified gravity models. This means that $f(R)$ gravity in addition to its applications in cosmology, can be a useful framework in wormhole physics.

Acknowledgment

H.G. would like to thank M.M. Sheikh-Jabbari, Ghadir Jafari and Davood Mahdavian Yekta for useful discussions.