A new class of f(R)-gravity model with wormhole solutions and cosmological properties

A spherically symmetric wormhole family of solutions, with null red-shift, in the context of f(R)-gravity is presented. The model depends on two parameters: m and β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document} and meets all requirements to be an asymptotically and traversable wormhole. To solve the field equations, an EoS is imposed: p⊥=-ρ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_{\perp }=-\rho $$\end{document}. It is found that for m=1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m=1$$\end{document} the solution satisfies the null energy condition, although F(R)<0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F(R)<0$$\end{document} everywhere. For m=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m=0$$\end{document}, the model satisfies the null energy condition away from the throat, where the function F(R) is everywhere positive and together with dF(R)/dR vanish at the throat of the wormhole. This fact is beyond the scope of the non-existence theorem. Furthermore, the cosmological viability of the model, to address the late – time accelerated epoch, is analyzed on the background of a flat FLRW space-time. The model satisfies consistency of local gravity tests, stability under cosmological perturbations, ghosts free and stability of the de Sitter point.


Introduction
The Einstein-Rosen bridge was developed in order to explain the particle problem in General Relativity, more precisely to try to explain fundamental particles such as electrons in terms of space-time tunnels threaded by electric lines of force [1]. However, a major problem is the instability of the geometric structure present in this kind of solutions within the arena of General Relativity. This means that the bridge is not able to remain open long enough for an object to pass through it (not even a photon). As far as this is concerned, Wheeler [2] tried to explain Einstein-Rosen bridge (today known as wormholes) in terms of topological entities called geons and, incidentally, provided the first (now familiar) diagram of a wormhole as a tunnel connecting two openings in different a e-mail: alvaro.restuccia@uantof.cl b e-mail: francisco.tello@ua.cl (corresponding author) regions of space-time. Basically, the geometry of a wormhole is formed by two mouths joined by a throat connecting two parallel Universes or two far away regions of the spacetime. Following Morris and Thorne pioneering work [3,4] the wormhole exists only if the wormhole throat is supported by the so-called exotic matter field. By exotic matter Morris and Thorne mean negative radial pressure p r , which entails the violation of null energy condition due to | p r | > |ρ| in the throat of the wormhole. What is more if an observer crosses the tunnel it will see a negative energy density ρ. Regarding this, many works available in the literature have addressed the energy violation problem at the wormhole throat, for example working out dynamical wormholes [5,6]. Moreover, theories containing high derivative terms in the gravitational sector through functions of the scalar curvature, enable to construct thin-shell structures driven by normal matter distributions [7,8]. Other interesting works concerning the study and analysis of wormholes in gravity theories include scalartensor gravity, such as Brans-Dicke theory, non-linear electrodynamics such as Born-Infield theory or higher dimensional theories like Kaluza-Klein, Einstein-Gauss-Bonnet or Einstein-Cartan theories [9][10][11][12][13][14][15][16][17][18][19][20][21][22].
Besides, modified gravity theories such as f (R) and f (R, T ) gravity are promised scenarios to investigate wormhole regions satisfying energy conditions. These theories attempt to explain the most recent data about the accelerated expansion of the Universe. This phenomenon has not been explained from the point of view of General Relativity and it is confirmed and supported by numerous observational data from Supernova type Ia [23], high Planck data [24] and large scale structure [25][26][27][28][29][30][31][32][33][34][35][36]. On the other hand, a possible explanation of this phenomenon is the so-called dark energy [37] a mysterious cosmic fluid, which has uniform density distribution and negative pressure.
So, as stated before, modified gravity theories are an alternative way to explain the above issue. f (R)-gravity [38] was the first modified gravity theory employed to explain the singularities introduced by the isotropic homogeneous cosmological model. After its presentation in the early 80's, Starobinsky [39] used this theory to face questions related to cosmic inflationary models. The main point of this theory is to introduce a function f (R) of the Ricci scalar R. Its consequence is the presence of high derivative terms, which could in principle explain the accelerated expansion and the existence of dark energy or dark matter without the addition of extra matter fields. Besides f (R, T ) gravity theory [40] can be seen as a generalization of f (R)-gravity. Particularly, f (R, T ) theory also modifies the material sector by including the trace T of the energy-momentum tensor, which breaks down the minimal-coupling matter principle between the gravitational sector and the material one. The study of wormhole solutions satisfying energy conditions in the arena of f (R, T ) theory was considered in [22,[41][42][43][44]. However, the study of wormhole geometries is more robust in the context of f (R)-gravity than it is performed in f (R, T ) theory. For example, it was found in [45][46][47] that high derivative terms coming from the f (R) function regularize energy conditions. In [48,49] the study of static wormhole solutions using power-law models R n was performed. Cosmological evolution and non-commutative techniques were treated in [50][51][52][53]. More recently the matching condition formalism to build thin-shell was done in [54][55][56] and the study of new wormhole solutions was performed in [57].
In a broader cosmological context, f (R)-gravity is a viable theory to address at the same footing the early-time inflation with the late-time acceleration process [58][59][60][61][62]. To face the ultra accelerated expansion or inflationary age of the Universe, the seminal work by Starobinsky [39] is considered as the building block in this direction. On the other hand, the present accelerated era of the Universe represents nowadays a great challenge to the astrophysical and cosmological community. Several works available in the literature have attempted to tackle this problem by proposing different f (R)-gravity models [63][64][65][66][67][68][69][70][71][72][73][74][75][76]. To be a viable models with a physical meaning, all these proposal should fulfill some general requirements [77]. What is more if the scalar curvature is convariantly constant ı.e, R ≡ R 0 where R 0 is representing the scalar curvature at the present epoch, the model contains a special kind of solutions: vacuum solutions, including the vacuum de Sitter space-time with cosmological constant (General Relativity solution) [78].
Motivated by these prior studies, we investigate the existence of wormhole solutions in the arena of f (R)-gravity. In order to achieve it we propose a new ansatz for the shape function b(r ) depending on two parameters m and β, and also impose an equation of state relating the components of the energy-momentum tensor threading the throat of the wormhole, specifically ρ ⊥ = −ρ. These components allow us to solve a first order differential equation for F(R(r )) and also find an explicit expression for the null energy condi-tion (NEC). To facilitate the mathematical treatment of the field equations we have set the red-shift function to be null. The resulting model satisfies all the basic requirements to be an asymptotically and traversable wormhole. Moreover, it is found that the matter content threading the throat satisfies the NEC in one case while in the second case this condition is partially violated near the throat, although it is satisfied away from it. To complete the study we emphasize that the main result on this scenario is the non-existence theorem proved in [79] and extended in [80,81]. It was proved in the Einstein frame in [79] and extended to scalar tensor theories and f (R)-gravity by performing a conformal map from the Jordan frame to the Einstein one [80,81]. The main point is that the wormhole topology together with the asymptotically flat spaces on both sides of the throat are preserved on the conformal map together with the NEC. In this concern, the feasibility of the existence of wormhole solutions in f (R) theory without ghost fields is subject to the violation of NEC at least on the throat of the wormhole. In this context one class of solutions we are going to present satisfies the null energy condition away from the throat, where the function F(R) is everywhere positive and together with d F(R)/d R vanish at the throat of the wormhole. This fact is beyond the scope of the nonexistence theorem, since the conformal map from the Jordan frame to the Einstein frame is not well defined at the throat. In this concern, the incompatibility between both frames was previously noted in the cosmological scenario, in the treatment of f (R) singularities (see [82] for a recent and detailed discussion about the correspondence of f (R)-gravity singularities in Jordan and Einstein frames in the cosmological context) and then translated into the wormhole studies. On the other hand, to check the feasibility of our model, we have explored the cosmological properties of the obtained f (R) models to address the latetime cosmic acceleration of the Universe (issue related with the dark energy problem). In this regard, the resulting f (R) Lagrangian satisfies the general requirements to face this point [77], such as: (i) ghost free consistency, (ii) consistency with local gravity and stability under cosmological perturbations, (iii) stability at late-time de Sitter point. On the other hand, we have extended our model to include the inflationary phase, adding a two step model [62], although it represents a toy inflationary-accelerated unification model to describe the early and late time epochs of our Universe. However, it shows the possibility of a natural extension of our model to address this problem. Furthermore, we have solved the field equations in the cosmological scenario to obtain the Hubble rate H (t) and scale factor a(t). To do it, we have imposed a flat Friedmann-Lemaitre-Robertson-Walker (FLRW) metric accompanied with an isotropic perfect fluid distribution. It is worth mentioning that f (R)-gravity offers in the cosmological framework a natural scenario to address the existence of an effective dark energy, without the need to introduce a negative pressure ideal fluid [59,83].
The article is organized as follows: Sect. 2 presents the wormhole anatomy. Section 3 develops the f (R) formalism, presenting the geometry and field equations. In Sect. 4 the field equations are solved and the model is studied in details. In Sect. 5 the main features of the models are discussed and in Sect. 6 some cosmological properties related with the accelerated expansion at late time of the Universe (dark energy) are discussed, and the Hubble rate and scale factor satisfying the field equations in the flat FLRW regime with an isotropic fluid are shown. Furthermore, the connection with the inflationary stage is also addressed. Finally Sect. 7 provides some conclusions.

