arxiv: v2 [gr-qc] 29 Oct 2018

Transcription

1 Confomally Fla Collapsing Sas in f (R) gaviy axiv: v [g-qc] 9 Oc 18 Soumya Chakabai Cene fo Theoeical Sudies, Indian Insiue of Technology, Khaagpu, Wes Bengal 71 3, India. Riupano Goswami, Sunil Mahaaj Asophysics and Cosmology Reseach Uni, School of Mahemaics, Saisics and Compue Science, Univesiy of KwaZulu-Naal, Pivae Bag X541, Duban 4, Souh Afica. Naayan Banejee Depamen of Physical Sciences, Indian Insiue of Science Educaion and Reseach, Kolkaa, Mohanpu Campus, Wes Bengal 74146, India. PACS Nos. 4..q; 4..Jb; 4.5.Kd; 4.7.Bw Keywods : gaviaional collapse, exac soluion, spheical symmey, f(r) gaviy Absac The pesen wok includes an analyical invesigaion of a collapsing spheical sa in f(r) gaviy. The ineio of he collapsing sa admis a confomal flaness. Infomaion egading he fae of he collapse is exaced fom he maching condiions of he exinsic cuvaue and he Ricci cuvaue scala acoss he bounday hypesuface of he sa. The adial disibuion of he physical quaniies such as densiy, anisoopic pessue and adial hea flux ae sudied. The inhomogeneiy of he collapsing ineio leads o a non-zeo acceleaion. The divegence of his acceleaion and he loss of enegy hough a hea conducion foces he ae of he collapse o die down and he fomaion of a zeo pope volume singulaiy is ealized only asympoically. 1 Inoducion A modificaion of Einsein s gaviy is one of he mehods employed o explain he wo phases of cosmic acceleaion. Modificaion of he mae seco (fo example, by including exoic scala fields) is ceainly he moe popula appoach o ealize an effecive negaive pessue and heefoe, he acceleaed expansion of he univese [1]. Modified gaviy on he ohe hand modifies he acion of Geneal Theoy of Relaiviy(GR) o wie down an effecive enegy-momenum enso wih a geomeical oigin. Pehaps he mos saighfowad modificaion of GR is undesood by wiing he Einsein-Hilbe lagangian as an analyic funcion f(r) of he Ricci scala R insead of R iself. The so-called f(r) soumya@cs.iikgp.ene.in Goswami@ukzn.ac.za Mahaaj@ukzn.ac.za naayan@iisekol.ac.in 1

2 heoy of gaviy can accoun fo boh he phases of cosmic acceleaion, ealy in he hisoy of he univese and he lae-ime acceleaed expansion. A bief idea, fom he moivaions of sudying he f(r) modificaions o he successes and shocomings of he heoy can be found in he eviews by Clifon, Feeia, Padilla and Skodis [], Soiiou and Faaoni [3], Nojii and Odinsov [4], Nojii, Odinsov and Oikinomou [5]. Many specific choices of f(r) have been pu fowad ove he yeas, o addess quesions posed by inceasingly advanced obsevaional cosmology, o simply fom a mahemaical cuiosiy fo a moe geneal heoy of gaviy. An f(r) model is expeced o descibe a coec cosmological dynamics, wihou compomising he sabiliy of is own. I should ace a pecise newonian and pos-newonian limi and cosmological peubaions well suied wih he cosmic micowave backgound and lage-scale sucue. Alhough hese esicions limi he feedom of choice of he funcional fom of f(r), hee ae viable models pesen in lieaue which aemp o descibe boh he phases of cosmic acceleaion. The specific case of f(r) = R+αR (α > ) caies addiional impoance since his was he vey fis model of inflaion poposed [6]. The pesen wok ages a diffeen poblem of gaviaional physics, ha of Gaviaional Collapse whee hese heoies do pose subsanial difficulies. Gaviaional collapse is defined as he implosion of a massive sella body which can no poduce sufficien inenal pessue o balance ou he inwad pull due o gaviy. Unless foced o sabilize in some kind of dynamical equilibium phase, he collapse can coninue unil a zeo pope volume is eached. A a zeo pope volume he cuvaue invaians, physical quaniies like densiy and pessue ae expeced o divege, poducing a space-ime singulaiy. The pedicabiliy of he collapse eaching a singulaiy a a finie fuue lies wihin he highly non-linea field equaions. Specific case sudies of gaviaional collapse indicae ha he ulimae singulaiy may o may no be hidden fom a disan obseve depending on iniial collapsing pofiles. One defines hese wo possible oucomes of a collapse as, a Black Hole whee he singulaiy is hidden behind a hoizon, o a Naked Singulaiy, whee he singulaiy can communicae wih an obseve hough exchange of phoon. Collapse of a sella objec consising of a fluid wih no pessue was addessed fo he vey fis ime by Da [7] and Oppenheime and Snyde [8]. Vaious aspecs and implicaions of gaviaional collapse was ecenly eviewed by Joshi [9, 1]. Howeve, he numbe of exising collapsing soluions in f(r) gaviy is ahe limied. The mos obvious eason fo his is he fac ha he field equaions of f(r) gaviy ae fouh ode in meic componens and an exac o appoximae soluion is no always guaaneed. The inees in f(r) collapse models has picked up quie ecenly. Possibiliy of ending up in a finie ime cuvaue singulaiy in a sella collapse pocess was invesigaed by Bamba, Nojii and Odinsov [11] and Abuzova and Dolgov [1] fo diffeen f(r) seup. A few ecen invesigaions by Boisov, Jain and Zhang [13] and Guo, Wang and Folov [14] demonsaed he uiliy of numeical simulaions in a sudy of f(r) collapse. The lae used a double-null fomalism in Einsein fame, educing he nonlinea cuvaue conibuion ino a nonminimally coupled scala field by a confomal ansfomaion. Kausa and Noueen [15] sudied he effec of pessue anisoopy and dissipaion on he collapse dynamics fo a specific f(r) model wien as f(r) R + αr n. A case of chaged and expansion-fee f(r) collapse was sudied by Shaif and Yousaf [16]. Vey ecenly, Cembanos, Cuz-Dombiz and Nunez [17] sudied unifomly collapsing cloud of self-gaviaing dus paicles in f(r) gaviy fo a geneal f(r). Fo a complee and consisen analysis of gaviaional collapse, he collapsing ineio mus be mached smoohly wih a pope exeio spaceime. This is a seco whee an f(r) collapse model can pose addiional difficuly. The fouh ode field equaions poduce exa consain equaions which mus be saisfied a he bounday hypesuface as well as he usual Isael-Damois [18, 19] condiions of Geneal Relaiviy. Fo a geneal scala enso heoy, he geneal juncion condiions wee eaed by Baabes and Bessange []. The addiional condiions in f(r) gaviy come fom he maching of Ricci

