arxiv: v1 [gr-qc] 8 Apr 2009

Transcription

1 On the Stbility of Sttic Ghost Cosmologies rxiv: v1 [gr-qc] 8 Apr 2009 John D. Brrow 1 nd Christos G. Tsgs 2 1 DAMTP, Centre for Mthemticl Sciences University of Cmbridge, Wilberforce Rod, Cmbridge CB3 0WA, UK 2 Section of Astrophysics, Astronomy nd Mechnics, Deprtment of Physics Aristotle University of Thessloniki, Thessloniki 54124, Greece Abstrct We consider the liner stbility of sttic universe filled with non-intercting mixture of isotropic rdition nd ghost sclr field. Unlike the conventionl Einstein sttic model, this cosmology is stble ginst homogeneous nd isotropic perturbtions. This is shown by mens of exct oscilltory solutions bout the originl sttic stte. Using covrint techniques, we lso exmine the liner response of the sttic ghost universe to ll types of inhomogeneous fluctutions, nmely density, vorticity nd grvittionl-wve perturbtions. The results show tht the sttic bckground constrins the liner evolution of these distortions, to the extent tht density perturbtions remin time invrint, vorticl distortions vnish nd grvittionl wves oscillte with constnt mplitude. We discuss potentil implictions of these results for pst-eternl initil sttes in clssicl generl reltivistic cosmology. 1 Introduction Cosmologists hve retined n interest in the consequences of dmitting ghost sclr fields into the Universe. These fields my hve negtive energy density [1], nd re fr from new, hving been incorported in the erly study of stedy-stte cosmology [2, 3], nd more recently in the investigtion of sclr fields [4], phntom drk mtter [5, 6], nd κ-essence [7]. They lso provide simple theoreticl lbortory for exploring the physicl consequences of cyclic closed universes which bounce t finite rdius. In prticulr, Brrow et l [8] exmined the unusul consequences tht rise for theories of vrying constnts in such oscillting universes. Ellis nd collbortors [9, 10] hve lso considered cosmologies tht emerge into expnsion from pst eternl Einstein sttic stte, reminiscent of the originl Eddington-Lemître cosmology fvoured by Eddington s n infinitely old universe tht cn nonetheless be thermodynmiclly young [11, 12]. One further consequence of these ghost cosmologies is the existence of sttic solutions with new properties. In this pper we will exmine the stbility of Einstein sttic universes in these theories. Two new fetures will be of interest. First, the existence of Einstein sttic solutions which re stble ginst conforml (i.e. homogeneous nd isotropic Friedmnn) perturbtions; nd second, the liner stbility of certin ghost-like sttic models under ll three types of inhomogeneous perturbtions (i.e. sclr, vector nd tensor modes). The inhomogeneous sitution shres some nlogies with the stbility shown by the Einstein sttic universe in the presence of conventionl (non-ghost) mtter fields, s discussed in [13]-[16]. 1 1 The stbility nlysis discussed there is lso rendered delicte becuse the Einstein sttic universe hs compct spce sections nd Killing vectors nd so is conicl point in the spce of ll solutions to the Einstein equtions 1