A quick review: WormHoles (WH) generalities
In curvature coordinates, the static and spherically symmetric line element describing a WH region is given by where d 2 ≡ sin 2 θ dφ 2 + dθ 2 . In Eq. (1) the metric functions depend only on the radial coordinate ı.e, = (r ) and = (r ). The function is referred as the red-shift function while e is related with the so called shape function b(r ), explicitly To be a traversable and asymptotically flat space-time we assume the WH geometry must satisfy the following general requirements [3,84] 1. The radius of the WH throat is defined as It is assumed that there is only one global minimum. Furthermore, the proper distant l is related with radial coordinate r by where l ∈ (−∞, +∞). The ± refer to both mouths of the WH. 2. For comparison with Schwarzschild solution the mass of the WH as seen for spatial infinity is given by b ± → 2G M ± .
3. There are two coordinate patches covering the range r ∈ [r 0 , +∞). Each patch covers an asymptotically flat region, and both match at r 0 , the throat of the WH. 4. The red-shift function must be finite everywhere, for r ≥ r 0 , in order to avoid an event horizon, so e > 0 for all r > r 0 . Moreover, where 0 a finite real number. 5. The shape function b(r ) must fulfill b(r 0 ) = r 0 , where the throat r 0 of the WH defines the spherical surface r = r 0 . 6. For all r > r 0 ⇒ b(r ) < r and b (r 0 ) < 1 or equivalently b (r ) < b(r )/r . This condition is known as the flare-out condition. Here b (r ) means db(r )/dr.
In the next section we discuss the explicit form of the field equations and the matter content surrounding the WH throat.

