Stability of spherically symmetric timelike thin-shells in general relativity with a variable equation of state S. Habib Mazharimousavi, M. Halilsoy, S. N. Hamad Amen Department of Physics, Eastern Mediterranean University, Gazima gusa, Turkey. Dated: August 16, 017) We study spherically symmetric timelike thin-shells in 3 1 dimensional bulk spacetime with a variable equation of state for the fluid presented on the shell. In such a fluid the angular pressure p is a function of both surface energy density σ the radius of the thin-shell. Explicit cases of the thin shells connecting two non-identical cloud of strings spacetimes a flat Minkowski spacetime to the Schwarzschild metric are investigated. PACS numbers: Keywords: Thin-shell; Spherically symmetric; Stability; Variable Equation of state arxiv:1708.04588v1 [gr-qc] 3 Aug 017 I. INTODUCTION Time-like thin-shells in spherically symmetric static spacetimes are the most interesting cosmological object which can be constructed in general relativity. Such models of cosmological objects have been used to analyze some astrophysical phenomena such as gravitational collapse supernova explosions. The seminal work of Israel in 1966 [1, ] provided a concrete formalism for constructing the time-like shells, in general, by gluing two different manifolds at the location of the thin-shell. As it was shown in [1], although the metric tensor of the shell which is induced by the bulk spacetime presented in both sides of the shell must be continuous, the extrinsic curvature across the shell is not continuous therefore matter has to be introduced on the shell. This formalism has been employed to study shells in general relativity by many authors for which a good review paper has been worked out by Kijowski et al in ef. [3]. In 1990 an exact solution for a static shell which surrounds a black hole was found by Frauendiener et al [4] its stability was also studied in efs. [5] [6]. In the formalism introduced by Israel there are some conditions which are called Israel junction conditions. These conditions provide a systematic method of finding the energy momentum tensor presented on the shell. In [7] a computer program was prepared to apply the junction conditions on the thin-shells in general relativity using computer algebra. Models of stars circumstellar shells in general relativity were studied in [8]. In [9, 10] the stability of spherically symmetric thin-shells was studied while the gravitational collapse of thin-shells was considered in [11] [1]. Thin-shells in Gauss-Bonnet theory of gravity has been studied in [13] while the rotating thin-shells has been introduced in ef. [14]. Stability of charged thin-shells was considered in [15] its collapse in isotropic coordinates was investigated in [16]. In [17] charge screening by thin-shells in a 1 dimensional regular black hole has been studied while the thermodynamics, entropy, stability of thin-shells in 1 dimensional flat spacetimes have been given in [18] [19]. ecently in [0] the stability of thin-shell interfaces inside compact stars has been studied by Pereira et al which is very interesting as they consider a compact star with the core the crust with different energy momentum tensors consequently with different metric tensors. Screening of the eissner-nordström charge by a thin-shell of dust matter has also been introduced recently in [1]. Finally one of the recent works published in this context is about thin-shells joining local cloud of strings geometries []. In this study first we give a brief review of the thin-shell formalism following that we employ a variable equation of state EoS) for the perfect fluid presented on the surface of the thin-shell of the form p = ψ σ, ) where p is the surface pressure, σ is the surface energy density is the radius of the shell. Let s add that the function ψ must be differentiable with respect to its arguments σ. This model of EoS has been introduced first in [3] was used recently in [4]. Following the generic stability analysis with a variable EoS we study two specific cases including a thin-shell connecting two spacetimes of different cloud of strings a vacuum flat Minkowski spacetime connected to a Schwarzschild metric. Electronic address: habib.mazhari@emu.edu.tr Electronic address: mustafa.halilsoy@emu.edu.tr Electronic address: sarbaz56@gmail.com
II. SPHEICALLY SYMMETIC TIMELIKE THIN-SHELLS In 3 1 dimensional spherically symmetric bulk spacetime a 1 dimensional timelike thin-shell divides the spacetime into two parts which we shall call as the inside outside of the shell. The line elements of the spacetimes in different sides must naturally be different otherwise the thin-shell becomes a trivial invisible object. Let s label the spacetime inside the shell as 1 outside the shell as. Hence the line element of each side may be written generically as ds a = f a r a ) dt a dr a f a r a ) r a dθ a sin θ a dϕ a) in which a = 1, for inside outside, respectively. We add that in general the coordinates i.e., t a, r a, θ a ϕ a need not be the same. In general a thin-shell is a constraint condition on the coordinates of the bulk spacetime but in our study the thin-shell is defined by F := r a τ) = 0 in which τ is the proper time measured by an observer on the shell such that on both sides we define f a ) ) dta 1 dτ f a ) 1) ) d = 1. ) dτ As one of the Israel junction condition, ds ts) from one side to the other of the thin-shell must be continuous. Hence, the coordinates on the shell i.e., τ, θ a ϕ a have to be identical on both sides so that we shall remove the sub index a. This results in a unique induced metric on the shell which is applicable to both sides expressed by ds ts) = dτ τ) dθ sin θdϕ ). 3) Before we proceed further let us note that although the proper time in different sides of the shell is common the coordinate times t a are different they are found from 3) as ṫ a = f a ) Ṙ f a ) in which a dot sts for the derivative with respect to the proper time τ. Let s also add that t 1 t refer to the coordinate times of the inner outer regions. For future use we set our coordinate systems of the bulk i.e., ds a = g µν a) dx a)µ dx a)ν the shell i.e., ds ts) = h ij dξ i dξ j as follow: for the bulk spacetime x a)µ = {t a, r a, θ a, ϕ a } ξ i = {τ, θ, ϕ} for the thin-shell. The second fundamental form or extrinsic curvature) tensor of the shell in each side can be found as K a) ij = n a) x a)γ γ ξ i ξ j x a)α x a)β ) Γa)γ αβ ξ i ξ j 5) in which n a) γ is the four-normal spacelike vector on each side of the thin-shell pointing outward given by ) n a) γ = Ṙ τ), ṫ a, 0, 0. 6) The nonzero components of the extrinsic curvature, therefore, are found to be K a) ττ = τ) f a ), 7) f a ) Ṙ 4) K a) θθ = τ) f a ) Ṙ 8) K ϕϕ a) = τ) f a ) Ṙ sin θ. 9) One observes that unlike the first fundamental form, the second fundamental form is not continuous in general. However, if f a ) f a ) are the same on both sides of the shell then K ij as well as h ij are both continuous. In case that K ij is not continuous it satisfies the other Israel junction condition [K j i ] δj i [K] = 8πGSj i 10)
3 in which [K j i ] = K)j i K 1)j i, [K] = trac[k j i ] = [Ki i ] S j i = diag σ, p, p) 11) is the energy-momentum tensor of the thin-shell. Herein, σ is the energy density p the angular pressure on the shell. We note that, as we have considered the bulk to be spherically symmetric, the pressures in θ ϕ directions are identical the energy-momentum tensor is of a perfect fluid type. The explicit form of the σ p are given by σ = 1 4πG f ) Ṙ f 1 ) Ṙ, 1) τ) p = 1 τ) f ) τ) f 1 ) 8πG f ) Ṙ f 1 ) Ṙ f ) Ṙ f 1 ) Ṙ. 13) τ) In the case of a static thin-shell one has to set τ) = which is a constant consequently σ 0 = 1 f ) ) f 1 ), 14) 4πG p 0 = 1 f ) 8πG f ) f 1 ) f 1 ) f ) ) f 1 ). 15) III. STABILITY Let s assume that our constructed thin-shell is at equilibrium at = which means Ṙ = = 0 therefore the energy density the pressures are given by Eqs. 14) 15). Any radial perturbation causes the radius of the shell to be changed in a dynamical sense. In other words after the radial perturbation becomes a function of proper time τ consequently the energy density the angular pressures are found to be Eqs. 1) 13) which also satisfy dσ d p σ) = 0. 16) This equation is the dynamic relation that connects p σ after the perturbation. Furthermore, any kind of fluid presented on the shell has to satisfy an equation of state which is nothing but a relation between p σ. This relation is traditionally expressed as but in our study we use a more general EoS given by [3, 4] A substitution in 16) yields p = p σ) 17) p = ψ, σ). 18) dσ d ψ, σ) σ) = 0 19) which is a principal equation that connects σ to after the perturbation. In addition to this, from the explicit form of σ in Eq. 1) we find Ṙ V, σ )) = 0 0)
4 in which V, σ )) = f 1 ) f ) f 1 ) f )) 8πGσ )) πgσ )). 1) This equation is a one dimensional equation of motion for the radius of the thin-shell after the perturbation. Together with Eq. 19) gives a clear picture of the motion of the thin-shell after the perturbation. More precisely, the solution of Eq. 19) is used in 0) the general motion of the radius of the thin-shell, in principle is found by solving Eq. 0). The nature of the motion after the perturbation depends on the form of the function ψ, σ) given by EoS the metric functions f 1 ) f ). We comment that the general one-dimensional equation of motion 0) is highly non-linear. In general we do not expect an exact, closed form solution for the radius of the thin-shell after the perturbation. A linearized version of this equation helps us to know the general behaviour of the motion of the thin-shell after the perturbation without going through the complete solution. As we have stated the thin-shell is in equilibrium at = so we exp V, σ )) about = keep it at the first nonzero term. This is called a linearized radial perturbation. The expansion of V, σ )) about = reads as V, σ) = V, σ 0 ) dv d ) 1 d V =0 d ) O ) 3) ) =0 where V, σ 0 ) dv d =0 both are zero identically i.e., the first because Ṙ 0 = 0 the second because = is the equilibrium radius in the sense that the force is zero there. Introducing x =, up to the second order we get ẋ ω x 0 3) in which Derivative with respect to τ implies ω = 1 d V d. 4) =0 ẍ ω x 0 5) which clearly for ω > 0 represents an oscillation about x = 0. This, however, means the radius of the thin-shell oscillates about the equilibrium radius =. This is what we mean by a stable condition. In other words if 1 d V d > 0 6) =0 the thin-shell oscillates remains stable. Unlike ω > 0 if ω < 0 then the motion of the radius of the thin-shell becomes an exponential form which yields an unstable thin-shell. To proceed further we need to calculate d V =0 d as V = V, σ) we shall need σ = dσ d σ = d σ d. Upon differentiation with respect to, 19) yields σ = ψ σ) ψ ) σ 3 ψ. 7) Next, we find V ) at = by applying σ σ whenever we need them which becomes V 0 = 16πGF 0H 0 ψ, [ ) )] H0 F 0 f 10 0 F0 H 0 f 0 0 F 0 H 0 F H) 0 ψ,σ [ ] 4F 4 0 f 10 f 10) F0 0f 10 H 3 0 [ 4H0 4 f 0 f 0) H0 0f F 0 H 0 ) F0 H 0 0 in which F 0 = f 10, H 0 = f 0, ψ, = ψ ψ,σ = ψ. =0 =0 σ 0] F 3 0 8)
5 IV. THIN-SHELL CONNECTING TWO SPACETIMES OF CLOUD OF STINGS In this section we consider f 1 = κ 1 f = κ in which κ 1 κ are two positive constants satisfying κ 1 > κ. Such spacetimes where κ 1, κ 1 represent the so called cloud of strings spacetime [, 5 8]. The energy momentum tensor at the equilibrium is given by S j i = κ1 κ diag 1, 1 ) 4πG, 1 9) which means σ 0 = κ1 κ 4πG 30) p 0 = 1 σ 0. 31) This energy momentum tensor satisfies the weak energy condition which states that σ 0 0, σ 0 p 0 0 therefore it is physical. To proceed with the stability analysis we must choose an EoS. For the first choice we set ψ = ω 1 3) ψ σ = ω 33) in which both ω 1 ω are constants. The general form of V ) given in Eq. 8) yields V ) = [ κ 1 κ 8Gπω1 0 ω 1) κ1 )] κ κ1 ). 34) κ In order to have V ) > 0 one must impose 0 8Gπω 1 0 ω 1) κ 1 κ ) < 0 35) which in turn implies that κ1 ) κ ω 1 < 8Gπ0 ω 1). 36) The case where ω 1 = 0 implies a linear perfect fluid which is stable for ω > 1 is unstable for ω < 1. For the case where ω 1 0 the situation becomes more complicated. In Fig. 1 we plot κ1 ) κ for κ 1 κ ) 8Gπ 0 ω 1 = 8Gπ 0 ω 1) 37) = 0.1, 0., 0.3 0.4. As it is imposed from the condition 35) the values of ω 1 ω under the lines for each specific choice of κ 1 κ ) 8Gπ 0 yields a region of stability while the opposite side above the lines) sts for the values of ω 1 ω that result in an unstable thin-shell. In ef. [] thin-shells joining local cloud of strings geometries has been studied by Eiroa et al. Although their line elements for the cloud of strings geometries was written in static axially symmetric form their equation of state was linear perfect fluid, their results is in agreement with what we found here in Eq. 36) provided ω 1 = 0. Let s add that ω 1 < 0 increases the region of stability in comparison with the linear perfect fluid with ω 1 = 0. This can be seen from Eq. 36) in Fig. 1.
6 FIG. 1: A plot of ω 1 with respect to ω for various values of κ 1 κ ) 8Gπ0 = 0.1, 0., 0.3 0.4. The arrows show the region of stability while the opposite side is the unstable zone for each case. V. THIN-SHELL CONNECTING VACUUM TO SCHWAZSCHILD In this part we consider the inner spacetime to be flat with f 1 = 1 the outer spacetime the Schwarzschild with f r ) = 1 m r. The closed forms of σ 0 p 0 are found to be 1 1 m σ 0 = 38) 4πG m 1 m p 0 =. 39) 8πG0 1 m These clearly satisfy the weak energy conditions i.e., σ 0 0 σ 0 p 0 0 provided > m. Furthermore, one finds in which Ω = which implies 1 m > 0, ψ, = ψ V 0 = 16GπΩ Ω 1 ψ, ψ,σ 1) Ω =0 0 40) ψ,σ = ψ σ. In this case also we set ψ = ω 1 ψ σ = ω =0 V 0 = 16GπΩ 1 Ω ω 1 ω 1) Ω. 41) 0
For the case ω 1 = 0 which corresponds to a linear perfect fluid one finds V 0 0 with ω 1 0 or equivalently ω 1. For the case ω 1 0 we have to work out the regions in the plane of ω 1 ω such that V 0 0. To find the region where V 0 0 we find ω in terms of ω 1 such that V 0 = 0. This amounts to in which Depending on the value of m > m, one finds ω = 1 χω 1 4) 4πG0 χ =. 43) 1 1 m 4πG 0 < χ <. 44) In Fig. we plot ω versus ω 1 for various values for χ = 0.1, 0., 0.3 0.4. Also the stability regions for each case is shown by an arrow indicator. We would like to add that such a thin-shell has been indirectly considered by Eiroa Simeone in [15] with a linear perfect fluid. In both Fig. 1 in [15], the upper left, represents the thin-shell connecting the vacuum to Schwarzschild with their parameter η. Presence of ω 1 changes the stability regime, especially if ω 1 < 0, 41) implies that the term due to ω 1 is positive which increases the stability region. 7 FIG. : A plot of ω with respect to ω 1 for various values of χ = 0.1, 0., 0.3 0.4. The arrows show the region of stability while the opposite side is the unstable zone for each case. Before completing this section we would like to find the explicit form of the energy density σ after the perturbation. In both examples we have worked out in this chapter we assumed ψ = ω 1 ψ σ = ω in which ω 1 ω are two constants. Integration with respect to σ results in ψ = ω 1 ω σ C 0 45)
in which C 0 is an integration constant. As p = ψ should give the equilibrium pressure at = we can find the value of C 0 as therefore the dynamic pressure becomes This EoS together with Eq. 19) gives the differential equation C 0 = p 0 ω 1 ω σ 0 46) p = ω 1 ) ω σ σ 0 ) p 0. 47) dσ d ω 1 ) ω σ σ 0 ) p 0 σ) = 0 48) which must be satisfied by σ. The solution of this equation is given by σ ) = ω σ 0 p 0 ω 1 ω 1 C 1 49) 1 ω 3 ω ω1) in which C 1 is an integration constant. Imposing σ ) = σ 0 yields C 1 = ω1) σ0 p 0 ω 1 0 1 ω Finally the closed form of the energy density is found to be σ ) = ω σ 0 p 0 ω 1 1 ω ω 1 3 ω 3 5ω ω ) ω1) 0 σ0 p 0 1 ω 8 ). 50) ω 1 3 5ω ω ). 51) We note that at =, σ ) reduces to σ 0 it is a function of as well as ω 1 ω. The case ω 1 = 0 admits σ ) = ω σ 0 p 0 1 ω 0 ) ω1) ) σ0 p 0 1 ω 5) while when ω = 0 we find σ ) = p 0 ω 1 ω 1 3 ) 0 σ 0 p 0 ω ) 1. 53) 3 In the case both ω 1 ω are set to zero the energy density becomes while σ ) = p 0 ) 0 σ 0 p 0 ) 54) p = p 0 55) even after the perturbation. According to the Fig., this is one of the case where the thin-shell is stable with V 0 = Ω 0 56) which is clearly positive. VI. CONCLUSION We studied the stability of the timelike thin-shells in spherically symmetric spacetimes with a variable EoS on the shell of the form p = ψ σ, ), [3, 4]. We recall that in previous studies this was taken simply as p = p σ). We analyzed the stability of the thin-shells against a radial perturbation by a linearized approximation we found a general condition to be satisfied in order to have a stable spherically symmetric thin-shell. We applied our results to
two explicit cases with the variable EoS on the shell numerically, as well as analytically we obtained the stability regions. Our first thin-shell joins two spacetimes of local cloud of strings, characterized by two distinct parameters κ 1 > 0 κ > 0 such that κ 1 > κ. The second thin-shell connects a vacuum flat spacetime to the Schwarzschild metric. In both examples the weak energy conditions are satisfied for which the stability regions are identified plotted. 9 [1] W. Israel. Nuovo Cimento B, 44:1, 1966). [] W. Israel. Nuovo Cimento B, 48:463, 1966). [3] D. Kijowski, J. Magli. Gen. elativ. Gravitation, 38:1697, 006). [4] C. Frauendiener, J. Hoenselaers W. Konrad. Classical Quantum Gravity, 7:583, 1990). [5] J. Brady, P.. Louko E. Poisson. Phys. ev. D, 44:1891, 1991). [6] B. G. Schmidt. Phys. ev. D, 59:04005, 1998). [7] P. Musgrave K. Lake. Class.Quant.Grav., 13:1885, 1996). [8] K. G. Zloshchastiev. Int.J.Mod.Phys. D, 8:549, 1999). [9] M. Ishak K. Lake. Phys. ev. D, 65:044011, 00). [10] F. S. N. Lobo P. Crawford. Class. Quant. Grav., :4869, 005). [11] J. Crisostomo. Olea. Phys. ev. D, 69:10403, 004). [1] S. del. Crisostomo, J. Campo J. Saavedra. Phys. ev. D, 70:064034, 004). [13] E. Gravanis S. Willison. Phys. ev. D, 75:08405, 007). [14]. J. Gleiser M. A. amirez. Class. Quant. Grav., 6:045006, 009). [15] E. F. Eiroa C. Simeone. Phys. ev. D, 83:104009, 011). [16] H. Beauchesne A. Edery. Phys. ev. D, 85:044056, 01). [17] S. H. Mazharimousavi M. Halilsoy. Eur. Phys. J. C., 73:57, 013). [18] J. P. S. Lemos G. M. Quinta. Phys. ev. D, 88:067501, 013). [19] J. P. S. Lemos G. M. Quinta. Phys. ev. D, 89:084051, 014). [0] J. G. Pereira, J. P. Coelho J. A. ueda. Phys. ev. D, 90:13011, 014). [1] S. H. Mazharimousavi M. Halilsoy. Eur. Phys. J. C., 75:334, 015). [] E. Eiroa, E. F. ubin de Celis C. Simeone. Eur. Phys. J. C., 76:546, 016). [3] N. M. Garcia, F. S. N. Lobo M. Visser, Phys. ev. D 86:04406, 01). [4] V. Varela Phys. ev. D 9:04400, 015). [5] P. S. Letelier, Phys. ev. D 0:194 1979). [6] S. G. Ghosh S. D. Maharaj, Phys. ev. D 89:08407 014). [7] S. G. Ghosh, U. Papnoi S. D. Maharaj, Phys. ev. D 90:044068 014). [8] S. H. Mazharimousavi M. Halilsoy, Eur. Phys. J. C 76:95 016).