3 3 scala and is nomal deivaive acoss he bounday. The impoance of hese addiional consains have been invesigaed only ecenly, by Deuelle, Sasaki and Sendouda [1], Senovilla [], Clifon [3], Clifon, Dunsby, Goswami and Nzioki [4], Ganguly, Gannouji, Goswami and Ray [5]. Vey ecenly Goswami, Nzioki, Mahaaj and Ghosh [6], sudied he dynamics of collapsing spheical sas in f(r) gaviy, fo a specific case of Saobinski model (f(r) = R + αr ). Consideing exisance of pessue anisoopy and hea flux, hey invesigaed he bounday maching wih a vacuum exeio soluion and obseved ha he collapsing soluion mimics he Lemaie-Tolman-Bondi soluion of GR [7, 8, 9]. Thei wok pu fowad he impossibiliy of maching a homogeneous collapsing sa in f(r) gaviy o a saic vacuum exeio spaceime acoss a fixed bounday suface. This makes exisence of an Oppenheime-Snyde-Da [7, 8] analogue of dus collapse model, i.e., a collapsing homogeneous dus ball wih a vacuum exeio, inadmissible in f(r) heoies. Simple examples of collapsing model in f(r) gaviy discussed by Chakabai and Banejee [3, 31] confimed his esul given by Goswami e. al. [6]. Theefoe, i emains vey impoan o sudy models of collapsing soluion in f(r) gaviy since many exising models become inconsisen due o he addiional esicions ove he Ricci Scala. The pesen wok deals wih he collapse of a spheical sa in f(r) gaviy addessing he addiional consains on he Ricci scala, wihou any a pioi specific choice of f(r) and a geneal enegy momenum disibuion. The soluion helps one o wie he expessions of densiy, anisoopic pessue and he adial hea flux of he collapsing fluid, explicily. Howeve, he non-lineaiy of he sysem is aken cae of by assuming he spaceime meic o be confomally fla. A vey apid collapsing naue of he spheical sa is obseved, which slows down wih ime. Fo diffeen choices of f(r), enegy condiions ae sudied as well. One salien feaue is ha in spie of his being a simple model, he allowed inhomogeneiy in g beges a nonivial acceleaion, which is no hee in he sandad meic ansaz, such as he Tolman-Bondi meic, which ae nomally aken up fo he collapse models. Secion and 3 defines he acion and he field equaions of f(r) gaviy unde confomal flaness. In secion 4, he ineio spaceime is mached wih a vacuum exeio soluion and he collapsing behavio is discussed. Secion 5 eas he possibiliy of he fomaion of a hoizon. In secion 6 and 7, he ime evoluion and adial disibuion of he physical quaniies ae pesened, along wih he ime evoluion of he enegy condiions, fo diffeen choices of f(r). We conclude he manuscip in secion 8. f(r) gaviy in meic fomalism The acion of f(r) gaviy can be wien as ( ) f(r) g A = 16πG +L m d 4 x, (1) whee L m is he Lagangian fo he mae disibuion. The vaiaion of he acion (1) wih espec o he meic enso leads o he following fouh ode paial diffeenial equaion, F(R)R µν 1 f(r)g µν µ ν F(R)+g µν F(R) = 8πGT m µν, () which defines he field equaions in f(r) gaviy, whee F(R) = df dr. Wiing his equaion in he fom of Einsein enso, one obains G µν = κ F (Tm µν +TC µν ), (3) whee Tµν C = 1 ( ) f(r) RF(R) g µν + µ ν F(R) g µν F(R). (4) κ

4 4 T C µν is echnically he cuvaue conibuion o he effecive enegy-momenum enso. The effecive enegy-momenum disibuion may o may no obey all he enegy condiions like an odinay fluid disibuion. 3 Confomal flaness and field equaions We assume ha he space-ime meic has a vanishing Weyl enso which indicaes ha i is confomally fla. Confomally fla adiaing fluid sphees and shea-fee adiaing sas ae vey well-sudied in lieaue [3]. The meic can be wien as ds = D(,) [ d d dω ], (5) whee D(, ) is he confomal faco and govens he evoluion of he sphee. The fluid inside he spheically symmeic body is assumed o be locally anisoopic along wih he pesence of hea flux. Thus he enegy-momenum enso is given by T µν = (ρ+p )u µ u ν p g µν +(p p )χ µ χ ν +q µ u ν +q ν u µ, (6) whee q µ = (,q,,) is he adially dieced hea flux veco, ρ is he enegy densiy, p he angenial pessue, p he adial pessue, u µ he fou-velociy of he fluid and χ µ is he uni space-like fou-veco along he adial diecion. The vecos ae nomalized as u µ u µ = 1, χ µ χ µ = 1, χ µ u µ =. (7) Inclusion of a non-zeo pessue anisoopy may play a cucial ole while dynamics of a massive collapsing sa is consideed. Fo discussions on hese issues, we efe he eade o he wok of Heea and Sanos [33], Heea e al. [34]. I was pecisely poved by Heea and Ponce de Leon ha unde he assumpion of confomal symmey, a smooh bounday maching wih he exeio geomey is only possible if he fluid has a pessue anisoopy [35]. A dissipaive pocess in he fom of hea conducion is hee in he sysem as well. Duing a gaviaional collapse, one expecs he collapsing sa o become inceasingly compac ove ime. As he dimension of he componen paicles and he mean fee pah becomes moe and moe compaable, a dissipaive pocess can help emi enegy o sele he sysem down o an equilibium, o slow down he collapse pocess such ha he zeo pope volume is eached only asympoically. Fo moe elevan deails we efe o he wok of Heea e. al. [34]. The fou-acceleaion is defined as a µ = u µ;ν u ν and he expansion scala as Θ = u β ;β. Fo a comoving line elemen (5), u µ = δµ D(,). (8) Using he equaions (5) and (6), he field equaions (, 3) can be wien fo he pesen case as 3Ḋ [ D D D 4D D + D ρ D = D F + f F R ] [ F + F + Ḋ F DF + ( + D D ) ] F, (9) F D [ D + Ḋ D +4D D +3D D = p D F f F + R ] [ F + F 3Ḋ F ] F DF 3D DF F, (1) F D [ D + Ḋ D +D D +D D D D = p D F f F + R ] [ F + F F F Ḋ F DF D F DF F F ], (11)