The f (R) formalism
In this section we present the f (R)-gravity formalism. The starting point is the modified Einstein-Hilbert action minimally coupled to a generic matter field given by where f (R) is an arbitrary smooth function of the Ricci's scalar R, g is the determinant of the metric tensor g μν (g μν its inverse), L m is the matter Lagrangian density, encodes the matter fields and κ = 8π G/c 41 . So, by taking variations with respect to g μν in the action (3) we arrive at the following general field equations where as always F ≡ d f/d R. Moreover, the right hand side of Eq. (4) represents the usual matter content. Considering the trace of the Eq. (4) we obtain where is the D'alembertian differential operator defined by The Eq. (5) shows that F(R) has a dynamical role in the theory ı.e, it is a fully dynamical degree of freedom. Furthermore, we have denoted the trace of the energy-momentum tensor by T (m) = T μ(m) μ . In the following we shall consider that the matter content threading the wormhole throat is taken to be an imperfect fluid distribution described by the following energymomentum tensor with ρ the energy density, p r and p ⊥ being the pressure waves in the principal directions ı.e, the radial and tangential ones respectively. The four-velocity of the above fluid distribution is characterized by the time-like vector u ν . Moreover χ ν is a unit space-like vector in the radial direction (orthogonal to u ν ) that is, Besides, throughout the study we shall assume a constant red-shift function, specifically = 0. This assumption implies that the gravitational redshift z is totally vanishing, of course Furthermore, from the mathematical point of view a vanishing red-shift function reduces the complexity of the field equations to be solved. So, taking into account the above considerations Eqs. (4) and (7) yield to the following field equations where F using Eqs. (1) and (6) is given by where primes denote differentiation with respect to the radial coordinate r . There are many ways to tackle the system (9)- (11). One possibility is to impose an adequate f (R) model supplemented by a suitable shape-function b(r ) or supplement the system with an equation of state relating the thermodynamic variables and impose a suitable shape function.
In this regard we follow the second approach. Nevertheless, among all the equations of state or possible equivalent relationships, only some are mathematically treatable, because the resulting expressions are highly complex to solve (at least analytically). So, following [47] we established p ⊥ = −ρ as a supplementary equation of state.

The Model
As was discussed in the previous section, in order to solve the f (R)-gravity system of equations (9)-(11), we impose the following shape function b(r ), where m, α and β are constant parameters. The parameter α and β have units of [length] and [length] 1+m respectively. To ensure the fulfillment of b(r 0 ) = r 0 and b (r 0 ) < 1 at the throat of the WH one has Furthermore, from the flare-out condition one obtains a bound on the constant parameter β as follows Besides, in to order to satisfy the condition b(r ) < r for all r > r 0 we have for m = 0 It follows that for all β > 0, Eq. (17) is always satisfied. We shall then restrict β to be strictly positive. Then the line element (1) representing the WH geometry is given by The above geometry reproduces the Minkowski spacetime geometry when r → ∞. Furthermore, it represents a family of WH depending on the choice of the m parameter. In this respect, power law shape functions have been considered by several authors in the arena of WH solutions.
Next, by using Eq. (15) and the equation of state p ⊥ = −ρ into the set of Eqs. (9)-(11) one arrives at the following first order differential equation for F(R(r )), The general solution of (19) is given by being F 0 an integration constant. Now we examined the output expressions obtained from Eq. (20) by taking different values for the running constant parameters m and β satisfying the mentioned requirements.
The election m = 1 yields from Eq. (20) to where A and B are defined by As it is observed from Eq. (21) the sign of the F(R(r )) function depends on the sign of the integration constant F 0 . If the sign of the F(R(r )) function is positive or negative is a very important issue regarding the existence of WH solutions in the arena of f (R)-gravity theory. As we will see later this fact is related with the fulfillment of energy conditions and ghost fields existence in the model. It should be noted that in order to have a dimensionless F(R(r )) function the dimensionality of the integration constant F 0 depends on the factors A and B.