2 Throughout this pper we ssume pseudo-riemnnin spcetime with Lorentzin metric of signture ( + + +). We lso dopt the 1+3 covrint pproch to generl reltivity nd cosmology, referring the reder to older nd more recent reviews of the formlism for further reding nd detils [20]-[22]. 2 Ghost sclr fields At the bsis of the 1+3 covrint formlism is the concept of the fundmentl observers, which when introduced llow for unique threding of the spcetime into time nd 3-dimensionl spce. The worldlines of these observers re tngent to timelike 4-velocity field, reltive to which physicl quntities nd equtions decompose into their timelike nd spcelike components. The formlism combines mthemticl compctness with physicl trnsprency nd hs been used to vriety of cosmologicl re generl reltivistic studies. Here we will use the covrint pproch to study cosmologicl models contining ghost sclr field (ψ) with zero potentil (i.e. V (ψ) = 0) nd Lgrngin of the form [8] L ψ = 1 2 ψ ψ g, (1) where g is the determinnt of the spcetime metric (g b ) nd is the stndrd covrint derivtive opertor. Then, the ssocited stress-energy tensor reds T (ψ) b = ψ b ψ cψ c ψg b. (2) To chieve n 1+3 decomposition of the ψ-field, we ssume tht ψ is timelike (nmely tht ψ ψ < 0) over given spcetime region. In tht cse ψ is norml to the spcelike hypersurfces ψ(x ) = constnt nd defines our fundmentl timelike 4-velocity field u = 1 ψ ψ, (3) with ψ = u ψ 0. This gurntees tht ψ 2 = ψ ψ > 0 nd tht u u = 1 s required. The bove 4-velocity lso defines the fundmentl time direction nd the tensor h b = g b + u u b tht projects into the 3-spce orthogonl to u. The projection tensor lso provides the covrint derivtive D = h b b, which opertes in the observers instntneous rest spce. In covrint terms, the kinemtics of the u -congruence (3) re determined by mens of the irreducible decomposition [11, 22] b u = 1 3 Θh b + σ b + ω b A u b, (4) where Θ, σ b nd ω b respectively describe the verge volume expnsion/contrction, sher distortions nd the vorticity of the flow. Following our 4-velocity choice, the lter vnishes identiclly (i.e. ω b = 0), mking the u -field irrottionl. 2 The 4-ccelertion vector (A ) which exhibit lineristion instbility bout it unless higher-order constrints re imposed [17]-[19]. 2 An direct consequence of our 4-velocity choice (see (3)) is tht D ψ = 0. This in turn gurntees tht ω b = D [b u ] = 0 nd the irrottionl nture of the flow (e.g. see [22]). 2

3 describes non-grvittionl/non-inertil forces nd is given by A = 1 ψ D ψ. (5) The bove, which my be seen s the momentum conservtion lw of the ψ-field (compre to Eq. (29b) in 4.2), implies tht ψ cts s n ccelertion potentil. Introducing the u -frme (3) lso fcilittes convenient fluid-like description of the ψ- field. In prticulr, with respect to the fundmentl observers, the energy-momentum tensor (2) recsts into T (ψ) b = ρ (ψ) u u b + p (ψ) h b, (6) where ρ (ψ) = p (ψ) = 1 2 ψ 2. (7) The right-hnd side of (7) is negtive by defult, ensuring tht ρ (ψ) < 0 nd consequently the ghost nture of the ψ-field. In the presence of non-vnishing potentil, sclr fields do not generlly behve like brotropic fluids. Here, however, the potentil hs been set to zero nd brotropic description of the ψ-field is therefore possible. According to (7), our sclr field corresponds to stiff medium with effective brotropic index w (ψ) = p (ψ) /ρ (ψ) = 1 nd corresponding sound-speed given by c 2 ψ = dp(ψ) /dρ (ψ) = 1. 3 Homogeneous ghost cosmologies 3.1 Sttic sptilly closed models The evolution of Friedmnn-Robertson-Wlker (FRW) cosmology with nonzero sptil curvture nd non-vnishing cosmologicl constnt is monitored by the set nd H 2 = 1 3 (ρ + Λ) K 2, Ḣ = H2 1 6 (ρ + 3p) Λ (8) ρ = 3H(ρ + p). (9) In the bove ρ nd p represent the totl energy density nd pressure of the mtter, K = 0, ±1 is the 3-curvture index nd Λ is the cosmologicl constnt. In ddition H = ȧ/ is the Hubble prmeter of the model nd is the cosmologicl scle fctor. To close this system one needs to supply n (effective) eqution of stte for the (totl) mtter sources. Consider sttic cosmology with zero cosmologicl constnt nd positive sptil curvture. Then, the cosmologicl scle fctor is time-independent ( = 0 = constnt) nd cn be seen s the rdius of the universe. Setting Λ = 0, K = +1 nd H = 0 = Ḣ, Eqs. (8)-(9) recst into the constrins ρ = 3 2, ρ + 3p = 0 nd ρ = 0, (10) 0 respectively. The former of the bove shows tht the size of our model is inversely proportionl to the energy density of its mtter content, nmely tht 0 = 3/ρ. The second gurntees tht the effective brotropic index of the system is w = p/ρ = 1/3. 3