5 5 and Ḋ D +4ḊD D = qd3 F + F F ḊF DF D F DF. (1) Thee ae fou independen equaions bu five unknowns, he componens of he enegy momenum enso ρ(,), p (,), p (,), q(,) and he confomal faco D(,). The saighfowad way of sudying soluion of he field equaions is o assume some fom of esicion ove he enegy momenum enso componens, fo example, pessue isoopy o vanishing pessue o a vanishing hea flux. Insead of specifying such an equaion of sae, we aemp o exac elevan infomaion egading he dynamics fom he maching condiions of he collapsing sphee wih a vacuum exeio descibed by he Schwazschild soluion. 4 Maching wih a vacuum exeio soluion and an exac soluion If one assumes ha he spheical collapsing sa is suounded by vacuum hen he exeio spaceime is well descibed by a Schwazschild geomey. Maching of he ineio wih a vacuum Schwazschild exeio a a bounday hypesuface equies he coninuiy of boh he meic and he exinsic cuvaue on he bounday hypesuface, discussed exensively by Sanos [36] and Chan[37]. Schwazschild soluion is given by ( ds = 1 M ) ( dt 1 M ) 1d R (dθ +sin θdφ ), (13) R R whee M is he oal mass conained by he ineio. Fo an ineio meic given by ds = S (,)d N (,)d C (,)(dθ +sin θdφ ), (14) he maching wih he meic (13) yields (using he maching of he second fundamenal fom o he exinsic cuvaue) ( C C ĊS CS NC ) = Σ N [ ( C NC S C Ṡ S Ċ ) ] Ċ + S [ (S ) C C N S + C C ( ) ] N C C, (15) C whee Σ is he bounday. Also hee we have used R = Σ C, Sd = Σ ( 1 M R ) dt (16) The Misne and Shap mass funcion [38], defined as m(,) = C (1+gµν C,µ C,ν ), (17) yields he mass conained by he suface defined by he adial coodinae. Fo a confomally fla meic defined by (5), he condiion (15) can be simplified ino D D +Ḋ D Ḋ D 3D D 4ḊD D 4D D =Σ. (18) Wiing he equaion (18) in ems of A(,) = 1 D(,), one ges Ä A +Ȧ A 3Ȧ A 3A A + A A =Σ. (19)

6 6 Le us assume ha he bounday of he collapsing sphee is a a adial disance given by = b. We also pu in he ansaz ha A(,) = η()+ξ(); η() govens he ime evoluion of he collapsing sphee and ξ() is a funcion of only, defined on he bounday hypesuface as ξ( b ) = ξ, () and ξ ( b ) = ξ 1, (1) ξ ( b ) = ξ. () Hee a pime denoes a deivaive wih espec o and a do epesens a deivaive wih espec o. Wih hese assumpions, he em Ȧ = and he equaion (19) simplifies consideably. A = b, i can be wien as η (η +ξ ) 3 η (η +ξ ) 3ξ 1 b (η +ξ ) + ξ1 =. (3) (η +ξ ) Equaion (3) can be inegaed o o yield a fis inegal as η (η +ξ ) 3 + 3ξ 1 b (η +ξ ) ξ 1 3 (η +ξ ) 3 = λ. (4) Hee, λ is a consan of inegaion. Fo all non-zeo values of λ, i is difficul o inegae he equaion (4)and wie η() in a closed fom. We discuss he soluion fo he simples case, i.e., fo λ =. Since we have assumed η() invesely popoional o he scale faco (Ḋ = Ȧ A ), we noe ha fo a collapsing soluion, η mus be posiive all he ime. Thus fo η >, equaion (4) is inegaed o yield an exac soluion wien as η() = ( ξ1 b 9 ξ ) 3ξ 1 8 b ( ), (5) whee is a consan of inegaion. Theefoe he confomal faco D(,) can be wien as D(,) = 1 [( ) ]. (6) ξ1 b 9 ξ +ξ() 3ξ 1 8 b ( ) Fom equaion (6), i appeas ha he collapsing sphee, eaches zeo pope volume only asympoically in comoving ime. Howeve he pope ime as calculaed by τ = D(,)d emains finie fo any obseve locaed a a given =. Invesigaing he coninuiy of he Ricci scala and is nomal deivaives acoss he bounday hypesuface, we can commen on he naue of ξ(). The maching cieion of he Ricci Scala was pesened in deails by Deuelle, Sasaki and Sendouda [1]. Fo moe ecen examples and discussions we efe o he woks of Clifon, Dunsby, Goswami and Nzioki [4], and Goswami, Nzioki, Mahaaj and Ghosh [6]. Fom he meic (5) and he exac soluion (6), we wie he Ricci scala as Hee, he em [ 9ξ 1ξ [ 9ξ1 ξ R = 81 8 b 3 8 b 81 ξ1 3 b ξ1 b + 3ξ 1ξ ][ ] ( )+Λ 1 (,ξ,ξ 1,ξ,, b ). (7) b + 3ξ 1ξ b ] govens he adial dependence of he scala and ime dependence is given by he em [( ) + Λ 1 (,ξ,ξ 1,ξ,, b )]. The em Λ 1 (,ξ,ξ 1,ξ,, b ) woks as a

7 7 paamee consising of ξ, ξ, ξ, ξ, ξ 1, ξ, and b which we wie in his fom fo he sake of beviy. Fo he coninuiy, he icci scala and i s spaial deivaive mus mach a he bounday suface i.e., Then, he Ricci scala can be wien in a fom like Compaing equaions (7) and (9), we wie R b = R b =. (8) R = ( b ) g(,). (9) 9ξ 1 ξ 81 8 b 3 ξ1 b + 3ξ 1ξ b = ( b ). (3) Wiing his equaion as a second ode non-linea diffeenial equaion of ξ(), we ge ξ = 8 b 9ξ 1 ( b ) 4ξ 3 + 9ξ 1 4 b. (31) The evoluion of ξ() can be wien fom his equaion in an explici fom as ξ() = 4 b 6 43 b ( 8 5 b 171ξ 1 39ξ ξ ) 1 3k k. (3) 9ξ 1 4 b Hee, k 1 and k ae abiay consans of inegaion. Diffeeniaing equaion (3) and puing = b gives a quadaic in ξ 1 as and heefoe ξ 1 can be wien as ξ 1 8k b ξ b =, (33) ξ 1 = 14k b ± To ensue ha ξ 1 is eal, one mus ensue ha k b 7 196k b b. (34) >, (35) which is valid fo any eal choice of k 1. Fom equaion (6), i is saighfowad o see ha fo a collapsing scenaio (Ḋ < ), ξ 1 mus be negaive. Fom equaion (34) i can be seen ha his pus addiional esicions ove k 1. We assume ha N 1 and N ae he negaive and he posiive oos of he em unde squae oo, especively. Theeafe he necessay condiions ha ensue a collapsing evoluion ae k 1 < N b, (36) and k 1 > N b. (37) The ime evoluion of he spheical disibuion is pesened gaphically in Figue 1 fo diffeen values of. I appeas ha, afe an iniial phase of expansion (depending on he consan of inegaion, which can be pu o wihou any loss of geneliy), he sphee eaches a maximum accessible volume. Theeafe a apid collapsing naue is pominen. Evenually he apidiy of he collapse diminishes wih