Solution #2: m = 0
This simple choice provides from Eq.
In this case, in order to satisfy the flare-out condition one can express β = γ r 0 , where γ is a constant parameter. To satisfy b (r 0 ) < 1, γ must be restricted to be γ > −1. Furthermore, from b(r ) < r we get The roots of the above inequality are r + = r 0 and r − = −γ r 0 . It is clear that r + > r − , then to cover the range being C ≡ (β −r 0 )/2(β +r 0 ). Again, the sign of the F(R(r )) function depends on the integration constant F 0 , however this time F 0 is dimensionless. It should be noted that the point r = r 0 makes the expressions (21) and (25) singular or null. This depends on the relation between r 0 and β which determine the sign of A, B and C.
Next, the f (R(r )) function can be obtained from where the general expression for the Ricci scalar corresponding to the line element (18) is given by To reduce the mathematical complications in the case m = 1 we shall assume β = r 2 0 . In this way the expression (21) becomes where the above choice leads to A = B = −1/6 and the resulting scalar curvature is So, combining Eqs. (26), (28) and (29) one gets  25) and (27) the following F(R(r )) function and Ricci scalar are obtained Thus, putting together Eqs. (26), (31) and (32) one arrives to The trend of the f (R(r )) function and the Ricci scalar R(r ) are depicted in Fig. 1. It is appreciated that R vanishes when r tends to infinity as expected (flat space-time).

Results and discussions
In this section, we discuss the behavior of the shape function b(r ), energy conditions and the presence of ghost fields of the WH models we have obtained.

Shape function behavior
As was pointed out in Sect. 2 the shape function b(r ) must respect some basic requirements in order to have an asymptotically and traversable WH space-time. So, taking into account these rules, the proposed ansatz given by (15) fulfills all the formalities. Furthermore, the above requirements are ensured only in the range m + 1 > 0 ⇒ m > −1. This is so because, if m is taken to be less than −1, then the shape function will not be finite when r goes to infinity. Moreover, the limit of b/r when r tends to infinity will be divergent. Besides, we discard the condition m + 1 = 0 ⇒ m = −1 because it does not lead to the Minkowski space-time in the mentioned limit, because the coefficient in front of dr 2 term in the line element will not be 1. In that case one needs to redefine the radial coordinate r but the solid angle is altered. As shown Fig. 2 for the chosen m the shape function b(r ) (blue line) remains finite everywhere. Moreover, the flareout condition is also satisfied, it is depicted by the green line, where b (r ) is always less than the straight line r = 1 at the neighborhood of the throat r 0 , and the blue line is representing the WH throat, as can be seen this curve cuts the radial axis just at r = 1 the throat of the WH geometry.

