arXiv:1905.08097v3 [gr-qc] 13 Aug 2019 yaptho h eSte erc hc smatched ul- stiff is a which with EOS filled metric, shell with Sitter a fluid to de trarelativistic layer the boundary a of through patch a by iigalgtmt lentv obakholes. black to alternative legitimate pro- the therefore a horizon, to viding event an shell non- without a the solution is singular gravastar matches The that geometry. Schwarzschild layer exterior boundary second a inwt qaino tt (EOS) re- a state interior of is an equation gravastar and with of a gion [15–24] consists form which instance, original configuration for its four-layer In (see, decade therein). last references the in- considerable in attracted has terest which 14], [13, Mottola and in holes black [12]. classical theory of Einstein paradoxes to the several a solutions the for as the up proposed [11] or clear been (see have 7–10] which (ECOs) [2, review), holes objects recent black compact to exotic alternatives holes so-called even black or non-Einsteinian to 6] due [5, as could interpreted GW. observations popu- be of the also most that detections found the future authors and are some However current stars of for astronomy neutron candidates the and assumed lar in commonly holes era is black new It that a objects. to spec- compact birth recent astrophysical gave electromagnetic as 4] the well [3, in as trum holes 2] black [1, of Collaborations observations Virgo and LIGO ∗ ‡ † [email protected] [email protected] [email protected] h su fsaiiyo rvsashsbe thoroughly been has gravastars of stability of issue The n fteeEO stegaatrpooe yMazur by proposed gravastar the is ECOs these of One the by (GW) waves gravitational of detection The 3 etod aeáia optçoeCgiã CC) Univer (CMCC), Cognição e Computação Matemática, de Centro ASnmes 04.50.Kd,04.70.Bw,04.30.-w,04.80.Cc numbers: PACS si algebraically the of rin in the class dominate of the tails end power-law to Schwarzschildian the belong times, at might echoes which of non mode absence against lating the stable h is black is Schwa feature a star the Schwarzschild peculiar of the mimicker to that good close prove very as files a star be Schwarzschild can the non Schwarzscmakes modes axial the quasinormal arbitrary that observable study shown the we been Here has perturbations. it symmetric Recently surface. pr positive the and at surface its outside geometry Schwarzschild h cwrshl tri nutaopc betbyn the beyond object ultracompact an is star Schwarzschild The .INTRODUCTION I. 1 nttt fPyisadRsac eteo hoeia Phy Theoretical of Centre Research and Physics of Institute .A Konoplya, A. R. tbeShazcidsasa lc-oemimickers black-hole as stars Schwarzschild Stable aut fPiooh n cec,Slsa nvriyi O in University Silesian Science, and Philosophy of Faculty ilkoMkaaSre,Mso 118 usa Federat Russian 117198, Moscow Street, Miklukho-Maklaya 6 2 epe redhpUiest fRsi RD University) (RUDN Russia of University Friendship Peoples p u blço E:02010 at nr,S,Brazil SP, André, Santo 09210-180, CEP: Abolição, Rua erčv á.1,C-4 1Oaa zc Republic Czech Opava, 01 CZ-746 13, nám. Bezručovo = p ǫ = ieie hr is there Likewise, . ,2, 1, − < ǫ ∗ .Posada, C. 0 described , 1, † .Stuchlík, Z. faShazcidbakhl.Smlrcnlsoswere conclusions Similar Pani those by hole. from drawn black different Schwarzschild is spec- gravastars a the of of that QNMs quasi- axial concluded of They their trum computed (QNMs). modes and Chirenti gravastars normal perturbations of paper, stability axial seminal the against investigated a [27] In Rezzolla and gravastar. pressure uous tr.Bcdh 3]soe hti n sue pos- a with assumes monotonically decreases one density which if constant density that of energy showed case itive [33] the Buchdahl config- to limiting restricted stars. this not of is existence uration The diverges. pressure nw htwe h aiso h cwrshl star Schwarzschild the bound of well Buchdahl radius is the the It reaches when density. spherically that energy a ex- uniform known for an with equations is mass Einstein’s which symmetric the inte- star, of Schwarzschild Schwarzschild solution or act the [32] reconsidered solution rior [31] Mottola and [30]. in consideration into taken and momentum not angular were (of resultant mass uncertainties the existing in large percents) ex- of a way, tens of This kind neu- [30]. rotating a star slowly adding the tron that, via of modes done After quasinormal was from case terms hole. non-rotating consider- the black of under Kerr trapolation object the central is from the follows ation which that value supposition that by discarded the of fixed are mass was gravastars the gravastar there the the though that observations, GW150914 argued the was by spacetime. it gravastar [30], the the In in of perturbations context Cardoso linear the recently, studied in More 19] instability. [18, ergoregion studied were perturbati gravastars Scalar of gravastars. thin-shell of oscillations ra- and the model, gravastar studied three-layer Horvat simple [25] a of Wiltshire stability and dial Visser investigated. oercnl,i oncinwt rvsas Mazur gravastars, with connection in recently, More suei h xenllyrwihvanishes which layer external the in essure gnal. peia etrain swl.Another well. as perturbations spherical tal. et peia etrain,adso that show and perturbations, spherical dw.Isedw bev nonoscil- a observe we Instead gdown. pca oe.A smttclylate asymptotically At modes. special idsa ssal gis spherically against stable is star hild 1 l.Tedcyn iedmi pro- time-domain decaying The ole. zcidlmta n ihs what wishes, one as limit rzschild n .Zhidenko A. and iaeFdrld B (UFABC), ABC do Federal sidade 2]cniee h tblt facontin- a of stability the considered [26] tal et isadAstrophysics, and sics uhallmt hc has which limit, Buchdahl 2]wocmue xa n polar and axial computed who [28] pava, ion , ,3, 1, R (9 = ‡ / 4) M h central the , tal et [29] ons 2 r, then there is a general upper mass bound M ≤ 4R/9 nearly Schwarzschild quasinormal modes end up not with which is independent of the relation between pressure and echoes, but with a new nonoscillating decaying mode, density. Configurations more compact than those within which, in its turn, goes over into the asymptotic power- the Buchdahl limit were believed to be unstable, as one law tails. The latter are identical to the Schwarzschild would have either assume the existence of the singularity ones. Thus, the l> 1 perturbations of the Schwarzschild inside the star or the pressure which is increasing out- star are indistinguishable from those of the black hole ward. The uniform density configurations were studied neither at the time of quasinormal ringing, nor at later in [34–37]. A model of black-hole mimicker in the form times, including the asymptotic tails. of a nonlocal star in the context of ghost-free, infinite This paper is organized as follows. In section II we derivative gravity was proposed in [38]. review the Schwarzschild star and the emergence of a Therefore, the existence of the Buchdahl limit together negative pressure interior when we consider ultracompact with the assumption of constant density, considered ‘un- configurations beyond the Buchdahl limit. In section III physical’ (see however [39]), has left the Schwarzschild we find quasinormal of axial gravitational perturbations star unstudied thoroughly (a notable exception can be via time-domain integration. Finally, in the conclusion found in [15]). Nevertheless, when one does consider we summarize the obtained results and mention some the constant density star below the Buchdahl limit, important open questions. 2M proper matching of the Schwarzschild star with the exte- rior Schwarzschild metric. 2 RSr RS 1 − − 1 − Note that in contrast to the model proposed originally R3 R in [13, 14], here the phase transition occurs first at the p(r)= ǫ · r r , (6) RS RSr2 center, once the Buchdahl bound is reached. The surface 3 1 − − 1 − R = RS is a null surface; however, there is no interior R R3 r r trapped surface and there is no formation of an event where ǫ denotes the mass-energy density, which is a con- horizon. The divergence in the pressure at R0 can be in- stant, tegrated through the Komar formula [44], if one assumes the presence of a transverse pressure at the hypersurface 3M 3RS R , such that = , r < R, 0 ǫ = 4πR3 8πR3 (7) ( 0, r ≥ R. f 8πǫ 8π r2(p − p)= R3δ(r − R ). (11) Notice that the pressure vanishes p(R)=0 at the sur- h ⊥ 3 0 0 face r = R. The Schwarzschild star is matched at the r boundary r = R, with the asymptotically flat exterior This δ function, associated to a surface layer of trans- Schwarzschild solution. Note that the central pressure verse stresses, provides a strict matching at the surface becomes infinite when the radius of the star reaches the as required by the Israel junction conditions [44, 45] on Buchdahl bound R = (9/8)RS [33]. However, a further null surfaces [46] (see [31] for details). analysis shows that (6) is regular except at some radius Note that in contrast with the model proposed in [13, R0 where the denominator in (6), 14] where the authors introduced a thin shell filled with a stiff fluid p = ǫ, to join the interior de Sitter with RS RS r2 the Schwarzschild exterior, in this new description there D ≡ 3 1 − − 1 − , (8) R R3 is no thin shell of matter. Instead, the strict matching r r between the ultracompact Schwarzschild star and the ex- vanishes in the range 0 III. GRAVITATIONAL PERTURBATIONS OF R/RS ωiRS ωlRS ULTRACOMPACT SCHWARZSCHILD STARS 1.00 0.7474 − 0.1779i 1.01 0.7486 − 0.1818i The axial gravitational perturbations for a relativistic 1.02 0.7430 − 0.1831i mass, with given mass-energy density ǫ and pressure p, 1.03 0.7398 − 0.1811i are determined by the master equation 1.04 0.7379 − 0.1777i 1.05 0.737 − 0.175i ∂2 ∂2 1.06 0.737 − 0.173i −0.11i − + Vl(r) Ψ(r, t)=0 (14) ∂t2 ∂r2 1.07 0.737 − 0.170i −0.10i ∗ 1.08 0.737 − 0.168i −0.08i where the effective potential Vl is given by [49, 50] 1.09 0.737 − 0.167i −0.070i 1.10 0.737 − 0.166i −0.057i l(l + 1) m(r) Vl(r)= f(r) − 6 − 4πp(r)+4πǫ , (15) 1.11 0.737 − 0.164i −0.042i r2 r3 1.12 0.737 − 0.163i −0.023i −9 where m(r) is the mass enclosed by r 1.13 0.218055 − 2.4 × 10 i 1.14 0.371288 − 1.2 × 10−6i r 1 15 0 470264 − 2 6 × 10−5 m(r)=4π ǫ · r2dr, (16) . . . i 0 Z TABLE I. Dominant quasinormal frequencies of axial gravita- and r∗ denotes the ‘tortoise’ coordinate defined in the tional perturbations (l = 2) obtained by fitting time-domain standard form by profiles. Initially the ringdown is governed by the oscillating mode ωi, at late times purely imaginary mode ωl dominates. dr The last three lines correspond to dominant frequencies of the dr∗ = . (17) constant density star above the Buchdahl limit, taken from f(r)h(r) [50]. Therefore, r∗ → ∞ correspondsp to spacial infinity and r∗ → −∞ when approaching the singular surface r → + where we have used the following definitions for the R0 . points: N = (u+∆, v+∆), W = (u+∆, v), E = (u, v+∆) Note that in the vacuum exterior (r ≥ R) ǫ = 0 and and S = (u, v). The initial data are specified on the two p(r)=0, hence the perturbation potential (15) reduces null surfaces u = u0 and v = v0. to the Regge-Wheeler [51] form From figure 3 we see that the ringdown of Schwarzschild stars is similar to the quasinormal ring- out RS l(l + 1) 3RS Vl (r)= 1 − − , r ≥ R. ing of the Schwarzschild black hole and becomes in- r r2 r3 distinguishable as the star’s radius approaches the (18) Schwarzschild radius. For larger Schwarzschild stars Here RS = 2M = 2m(R), where M corresponds to the (R . 9RS/8), the purely imaginary mode appears at total mass of the star. late time. A similar phenomenon was observed in the In figure 1 we show the behavior of the effective po- Einstein-Gauss-Bonnet theory, where a purely imaginary tential as a function of the tortoise coordinate: at the nonperturbative mode appears in the spectrum [55], lead- star’s surface it has a discontinuity due to discontinuity ing to instability for sufficiently large coupling [56, 57]. of ǫ. Since p(r) grows unboundedly as r → R0 the po- Although the effective potential for the gravitational per- tential has a negative gap from the righthand side of the turbations of Schwarzschild stars also has a negative gap, − peak, approaching 0 as r∗ → −∞. As the radius of the it does not become deeper as l grows, so that we do not star approaches the Schwarzschild radius, the disconti- have any indications of the eikonal instability [58]. In nuity and the negative gap become smaller and can be addition, we have observed that for l > 2 the imaginary neglected, if R is sufficiently close to RS (see Fig. 2). mode does not dominate at late time. In order to produce the time-domain profiles, we in- Note, that the quasinormal spectrum of the tegrate the wavelike equation (14) rewritten in terms of Schwarzschild black hole has also the purely imaginary the light-cone variables u = t − r∗ and v = t + r∗. The mode called algebraically special. However, the alge- discretization scheme was suggested in [52] and used in a braically special mode of the Schwarzschild black hole number of papers (see for instance [53, 54] and references corresponds to high imaginary part, so that it does not therein) showing very good accuracy. It has the following dominate in the signal and, consequently, we cannot see form: it in the time-domain profile. The algebraically special mode for the Schwarzschild solution has the form [59] Ψ(N)=Ψ(W )+Ψ(E) − Ψ(S) (19) V (W )Ψ(W )+ V (E)Ψ(E) (l − 1)l(l + 1)(l + 2) −∆2 + O(∆4) , ωaRS = −i , (20) 8 6RS 5 2 VHr*LRS 2 VHr*LRS 0.6 0.006 0.5 0.4 0.004 0.3 0.002 0.2 0.1 r*RS -35 -30 -25 -20 -15 -10 r*RS -5 5 FIG. 1. Effective potential for the axial gravitational perturbations (l = 2) of Schwarzschild star R = 1.1RS . Left panel: negative gap when approaching the singular point r → R0. Right panel: discontinuity at r = R. 2 VHrLRS so that the algebraically special mode for l = 2 is ω = −4i/RS what does not exclude possibility that the 0.6 purely imaginary mode may go over into the known al- 0.5 gebraically special mode in the Schwarzschild limit. If so, we should assume that the purely imaginary mode, 0.4 which we observe here for l = 2, exists also for higher l, but it does not dominate at late times for l > 2 even 0.3 when R → 9RS/8, because at higher l the algebraically special mode (20) is strongly damped. 0.2 After the ringdown phase and the transition consisting of exponential decay we see the asymptotic power-law 0.1 tails, which are identical to those for the Schwarzschild black hole [60, 61], rRS 1.0 1.5 2.0 2.5 |Ψ|∼ t−(2l+3). (21) FIG. 2. Effective potentials for the axial gravitational pertur- bations (l = 2) of the Schwarzschild black hole (black) and for We also notice that the ringdown of the Schwarzschild the interior (R0 10-6 while the asymptotic tails are identical to those for the Schwarzschild ones. 10-9 IV. CONCLUSIONS 10-12 We have shown that the Schwarzschild stars are stable tRS 50 100 150 200 250 300 350 not only against spherical perturbations, but also against axial nonspherical ones, which is confirmed by the de- FIG. 3. Ringdown of the axial gravitational perturbations cayed time-domain profiles for the corresponding pertur- (l = 2) for the Schwarzschild black hole (black) and for the bation equations. This is one more step to establishing Schwarzschild stars: R = 1.04RS (blue, lower), R = 1.08RS = 1 1 the viability of the Schwarzschild star model. At the (red), R . RS (magenta, upper). same time, the Schwarzschild star can be a very good 6 ‘mimicker’ of the gravitational wave response of a black is whether the purely imaginary modes which we found hole at l = 2 and higher multipoles, because it can ap- belong to the class of algebraically special modes. The proach the Schwarzschild spectrum as closely as neces- latter can be checked by numerical calculations of quasi- sary. The latter is possible because the internal surface normal modes in the frequency domain, for example, with with an infinite redshift provides the same boundary con- the help of the Frobenius method and accurate extrapo- ditions for quasinormal mode problem as a black hole lation to the Schwarzschild limit. Finally, an important does. Moreover, at later times, when gravitational echoes question is whether the lower multipoles (l =0, 1) could should show whether we have a gravastar or a classical be seen in the gravitational wave signal we observe nowa- black hole, the purely imaginary damping mode is domi- days. If so, then this could be the way to distinguish the nating, instead of echoes, which goes over into the usual Schwarzschild star from a black hole. Schwarzschild asymptotic power-law tails. The present study raises a number of interesting ques- tions. First of all, whether the polar gravitational pertur- ACKNOWLEDGMENTS bations are stable as well. In case of complete stability, the rotating model must be constructed beyond the linear This paper was supported by the grant 19-03950S of order in rotation parameter and its stability, quasinor- Czech Science Foundation (GAČR) and the “RUDN Uni- mal modes and shadows investigated. Another question versity Program 5-100”. [1] Abbott B P et al. (LIGO Scientific, Virgo) 2016 Phys. 104013 (Preprint 1408.6929) Rev. Lett. 116 241103 (Preprint 1606.04855) [21] Pani P 2015 Phys. Rev. D92 124030 [Erratum: Phys. [2] Abbott B P et al. (LIGO Scientific, Virgo) 2016 Phys. Rev.D95,no.4,049902(2017)] (Preprint 1506.06050) Rev. Lett. 116 061102 (Preprint 1602.03837) [22] Völkel S H and Kokkotas K D 2017 Class. Quant. Grav. [3] Akiyama K et al. (Event Horizon Telescope) 2019 Astro- 34 125006 (Preprint 1703.08156) phys. J. 875 L1 [23] Nakao K i, Yoo C M and Harada T 2019 Phys. Rev. D99 [4] Goddi C et al. 2016 Int. J. Mod. Phys. D26 1730001 044027 (Preprint 1809.00124) [1,863(2017)] (Preprint 1606.08879) [24] Wang Y T, Zhang J and Piao Y S 2019 Phys. Lett. B795 [5] Yunes N, Yagi K and Pretorius F 2016 Phys. Rev. D94 314–318 (Preprint 1810.04885) 084002 (Preprint 1603.08955) [25] Visser M and Wiltshire D L 2004 Class. Quant. Grav. 21 [6] Konoplya R and Zhidenko A 2016 Phys. Lett. B756 350– 1135–1152 (Preprint gr-qc/0310107) 353 (Preprint 1602.04738) [26] Horvat D, Ilijić S and Marunović A 2011 Class. Quant. [7] Cardoso V, Franzin E and Pani P 2016 Phys. Grav. 28 195008 (Preprint 1104.3537) Rev. Lett. 116 171101 [Erratum: Phys. Rev. [27] Chirenti C B M H and Rezzolla L 2007 Class. Quant. Lett.117,no.8,089902(2016)] (Preprint 1602.07309) Grav. 24 4191–4206 (Preprint 0706.1513) [8] Cardoso V, Hopper S, Macedo C F B, Palenzuela C [28] Pani P, Berti E, Cardoso V, Chen Y and Norte R 2009 and Pani P 2016 Phys. Rev. D94 084031 (Preprint Phys. Rev. D80 124047 (Preprint 0909.0287) 1608.08637) [29] Cardoso V, Crispino L C B, Macedo C F B, Okawa H and [9] Abedi J, Dykaar H and Afshordi N 2017 Phys. Rev. D96 Pani P 2014 Phys. Rev. D90 044069 (Preprint 1406.5510) 082004 (Preprint 1612.00266) [30] Chirenti C and Rezzolla L 2016 Phys. Rev. D94 084016 [10] Konoplya R A and Zhidenko A 2016 JCAP 1612 043 (Preprint 1602.08759) (Preprint 1606.00517) [31] Mazur P O and Mottola E 2015 Class. Quant. Grav. 32 [11] Cardoso V and Pani P 2019 Living Rev. Rel. 22 4 215024 (Preprint 1501.03806) (Preprint 1904.05363) [32] Schwarzschild K 1916 Sitzungsber. Preuss. Akad. Wiss. [12] Wald R M 2001 Living Rev. Rel. 4 6 (Preprint gr- Berlin (Math. Phys.) 1916 424 qc/9912119) [33] Buchdahl H A 1959 Phys. Rev. 116 1027 [13] Mazur P O and Mottola E 2001 (Preprint gr-qc/0109035) [34] Boehmer C G 2004 Gen. Rel. Grav. 36 1039–1054 [14] Mazur P O and Mottola E 2004 Proc. Nat. Acad. Sci. (Preprint gr-qc/0312027) 101 9545–9550 (Preprint gr-qc/0407075) [35] Stuchlik Z 2000 Acta Phys. Slov. 50 219–228 (Preprint [15] Cattoen C, Faber T and Visser M 2005 Class. Quant. 0803.2530) Grav. 22 4189–4202 (Preprint gr-qc/0505137) [36] Stuchlík Z, Hledík S and Novotný J 2016 Phys. Rev. D94 [16] DeBenedictis A, Horvat D, Ilijic S, Kloster S and 103513 (Preprint 1611.05327) Viswanathan K S 2006 Class. Quant. Grav. 23 2303– [37] Stuchlík Z, Schee J, Toshmatov B, Hladík J and Novotný 2316 (Preprint gr-qc/0511097) J 2017 JCAP 1706 056 (Preprint 1704.07713) [17] Lobo F S N and Arellano A V B 2007 Class. Quant. Grav. [38] Buoninfante L and Mazumdar A 2019 Phys. Rev. D100 24 1069–1088 (Preprint gr-qc/0611083) 024031 (Preprint 1903.01542) [18] Cardoso V, Pani P, Cadoni M and Cavaglia M 2008 Phys. [39] Misner C W, Thorne K S and Wheeler J A 1973 Gravi- Rev. D77 124044 (Preprint 0709.0532) tation (San Francisco: Freeman and Co.) [19] Chirenti C B M H and Rezzolla L 2008 Phys. Rev. D78 [40] Posada C 2017 Mon. Not. Roy. Astron. Soc. 468 2128 084011 (Preprint 0808.4080) [41] Posada C and Chirenti C 2019 Class. Quant. Grav. 36 [20] Sakai N, Saida H and Tamaki T 2014 Phys. Rev. D90 065004 (Preprint 1811.09589) 7 [42] Chandrasekhar S 1964 Astrophys. J. 140 417–433 [Erra- [53] Konoplya R A, Zinhailo A F and Stuchlík Z 2019 Phys. tum: Astrophys. J.140,1342(1964)] Rev. D99 124042 (Preprint 1903.03483) [43] Konoplya R A, Stuchlík Z and Zhidenko A 2019 Phys. [54] Konoplya R A and Zhidenko A 2011 Rev. Mod. Phys. 83 Rev. D99 024007 (Preprint 1810.01295) 793–836 (Preprint 1102.4014) [44] Poisson E 2004 A Relativist’s Toolkit (Cambridge, UK: [55] Konoplya R A and Zhidenko A 2008 Phys. Rev. D77 Cambridge University Press) 104004 (Preprint 0802.0267) [45] Israel W 1966 Nuovo Cim. B44S10 1 [56] González P A, Konoplya R A and Vásquez Y 2017 Phys. [46] Barrabes C and Israel W 1991 Phys. Rev. D43 1129–1142 Rev. D95 124012 (Preprint 1703.06215) [47] Bekenstein J D 1973 Phys. Rev. D7 2333 [57] Konoplya R A and Zhidenko A 2017 JHEP 09 139 [48] Bardeen J M, Carter B and Hawking S W 1973 Commun. (Preprint 1705.07732) Math. Phys. 31 161 [58] Cuyubamba M A, Konoplya R A and Zhidenko A 2016 [49] Chandrasekhar S and Ferrari V 1991 Proc. Roy. Soc. Phys. Rev. D93 104053 (Preprint 1604.03604) Lond. A432 247–279 [,1025(1991)] [59] Chandrasekhar S and Detweiler S L 1976 Proc. Roy. Soc. [50] Chandrasekhar S and Ferrari V 1991 Proc. Roy. Soc. Lond. A350 165–174 Lond. A434 449 [60] Price R H 1972 Phys. Rev. D5 2419–2438 [51] Regge T and Wheeler J A 1957 Phys. Rev. 108 1063 [61] Price R H 1972 Phys. Rev. D5 2439–2454 [52] Gundlach C, Price R H and Pullin J 1994 Phys. Rev. [62] Kokkotas K D 1994 Mon. Not. Roy. Astron. Soc. 268 D49 883–889 (Preprint gr-qc/9307009) 1015–1018 ISSN 0035-8711