8 A(, ) 1. 1 A(, ) Figue 1: Time-evoluion of he scale faco a diffeen adial disance. The yellow gaph shows he evoluion a = b. The ed gaph shows he evoluion a = b and he blue gaph = b 1. The gaph on he op indicaes an iniial bief peiod of expansion befoe he collapse begins. ime and a zeo pope volume is eached only asympoically. Diffeen cuves show he naue of he collapse fo diffeen shells labelled by diffeen values of he adial coodinae. The shells sa collapsing a a diffeen ae, as in he innemos shell sas collapsing moe apidly, followed by subsequen shells of he fluid. Howeve, wih ime he shells come close and close ogehe and he cuves denoing hei collapse follow a moe-o-less simila ajecoy. The yellow gaph (lowemos) shows he evoluion a = b. The puple gaph (middle) shows he evoluion a = b while he blue gaph (op) is fo = b 1. and Now, using (8) and he meic (5), he fou-acceleaion and is magniude(scala) can be wien as a 1 = D (,) D(,) = ( ) ( ( ) ), (38) a µ a µ = D (,) D(,) 4 = ( ). (39) We noe ha he magniude of fou-acceleaion is a funcion of adial disance fom he cene only. We plo a µ a µ as a funcion of in Figue. 5 Evoluion of Densiy, Radial Pessue, Tangenial Pessue and Hea Flux Since we have no assumed any specific equaion of sae, and he exac soluion is obained saighaway fom a smooh bounday maching, he geneal expession of physical quaniies like densiy, adial and angenial pessue, hea flux, can be wien in ems of he meic componens fom he field equaions (9), (1), (11) and (1). These feaues ae all quie geneic, as no specific funcional fom of f(r) is assumed as ye. In his secion, we pesen a gaphical analysis of physical quaniies duing he collapse fo hee popula and cosmologically impoan classes of f(r) models; f(r) = R, f(r) = e αr and f(r) = R+αR, fo diffeen choices of he iniial paamees. The value of he paamees govens naue of he evoluion, howeve, he main qualiaive naue emains simila. Fo f(r) = R+αR, i is necessay

9 9 1 1 a Α aα Figue : Evoluion of he magniude of fou-acceleaion as a funcion of o assume α o be small, ohewise he collapsing disibuion doesn mach popely wih a Scwazschild soluion a he bounday. We have woked wih a fixed bounday hypesuface a b = 1 and α 1 3 o lowe, fo f(r) = R+αR, and wih α.16 fo f(r) = e αr. We keep ohe paamees fixed and plo he physical quaniies as a funcion of fo diffeen values of ime. Saing fom = 1., we plo upo = 11. a equal inevals of. unis of ime. We wie he expession of densiy as ρ = [( 3 Ḋ D D D 4D D + D D F F Ḋ F DF F F D F DF ) ( 1 f D F R )] F. (4) In Figue 3, we plo he disibuion of densiy of he collapsing sphee as a funcion of he adial disance fo diffeen epoch. The lowemos cuve shows he disibuion when = 1. and he uppemos one is fo = 11.. Gadually, he densiy inceases wih adial disance deceasing, being he maximum a he cene ( = ). Theefoe, i is clea ha he collapsing mae keeps on piling up a he cene and he coe densiy inceases wih ime. The densiy is sicly posiive houghou he collapse and i becomes zeo aound he bounday hypesuface. Radial pessue is govened by he equaion D D + Ḋ D +4D D +3D D = D [ p F f F + R ] [ F + F 3Ḋ F ] F DF 3D DF F. (41) F We sudy he disibuion of adial pessue of he collapsing sphee as a funcion of he adial disance in Figue 4. The adial pessue is sicly posiive houghou he collapse. Gadually, he densiy inceases wih adial disance deceasing, being he maximum a he cene ( = ). I becomes zeo a he bounday of he sa and inceases a he cene, i.e. = as ime goes. Diffeen cuves in a single gaph denoes evoluion fo diffeen values of ime. The lowemos cuve shows he disibuion when = 1. and he uppemos one is fo = 11.. One impoan poin o noe hee, a, he adial pessue aains a vey high bu a finie value, which is ou of ange ( 1 1 o highe) of he gaphs in Figue 4. The evoluion of angenial pessue is govened by he equaion D D + Ḋ D +D D +D D D D = D [ p F f F + R ] + [ F F F F Ḋ F ] DF D F DF F, (4) F We plo he angenial pessue as a funcion of in Figue 5. The angenial pessue is negaive houghou he collapse fo R + αr and R. Diffeen cuves in a single gaph denoes evoluion fo

10 ρ 3 ρ ρ Figue 3: Evoluion of he densiy wih adial disance; = o = 1 fo diffeen values of ime fo f(r) = R+αR (op lef), f(r) = R (op igh) and fo f(r) = e αr (boom). Diffeen cuves in a single gaph denoes evoluion fo diffeen values of ime; he lowemos cuve denoes evoluion fo = 1., he cuve second fom he boom denoes evoluion fo = 1. and so on. The opmos cuve denoes evoluion fo = 11.. diffeen values of ime. The uppemos cuve shows he disibuion when = 1. and he lowemos one is fo = 11.. The angenial pessue becomes zeo a he bounday of he sa and has he maximum value as one appoaches he cene, i.e. =. Fo e αr, as he lowe gaph suggess, he evoluion wih adial disance is negaive nea he cene of he collapsing sa and zeo a he bounday. Howeve, as one appoaches he bounday hypesuface fom he cene, he evoluion becomes posiive fom a negaive pofile. In his case as well, aound, he angenial pessue aains a vey high bu a finie value, ou of ange of he gaphs, as is clea fom he lowemos gaph of in Figue 5. Ḋ D +4ḊD D = qd3 F + F F ḊF DF D F DF. (43) In Figue 6, he evoluion of hea flux is ploed as a funcion of fo diffeen values of ime. The hea flux mainains a negaive pofile houghou he collapse. I vanishes boh a he cene and a he bounday hypesuface. The adial disance fom he cene a which he hea flux is maximum depends on he choice of f(r). Diffeen cuves in a single gaph denoes evoluion fo diffeen values of ime. The uppemos cuve shows he disibuion when = 1. and he lowemos one is fo = Fomaion of a hoizon The condiion of he fomaion of an appaen hoizon is given by he equaion g µν Y,µ Y,ν =, (44)