Thermodynamic variables and energy conditions
It is well known that the matter distribution surrounding the space-time can be composed of a large number of material fields. While one has an idea of how the matter content described by the energy-momentum tensor behaves, obtaining an exact description of it can be a very complex task. As far as this is concerned, the so-called energy conditions constitute a way of testing the positivity of energy density ρ, as is required by a real and well behaved matter content.
Despite energy conditions are not part of fundamental physics, they are useful to characterize the type of fluid distribution one deals with [84]. As it is well known in the framework of General Relativity, the so called null energy condition (NEC): ρ + p r ≥ 0 and ρ + p ⊥ ≥ 0 (for anisotropic matter distributions) is violated in order to have a spherically and static traversable WH space-time. Since, in that case the radial pressure p r is negative and greater in magnitude than the energy density ρ. It may also occur that the energy density ρ becomes not strictly positive at all points. Besides, in General Relativity the energy conditions establish that gravity should be attractive and they can be derived from the Raychaudhuri equation which is given by where θ , σ μν , ω μν are respectively the expansion, shear and rotation associated with the congruence of time-like geodesics specified by the vector field u ν . So, the Raychaudhuri equation for the congruence of null geodesics (defined by the null vector k μ ) is as follows So, to satisfy the attractive gravitational force condition one needs R μν u μ u ν ≥ 0 and R μν k μ k ν ≥ 0 [85]. It should be noted that the Raychaudhari equations are purely geometric description. Since Einstein fields equations relate the geometry of the space-time R μν with the matter content T μν , Raychaudhuri equations can be combined with the General Relativity ones to impose certain conditions on the energymomentum tensor for its physical viability. Hence, in Einstein theory one requires R μν k μ k ν ≥ 0, which implies T μν k μ k ν ≥ 0 (NEC). Nevertheless, in the arena of modified gravity theories, such constraints are not straightforward.
Specifically, in f (R) gravity the geometric sector is too complicated due to the introduction of high power terms of the Ricci scalar R. Therefore, to deal with energy conditions in this case one takes the total or effective thermodynamic quantities, that is those formed by the combination of the normal matter sector described by T m μν and the curvature energymomentum tensor T c μν . It should be noted that violations of the energy conditions have sometimes been presented as only being produced by nonphysical stress energy tensors. Nevertheless, they can be violated in many cases, for example the minimally coupled scalar field and curvature-coupled scalar field theories. As was discussed in [81] the possibilities to obtain WH solutions in f (R)-gravity theory can be analyzed by performing a conformal map from the Jordan frame to the Einstein one and apply there the known results given in [79]. Nevertheless, there are possibilities of getting geometries connecting two far away flat or infinity regions that is, wormholes, driven by normal matter that is, matter satisfying the NEC [57]. However, the price to pay is the presence of ghost fields [80,81].
So, to check the energy conditions it is necessary to know the behavior of the principal thermodynamic variables that characterize the fluid distribution, namely the density ρ, the radial p r and tangential p ⊥ pressures. These physical quantities are given by Eqs. (9)- (11), where the f (R(r )) function for each model are expressed by Eqs. (30) and (33).
For m = 1 we have and for m = 0 one arrives to In both cases the tangential pressure is determined by the condition p ⊥ = −ρ. As can be seen the behavior of the principal thermodynamic variables depends on the integration constants F 0 and f 0 . The Fig. 3 displays the behavior of the main physical quantities for both models. As can be seen in the upper left panel for m = 0 the density ρ is positive defined everywhere for all r ∈ [r 0 , +∞), while the upper right panel shows the trend of the radial p r and tangential p ⊥ pressures, being both quantities negative everywhere. As it is well known a negative radial pressure causes a repulsive gravitational force maintaining the throat of the WH open. In comparing the density and radial pressure in magnitude it is observed that the former is greater than the second one, namely |ρ| > | p r |. This fact as we will see soon ensure the satisfaction of the null and weak energy conditions at least at the throat of the WH and its neighborhood. On the other hand for the case m = 1 the situation concerning the density is similar to the previous case, that is, positive defined everywhere. However, the radial pressure is positive at the throat but takes negative values away of it. Now, we proceed to analyze the behavior of the NEC and WEC. So, combining Eqs. (5) and (9) we get Next, adding Eqs. (9) and (40) we obtain The above equation is just the NEC. As can be seen this expression is independent of f (R(r )).
It is clear from the previous expressions (42) and (43) that ρ+ p r > 0 for F 0 < 0. In this case to contrast the implications of having a WH solution that satisfies NEC and another that violates it, we consider F 0 = −1 for both cases ı.e, m = 1 and m = 0. It is worth mentioning that to satisfy the positiveness of ρ + p r in the case m = 0 there are several more restrictions than just considering F 0 < 0. Due to the r > r 0 for all r ∈ [r 0 , +∞) the denominator of Eq. (43) is always negative, but this negative sign is cancel out by the global factor β taking into account that β < 0. Nevertheless, on the numerator appears the term (r 0 − β), thus to assure a positive ρ + p r quantity r 0 must be greater than β in modulus. Therefore, our choice β = −r 0 /2 satisfies the above requirement. Finally, the term −4r 2 dominates over the rest part, then the previous choice on F 0 is necessary to satisfy the mentioned condition.  Figure 4 shows the trend of the NEC for both cases. The left panel corresponds to m = 1, as illustrated the NEC is positive for all r belonging to [r 0 , +∞). Moreover, at large distance the NEC is saturated ı.e, ρ + p r = 0. On the other hand the right panel displays a different situation for the m = 0 case. In this opportunity the NEC is satisfied at the throat of the WH (see Eq. (43)) and its neighborhood, after that takes negative values, this means violation of the NEC. However, far away from the throat is saturated. Furthermore, as was pointed out earlier, in both cases ρ > 0, hence the WEC (ρ ≥ 0 & ρ + p r ≥ 0) is also satisfied everywhere for m = 1 and partially for m = 0, because the NEC is violated in some regions. The fulfilment of NEC and WEC in the model m = 1 implies that the dominant energy condition (DEC) is also satisfied, that is, ρ ≥ |p r |. In considering the NEC, WEC and DEC in the tangential direction all are satisfied for m = 1 and m = 0, while the strong energy condition (SEC) ı.e, ρ + p r + 2 p ⊥ ≥ 0 is violated for each model due to the condition p ⊥ = −ρ.

