Wormholes admitting conformal Killing vectors and supported by generalized Chaplygin gas

When Morris and Thorne first proposed that traversable wormholes may be actual physical objects, they concentrated on the geometry by specifying the shape and redshift functions. This mathematical approach necessarily raises questins regarding the determination of the required stress-energy tensor. This paper discusses a natural way to obtain a complete wormhole solution by assuming that the wormhole (1) is supported by generalized Chaplygin gas and (2) admits conformal Killing vectors.


Introduction
Wormholes are handles or tunnels in spacetime connecting different regions of our Universe or different universes altogether. That wormholes could be actual structures suitable for interstellar travel was first proposed by Morris and Thorne [1]. Such a wormhole could be described by the static spherically symmetric line element [1] ds 2 = −e 2Φ(r) dt 2 + dr 2 1 − b(r)/r + r 2 (dθ 2 + sin 2 θ dφ 2 ), using units in which c = G = 1. Here Φ = Φ(r) is called the redshift function, which must be everywhere finite to avoid an event horizon. The function b = b(r) helps determine the spatial shape of the wormhole and is therefore called the shape function. The spherical surface r = r 0 is the throat of the wormhole and must satisfy the following conditions: b(r 0 ) = r 0 , b(r) < r for r > r 0 , and b ′ (r 0 ) < 1, now usually called the flare-out condition. This condition refers to the flaring out of the embedding diagram pictured in Ref. [1]. The flare-out condition can only be satisfied by violating the null energy condition. The Einstein field equations in the orthonormal frame, Gμν = 8πTμν, yield the following simple interpretation for the components of the stress-energy tensor: Ttt = ρ(r), the energy density, Trr = p r , the radial pressure, and Tθθ = Tφφ = p t , the lateral pressure. For the theoretical construction of the wormhole, Morris and Thorne then proposed the following strategy: retain complete control over the geometry by specifying the functions b(r) and Φ(r) to obtain the desired properties of the wormhole. The problem with this strategy is that it relies on the engineering team to manufacture or search for those materials or fields that yield the required stress-energy tensor.
This purely geometric approach can be supplemented by an appropriate equation of state, an example of which is p = ωρ, ω < −1, representing phantom dark energy [2,3,4,5]. The energy density may also be known for physical reasons, as in the case of the Navarro-Frenk-White density profile for dark matter [6,7] where r s is the characteristic scale radius and ρ s is the corresponding density.
In this paper we consider another equation of state, p = −A/ρ, A > 0, representing Chaplygin gas or, more generally, p = −A/ρ α , 0 < α ≤ 1, called generalized Chaplygin gas [8,9,10,11,12]. Cosmologists became interested in this model when it was shown that Chaplygin gas is a candidate for unifying dark matter and dark energy. To see this, consider the energy conservation equationρ = −3ȧ(ρ + p)/a in a flat FLRW spacetime and then substitute p = −A/ρ α . The result is where B is an integration constant. It is seen that ρ ∼ a −3 at early times, implying that ρ behaves like matter, while in later times it behaves like a cosmological constant (ρ ≡ constant).
In a cosmological setting one would normally assume a homogeneous distribution of matter. It was pointed out in Ref. [13], however, that the generalized Chaplygin equation of state is that of a polytropic gas with a negative polytropic index. Thus inhomogeneous structures may arise from a density fluctuation in the cosmological background. Also, a Born-Infeld phantom gravastar has been constructed by replacing the interior de Sitter solution with the Chaplygin gas equation of state [14]. A similar problem arises with another type of dark energy, phantom dark energy. Here Sushkov and Kim [15] have shown that away from the throat the radial and transverse pressures converge fairly quickly. Accordingly, we will follow Refs. [10,11] and assume that the equation of state is The tangential pressure can be determined from the Einstein field equations.
To obtain a complete wormhole solution, we are going to make the additional assumption that our spacetime admits a one-parameter group of conformal motions, i.e., conformal Killing vectors, discussed next.

Conformal Killing vectors
As noted above, we assume that our spacetime admits a one-parameter group of conformal motions, which are motions along which the metric tensor of a spacetime remains invariant up to a scale factor. In other words, there exists a set of conformal Killing vectors such that where the left-hand side is the Lie derivative of the metric tensor and ψ(r) is the conformal factor. The vector ξ characterizes the conformal symmetry since the metric tensor g µν is conformally mapped into itself along ξ. The assumption of conformal symmetry has led to numerous new solutions, as well as new geometric and kinematical insights [16,17,18,19,20,21,22]. Given a noncommutative-geometry background, exact solutions of traversable wormholes admitting conformal motions are discussed in Ref. [23]. Two earlier studies assumed a non-static conformal symmetry [24,25].
As in Ref. [23], we use the following form of the metric to discuss conformal symmetry: The Einstein field equations are: e −λ 1 and Eq. (9) can actually be obtained from the conservation of the stress-energy tensor, i.e., T µν ;ν = 0. So we need to use only Eqs. (7) and (8).
To discuss the assumption of conformal symmetry in Eq. (5), we follow Herrera and Ponce de León [16] and restrict the vector field by requiring that ξ α U α = 0, where U α is the four-velocity of the perfect fluid distribution. The assumption of spherical symmetry then yields ξ 0 = ξ 2 = ξ 3 = 0 [16]. Eq. (5) now produces the following results: and These equations, in turn, yield e ν = Cr 2 and where C and a are integration constants. In order to make use of these equations, it is convenient to write Eqs. (7) and (8) in the following forms:

The solution
To obtain a wormhole solution, we start with the equation of state (4), p r = −A/ρ α , and substitute Eqs. (15) and (16) to obtain 1 8π This equation can be put into a more transparent form by noting that 2ψψ ′ = (ψ 2 ) ′ and eliminating the negative sign: Now observe that since A and ρ α are positive, the left side is also positive, allowing us to raise each side to the power 1/α. We can thereby obtain Before considering a possible solution, observe that by comparing Eqs. (1) and (6), we have in view of Eq. (14), So to obtain b(r), Eq. (19) must be solved for ψ 2 (r). To satisfy the condition b(r 0 ) = r 0 , we must have ψ 2 (r 0 ) = 0, which becomes the initial condition for Eq. (19). Moreover, since 1 − 3ψ 2 /a 2 > 0 by Eq. (18), we also have 1 − ψ 2 /a 2 > 0, so that b(r) > 0.
To check the flare-out condition b ′ (r 0 ) < 1, we obtain from Eqs. (20) and (19), So to meet the flare-out condition, the constant A from the Chaplygin equation of state must satisfy the inequality  which agrees with Refs. [10,11]. Now we turn our attention to Eq. (19), recalling the initial condition ψ 2 (r 0 ) = 0. Since this equation does not have a closed-form solution, we will use a numerical approach. To do so, we choose an arbitrary value for r = r 0 and some values of A and α that satisfy the above conditions. It becomes apparent immediately that the solution for b(r) is independent of the integration constant a. The reason for this can be seen by writing Eqs. (19) and (20) in the following respective forms: and So the solutions of Eqs. (19) and (22) have the same qualitative forms. (In other words, by rescaling ψ, the constant a could be eliminated.) The plots for ψ 2 and b(r), using some typical values of the parameters, are shown in Fig. 1. For any particular choice of A, the function ψ 2 becomes zero for some r = r 1 and hence negative for r > r 1 . Since the corresponding function ψ is now imaginary, b = b(r) is undefined for r > r 1 . (Observe that b(r 1 ) = r 1 , as shown in Fig. 1; also, b(r) < r for r 0 < r < r 1 .) We conclude that the wormhole material is confined to some interval [r 0 , r 1 ]. This conclusion is consistent with Refs. [10,11], which state that the dimensions of the wormhole cannot be arbitrarily large.
As a final comment, while the length of the interval [r 0 , r 1 ] is independent of a, it does depend on A: the closer A is to the upper limit 1/(8πr 2 0 ) α+1 , the closer b ′ (r 0 ) is to unity and hence the smaller the interval [r 0 , r 1 ]; this behavior can also be seen from Fig. 1.

The redshift function
We see from Eq. (13), e ν = Cr 2 , that the wormhole spacetime cannot be asymptotically flat. Normally one would now introduce an extra requirement by stating that the wormhole material has to be cut off at some r and joined to an exterior Schwarzschild solution Since the need for a cut-off is already known from the previous section, we are not dealing with a new requirement. Unfortunately, the previous cut-off at r = r 1 cannot be used here because b(r 1 ) = r 1 . Instead, we need to choose some r = r 2 < r 1 . We can therefore complete the wormhole solution by noting that So for e ν = Cr 2 2 , we have Cr 2 2 = 1 − 2M/r 2 and the integration constant becomes The need for a cut-off at r = r 2 < r 1 , due to the coordinate singularity at r = r 1 , is consistent with Refs. [26,27], which consist of a detailed investigation of the Tolman-Oppenheimer-Volkoff equations for Chaplygin and generalized Chaplygin gas, respectively. The main conclusion in these studies is much more general in the sense that the scalar curvature was found to become singular at some finite distance; moreover, the Universe is not asymptotically flat. That our wormhole spacetime cannot be arbitrarily large was also shown in Ref. [10] without relying on conformal symmetry. However, in the present study, conformal symmetry is needed to determine the redshift function (to obtain a complete wormhole solution) but not for the wormhole solutions in Refs. [26,27].

Conclusion
For the theoretical construction of a traversable wormhole, Morris and Thorne proposed the following strategy: retain complete control over the geometry by choosing functions Φ(r) and b(r) that yield the desired properties, thereby leaving the components of the stress-energy tensor unspecified. In this paper we obtained a complete wormhole solution by (1) adopting the equation of state p r = −A/ρ α representing generalized Chaplygin gas and (2) assuming that the wormhole admits a one-parameter group of conformal motions. These two assumptions complement each other via the requirement that the wormhole spacetime must be cut off at some r = r 2 and joined to an exterior vacuum solution.