11 p_ p_ p_ Figue 4: Evoluion of adial pessue wih adial disance; = o = 1 fo diffeen values of ime fo f(r) = R + αr (op lef), f(r) = R (op igh) and fo f(r) = e αr (boom). Diffeen cuves in a single gaph denoes evoluion fo diffeen values of ime; he lowemos cuve denoes evoluion fo = 1., he cuve second fom he boom denoes evoluion fo = 1. and so on. The opmos cuve denoes evoluion fo = 11.. whee Y(,) is he adius of he wo-sphee, given by D(,) in he pesen case. Fom he evoluion of he scale faco given by equaion (6) one can ewie his condiion as The ime of fomaion of an appaen hoizon is hen given by η (ξ +η) ξ +ξ (ξ +η) =. (45) [ 8b ( ap = ± ξ ξ + ξ )]1 1 b ξ. (46) 3ξ 1 9 Since a collapsing scenaio (Ḋ < ) is only ensued when ξ 1 is negaive, and b is always posiive, a fomaion of hoizon is possible if and only if (ξ ξ + ξ 1 b 9 ξ ) >. ξ is given by ξ() = 4 b 6 43 b ( 8 5 b 171ξ 1 39ξ ξ ) 1 3k k. (47) 9ξ 1 4 b If his condiion is saisfied, we acually have an example of a hoizon wihou any singulaiy o cove, since he zeo pope volume singulaiy is eached only asympoically in comoving ime, wheeas he hoizon foms in a finie comoving ime. 7 Evoluion of he enegy condiions In ode ha a collapsing fluid is physical, i mus obey some enegy condiions which we biefly discuss below.

12 p_ -15 p_ p_ Figue 5: Evoluion of angenial pessue wih adial disance = o = 1 fo diffeen values of ime fo f(r) = R + αr (op lef), f(r) = R (op igh) and fo f(r) = e αr (boom). Diffeen cuves in a single gaph denoes evoluion fo diffeen values of ime; he uppemos cuve denoes evoluion fo = 1., he cuve second fom he op denoes evoluion fo = 1. and so on. The lowemos cuve denoes evoluion fo = The Null Enegy Condiion holds if fo all null vecos T µν k µ k ν (48). The weak enegy condiion holds if fo any non-spacelike veco w α, T αβ w α w β (49) This is basically equivalen o saying ha he enegy densiy of he fluid is non-negaive. 3. The dominan enegy condiion holds good if fo a imelike veco w α, T αβ w β is imelike o null. This can be inepeed as a condiion such ha he speed of enegy flow of mae is less han he speed of ligh fo any obseve. 4. The song enegy condiion is saisfied if fo any imelike uni veco w α T αβ w α w β +T, (5) whee T is he ace of he enegy momenum enso. The song enegy condiion can only be violaed if he enegy densiy is negaive o if hee exiss a lage negaive pessue componen of he enegy momenum enso. The enegy condiions wee invesigaed fo an impefec fluid in igoous deails by Kolassis, Sanos and Tsoubelis [39], and Pimenel, Loa-Clavijo and Gonzalez [4] quie ecenly. Algebaically, he enegy condiions o be saisfied ae as follows,

13 q q q Figue 6: Evoluion of hea flux wih adial disance = o = 1 fo diffeen values of ime fo f(r) = R+αR (op lef), f(r) = R (op igh) and fo f(r) = e αr (boom). Diffeen cuves in a single gaph denoes evoluion fo diffeen values of ime; he uppemos cuve denoes evoluion fo = 1., he gaph second fom he op denoes evoluion fo = 1. and so on. The lowemos gaph denoes evoluion fo = (null enegy condiions; N EC1 and N EC) ρ+p q, (51) ρ p +p +, (5). weak enegy condiions (WEC) 3. dominan enegy condiions (DEC1 and DEC) ρ p +, (53) ρ p, (54) ρ p p +, (55) 4. song enegy condiions (SEC) p +, (56) whee = (ρ+p ) 4q. We plo he evoluion of he enegy condiion wih espec o ime fo diffeen values of adial coodinae, anging fom he cene a = o he bounday a = 1.

14 {Δ-p_+p_+ρ} {Δ-p_+ρ} {ρ-p_} {Δ-p_-p_+ρ} {Δ+p_} {p_-q+ρ} Figue 7: Evoluion of he enegy condiions wih ime, wih diffeen cuves in a single gaph denoing evoluion fo diffeen adial disance fom he cene fo f(r) = R+αR. Evoluion of NEC (op lef gaph), WEC (op igh gaph), DEC1 (middle lef gaph), DEC (middle igh gaph) and SEC (boom lef gaph) suggess ha hese enegy condiions ae valid houghou he collapse. Fo each of hese gaphs, he lowemos cuve denoes he evoluion when and he opmos cuve denoes he evoluion when b. The NEC1 is given in he boom igh gaph as a funcion of adial disance which is clealy violaed houghou he couse of he collapse. We now invesigae he evoluion of enegy condiions fo a few choice of f(r). We include a model wih f(r) = R+αR in ou sudy. As discussed by Clifon, Feeia, Padilla and Skodis, such an f(r) model descibes well he empeaue anisoopies in Cosmic Micowave Backgound [41]. Pesence of a em αr seves even moe pupose, fo insance, mae insabiliy and difficulies egading local gaviy consains in models like f(r) R α R n can be solved by inclusion of an R em, as discussed by Paul, Debnah and Ghoshe [4]. Addiion of an R em ino he acion can also lead o consisen f(r) model ha passes he sola sysem ess, bypassing insabiliy poblems as shown by Nojii and Odinsov [43, 44]. Since we ae sudying a scenaio of gaviaional collapse ahe han a cosmological one, we also sudy he evoluions fo a couple of f(r) models which gow apidly wih R a high cuvaue (f(r) = R and f(r) = e αr ). We also sudy a case whee f(r) R 1+δ, whee δ is a peassigned small quaniy. Clifon and Baow sudied he weak-field popeies of such a model and found ha i is compaible wih local asonomical ess (e.g. peihelion shif), iff he f(r) lagangian is vey close o Geneal Relaiviy (in f(r) R (1+δ), δ mus be of he ode of 1 19 o less [45]). We make a noe hee ha a choice of f(r) R is aguable since his heoy does no have a pope Newonian limi as discussed by Pechlane and Sexl [46] in hei discussion on quadaic lagangians