Ghost fields
The existence of WH regions, at least from the theoretical point of view in the context of f (R)-gravity theory should be supplemented with an additional analysis. As was pointed out by Bronnikov and Starobinsky [80,81] a WH space-time with matter fields satisfying the NEC cannot exist if This implies that the existence of such solutions violates the NEC (like in General Relativity). The first statement of (44) assures that ghost fields (gravitons with negative kinetic energy) are absent in the solution. Hence, if F(R) < 0 then the theory propagates ghost fields (for a detailed discussion about this subject in modify gravity theories see [58] and Appendix A). The second statement of (44) is refereed to the stability of the WH throat. 2 So, to have a better understanding of how (44) works we will re-express Eqs. (21) and (25) as functions of the Ricci scalar and its value at the throat. So, from the expressions (29) and (32) where in order to avoid complex numbers R 0 < 0. For m = 0 we have So, using Eqs. (28), (30), (31), (33), (45) and (46) one arrives to for m = 1 and for m = 0. As it is depicted in Fig. 5 (right upper panel) in the m = 1 case F(R) is negative for all R. This means that the first statement of (44) is not met. A negative F(R) implies the present of ghost fields in the solution. However, as was stated in [80,81] it is not possible to obtain a free ghost field WH solution satisfying NEC. Of, course as pointed out earlier, the case m = 1 satisfies the NEC but contains states with negative kinetic energy. The second requirement of (44) is satisfied for this model as shows the left panel of Fig.  6. The case m = 0 entails a more interesting situation. As it is appreciated from expression (50) the F(R) function is everywhere positive and zero at R = R 0 = 1 (or equivalently at r = r 0 , since there is a one to one relation between r and R). This implies that at the throat the gravitational forces become infinite. Moreover, from the right panel in Fig. 6 is evident that d F/d R = 0 at R = R 0 = 1. This particular model is not within the scope of the analysis given in [80,81], because in these works only the case F(R) = 0 and d F/d R = 0 is studied. In general, for the case m = 0 the expression of (51) as can be seen the numerator of the above expression vanishes for every R = R 0 . Besides, this solution constrains ghost fields as Fig. 5 illustrates (right panel in the lower row) and partially violates the NEC as can be seen from the left panel in Fig. 4. In this case, ghosts appear because the NEC is satisfied in some regions. At this point, it is worth mentioning that the general solution of (20) with m = 0 and without any specification about the values of γ and r 0 is given by so regardless of which values are chosen for the parameters γ and r 0 this particular model always leads to F(R(r )) = 0 at the throat, or using (50) in terms of the scalar curvature R, one obtains at R 0 : F(R) = 0 and d F(R)/d R = 0.

Cosmological properties
In this section we discuss some cosmological aspects associated to f (R)-gravity. As it is well known this modify gravity theory emerges as an alternative for a unified description of the early-time inflation with late-time cosmic acceleration, without adding unknown forms of dark components ı.e, dark energy and dark matter [59]. Although several models of gravity f (R) have been proposed to unify the aforementioned problems [60][61][62], there are models that only address either early or late stage of the Universe. Considering, the latter a wide range of works available in the literature have faced the accelerated expansion of the Universe to give an explanation to the existence of dark energy from the perspective of f (R)-gravity [63][64][65][66][67][68][69][70][71][72][73][74][75][76].

Early and late time phases
To tackle the dark energy issue from the arena of f (R)gravity there some general requirement that any viable f (R) model should satisfy [77] 1. F(R) > 0 for R ≥R (R > 0), whereR is the Ricci scalar at the present epoch.

d F(R)
d R > 0 for R ≥R. This is required for consistency with local gravity tests, for the presence of the matterdominated epoch and for the stability of cosmological perturbations. 3. f (R) → R for R >>R. This is required for consistency with local gravity tests and for the presence of the matterdominated epoch.
This is required for the stability of the late-time de Sitter point.
To check the feasibility and viability of our model, we have analyzed in details the above criteria for the case given by m = 0 (the case m = 1 is not analyzed due to the general F(R) function without any specification about β and r 0 cannot be integrated to obtain the f (R) function). So, from expression (52) one has where c is a constant parameter replacing r 0 . Integration of (53) with respect to R leads to Next, from (53) the condition 1 is fulfilled if for all R ≥R > 0 iff F 0 > 0 and cβ < 0 where β > 0 and c < 0, thus F(R) > 0. To check the second condition from Eq. (53) one gets As before, this condition is satisfied for all β > 0 and c < 0. Next, when R >>R the special function tends to a real positive number, namely K . Then from (54) we have Now where at large scalar curvature R the constants f 0 can be neglected, or it can be interpreted as f 0 = −2 , thus one ends with Einstein theory including the cosmological constant term. Then condition 3 is satisfied. As can be seen from (57) in the limit R >>R Einstein's gravity is recovered as expected for viable and realistic f (R) models going from a matter-dominated epoch to a dark energy Universe [86]. The last condition entails a special kind of solutions within the f (R)-gravity background, those are the so-called de Sitter point class of solutions. This class of solutions satisfy the following constraint obtained from the trace of the equations of motion (5), which is trivial in the Einstein theory but gives precious dynamical information in the modified gravitational models [62]. So, this condition reads This constraint reveals the existence of maximally symmetric vacuum solutions in the theory [78]. So, if a covariantly constant scalar curvature R ≡R = constant satisfies (58) given any f (R)-gravity model, then the theory contains the General Relativity de Sitter solution with constant curvaturẽ R. Therefore, by replacing in Eq. (58) one arrives to It is obvious that the previous equation has always a solution which depends on three parameters, namely {c, β,R}. Next, the stability condition 4: 0 < R F d F d R < 1 is always satisfied. In the present case this condition reads Indeed, from (60) the denominator is always greater than the numerator, hence the expression is always positive and less than 1 for every c and β. It is clear from (54) that the model reproduces the dominated matter epoch expressed by (57). As it is well-known, this stage corresponds to early times or equivalently large curvature [86]. On the other hand, the late time stage (small curvature) described by (54) can be stated as follows 3 As can be seen at late times the asymptotic behavior of the present model leads to a model with positive powers of curvature [86]. At this point, it is worth mentioning that the present model does not unify the inflationary stage with the accelerated expansion. Notwithstanding, without loss of generality Eq. (53) can be generalized as follows being δ the delta Dirac function and R 0 and R I (R I >> R 0 ) are representing transition scalar curvatures (for further details see [62]). The above extension is completely plausible, since the functional given by (62) is valid for all R = R 0 and R = R I , being δ(R − R 0 ) and δ(R − R I ) vanishing. So, integrating the above expression one gets where is the Heaviside step function. As stated in [59,62] any viable f (R) model unifying inflationary phase with accelerated stage must satisfy the following requirements 3 For further details see the Appendix B f (0) = 0 and lim withR >> R I >> R 0 . Of course when R = 0 the special function tends to 1 and (63) provides f (0) = 0, ensuring the disappearance of the cosmological constant in the limit of flat space-time [62]. On the other hand, in the limit R →R Eq. (63) yields to −2 ( 0 + I ), where I is the inflation cosmological constant andR being the corresponding transition large scalar curvature. Despite the coupling between our model and the two step model provides a well posed scenario to unify the early and late time stages of the Universe, the introduction of Heaviside and Dirac distributions makes it a simply toy model [62]. However it shows a natural way to extend our model to address this point.