4 Suppose now tht mtter forms non-intercting mixture of isotropic rdition, with p (γ) = ρ (γ) /3, supplemented by ghost sclr field (ψ) with zero potentil nd n effective eqution of stte given by (7). When both components re comoving, their totl energy density nd pressure, reltive to fundmentl observer, re ρ = ρ (γ) + ρ (ψ) nd p = p (γ) + p (ψ) respectively. In tht cse expressions (10) reduce to nd ρ (γ) + ρ (ψ) = 3 2, ρ(γ) + 2ρ (ψ) = 0 (11) ρ (γ) + ρ (ψ) = 0. (12) Of the bove constrints, (11b) combines with (11) to ensure tht the energy densities of the two species re directly relted, ccording to This result then, together with (11), leds to ρ (γ) = = ρ (ψ) = 1 2 ρ(γ) = ρ. (13) 6 ρ (γ) nd ρ (ψ) = = 3, (14) ρ (ψ) expressing the rdius of our sttic model in terms of the energy densities of the individul mtter components. We finlly note tht, for non-intercting fluids, the conservtion lw (12) is lso seprtely stisfied (i.e. ρ (γ) = 0 = ρ (ψ) ). 3.2 Almost-sttic oscilltory models Before looking into inhomogeneous cosmologies with ghost-like sclr fields, we will first consider n exctly soluble lmost-sttic homogeneous model. The Friedmnn eqution for the scle fctor, (t), of closed universe contining ghost sclr field with density ρ (ψ) = ψ 2 /2 = Ψ 6 < 0, nd rdition with density ρ (γ) = Γ 4, reds ȧ 2 2 = Ψ 6 + Γ 4 1 2, (15) where Ψ nd Γ re positive constnts. Expressed in terms of conforml time (η, with dη = 1 dt), the bove integrtes to give 2 (η) = 1 2 { Γ + } Γ 2 4Ψ sin[2(η + η 0 )], (16) with η 0 constnt when Γ 2 4Ψ. Identifying the globl mximum nd the globl minimum of the expnsion with mx = [Γ + Γ 2 4Ψ]/2 nd min = [Γ Γ 2 4Ψ]/2 respectively, we find tht the scle fctor evolution is given by 2 (η) = 1 2 { } 2 mx + 2 min + (2 mx 2 min )sin [2(η + η 0)]. (17) 4

5 This solution oscilltes bout finite rdius between successive mxim nd minim. 3. If we hd chosen K = 0 then there would hve been single minimum (nd no mximum) with collpse for < min nd expnsion for > min. Since ψ = 2Ψ 3, the sclr field evolution given by [8] [ Γ tn (η + ψ = ±2tn 1 η0 ) + ] Γ 2 4Ψ 2. (18) Ψ With mx min, we hve min = Ψ/Γ nd mx = Γ. We cn then see tht the bounce durtion is t 2 min / mx. Since ψ 6Γ 3/2 /Ψ ner the bounce, we find tht ψ 6, independently of initil conditions, during ech bounce. The extreme cse Γ 2 = 4Ψ is sttic universe with = min = mx = Γ/2 = Ψ. Setting ȧ = 0 nd ä = 0, we cn see tht this cse is relized when ρ ψ = ρ r /2, giving = 6/ ρ r. We cn see explicitly tht the solution = Γ/2 is stble ginst homogeneous nd isotropic perturbtions becuse they led to universe with regulr sinusoidl oscilltions of rbitrrily smll mplitude round = Γ/2. These solutions cn be interpreted s lmost sttic oscilltory universes or cyclic universes of smll mplitude. Such solutions re described by (17) in the cse where mx min. This sitution differs from tht found in generl reltivity in the bsence of ghost fields [13]-[16]. We turn now to consider the stbility of these sttic solutions ginst inhomogeneous density perturbtions using guge-invrint formlism. 4 Perturbed sttic ghost cosmologies 4.1 Describing inhomogeneity In ccord with the 1+3 covrint pproch to cosmology, sptil inhomogeneities in the density distribution of the mtter fields re described vi the ssocited comoving frctionl density grdients. For the totl fluid, these re given by the dimensionless rtio [23] = ρ D ρ, (19) with D = h b b representing the covrint derivtive opertor in the observers rest spce. Similrly, the grdients (γ) = ρ (γ) D ρ (γ) nd (ψ) = ρ (ψ) D ρ (ψ), (20) monitor inhomogeneities in the densities of the two individul mtter components. When the bckground model is sptilly homogeneous ll of the bove grdients vnish nd therefore describe liner density perturbtions in guge invrint mnner [24]. Also, definitions (19) nd (20) imply tht = 1 ρ ( ρ (γ) (γ) ) + ρ (ψ) (ψ). (21) The density grdients re supplemented by dditionl inhomogeneity vribles, of which the most importnt describes fluctutions in the volume expnsion (Θ). The ltter is represented 3 For mthemticl simplicity we hve specilised the mtter content of the model to be rdition plus the ghost field. The qulittive behvior of the Friedmnn eqution would be similr if the rdition ws replced by ny perfect fluid with ρ > p > ρ/3. 5