15 {Δ-p_+p_+ρ} {Δ-p_+ρ} {ρ-p_} {Δ-p_-p_+ρ} {Δ+p_} {p_-q+ρ} Figue 8: Evoluion of he enegy condiions wih ime, wih diffeen cuves in a single gaph denoing evoluion fo diffeen adial disance fom he cene fo f(r) = R. Evoluion of NEC (op lef gaph), WEC (op igh gaph), DEC1 (middle lef gaph), DEC (middle igh gaph) and SEC (boom lef gaph) suggess ha hese enegy condiions ae valid houghou he collapse. Fo each of hese gaphs, he lowemos cuve denoes he evoluion when and he opmos cuve denoes he evoluion when b. The NEC1 is given in he boom igh gaph as a funcion of adial disance which is clealy violaed houghou he couse of he collapse. in geneal elaiviy. Moeove, an f(r) model is compaible wih local asonomical ess only if he lagangian is vey close o geneal elaiviy (fo f(r) R (1+δ), δ can be amos of he ode of 1 19 [45]). Howeve, we keep he case of f(r) R aleas fo he sake of illusaion. 1.f(R) = R+αR Evoluion of he enegy condiions wih ime is shown in figue 7 fo f(r) = R+αR wih α 1 3. Diffeen cuves in a single gaph denoes he evoluion fo diffeen adial disance fom he cene. Evoluion of NEC (op lef gaph), WEC (op igh gaph), DEC1 (middle lef gaph), DEC (middle igh gaph) and SEC (boom lef gaph) shows ha hey ae posiive houghou he collapse and heefoe, holds ue. Fo each of hese gaphs, he lowemos cuve denoes he evoluion nea he cene of he collapsing sa ( ) and he opmos cuve denoes he evoluion close o he bounday ( b ). The NEC1 is given in he boom igh gaph as a funcion of adial disance which is negaive and is clealy violaed houghou he couse of he collapse..f(r) = R

16 {Δ-p_+p_+ρ} {Δ-p_+ρ} {Δ-p_-p_+ρ} {Δ+p_} Figue 9: Evoluion of he enegy condiions wih ime, wih diffeen cuves in a single gaph denoing evoluion fo diffeen adial disance fom he cene fo f(r) = e αr. Evoluion of NEC (op lef gaph), WEC (op igh gaph), DEC (boom lef gaph) and SEC (boom igh gaph) suggess ha hese enegy condiions ae valid houghou he collapse. Fo each of hese gaphs, he lowemos cuve denoes he evoluion when and he opmos cuve denoes he evoluion when b. Evoluion of he enegy condiions wih ime is shown in figue 8 fo f(r) = R. Diffeen cuves in a single gaph denoes he evoluion fo diffeen adial disance fom he cene. Evoluion of N EC (op lef gaph), W EC (op igh gaph), DEC1 (middle lef gaph), DEC (middle igh gaph) and SEC (boom lef gaph) shows ha hey ae posiive houghou he collapse and heefoe, holds ue. Fo each of hese gaphs, he lowemos cuve denoes he evoluion nea he cene of he collapsing sa ( ) and he opmos cuve denoes he evoluion close o he bounday ( b ). The NEC1 is given in he boom igh gaph as a funcion of adial disance which is negaive and is clealy violaed houghou he couse of he collapse. 3.f(R) = e αr Evoluion of he enegy condiions wih ime fo f(r) = e αr ae ploed in Figue 9 and 1. Diffeen cuves in a single gaph denoes he ime evoluion fo diffeen adial disance fom he cene. Evoluion of NEC (op lef gaph of Figue 9), WEC (op igh gaph of Figue 9), DEC (boom lef gaph of Figue 9) and SEC (boom igh gaph of Figue 9) suggess ha hese enegy condiions ae sicly posiive and valid houghou he collapse. The lowemos cuve fo each of hese gaphs denoes he evoluion when and he opmos cuve denoes he evoluion when b. Fo f(r) = e αr, he evoluion of NEC1 is shown on he lef gaph of Figue 1. The lowemos cuve denoes he evoluion when, i.e., close o he cene, whee he enegy condiion is clealy violaed. Howeve, he evoluion vey close o he bounday is denoed by he opmos cuve which mainains a posiive pofile houghou he collapse, i.e., he enegy condiion is valid close o he bounday hypesuface b. A simila evoluion is seen fo he evoluion of DEC1 which is ploed on he igh gaph of Figue 1. One can infe fom his esul ha even if all he enegy condiions ae well-saisfied a he bounday suface of he dying sa, i does no always guaanee he validiy of all he enegy

17 17 condiions in as one appoaches he cene {p_-q+ρ} {ρ-p_} Figue 1: Evoluion of he enegy condiions wih ime, wih diffeen cuves in a single gaph denoing evoluion fo diffeen adial disance fom he cene fo f(r) = e αr. Evoluion of NEC1 is shown on he lef gaph. The lowemos cuve denoes he evoluion when, i.e., close o he cene, whee he enegy condiion is clealy violaed. Howeve, he evoluion vey close o he bounday is denoed by he opmos cuve which mainains a posiiove pofile houghou he collapse, i.e., he enegy condiion is valid close o he bounday hypesuface b. A simila evoluion is seen fo he evoluion of DEC1 which is ploed on he igh gaph. 4.f(R) = R (1+δ) Evoluion of he enegy condiions wih ime is shown in figue 11 fo an f(r) model fo which he depaue of he f(r) modificaion fom Einsein s gaviy is vey small. This is given by f(r) = R (1+δ), whee δ is a vey small posiive numbe, 1 19 o lowe. Diffeen cuves in a single gaph denoes he evoluion fo diffeen adial disance fom he cene. Evoluion of NEC (op lef gaph), WEC (op igh gaph), DEC1 (middle lef gaph), DEC (middle igh gaph) and SEC (boom lef gaph) shows ha hey ae posiive houghou he collapse and heefoe, holds ue. Fo each of hese gaphs, he lowemos cuve denoes he evoluion nea he cene of he collapsing sa ( ) and he opmos cuve denoes he evoluion close o he bounday ( b ). The NEC1 is given in he boom igh gaph as a funcion of adial disance which is negaive and is clealy violaed houghou he couse of he collapse. These behavio is in fac, qualiaively simila as ha of an R model. The enegy condiions ae acually some consains o be followed such ha he noion of locally posiive enegy densiy holds eveywhee. In he pesen case, some ae indeed violaed. Bu one mus ealize ha his is no eally a seious infingemen on he physical plausibiliy of he sysem. Fis of all, he posiiviy of enegy is neve violaed, so he sysem is quie well posed fo a quanizaion scheme. Secondly, he enegy condiions ha we deived ae no fo he fluid alone, i akes cae of he conibuion of cuvaue o he effecive sess-enegy enso as well. The violaion of some of he enegy condiions can be aibued o he conibuion fom cuvaue, which does no eally have he onus of saisfying he enegy condiions. Thee ae los of examples of quie easonable classical sysems, which can violae some of he enegy condiions. The simples example is a massive scala field non-minimally coupled o gaviy, as discussed by Visse and Bacelo [47, 48]. One ineesing oucome of he soluions of gaviaional field equaions ha disobey enegy condiions is ha hey can give ise o Loenzian womholes as discussed by Bacelo [48], Flanagan and Wald [49].