Late time accelerated phase in the flat FLRW space-time
To explore the late time accelerated phase of our Universe, we shall consider that the geometric background is described by a flat Friedmann-Lemaitre-Robertson-Walker (FLRW) metric [59,83], with line element given by, with a(t) denoting as usual the scale factor. Next, taking into account a perfect fluid matter distribution the field equation (4) adopt the following form where dots represent differentiation with respect to the cosmic time t and being H the Hubble rate defined as H (t) ≡ a(t)/a(t). As usual we shall assume that ρ and p are related by the following equation of state being ω a constant parameter, namely the equation of state parameter. So, plugging the model expressed by (61) 4 with the previous ingredients into the field equation (67) one gets where D ≡ 2F 0 (β − c) /5 (−2βc) 1/4 . Besides, for the sake of simplicity we have fixed f 0 = 0. The FLRW line element given by (65) leads to So, Eq. (70) can be written as The above equation admits the following solution in inverse powers of t for the Hubble rate H (t) as follows where I and J are given by then the scale factor a(t) is being .
It should be noted that the first terms of (73) and (75) yield to the well-known results provided by General Relativity [87], hence next terms can be seen as corrections introduced by f (R)-gravity. Now, from the field equations given above (67)-(68) the acceleration equation reads being ρ eff ≡ ρ + ρ R and p eff ≡ p + p R , where and From Eqs. (77)-(79), the acceleration condition, for a dust dominated model ( p = 0 ⇒ ω = 0 ), leads to [88] the above condition can be cast in terms of the curvature contributions as follows In Fig. 7 are displayed the evolution of the density ρ(t) (left panel) against the cosmic time t and the trend of the scale factor a(t) (right panel). On the other hand Fig. 8 shows the behavior of the effective equation of state parameter ω eff . As can be seen the condition (80) is satisfied. Then the model describe a Universe in an accelerated phase as desired [88].