6 by the orthogonlly projected grdient Z = D Θ, (22) which is lso independent of the guge choice. Consider the bckground model described in 3.1 nd perturb it by introducing wek inhomogeneities nd nisotropies. To liner order, definitions (19), (20) nd (22) red nd = 0 ρ D ρ, (γ) = 0 ρ (γ) D ρ (γ), respectively. Also to first order, reltion (21) reduces to (ψ) = 0 ρ (ψ) D ρ (ψ) (23) Z = 0 D Θ, (24) = 2 (γ) since ρ (ψ) = ρ (γ) /2 = ρ in the bckground (see conditions (11)). 4.2 Liner density inhomogeneities (ψ), (25) In non-intercting multi-component systems, the individul members observe seprte conservtion lws. Since the ψ-field does not interct with the rditive component, their ssocited liner momentum conservtion formule re ( w (i)) A = c 2(i) s (i), (26) with i = γ, ψ nd A being the common 4-ccelertion vector (e.g. see [22]). In the bsence of interctions, liner density inhomogeneities in ech one of the species evolve ccording to ( = 1 + w (i)) Z, (27) where (i) Z = 1 ) (1 2 ρ + 3c 2 s 0 2 ρ(1 + 3w)A + 0 D D b A b, (28) to liner order. Note tht ρ, p, w nd c 2 s correspond to the totl fluid, nd is given by (19). When deling with perturbed sttic universe filled with non-intercting mixture of rdition nd ghost sclr field, we hve w (γ) = 1/3 = c 2(γ) s nd w (ψ) = 1 = c 2(ψ) s t the zero perturbtive level. Hence, the liner expressions (27)-(26) led to 4 0 A = (γ), 2 0 A = (ψ) (29) nd (γ) = 4 3 Z, (ψ) = 2Z. (30) The first of these sets implies tht (γ) = 2 (ψ), while second ensures tht (γ) (with C being time-independent). Therefore, to liner order (γ) = 2 (ψ) /3 + C = 2 (ψ) = 3 2 C. (31) 6

7 In other words, liner density inhomogeneities (both in the individul species nd the totl fluid see Eq. (25)) remin constnt in time. This result seems to reflect the high symmetry of the unperturbed model nd the properties of the two mtter fields. In prticulr, the sttic bckground, the bsence of cosmologicl constnt nd the non-intercting nture of the mteril mixture men tht the energy densities of the constituent species re directly relted (see constrints (11)-(14)). When the two mtter fields lso shre the sme 4-velocity, s it hppens here, the forementioned constrints led to the momentum conservtion lws (29). The ltter, together with the propgtion formule (30), gurntee tht ny liner density inhomogeneities tht might exist will not chnge in time. 4 Result (31) lso gurntees the time independence of liner density perturbtions (described vi the sclr = 0 D ) s well s of shpe distortions (monitored by mens of the trce-free tensor b = 0 D b ). On the other hnd, vorticl perturbtion in the density distribution of the species (described vi the ntisymmetric tensor [b] = 0 D [b ] ) vnish identiclly. To verify this recll tht (e.g. see Appendix A.3 in [22]) D [b D ] ρ = ρω b = 0, (32) with the null result gurnteed on ll perturbtive levels by the irrottionl nture of the u - frme (see 2 erlier). 4.3 Liner grvittionl-wve perturbtions So fr, our nlysis hs shown tht sclr (density) perturbtions re time independent, while vector (rottionl) distortions vnish. The remining third type, concerns pure tensor fluctutions. Grvittionl wves re covrintly described by the trnsverse prts of the electric (E b ) nd the mgnetic (H b ) components of the Weyl tensor. To liner order, these re isolted by switching off both the sclr nd the vector perturbtions (to gurntee tht D b E b = 0 = D b H b e.g. see [22]). On sttic, FRW-type bckground, the liner propgtion of these two fields is monitored by the set nd Ė b = 1 2 ρ(1 + w)σ b + curlh b (33) Ḣ b = curle b. (34) where curle b = ε cd D c E b d for ny trceless nd trnsverse tensor E b. Also to first order, the mgnetic prt of the Weyl field is relted to the sher by mens of curlh b = 3 2 σ b D 2 σ b, (35) 0 with D 2 = D D representing the orthogonlly projected covrint Lplcin nd 0 the rdius of the bckground 3-D hypersurfces. The bove ensures the direct interconnection between the mgnetic Weyl field nd the sher t the liner perturbtive level. This in turn llows us to replce Eq. (34) with the sher propgtion formul σ b = E b. (36) 4 The time independence of (γ) nd (ψ) implies tht the liner expnsion inhomogeneities must vnish t ll times (see Eqs. (30)). This imposes strict condition to the right-hnd side of (28), which must remin zero t ll times. Checking the consistency of the liner condition goes beyond the scope of this rticle. 7

8 The ltter combines with (33) to give the wve-eqution σ b = 1 2 ρ(1 + w)σ b 3 2 σ b + D 2 σ b, (37) 0 of the trnsverse prt of the sher tensor. When the sttic bckground contins mixture of non-intercting rdition nd ghost sclr field, w = 1/3 nd ρ = 3/ 2 0 (see 3.1). Applied to this environment, Eq. (37) reduces to σ b = 2 2 σ b + D 2 σ b. (38) 0 Our next step is to Fourier decompose the lst expression. Introducing the stndrd tensor hrmonic functions, we my write σ b = n σ (n)q (n) b, where Q(n) b re the forementioned hrmonics with Q (n) b = Q (n) b, Q (n) b = 0 = Db Q (n) b nd 5 ( ) D 2 Q (n) n 2 b = Q (n) b. (39) 0 Substituting decomposition σ b = n σ (n) Q (n) b into Eq. (38), the hrmonics decouple nd we rrive t the wve eqution of the n-th sher mode ( n 2 ) + 2 σ (n) = σ (n), (40) where n 2 3 due to the compctness of the unperturbed sptil sections. Expression (40) governs the liner evolution of grvittionl wves on sttic, FRW-type bckground tht hs sphericl sptil geometry nd contins mixture of non-intercting blck-body rdition with ghost sclr field. The solution of (40) hs the simple oscilltory form [( ) ] [( ) ] n σ (n) = C 1 sin n t + C 2 cos t, (41) 2 0 which ensures tht grvittionl-wve perturbtion oscillte with time-independent mplitude. In other words, the pure tensor modes re neutrlly stble t the liner level on ll scles. 5 Non-liner grvittionl wve nisotropies We cn see tht the ddition of lrge mplitude grvittionl wve nisotropies cn render the sttic universe unstble. This cn be seen by exmining Eq. (15) for the isotropic nd homogeneous bckground. Suppose tht we dd simple nisotropy to the expnsion, with isotropic 3-curvture. This is of the Binchi I type nd dds new nisotropy energy density term for the sher, σ 2 = Σ 2 / 6 with Σ constnt, to the Friedmnn eqution, where (t) is now the geometric men scle fctor nd the Friedmnn eqution becomes (ȧ ) 2 = Σ2 6 Ψ 6 + Γ (42) 5 Angled brckets denote the symmetric nd trce-free prt of orthogonlly projected second rnk tensors. 8

9 Clerly, if the nisotropy is sufficiently lrge, so tht Σ 2 > Ψ, then neither the sttic solution nor the exct solution with finite non-singulr oscilltions round it (16)-(17) ny longer exist. This sitution requires nisotropic perturbtions to the FRW bckground tht re fr from liner. We cn envisge n irregulr erly stte in which the ghost field, rdition, nisotropy nd 3- curvture index, K, ll vry in spce. The effective Friedmnn eqution for the sptilly-vrying scle fctor ( x,t) would hve the form (ȧ ) 2 = Σ2 ( x) 6 Ψ( x) 6 + Γ( x) 4 K( x) 2. (43) Regions where Σ 2 ( x) > Ψ( x) would not oscillte through successive non-singulr cycles; regions where k( x) 0 would bounce only once nd then expnd forever. Only those regions with low nisotropy nd positive curvture would evolve through succession of periodic cycles. If there ws production of rdition entropy from cycle to cycle, in ccord with the Second Lw of thermodynmics, then this would be equivlent to n increse in the vlue of Γ from cycle to cycle [25] nd corresponding increse in the mplitude of the oscilltions round the sttic universe in Eq. (16), since they scle with Γ. They would pproch the zero curvture expnsion from below. If positive cosmologicl constnt were dded to the right-hnd side of Eq. (43), then, no mtter how smll its mgnitude, it would eventully come to control the dynmics if successive cycles grow in mplitude. Therefter, the oscilltions would cese nd the expnsion would symptote towrds n ever-expnding sptilly-flt de Sitter stte, s shown in [25]. More complicted forms of nisotropy, including 3-curvture nisotropies, would disply similr instbility becuse t erly times they include lso the simple nisotropic stress considered here. 6 Discussion We hve shown tht in simple closed Friedmnnin cosmologies contining blck-body rdition nd ghost sclr field, there exist exct sttic solutions, which re stble ginst conforml homogeneous nd isotropic perturbtions (tht need not be smll). This is in shrp contrst to the instbility shown by the Einstein-sttic cosmology, when conventionl mtter is used [11]. Our sttic universe lso ppers stble ginst ll three type of liner inhomogeneous perturbtions, in wy nlogous to tht seen in [13]-[16]. However, if lrge-mplitude homogeneous nd nisotropic perturbtions re introduced, the model becomes unstble (it no longer exists), with its pst nd future evolution moving towrds strong curvture singulrity. In the bsence of these lrge perturbtions the sttic solutions with stble bounded oscilltions could provide vible pst-eternl stte for the Universe, lthough it would be subject to evolution if entropy ws required to increse with time. Our nlysis is of interest for studies of the very erly universe incorporting theories of prticle physics with new types of ghost mtter field. It shows tht it is possible for universes to get trpped in erly-time evolutionry trcks tht exhibit bounded oscilltions bout non-singulr sttic stte. Similr behviour might lso rise in higher-order grvity theories [26]-[28], which introduce stresses tht mimic the presence of ghost fields nd would prevent some domins from prticipting in infltion. References [1] Gibbons G W, hep-th/

10 [2] Hoyle F 1948 Mon. Not. R. Astron. Soc [3] McCre W H 1951 Proc. Roy. Soc. A [4] Burd A B nd Brrow J D 1988 Nucl. Phys. B [5] Cldwell R 2002 Phys. Lett. B [6] Brown M G, Freese K nd Kinney W H 2008 JCAP [7] Armendriz-Picon C, Mukhnov V nd Steinhrdt P 2001 Phys. Rev. D [8] Brrow J D, Kimberly D nd Mgueijo J 2004 Clss. Quntum Grv [9] Ellis G F R nd Mrtens R 2004 Clss. Quntum Grv [10] Ellis G F R, Murugn J nd Tsgs C G 2004 Clss. Quntum Grv [11] Eddington A S 1930 Mon. Not. R. Astron. Soc [12] Brrow J D nd Mestel L 2004 Astronomy nd Geophysics [13] Hrrison E R 1967 Rev. Mod. Phys [14] Gibbons G W 1987 Nucl. Phys. B [15] Gibbons G W 1988 Nuc. Phys. B [16] Brrow J D, Ellis G F R, Mrtens R nd Tsgs C G 2003 Clss. Quntum Grv. 20 L155 [17] Brill D nd Deser S 1973 Commun. Mth. Phys [18] Brrow J D nd Tipler F J 1979 Phys. Rep [19] Losic B nd Unruh W 2005 Phys. Rev. D [20] Ellis G F R 1973 in Crgèse Lectures in Physics, ed. E Schtzmn (New York: Gordon nd Brech) 1 [21] Ellis G F R nd vn Elst H 1999 in Theoreticl nd Observtionl Cosmology, ed M Lchièze- Rey (Dordrecht: Kluwer) 1 [22] Tsgs C G, Chllinor A nd Mrtens R 2008 Phys. Rep [23] Ellis G F R nd Bruni M 1989 Phys. Rev. D [24] Stewrt J M nd Wlker M 1974 Proc. R. Soc. A [25] Brrow J D nd D browski M 1995 Mon. Not. R. Astron. Soc [26] Brrow J D nd Hervik S 2006 Phys. Rev. D [27] Brrow J D nd Hervik S 2006 Phys. Rev. D [28] Kloper N 1991 Phys. Rev. D