18 {Δ-p_+p_+ρ} {Δ-p_+ρ} {ρ-p_} {Δ-p_-p_+ρ} {Δ+p_} {p_-q+ρ} Figue 11: Evoluion of he enegy condiions wih ime, wih diffeen cuves in a single gaph denoing evoluion fo diffeen adial disance fom he cene fo f(r) = R (1+δ). Evoluion of NEC (op lef gaph), WEC (op igh gaph), DEC1 (middle lef gaph), DEC (middle igh gaph) and SEC (boom lef gaph) suggess ha hese enegy condiions ae valid houghou he collapse. Fo each of hese gaphs, he lowemos cuve denoes he evoluion when and he opmos cuve denoes he evoluion when b. The NEC1 is given in he boom igh gaph as a funcion of adial disance which is clealy violaed houghou he couse of he collapse. 8 Conclusion We invesigaed a collapsing spheical sa in a vey geneal seup of f(r) modificaions of gaviy, fo a confomally fla spaceime. I was discussed quie ecenly ha he smooh maching of a collapsing ineio wih a vacuum exeio is no a all a saighfowad ask in he conex of f(r) heoies [6]. The exa maching condiions on he Ricci scala enfoce song esicions on he allowed soluions. The pesen wok uilizes hese condiions equied o close he sysem of equaions. I deseves menion ha some classes of collapsing models allowed in geneal elaiviy ae in fac uled ou in f(r) gaviy models [6]. The soluions pesened hee saisfy all he bounday maching condiions and hey descibe he ineio of a collapsing sa wih pessue anisoopy and hea flux, wihou any assumpion of equaion of sae o any paicula choice of f(r) a he ouse. I is quie ineesing o noe ha he ineio disibuion needs o be inhomogeneous fo he pope maching a he bounday. Howeve, whehe he mae disibuion of he sa saisfies all he easonable enegy condiions o no, depends on he choice of funcional fom of f(r). We noe hee ha he sysem unde consideaion in he manuscip is chosen o be confomally fla

19 19 a he ouse. This is indeed a vey special siuaion, in he sense ha hee emains only one scala gaviaional degee of feedom, namely he confomal faco D. This migh have been esponsible fo a loss of genealiy in he collapsing evoluion. Howeve, he high non-lineaiy of he fouh ode field equaions infac compels one o eso o a simplified sysem which can descibe he oveall physics noneheless. The esuls in he manuscip heefoe may no epesen a complee geneal spheical gaviaional collapse, bu ae physically elevan and ineesing as an example iself. We also noe ha a Schwazschild soluion can be an exac soluion fo f(r) heoies only when he cuvaue becomes consan. Fom he consan cuvae condiion Rf (R) = f(r), one can ague ha he exponenial f(r) cases may no poduce a consan cuvaue soluion a all. This is no oo unexpeced eihe, since in he pesence of a highe cuvaue em in he lagangian of gaviy he smooh maching of he ineio wih an exeio Schwazschild meic may be spoil someimes, as discussed by Banejee and Paul [5], Chakabai vey ecenly [51]. The collapsing sa eaches a zeo pope volume only asympoically, since he ae of collapse dies down wih ime. Theefoe, whehe he collapse ends up in a singulaiy o no becomes acually ielevan. Recen sudies of confomally fla collapsing soluions have shed ligh on many diffeen possible oucomes of a gaviaional collapse. Fo insance, a confomally fla collapsing pefec fluid was sudied by Hamid, Goswami and Mahaaj [5] quie ecenly. I was poved ha if he song enegy condiion is obeyed, he collapse necessaily leads o a black hole end-sae. A spaially homogeneous massive scala field collapse also leads o a black hole end-sae unde confomal flaness, as discussed by Chakabai and Banejee [53]. I was also poved ha in he pesence of pessue anisoopy and a dissipaion, he spaceime singulaiy can emain naked [54]. In he pesen case, howeve, i is obvious ha fo a simila collapsing body in f(r) gaviy, he end sae may no end up in eihe a black hole o a naked singulaiy afeall. Thee ae quie a few examples of specific f(r) gaviy models whee a spheical collapse o a singulaiy a a finie fuue is possible [3, 31]. In he pesen case ineio inhomogeneiy leads o an acceleaion unlike he sandad Tolman-Bondi collapsing soluions in geneal elaiviy. The divegence of his acceleaion and he dissipaive effecs like he hea conducion should be he key o his infinie delay in he fomaion of he singulaiy as expeced fom he Raychaudhui equaion [55]. 9 Acknowledgemens N.B. wishes o hank he Asophysics and Cosmology Reseach Uni a Univesiy of Kwazulu-Naal fo a wam hospialiy, whee he wok was iniiaed. S.C. wishes o hank Pof. Sayan Ka fo useful discussions. S.C. was suppoed by he Naional Pos-Docoal Fellowship (file numbe: PDF/17/75) fom he Science and Engineeing Reseach Boad (SERB), Govenmen of India. Refeences [1] A. A. Saobinsky, Phys. Le. B. 91, 99 (198); D. Kazanas, Asophys. J. 41, L59 (198); A. H. Guh, Phys. Rev. D. 3, 347 (1981); K. Sao, Mon. No. Royal Ason. Soc. 195, 467 (1981); D. H. Lyh and A. Rioo, Phys. Rep. 314, 1 (1999); B. A. Basse, S. Tsujikawa, and D. Wands, Rev. Mod. Phys. 78, 537 (6). [] T. Clifon, P. G. Feeia, A. Padilla, and C. Skodis, Phys. Rep. 513, 1 (1). [3] T. P. Soiiou and V. Faaoni, Rev. Mod. Phys. 8, 451 (1). [4] S. Nojii and S. D. Odinsov, Phys. Rep (11).

