Wormhole solutions with a polynomial equation-of-state and minimal violation of the null energy condition

This paper discusses wormholes supported by general equation-of-state , resulting in a significant combination of the linear equation-of-state and some other models. Wormhole with a quadratic equation-of-state is studied as a particular example. It is shown that the violation of null energy condition is restricted to some regions in the vicinity of the throat. The combination of barotropic and polytropic equation-of-state has been studied. We consider fluid near the wormhole throat in an exotic regime which at some r=r1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r=r_{1}$$\end{document}, the exotic regime is connected to a distribution of asymptotically dark energy regime with -1<ω<-1/3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-1<\omega <-1/3$$\end{document}. We have presented wormhole solutions with small amount of exotic matter. We have shown that using different forms of equation-of-state has a considerable effect on the minimizing violation of the null energy condition. The effect of many parameters such as redshift as detected by a distant observer and energy density at the throat on the r1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r_1$$\end{document} is investigated. The solutions are asymptotically flat and compatible with presently available observational data at the large cosmic scale.


Introduction
The significant development in wormhole physics has been started by Morris and Thorne [1]. Wormhole is an exact solution of the Einstein field equations which can be used as a hypothetical shortcut between points in a universe or between two different universes. Wormholes are not famous as black holes, but the possibility that wormholes and black holes are in fact very similar object has been presented [2][3][4][5]. Wormholes have not been observed experimentally but there is not any observational reason to rule out this theory. Wormhole is a good candidate which would provide us with a practically unlimited possibility for interstellar travel. In wormhole study, a description of the matter content of the a e-mail: fparsaei@gmail.com (corresponding author) b e-mail: rastgoo@sirjantech.ac.ir wormhole is essential to present wormhole physics. Morris-Thorne wormholes violate classical energy conditions [6]. Null energy condition (NEC) specified by T μν k μ k ν ≥ 0, in which k μ is any null vector and T μν stress-energy tensor. The matter which violates the NEC is called exotic. Exotic matter is the main ingredient in wormhole theory. It was shown that violation of the NEC is a generic and universal feature of the traversable wormhole in General Relativity (GR) [7,8]. Of course, wormholes with NUT parameter [9][10][11] may be considered as traversable wormhole without exotic matter but, one should keep in mind that cosmological constant is exotic in this realm. Some researchers try to solve the problem of exotic matter by investigating wormhole in modified theories of gravity. Wormhole in modified gravitational theories such as Brans-Dicke [12][13][14], f (R) gravity [13,[15][16][17][18][19], curvature matter coupling [20,21], and brane-world [22][23][24][25][26][27] has been studied. For a review on traversable wormhole solutions in modified theories of gravity, see [28] and references therein. In most of these theories, an effective stress-energy tensor which contains the higher-order curvature has appeared in the right side of the Einstein field equations. So the effective stress-energy tensor violates NEC instead of ordinary matter. Investigation of wormholes in higher-dimensional spacetime was done [29,30]. Wormholes in the framework of Lovelock [31][32][33] which is considered as the most theory of gravitation in n dimensions have been studied. Arkani-Hamed et al. have studied Euclidean wormholes in string theory [34].
Application of quantum inequalities imposed a bound on exotic energy densities [35]. The quantum effects allow local violation of energy conditions, but this bound implies that any inertial observer in flat spacetime cannot see an arbitrarily large negative energy density which lasts for an arbitrarily long period of time [36]. Kuhfittig has shown that a wormhole supported by only small amounts of the exotic matter really can be traversable [36]. Wormholes with small amount of violation of the energy conditions in the framework of GR have been studied in the literature. Somebody have used the cut and paste method to minimize the violation of NEC [38][39][40][41][42]. In the cut and paste method, the single manifold M is the combination of two spherically symmetric spacetime M int and M ext which the interior space is related to the wormhole and the exterior one is usually the vacuum space. The junction conditions originated from Israel junction condition, should be considered in this surgery. Since the wormhole geometry is confined to a local spacetime in the cut and paste method, the violation of NEC is due to this part and can be arbitrarily small. We believe that thin shells are a mathematical abstraction, for physical reasons, it is better to minimize the usage of the cut and paste method.
Theoretical physicists have used the different equationof-state (EoS) to investigate cosmological as well as astrophysical phenomena. The EoS of fluid supporting wormhole geometry is an essential equation in studying wormhole theory. Recent observations indicate that the Universe is spatially flat with low matter density and expanding with acceleration [43][44][45][46]. Fluid with an EoS, p = ωρ, and positive energy density is the most considerable one to explain the evolution of the cosmos. The regime with −1 < ω ≤ 0 is called dark energy and the special case with ω < −1 is called phantom. Accelerated expansion of the Universe is reachable for fluid with ω ≤ − 1 3 . The phantom fluid is found to be compatible with most of classical test of cosmology such as CMB (cosmic microwave background), anisotropy, and mass power spectrum. Although the exotic nature of phantom energy leads to some unusual cosmological consequences such as big rip, phantom fluid is a good candidate to support wormhole solutions. It violates the NEC and is speculated to be a possible driving late time cosmic acceleration. Many asymptotically flat [47][48][49][50] or not asymptotically flat [49,[51][52][53][54] wormhole solutions with phantom fluid source have been presented. Wormholes with other forms of fluid, like Chaplygin EoS [55][56][57][58], polytropic EoS [59], and modified Chaplygin EoS [60] have been studied in the literature. Jamill et al. have studied wormhole supported polytropic phantom energy which is a generalization of phantom energy and in some cases Chaplygin-gas models [61]. In all of these papers, the violation of the NEC is inevitable in the whole of wormhole spacetime.
Some authors have studied the wormhole with a variable EoS parameter in which ω = p ρ is a function of radial coordinate [62][63][64][65][66][67]. In some of these studies, the violation of energy condition is inevitable in the whole of spacetime [62]. In [63], the cut and paste is used to find wormhole solutions of the finite size which minimize the violation of energy conditions. Lopez et al. have assumed an EoS in which the sum of the energy density and radial pressure is proportional to a constant with a value smaller than that of the inverse area characterizing the system [64]. They have found solutions which are not asymptotically flat. Azreg-Aïnou by consid-ering two barotropic equations of state for lateral and radial pressures has presented wormholes with no gluing effects [65]. In another method, Remo and Lobo have proposed solutions in which phantom fluid is concentrated in the neighborhood of the throat to ensure the flaring out condition [66]. Some of their solutions are asymptotically flat. They have investigated the possibility that these phantom wormholes be sustained by their own quantum fluctuations. In [67], a general formalism to find asymptotically flat wormholes with variable EoS parameter has been presented. Also, the physical difference between the cut and paste method and the intrinsically asymptotically flat wormholes which minimize the NEC has been discussed. Note that, it is not a necessary condition that the EoS of wormhole or any astrophysical object should be the same as the EoS of the Universe. The only condition is that the EoS of wormhole should be fitted by EoS of the Universe in the large cosmological scale. This motivation helps us to find wormhole solutions with EoS that is asymptotically linear but have additional terms. So this work is focussed on the exploration of wormhole solutions using the different form of EoS. We will seek solutions which restrict the violation of the NEC to some regions in the vicinity of the throat. The theory is GR but the matter content of wormhole is investigated in more details. Actually, we considered an EoS which has a linear term and some extra terms. This method can be used to sustain wormhole solutions with minimum violation of energy conditions. We study some new wormhole solutions which are asymptotically flat. We will present solutions by considering a special shape function.
The paper is organized as follows: In the next section, the basic structure of wormhole theory is presented. In Sect 3, by considering a power-law shape function, wormhole with quadratic EoS is studied and solutions with minimum violation of the NEC are presented. Wormholes with a combination of barotropic and polytropic EoS are investigated in Sect. 4. Discussions and Concluding remarks are presented in the last section.

The general model of wormhole theory
Usually, the metric of a static and spherically symmetric wormhole is given by [1] where d 2 = (dθ 2 + sin 2 θ dφ 2 ). At a minimum radial coordinate, r 0 , with b(r 0 ) = r 0 the wormhole connects two different worlds or two distant parts of the same universe. Here b(r ) is called the shape or form function and r o is the throat of the wormhole. The function φ(r ) is called redshift function because the gravitational redshift as measured by a distant observer is The redshift function should be finite everywhere to avoid the existence of the horizon. The conditions and are necessary to support a traversable wormhole. Note that the prime denotes the derivative d dr . In this paper, we investigate asymptotically flat wormhole solutions so the condition is also imposed. Considering an anisotropic fluid in the form T μ ν = diag(−ρ, p, p t , p t ) and using the Einstein field equations, the following distribution of matter (with 8π G = c = 1) are obtained, where ρ(r ) is the energy density, p(r ) is the radial pressure and p t (r ) is the lateral pressure. Equation (8) can be considered as a consequence of the conservation of the stress-energy tensor, T μν ;μ = 0. Although the wormhole solutions are anisotropic in general form, we can use an EoS for radial pressure and find lateral pressure through the Einstein equations. This is motivated by the discussion of the inhomogeneities that may arise due to gravitational instabilities. Indeed, as the dark energy EoS represents a spatially homogeneous cosmic fluid and is assumed not to cluster. It is also possible that inhomogeneities may arise due to gravitational instabilities. Thus anisotropic wormholes may have originated from density fluctuations in the cosmological background, resulting in nucleation through the respective density perturbations. One can also deduce the possibility that these structures are sustained by their own quantum fluctuations [66]. Because of the aforementioned reason, the pressure in the EoS may be regarded as a radial pressure, and the tangential pressure may be deduced through the Einstein field equations. This method has been used extensively in the literature, i.e. phantom wormholes [47][48][49][50][51][52][53][54]. Also Sushkov and Kim have studied a time-dependent solution describing a spherically symmetric wormhole in a cosmological setting with a ghost scalar field [68]. More specifically, they have shown that the radial pressure is negative throughout the spacetime, and for large values of the radial coordinate, equals to the lateral pressure, which demonstrates that the ghost scalar field behaves essentially as dark energy. There are many algorithms to construct wormhole geometry theoretically. Assuming desired forms of the redshift and shape functions and then the corresponding matter field is determined by the Einstein field equations. In another method, the redshift or shape functions are obtained as solutions by solving the Einstein field equations with prescribed matter field configuration. The energy density and pressure in the first method usually are not realistic. In the second algorithm, one may consider an EoS which has suitable properties to describe the physical matter. The second algorithm has been used by many authors in the literature to find a wormhole solution with different forms of EoS [47][48][49][50][51][52][53][54][55][56][57][58][59][60][61][62][63]. In the present article, we used the same algorithm to find asymptotically flat wormhole exact solutions with some other forms of EoS.
Generally, an EoS is added to set of equations and then some strategies have been used to solve the equations and find unknown functions. The EoS is very important because it explains the physical fluid which is essential to sustain the wormhole. Commonly, the essential fluid need to construct the wormhole geometry is not isotropic. We can use the radial pressure in the EoS, p = f (ρ), which was first presented in the study of phantom wormhole solutions [51]. In recent studies, fluid with a linear EoS manifests itself as the source causing a rapid accelerated expansion of the Universe on a large cosmic scale and EoS is the most famous in studying cosmos. As it was mentioned the accelerated expansion of the Universe is reachable for fluid with ω ≤ −1/3. Fluid with −1 < ω ≤ −1/3 and ρ > 0 satisfies the NEC which can cause accelerated expansion of the Universe. Phantom fluid violates the NEC every where in the spacetime but a linear EoS with ω > −1 and an extra term may violate the NEC only in some region of the spacetime. So we try to find some solutions with barotropic EoS and a mixed energy density. We consider an EoS as follows: At large distance from the throat, the physics of wormhole should be compatible with the cosmos. So the condition leads to an asymptotically linear EoS on the large scale. Wormhole models building in cosmology based on two main ingredients: a theory of gravity and a description of the matter content of the Universe. In this paper, the gravity sector of the theory is not completely fixed. There are some free parameters. The matter sector is represented in the field equations by the energy-momentum tensor and some EoS rather than conventional ones have been considered. The conventional EoSs like linear, polytropic, Chaplygin, power-law and so on are models to describe cosmology. The wormhole theory has been explored in those models to achieve a good understanding of our Universe. Since there are not any observational data about wormhole, it seems that the EoS of wormhole hasn't a specific form. So, any EoS that has the necessary conditions in the local view and is compatible with cosmological known EoS in global view, can be considered as a candidate. It is clear that wormhole geometries admit an EoS in the general form p = f (ρ). Mathematically, this EoS can be expressed in a Taylor series so any general EoS of the wormhole can be expressed in the form of Eq. (10). The EoS (10) helps us to investigate a large class of wormhole solutions which violates the NEC only in some regions near the wormhole throat. Phantom wormholes violate the NEC in the whole of spacetime but the additional term g(ρ) in Eq. (10) provides a way to construct wormhole solutions with the minimum violation of energy condition and considering ω > −1 instead of phantom EoS in large scale. In asymptotically flat spacetime at large radial coordinate, ρ(r ) and p(r ) should tend to zero. So to have compatibility with cosmos in large scale, we have In the following, we consider some strategy to find exact wormhole solutions with different forms of g(ρ).

Wormhole with a quadratic equation of state
Ananda and Bruni have discussed cosmological dynamics and dark energy with a quadratic model of EoS [69,70]. They have shown this model evolves from a phantom phase asymptotically approaching a de Sitter phase instead of evolving to a big rip. In [69,70] and references therein, it has been addressed the motivations behind considering a quadratic EoS. For example, the loop quantum gravity corrections which result in a modified Friedmann equation provides a quadratic EoS. The modification is appearing as a negative term which is quadratic in the energy density. Also the studies of k-essence fields as unified dark matter (UDM) models indicated that a fluid with a closed-form barotropic EoS is a good candidate to describe the general k-essence field. Now, we consider the same EoS in the following form It should be noted that the correct form of EoS is as follows but for the sake of simplicity, we set ρ c = 1. This EoS is suitable to describe a standard fluid at high pressure and high energy density to vanishing presumer at low energy density(lim ρ→0 p(ρ) = 0). Wormholes with a quadratic EoS has been studied by Rahaman et al. [72]. They have considered the shape function as a series and have found the form of coefficients of the series solution. They have studied the possibility of finding solutions for this particular shape function. They have not reached in a reasonable proof for the condition (4) in their research. They only discussed some general mathematical possibility. Let us study wormhole with a vanishing redshift function, putting φ = 0 in Eq. (7), leads to So by considering (13), one can get Finding exact solution for this equation is very difficult. Thus we try to find solutions with nonconstant redshift function.
For the sake of simplicity, we set r 0 = 1 in the recent part of this paper. Several authors have investigated wormhole by considering power-law shape function [49][50][51][52][53][54][55][56][57][58][59]61]. We choose a special shape function in the following form [50] b The mass function related to this shape function is The condition 1 > α should be imposed to satisfy asymptotically flat condition (5). It is clear that for α > 0 the mass is unbounded and for 0 > α mass is bounded. This shape function shows an energy density in the form which is a positive smooth function of r for A α > 0 . This energy density function has a maximum at the throat and tends to zero at large infinity distance from the throat. Equation (20) relates the constant A and α to physical parameter ρ 0 . Someone may believe that specifying a shape function lacks physical justification and motivation for the stress-energy tensor. But we can consider another alternative approach that has been used by many authors. As it was done in [50,51], one can consider a specific energy density profile instead of considering a known shape function. In this context, the energy density profile (19) can be considered known which leads to the shape function (17). We think that these two approaches are intrinsically the same. So we consider only the known shape function instead of considering a smooth energy density distribution. The presented algorithm could be used for other forms of energy density distribution. Because of the aforementioned properties, the defined shape function is a good candidate to find wormhole exact solutions. Now, we try to find solutions for φ(r ) by using the Einstein field equation. Using Eq. (7) leads to Putting p(r ) from (13) and (20) in this equation yields Calculating this integral for a general form of A and α is difficult, so we analyzed some special cases. As the first example, consider α = −1 which leads to with To avoid from horizon, we should set c 5 = c 6 = c 7 = 0 which leads to Equation (19) implies that for a positive energy density, A should be negative. Defining the function helps us to investigate the violation of the NEC. The function F(r, A) is depicted as the surface in Fig. 1, which is negative through the spacetime so the NEC is violated everywhere. Let us try to find solutions which satisfy the NEC in some regions of spacetimes. We consider A = 1 and α = 1/2. Equation (22) provides The condition of an event-horizon-free spacetime requires that φ(r ) be finite everywhere. Therefore the coefficient of ln(r ) should be equal to zero. This yields Using Eqs. (2) and (27), one can easily show that It relates the redshift as detected by a distant observer and ω 1 . Since −1 < ω < −1/3 is related to dark energy regime, it is easy to show that is acceptable in this regime. Now, we devote some words to the possibility of the violation of the NEC. We define the function The sign of this function indicates the violation of the NEC. The negative sign corresponds to exotic matter. We have plotted F(r, ω) against r in Fig. 2 for some different ω. It is clear that the violation of the NEC is restricted to some Here r 1 is the radial where the sign of F(r, ω) changes. We have plotted r 1 as a function of ω in Fig. 3. This figure shows that as ω increases r 1 decreases. Since ω < −1/3 leads to accelerated expansion of the Universe, the minimum of r 0 is given by It is also obvious that lim ω−→−1 r 1 (ω) −→ ∞. Although the violation of the NEC is resulted from the flaring out condition which appears directly in ρ + p, it is necessary to check the other NEC ρ + p t ≥ 0. So we have defined the function We have plotted F 1 (r, ω) versus r and ω in Fig. 4. Since  F 1 (r, ω) is positive through the entire range of r and ω, we can conclude that ρ + p t ≥ 0 holds everywhere. In general, we can conclude that the selected shape function (17) can be used to construct wormhole geometries with a quadratic EoS fluid in which the violation of the NEC is considerable only near the wormhole throat.  In the next sections, we will try to find solutions with a more general EoS. The class of solutions has been presented in this section seems to be a special case of the solutions in the next section.

Combination of barotropic and polytropic equation of state
Wormhole with a polytropic EoS has been studied in the literature, kjuy [59,61]. On the other side, many papers have discussed wormhole with a barotropic EoS [49][50][51][52][53][54]71]. Now, we consider wormhole with a combination of barotropic and polytropic EoS: An important feature of this EoS is the reduction of exotic matter in constructing wormhole theory. We use the method of the previous section to find exact wormhole solutions. If we let A = 1 in the shape function (17) and use Eq. (35), the form of integral (22) becomes as follows where Considering the derivative of redshift function finite at the throat (the term behind integral (21) should be finite at the throat) leads to The solutions of φ α m for a general α and m appear in the hypergeomtric function form, so we will solve some special cases. An interesting solution corresponds to α = 1 2 and For these values, we can verify that n = 3+s 5 and Finally, from Eqs. (36)(37)(38)(39)(40) for n = 3+s 5 and α = 1 2 one can get By taking into account Eqs. (2) and (41), one can verify that where H 5n−5 is the (5n − 5)-th harmonic number. Equation (42) shows the relation between z and n. To check the NEC, we use the function It is clear that F(n, r, ω) depends on n, ω and r . Since ρ > 0, one can deduce that Here r 1 is the radial where the sign of F(n, r, ω) changes. We have plotted r 1 (n, ω, α) as a function of ω and α for n = 3 in Fig. 5. This figure shows that the value of r 1 increases as α or ω decreases, so in the solutions with larger ω or α, the violation of the NEC is limited to a smaller space near the throat. We can use the same analyzes for the dependence of r (n, ω, α) on n and α or n and ω when the third parameter is fixed. Figures 6 and 7 describe these dependencies. Figure 6 indicates that r 1 is an ascending function as α or n decreases. Also, Fig. 7 explains that r 1 is an ascending function as ω or n decreases which shows that all of three parameters generally have the same behaviour. Observe that as n −→ ∞ then r 1 −→ 1 which means that for larger n the violation of the NEC is limited to a very thin radial in the vicinity of the throat.
For this class of solutions, investigating the condition ρ + p t ≥ 0, for a general n is too complicated. So, we will    Fig. 8 The plot depicts the function F 3 1 (r, ω) as a function of ω and r for. It is clear that F 3 1 (r, ω) is positive through the entire range of ω and r . This shows that condition ρ + p t ≥ 0 holds everywhere. See the text for details for each n then we can plot F n 1 (r, ω versus r and ω. As an example, we have plotted F 3 1 (r, ω) against r and ω in Fig. 8. This figure guarantees that ρ + p t ≥ 0 holds everywhere. Let us seek a special example for wormhole solutions violating the NEC only in the vicinity of the wormhole throat in detail. If we set α = 1 2 , ω = − 2 3 and n = 3 , leading to the line element The stress-energy tensor components are as follows p t (r ) = 1 36 r 13 (3r 21/2 + r 10 + r 19/2 + r 9 + r 17/2 +r 8 + r 15/2 + r 7 + r 13/2 + r 6 + 133r 11/2 −4r 5 − 4r 9/2 − 4r 4 − 4r 7/2 − 4r 3 − 4r 5/2 Then r 1 = 4 1/5 1.32 and the total amount of exotic matter can be measured by [74] I = 8π We have plotted ρ + p and ρ + p t as a function of r in Fig.  9 which implies that the NEC is violated only in the interval r 0 ≤ r < r 1 with the total amount of the NEC violation I −8π × 0.056 and z = e 3/4h 10 10.58. To summarize, Fig. 9 The plot depicts the function p + ρ (solid line) and p t + ρ (dotted line) as a function of r where n = 3, ω = −2/3 and r 1 1.32. It is clear that for r > r 1 the NEC is satisfied. See the tex for details the line element (46) is a special example of wormholes with an unbounded mass which has been constructed by a fluid in the form of Eq. (35). One can construct many different wormholes solutions by considering other forms of b(r ) and p(ρ) with the minimum violation of the NEC by carefully fine-tuning the parameters to find φ(r ).
The same analyse can be used for a general EoS in the form (51) where n i > 1. The physical properties of this EoS when n i is an integer has been investigated in [73]. The condition, avoiding horizon at the throat, gives It is easy to show that where As it was mentioned before, calculating φ α m i for general α and m i leads to hypergeometrical functions. So one can study the special cases by choosing the suitable values for α and m i . Again, we put α = 1 2 and m i = s 2 , s = 3, 4, 5, ... or Fig. 10 The plot depicts the function ρ − (ω 2 ) as a function of ω 2 . It is clear that ρ − is in the possible range 0 ≤ ρ ≤ 1 2 so ρ − is acceptable. See the tex for details n i = 3+s 5 , then For example, by using Eq. (55) for the case where the only non vanishing n i are n 1 = 2 and n 2 = 3 with then for ω = −α = − 1 2 , Eq. (52) yields This implies that we can not choose all of ω i in Eq. (52) freely. It is easy to show that p + ρ changes sign in We have plotted ρ − (ω 2 ) and ρ + (ω 2 ) as a function of ω 2 in Figs. 10 and 11 respectively. Since for α = 1 2 , the possible range for energy density is 0 ≤ ρ ≤ 1 2 , one can deduce from Fig. 11 The plot depicts the function ρ + (ω 2 ) as a function of ω 2 . It is clear that ρ + is not in the possible range 0 ≤ ρ ≤ 1 2 so ρ + is not acceptable. See the tex for details these figures that ρ + is not acceptable. Now, it is easy to show that r 1 (ω 2 ) has been plotted against ω 2 in Fig. 12. This plot shows that r 1 is an ascending function of ω 2 . It means that the smaller ω 2 will present a wormhole with smaller amount of the NEC violation. We analyze three cases ω 2 = −10 (case 1), ω 2 = −2 (case 2) and ω 2 = 0 (case 3) which are related to following EoS's respectively.
We have plotted p(ρ) as a function of ρ for three cases and compared with p = ωρ in Fig. 13. This figure shows that p 1 have more consistency with linear EoS. It should be noted from Fig. 12 that the NEC is violated for ω 2 = −10 (case 1) in a smaller region in comparison to the other cases. It is also of particular interest that p 3 is a quadratic EoS which is achieved by putting ω 2 = 0 in Eq. (56). So one should note that if EoS has more terms, the opportunity to fit the EoS with linear EoS is more accessible and the NEC violation is Fig. 12 The plot depicts the function r 1 (ω 2 ) as a function of ω 2 where ω = − 1 2 . It is evident that r 1 increases as ω 2 increases which shows that the NEC violation is restricted in a smaller region for smaller ω 2 . See the tex for details restricted to a smaller region. Also, EoS with few terms can be considered as a special case of a more general EoS. This result may be considered as the one of the main points of our work. It should be mentioned that the conditions (3) and (4) are necessary but not sufficient to check wether wormhole is traversable or not. We should checked the parameters where a is the radial acceleration that traveler feels and a t is the maximum tidal acceleration felt by traveler [1]. As an example, for metric (46) we have calculated the maximum value of |a(r )| in the range 1 ≤ r < ∞ which is equal to a max = 0.52. So, the condition is satisfied. Also, the maximum value of a t is ( a t ) max = 0.25. Therefor, the condition is satisfied. Note that we have considered 8π G = C = 1. In [1], it was mentioned that are the necessary conditions in a journey. Checking these conditions analytically for our solutions is so difficult but one can consider some value for l 1 and l 2 then try to check these conditions numerically. As an example, if we set l 1 = −l 2 = −3 × 10 8 for metric (46) then we can find τ 1 19.89 and t 1 20.002 . So total time to travel conditions are satisfied. During this calculation, we consider γ = 1/ √ 0.99. Reader can used the same procedure for other solutions.

Discussion and concluding remarks
Wormholes have attracted a lot of attention in the literature and a large number of models have been proposed in this realm. Many of these studies are presented to answer to the basic question, how to avoid or at least minimize the requirement of the exotic NEC violating matter? Minimizing the usage of exotic matter for the physical viability of wormhole has received considerable attention. We have shown the fact that the choice of EoS for the description of matter presenting the wormhole has great relevance in the violation of the NEC in the wormhole static solutions. In the following discussion, we will see how our method is interpreted as a solution for wormhole and how the definition of EoS affects the violation of energy conditions. We also discuss the priority of this method than the cut and paste method.
Actually, we can divide the geometry of the wormhole into two regimes; high energy density regime and low energy density regime. If we consider ρ as a steady function of r with a minimum and a maximum, then low regime is related to the regime where ρ tends to a minimum when the high regime is relevance to the vicinity of the maximum. This motivation helps us to consider an EoS with two terms. The first term is the dominant term in the low regime. The violation of the NEC can be caused by the combination of these two terms. On the other hand, it supports the possibility of confining the exotic matter in some finite regions of the spacetime. The shape function (17) present a smooth energy density function which has a maximum at the throat and tends to zero at large distance from the throat. It seems that this is a suitable choice to present the desired energy density. It should be noted that this algorithm can be used for other forms of shape functions or energy density. Any physical energy density which has a maximum near the throat and tend to zero at large distance may be considered as a candidate. The mathematical difficulty of solving equations to find redshift function is different for different choices of shape function.
Generally, in our method EoS is as follows The former function, L, in the recent equation gives information about the behavior of the fluid supporting wormhole far from the throat while combination of L and H specifies the behavior of the fluid in the vicinity of the wormhole throat. It should be noted that the portion of the H is more than L near the wormhole throat. The parameter (r ) = H L can be used as a measure to quantify the portion of each term in different regions. It is clear that lim r →∞ → 0. Since p(ρ) = ωρ is the most usual EoS which has been used in the literature to describe the cosmos, we set L(ρ) = ωρ to have consistency with physics of cosmos in large scale. Note that Eq. (11) provides the difference between the low regime and high regime in the EoS. This is the essential condition for the function H (ρ). As an example, wormhole with a modified Chaplygin EoS [60] with α > 0 is not a suitable case in our method. Because Eq. (11) is not satisfied. The use of the mixed EoS (68) offers two important advantage over p(ρ) = ωρ: First, in the wormhole theory, violation of the NEC is an essential condition at the throat, so one should assume ω < −1 for satisfying flaring out condition. This directed us toward phantom wormhole, but phantom models have problems such as Big Rip singularity or problem of instability in quantum field theory. In phantom model, the violation of the NEC is inevitable in the entire range of spacetime. Mathematically, this assumption is not necessary when the EoS is considered as a mixed model, thereby producing a smaller violation of energy conditions. Second, in the cut and paste method the violation of the NEC can be confined to a small region near the throat. The wormhole metric was matched to an exterior vacuum metric to keep the exotic matter within a finite region of space.
We have found an asymptotically flat spacetime without having to go through acrobatics of cutting and pasting. Our method seems to be more physically than this method. In the cut and paste method, the geometry is divided into two different parts the interior part is the wormhole line element and the exterior part is the geometry of the general cosmos spacetime or at least a geometry which asymptotically is compatible with the metric of cosmos. The surgery may need a surface stress-energy tensor in the boundary of the cut and paste. In our method, the geometry is single and there is not any need to surface stress-energy tensor. The EoS is compatible with cosmos EoS in large scale intrinsically. We have utilized a specific shape function which describes a smooth energy density function with a maximum near the throat and tends to zero at large distance from the throat. Although some classes of EoS are investigated in this paper, this method can be used with the other EoS or shape function to find new solutions. For solutions with quadratic EoS, avoiding horizon at the throat has led to a relation between ω and ω 1 . The same condition is imposed between the ω i parameters on the other EoS. So we can not choose EoS without considering this constraint. We have shown that for a combination of barotropic and polytropic EoS r 1 is dependent on the parameter ω, n and α in which the larger n proposed a fluid with a smaller amount of exotic matter. The general dependence of r 1 has been discussed in detail. We have presented solutions with a more general EoS which are a polynomial function of ρ. The previous solutions can be considered as a special case of these solutions. We have shown that for three individual cases (Eq. (61)), if we choose an EoS with more additional terms, the opportunity to have more consistency with linear EoS will increase. Also, one can conclude that a more general EoS instead of a linear one may present solutions with fewer violation of the NEC. As we know, different astrophysical objects may have different EoS in the local view in comparison with global view. In the other word, the EoS of fluid supporting wormhole near the throat has a behavior different from the behavior of the EoS of the Universe in large scale. This local behaviour of fluid near the throat will tend to a global EoS in the large scale. This method seems to be more physically than considering a unique EoS in all regions of the spacetime. This is the main point of our work which leads to the minimum violation of the NEC.