Concluding remarks
In this work a spherically and static wormhole family with vanishing red-shift function was studied in the arena of f (R)gravity theory. The field equations were solved by imposing a suitable shape function satisfying all the requirements to describe a traversable wormhole and asymptotically flat space-time at infinity and an equation of state, specifically p ⊥ = −ρ. The imposition of this relation leads to a first order differential equation for F(R(r )) function. Once the F(R(r )) function is obtained it can be expressed in the form F(R) by inverting the radial coordinate r from the Ricci scalar expression.
As it is well known in the context of General Relativity a wormhole space-time is only possible if the NEC is not satisfied at the throat and its neighborhood [79]. The same argument is translated into the f (R)-gravity theory in the absence of ghost fields. Moreover, under some special con-ditions wormhole geometries respecting NEC are possible in the background of f (R)-gravity. Regarding the present study we can summarize the main results for each model as follows: m = 0 1. For the considered space parameter, that is {r 0 , β, F 0 } = {1, − 1 2 , −1}, it can be seen from Fig. 3 (upper panels) that the density ρ is positive defined everywhere, contrarily to what happens with the radial pressure p r which is negative for all r ∈ [r 0 , +∞), so the WH geometry is supported by a repulsive gravitational force introduced by this negative radial pressure. However, it should be noted that there are some points where | p r | > |ρ|, then the NEC and WEC are locally violated as it is illustrated in Fig. 4 (left panel). 2. As Fig. 5 shows (right panel in the lower row ) F(R) is strictly negative everywhere. One should take into account that as stated in [80,81] in the framework of f (R)-gravity theory it is not possible to obtain a wormhole solution respecting energy conditions (not supported by exotic matter) without the presence of ghost fields. m = 1 1. In this case the space parameter was fixed to {r 0 , β, F 0 } = {1, 1, −1} and in distinction with the previous case the the WEC and NEC are satisfied everywhere. This is so because the density ρ is positive defined at all points r ∈ [r 0 , +∞) and greater in magnitude than the radial pressure p r , which is negative everywhere. The fact of fulfilling the energy conditions (WEC and NEC) makes the f (R)-gravity arena an interesting and promising field to study the main properties of wormhole structures without any exotic matter distribution as occurs in General Relativity. 2. Although this case satisfies the WEC and NEC everywhere, as stated before then the theory can not have absence of ghosts. In fact, if the wormhole structure is threading at its throat by normal matter distribution respecting the mentioned energy conditions, then the solution must have ghost fields. The main problem with the ghost fields is that the f (R)-gravity model associated with the wormhole geometry can not describe or explain cosmological issues, since this requires a positive defined energy.
To close the wormhole scenario, we remark that it is possible to build viable wormholes space-time in the f (R)-gravity context respecting all the general requirements to be a traversable region connecting two asymptotically Minkowskian or infinite spaces, either violating or satisfying NEC. However, in the latter case, the conformal map [80,81] must be ill defined at some points since in the Einstein frame there are not wormholes satisfying NEC [79].
An important property of one of our models is the cosmological viability to address the late-time accelerated epoch of our Universe. This point was analyzed by checking that the f (R)-gravity models satisfy the general requirements given in [77]. It is worth mentioning that to study the cosmological implications described by the resulting models with m = 0 and m = 1, the space parameter {r 0 , β, F 0 } should be reset, because the previous assignation was to describe an asymptotically and traversable wormhole geometry satisfying the energy conditions (partially or totally). For the case m = 1 it is not possible in general to analyze its cosmological viability, because the general expression given by (21) is not tractable mathematically. in distinction, the case m = 0, as was shown in Sect. 6 has good properties to address the late-time acceleration era of the Universe. In this concern, we have checked the general criteria [77] to confirm these good properties. Among these features, the model f (R) corresponding to m = 0 respects: (i) consistency with local gravity tests, (ii) stability under cosmological perturbations, (iii) ghosts free and (iv) stability of the de Sitter point. In considering the last point, it is remarkably that this model contains vacuum solutions ı.e, solutions with constant curvature R 0 and null energymomentum tensor T μν = 0, which also include the de Sitter's solution to the vacuum Einstein field equations with cosmological constant [78]. Furthermore, to check the accelerated phase of our Universe at late time for the present model, we have explored the behaviour of the scale factor a(t) and effective equation of state parameter ω eff on the background of a flat Friedmann-Lemaitre-Robertson-Walker space-time. As it is corroborated in Fig. 8 the equation of state parameter satisfies the condition ω eff < −1/3 as desired for an accelerated expanding Universe [88]. At this point it should be noted that the present model is in accordance with some wellknown recognize models previously reported [59,62,86] in the arena of f (R)-gravity theory to deal with cosmological open issues, facing both: the early and late time epochs of our evolving Universe. In considering the former, our model does not include early time phase (inflationary phase), since the model is not capable to produces the minimal ingredients to describe this stage of the Universe [59]. Nevertheless, as was discussed in Sect. 6, the model can be consistently coupled with a two step model [62] to unify both stages, although the output is only a toy model it shows a natural way to extend our model to address this point. For these reasons we only explore in some details the late time accelerated phase, which in the limit of high curvature reproduces Einstein gravity theory. Finally, we want to highlight that the model obtained in this work is able to describe for certain values of the parameters {m, β, r 0 , F 0 } a wormhole structure without ghost, partially satisfying the NEC, while for other values of aforementioned parameters, describes a consistent and viable cosmological model to explain late stages of the cosmic evolution of the Universe in a flat FRLW background. new metric, in the Einstein frame, is 2L E = Sgn( ) R +g μν ∂ μ ϕ∂ ν ϕ + W (ϕ), where The condition of absence of ghosts is If F(R) < 0 the graviton is a ghost field. There is no way to avoid this behavior. Besides, scalar-tensor theories even when the graviton has positive energy may have ghosts on the scalar sector. In [58,59] a procedure to avoid them has been presented.

Appendix B: Series representation of Appell functions
The AppellF 1 series is defined for |x| < 1 and |y| < 1 by the double series where (q) i is the Pochhammer's symbol. So, when x → 0 and y → 0 the series expansion is given by