20 [5] S. Nojii, S. D. Odinsov and V. K. Oikonomou, Phys. Rep (17). [6] A. A. Saobinsky, Phys. Le. B. 91, 99 (198). [7] B. Da, Z. Phys. 18, 314 (1938). [8] J. R. Oppenheime and H. S. Snyde, Phys. Rev. 56, 455 (1939). [9] P. S. Joshi, Global Aspecs in Gaviaion and Cosmology, Claendon Pess, Oxfod Univesiy Pess (1993). [1] P. S. Joshi, Pamana 55, 59 (). [11] K. Bamba, S. Nojii, and S. D. Odinsov, Phys. Le. B. 698, 451 (11). [1] E. V. Abuzova and A. D. Dolgov, Phys. Le. B. 7, 89 (11). [13] A. Boisov, B. Jain, and P. Zhang, Phys. Rev. D. 85, (1). [14] J. Guo, D. Wang, and A. V. Folov, Phys. Rev. D. 9, 417 (14). [15] H. R. Kausa and I. Noueen, Eu. Phys. J. C. 74, 76 (14). [16] M. Shaif and Z. Yousaf, Eu. Phys. J. C. 73, 633 (13); M. Shaif and Z. Yousaf, Phys. Rev. D 88, 4 (13). [17] J. A. R. Cembanos, A. d. l. Cuz-Dombiz, and B. Mones-Nunez, JCAP 14, 1 (1). [18] W. Isael, Nuovo Cim. B. 44, 1 (1966). [19] J. G. Damois, Memoial des Sciences Mahemaiques, Gauhie-Villas, Pais 5, 1 (197). [] C. Baabes and G. F. Bessange, Class. Quan. Gav. 14 : 85 (1997). [1] N. Deuelle, M. Sasaki, and Y. Sendouda, Pog. Theo. Phys. 119, 37 (8). [] J. M. M. Senovilla, Phys. Rev. D. 88, 6415 (13). [3] T. Clifon, Class. Quan. Gav. 3, 7445 (6). [4] T. Clifon, P. K. S. Dunsby, R. Goswami, and A. M. Nzioki, Phys. Rev. D. 87, (13). [5] A. Ganguly, R. Gannouji, R. Goswami, and S. Ray, Phys. Rev. D. 89, 6419 (14). [6] R. Goswami, A. M. Nzioki, S. D. Mahaaj, and S. G. Ghosh, Phys. Rev. D. 9, 8411 (14). [7] G. Lemaie, Ann. Soc. Sci. Buxelles p. 51 (1933). [8] R. C. Tolman, Poc. Nal. Acad. Sci. USA., 169 (1934). [9] H. Bondi, Mon. No. R. Ason. Soc. 17, 41 (1948). [3] S. Chakabai and N. Banejee, Gen. Relaiv. Gavi. 48 : 57 (16). [31] S. Chakabai and N. Banejee, Eu. Phys. J. Plus. 131 : 144 (16).

21 1 [3] M. M. Som and N. O. Sanos, Phys. Le. A. 87, 89 (1981); S. R. Maii, Phys. Rev. D. 5, 518 (198); B. Modak, J. Asophys. Ason. 5, 317 (1984); A. Banejee, S. Choudhuy and B. Bhui, Phy. Rev. D. 4, 67 (1989); L. K. Pael and R. Tikeka, Mahemaics Today IX, 19 (1991); D. Schafe and H. F. Goenne, Gen. Relaiv. Gav. 4, 119 (); B. V. Ivanov, Gen. Relaiv. Gav. 44, 1835 (1); L. Heea, G. Le Denma, N. O. Sanos and A. Wang, In. J. Mod. Phys. D. 13, 583 (4); S. Chakabai and N. Banejee, Eu. Phys. J. C. 77 : 166 (17); N. Banejee and S. Chakabai, Phys. Rev. D. 95, 415 (17). [33] L. Heea and N. O. Sanos, Phys. Rep. 86, 53 (1997). [34] L. Heea, A. Di Pisco, J. Main, J. Ospino, N. O. Sanos, and O. Toconis, Phys. Rev. D. 69, 846 (4). [35] L. Heea and J. Ponce de Leon, J. Mah. Phys. (N.Y.) 6, 18 (1985). [36] N. O. Sanos, Mon. No. R. Ason. Soc. 16, 43 (1985). [37] R. Chan, Mon. No. R. Ason. Soc. 316, 588 (). [38] C. W. Misne and D. Shap, Phys. Rev. 136, B571 (1964). [39] C. A. Kolassis, N. O. Sanos and D. Tsoubelis, Class. Quanum Gav (1988). [4] O. M. Pimenel, F. D. Loa-Clavijo, and G. A. Gonzalez, Gen. Relaiv. Gavi. 48, 14 (16). [41] T. Clifon, P. G. Feeia, A. Padilla, and C. Skodis, Phys. Rep. 513, 1 (1). [4] B. C. Paul, P. S. Debnah, and S. Ghose, Phys. Rev. D. 79 : 83534, (9). [43] S. Nojii and S. D. Odinsov, Phys. Rev. D. 68, 1351 (3). [44] S. Nojii and S. D. Odinsov, Gen. Rel. Gav. 36, 1765 (4). [45] T. Clifon and J. D. Baow, Phys. Rev. D. 7, 135 (5). [46] E. Pechlane and R. Sexl, Commun. Mah. Phys., 165 (1966). [47] M. Visse and C. Bacelo, COSMO-99: pp. 98 (). [48] C. Bacelo and M. Visse, Phys. Le. B. 466, 17 (1999). [49] E. Flanagan and R. Wald, Phys. Rev. D. 54 : 633, (1996). [5] N. Banejee and T. Paul, Eu. Phys. J. C. 78 : 13 (18). [51] S. Chakabai, Eu. Phys. J. C. 78 : 96 (18). [5] A. I. M. Hamid, R. Goswami, and S. D. Mahaaj, Class. Quan. Gav. 31, 1351 (14). [53] S. Chakabai and N. Banejee, Eu. Phys. J. C. 77 : 166 (17). [54] S. Chakabai and N. Banejee, Phys. Rev. D. 95, 415 (17). [55] A.K. Raychaudhui, Phys. Rev., 98, 113 (1955).