arXiv:1607.03593v1 [gr-qc] 13 Jul 2016 hc ssilukona ula n h super-nuclear the to and core, nuclear used densities. neutron-star at be unknown the still of can is state which moments of equation multipole measure- the first Thus, constrain the structure. of internal ments its on ways similar e.g. (see, distribution 2]). matter re- [1, the Refs. is of which stiffness body, the rotating to a the lated of instance, deformability the For about information encodes moments. moment quadrupole multipole spin-induced its non- in a th encoded of are matter properties of The distribution self-gravitating body. spherical, the gravita- outside the fields of tional moments multipole be the may by symmetry characterized spherical from with Deviations us inside object. distribution the provides matter forces. the deformations on tidal information such important and/or to centrifugal response to The due bi- symmetry in spher- from ical are deviates and/or shape their rotate however, systems, they nary If symmetric. spherically ∗ ‡ † [email protected] [email protected] [email protected] h utpl oet fanurnsa eedin depend star neutron a of moments multipole The roughly are objects astrophysical Self-gravitating ia eomblt n -oeQrltosfrgravastars for relations I-Love-Q and deformability Tidal 4 iatmnod iia Spez”Uiestad oa&Se & Roma Universit`a di ”Sapienza” Fisica, di Dipartimento nttt ueirTenc S,Uiesdd eLso - Lisboa de Universidade IST, T´ecnico — Superior Instituto ASnmes 42.q 43.w 47.w 04.70.-s. 04.70.Bw, 04.30.-w, 04.20.-q, numbers: PACS gravastars. like objects compact Final exotic state. detectio constrain of gravitational-wave the equations neutro from of the deformability class tidal to limited a similarly universal for that, no emerge find possibility might and the shell, deny matter cannot the we for l state negativ black-hole of are the equations moment approach quadrupole continuously and relations number I-Love-Q Love The drasti are star. gravastar tron thin-shell of a number of relations Love I-Love-Q calculat tidal n The we electric) horizon (quadrupolar expansion, the event field and an small-tidal moment, cl possess and large small-spin not a a do in of that properties black-ho objects the compact deformation approach exotic the they investigate how we and here objects relation compact “I-Love-Q” other interesting These for An astronomy. state. gravitational-wave and of som physics equation through stellar related the are on models quark-star and neutron-star h oeto nri,tesi-nue uduoemoment, quadrupole spin-induced the inertia, of moment The 1 .INTRODUCTION I. eateto hsc,Rky nvriy Nishi-ikebukur University, Rikkyo Physics, of Department 3 srnmclIsiue oouUiest,Aramaki-Aoba University, Tohoku Institute, Astronomical 5 etoMliicpia eAto´sc ETA Departa CENTRA, — Astrof´ısica de Multidisciplinar Centro 2 eateto ahmtc n hsc,Gaut colo S of School Graduate Physics, and Mathematics of Department aiUchikata Nami sk iyUiest,Smysik,Oaa5888,Japa 558-8585, Osaka Sumiyoshi-ku, University, City Osaka 1,2 , ∗ Dtd uy1,2016) 14, July (Dated: hjnYoshida Shijun en , in NNRm1 izaeAd oo5 08 oa Italy Roma, 00185 5, Moro Aldo Piazzale Roma1, INFN zione tvt se .. e.[]frarve) nti context, this mass In the namely review). rel- moments, a general multipole for of lowest-order [7] the and null- Ref. metric e.g., provide Kerr (see, to ativity the black or of objects distinguish tests compact to hypothesis exotic used three from least be holes at can of moments [3– measurements quantities multipole independent two Thus, these by 6]. only determined uniquely ul hrceie yol hi mass their are only holes by black Astrophysical characterized black- the fully of theorems. tests to no-hair related hole is objects compact of ments a momentum lar needn fteeuto fsaet ihnafwper- few a within to state nearly of are equation that the of “I-Love-Q”) independent as (dubbed relations versal nue uduoemoment quadrupole induced deformability rnsas h oeto inertia of re- moment More neu- the deformed slightly stars, 9]). for tron [8, that, found Refs. been see has it (e.g., mo- cently, found multipole been lowest-order have the ments and inertia of the neutron-star between moment a the relations considering for universal By state approximately of interior, equations state. of of set reasonable equation mildly only the are on which dependent stars compact of moments mul- tipole lowest-order the among relations approximate ous role. prime a play would M L v oic as1 0901Lso,Portugal Lisboa, 1049-001 1, Pais Rovisco Av. UL, nte neetn plcto ftemlioemo- multipole the of application interesting Another thsbe eetypitdotta hr xs vari- exist there that out pointed recently been has It nua momentum angular , facmatbnr nprlcnb sdto used be can inspiral compact-binary a of n al ifrn rmtoeo nodnr neu- ordinary an of those from different cally 3 , mt ecnie ait fpolytropic of variety a consider We imit. rbe swehrsmlrrltoshold relations similar whether is problem † y edsushwamaueeto the of measurement a how discuss we ly, h oeto nri,tequadrupole the inertia, of moment the e astro- for implications important have s -trcs,a prxmt universality approximate an case, n-star ,Tsiak,Tko1180,Japan 171-8501, Tokyo Toshima-ku, o, t nteILv- eain.However, relations. I-Love-Q the in ity n al Pani Paolo and s fti-hl rvsas hc are which gravastars, thin-shell of ass oak,Sna 8-58 Japan 980-8578, Sendai Aoba-ku, , rvsaswt oyrpcti shell. thin polytropic a with gravastars mildly only depend which relations e rasaeiesnuaiy Working singularity. spacetime a or elmt oase hs questions, these answer To limit. le o escmatmdl n the and models compact less for e λ n h ia oenme of number Love tidal the and J tetdlLv ubr,adtespin- the and number), Love tidal (the ihpltoi hnshells thin polytropic with l ihrmlioemmnsbeing moments multipole higher all , et eF´ısica, de mento cience, n 4,5‡ J n uduoemoment quadrupole and , Q lostsysm uni- some satisfy also , I h ia quadrupole tidal the , M n angu- and Q , 2 cent level [10, 11]. Here, the tidal Love number λ is stars. These peculiar models with strongly anisotropic defined as the ratio of the tidally-induced quadrupole pressure have negative Love number and quadrupole moment of a nonspinning star to the strength of the per- moment, in contrast with the neutron-star case. A nega- turbing quadrupolar tidal field. tive value of the spin-induced quadrupole moment ba- These nearly universal relations will be helpful to sically indicates that the matter distribution of the object extract physically important quantities from observa- is prolate. Indeed, Yagi and Yunes observed that the ro- tional data if they exist within the required accuracy. tating stars with strongly anisotropic pressure become For example, the I-Love-Q relations can break degen- prolate despite of the centrifugal force, due to a strongly eracy in the models used for X-ray and gravitational- anisotropic pressure [20]. wave observations of neutron stars. For these reasons, Similar results have been also obtained for another pe- nearly-universal relations of compact objects have at- culiar compact star model, namely for thin-shell gravas- tracted much attention in the last few years. Although tars. The original gravastar model has been proposed the I-Love-Q relations were originally found [10, 11] by Mazur and Mottola [32] as an alternative to the fi- for slightly deformed, isolated, nonmagnetized com- nal state of the stellar evolution of very massive stars, pact stars, they have been extended to more general whose collapse would form a black hole in the standard cases, namely dynamical configurations [12], rapid rota- scenario. The thin-shell gravastar model is a simpli- tion [13–16], nonbarotropic [17] and anisotropic [18–20] fied version of the original gravastar and is composed fluids, strong magnetic fields [21], and also deviations of a vacuum core with a positive cosmological constant from general relativity [10, 11, 22–24]. Possible expla- (de Sitter core) surrounded by an infinitesimally thin nations for the emergence of this approximate univer- shell, which is required to match the interior core to the sality are given in Refs. [25, 26]. Similar universal rela- Schwarzschild exterior metric [33]. Although this sce- tions have been also found among higher-order multi- nario lacks a precise formation mechanism, the de Sit- pole moments induced by tidal effects for nonrotating ter core is assumed to appear through a quantum phase neutron stars [27]. The inclusion of the rotational Love transition in the vicinity of the would-be event horizon numbers in the I-Love-Q relations has been also argued during the gravitational collapse of very massive ob- in Ref. [28]. jects. Gravastars are a hypothetical exotic compact ob- These universal relations among higher-order multi- jects which can be as compact as black holes but are free pole moments of a neutron star are reminiscent of the from the theoretical problems associated with an event black-hole no-hair theorems, namely they show that horizon and a spacetime singularity [32]. only a handful of low-order multipole moments can The possibility of testing the gravastar scenario with (approximately) characterize the gravitational field of a electromagnetic and gravitational-wave observations compact object. has been recently argued in Refs. [34–37]. Some mod- To investigate this problem and to elucidate an ori- els of thin-shell gravastars are stable against small radial gin of the I-Love-Q relation, it is interesting to under- disturbances [33] although nonlinear instability might stand how the I-Love-Q relations behave in the black- occur for those models whose compactness is larger hole limit. Since the ratio of the mass M to the ra- than 1/3 (and therefore possess a light ring) [38, 39] and dius R (i.e., the compactness) of any standard spheri- for highly-spinning compact models due to the ergore- cal perfect-fluid star has the upper limit M/R < 4/9 gion instability [40]. (in geometrized units, see, e.g., Refs. [29, 30]), the com- Recently, one of us investigated the I-Love-Q relations pactness of a perfect fluid star is disconnected from that for a particular model of thin-shell gravastar in which of a black hole, M/R = 1/2. To examine the black-hole the energy density of the thin shell vanishes [41]. Simi- limit of the I-Love-Q relations, therefore, we have to con- larly to the case of fluid stars with maximum anisotropic sider somewhat peculiar compact objects that can sus- pressure, in this model, I, λ and Q smoothly connect tain higher compactness. Yagi and Yunes [18–20] have to their black-hole values as the compactness increases. studied this problem using fluid stars with anisotropic Furthermore, low-compactness gravastars have nega- pressure. In the models they employed [31], stars with tive Love number and quadrupole moment and become the maximally anisotropic pressure can approach the prolate shaped when rotating. At the same time, by as- black-hole limit continuously, i.e. their compactness can suming that the equation of state for the thin shell is be as high as M/R 1/2 in the nonrotating case. They given by the same form as that given by equilibrium se- confirmed that the universality∼ of the I-Love-Q relations quences of spherical solutions with fixed values of the hold for stars with weakly anisotropic pressure and that gravitational mass, two of us [42] have studied rota- the tidal Love number and the quadrupole moment of tional effects on a thin-shell gravastar in the absence of anisotropic fluid stars continuously approach their cor- an external tidal field. Similarly to the case of gravas- responding black hole’s values as the compactness in- tars with zero thin-shell energy, we obtained prolate creases. On the other hand, they observed that, for shaped rotating gravastars and showed that some solu- strongly anisotropic pressure models, the behavior of tions may have the same quadrupole moment of a black the Love number and quadrupole moment are quite dif- hole with same mass and spin. This implies a confu- ferent from those of standard neutron stars and quark sion problem, namely these particular solutions cannot 3 be distinguished from a black hole through independent II. FORMULATION measurements of their mass, spin and quadrupole mo- ment. Thin-shell gravastars can be constructed by matching Given the plethora of applications related to the mul- two distinct spacetimes at the shell location. As pre- tipole moments and the tidal deformability of relativis- viously discussed, we focus on small deviations from tic compact objects, the motivation for the present study spherically symmetry, which can be investigated by us- is manifold. On the one hand, we extend our previ- ing Hartle and Thorne’s formalism [43, 44]. As unper- ous work [41, 42] to investigate the universality and turbed solutions, we consider a spherically-symmetric the black-hole limit of the I-Love-Q relations for vari- thin-shell gravastar [33] in which the de Sitter and ous thin-shell gravastars with a generic polytropic equa- Schwarzschild spacetimes are patched together. We are tion of state for the shell. We use the Hartle-Thorne for- interested in the dominant quadrupole deformations malism [43, 44], assuming the tidal and rotational ef- and, therefore, we do not consider in detail spherically- fects on the gravastar are sufficiently small, which is a symmetric perturbations which are also induced by ro- standard assumption to study the tidal Love numbers tation and by the tidal field. and the spin-induced multipole moment for compact As shown later, calculations of physical quantities re- objects [9, 45–48]. Since the thin-shell gravastar is com- lated to the rotational and tidal deformations of the posed of two distinct spacetimes, namely a de Sitter core gravastar are very similar to each other. The differences and a Schwarzschild exterior in the spherically symmet- appear only in the boundary conditions at infinity and ric case, we need to take account of the junction condi- in the time-reversal symmetry. Since the formulation tions of spacetime at the thin-shell position, where the for spin-induced deformations may be reduced to that two spacetimes are matched together [49, 50]. for the tidal deformations, we first give a concise sum- On the other hand, our novel results for the tidal de- mary of the master equations governing slowly rotating formability of thin-shell gravastars allow us to investi- thin-shell gravastars, as derived in Ref. [42]. We then gate the extend to which gravitational-wave observa- briefly discuss the prescription for obtaining tidally de- tions [51, 52] can constrain gravastar models. Tidal ef- formed nonrotating gravastars in the same framework. fects enter the two-body inspiral gravitational wave- In the following, we basically follow the treatment of forms at high post-Newtonian order [47, 48] (cf., e.g., Ref. [42] (see, also, Ref. [57], in which slowly rotating Ref. [53] for a review). Because the tidal Love num- regular black holes are studied with a similar perturba- bers of static [54, 55] (and, presumably, also rotating [28, tive approach). 56]) black holes are identically zero, any gravitational- wave measurement of a nonvanishing tidal deforma- bility would imply that one of the two objects is not a A. Interior and exterior spacetimes of a gravastar black hole. Conversely, gravitational-wave observations of compact-binary inspirals [51, 52] may be used to put Following Hartle and Thorne’s treatment [43, 44], we upper bounds on the tidal Love number of the two bod- assume the spacetimes outside and inside the slowly ro- ies, thus constraining exotic alternatives. tating thin-shell gravastar to be described by a common The plan of this paper is the following. In Sec. II, we form of the metric tensor, gαβ, namely briefly give the formulation for constructing models of ds2 = g dxαdxβ distorted thin-shell gravastars whose small deformation αβ is caused by the centrifugal and the tidal forces, respec- = f (r)(1 + 2ǫ2h(r, θ))dt2 tively. In Sec. III, we present numerical results for sta- − 1 2ǫ2m(r, θ) tionary and static axisymmetric models. We first con- + 1 + dr2 (2.1) f (r) r f (r) struct spherically symmetric, thin-shell gravastars and examine their radial stability, which depends on the + r2(1 + 2ǫ2k(r, θ)) dθ2 + sin2 θ dφ ǫω(r)dt 2 equation of state for the thin-shell matter. Next, we { − } present the results for the I-Love-Q relations of thin- + (ǫ3), h i shell gravastars with polytropic equation of state and O where ǫ stands for the smallness parameter for pertur- investigate their black-hole limit. Finally, in Sec. IV bations, and we have used the Greek letters (α, β, γ,...) we discuss how a measurement of the tidal deformabil- to denote spacetime indices. Here, the coordinate func- ity from the gravitational-wave detection of a compact- tions xα and the metric functions appearing in Eq. (2.1) binary inspiral can be used to constrain models of thin- are given by shell gravastars. We conclude in Sec. V. In this study, + µ + + + + we use geometrized units in which G = vc = 1, where (x ) =(t , r , θ , φ ), G and vc are the gravitational constant and the speed + + 2M of light, respectively. Note that, as mentioned later, we f (r) = f (r ) = 1 + , employ the symbol “c” to denote the compactness of the − r (2.2) object, i.e., c := M/R with M and R being the mass and (h(r, θ), k(r, θ), m(r, θ), ω(r)) = radius of the object with spherical symmetry. (h+(r+, θ+), k+(r+, θ+), m+(r+, θ+), ω+(r+)), 4 for the spacetime outside the thin shell, and the Legendre function of the second kind. The explicit 2 1 µ forms of Q2 and Q2 read (x−) =(t−, r−, θ−, φ−), 2 x(5 3x2) 3(x2 1) x + 1 (r−) Q 2(x) = + log , (2.12) f (r) = f −(r−) = 1 , 2 2− − − L2 (2.3) x 1 2 x 1 − 2 − (h(r, θ), k(r, θ), m(r, θ), ω(r)) = 1 2 2 3x 3x x + 1 Q2 (x) = x 1 −2 + log . (2.13) (h (r , θ ), k (r , θ ), m (r , θ ), ω (r )), − x 1 2 x 1 − − − − − − − − − − − p − − The interior solutions that are regular at the center of the for the spacetime inside the thin shell, where M and L star are given by (see, e.g., Refs. [41, 42, 57]) are, respectively, the mass of the unperturbed spheri- cal gravastar and the de Sitter horizon radius, defined ω− = C1, (2.14) by L = √Λ/3, with Λ being the positive cosmologi- = cal constant. Here and henceforth, the superscripts “+” m0− 0, (2.15) and “ ” indicate quantities defined outside and inside h− = C , (2.16) − 0 2 the thin-shell gravastar, respectively. When the distinc- 2 2 C3 3L + 5r 3L f −(r)Arctanh(r/L) tion is not necessary, however, the superscripts are fre- h− = − + , 2 8r2 L2 f (r) r quently omitted. To achieve separation of variables, − we expand the perturbation functions h(r), k(r), m(r) in (2.17) 2 2 2 2 Legendre polynomials Pl(cos θ) as follows [43, 44] C3 3L + 4r 3(L + r )Arctanh(r/L) k− = , 2 8r2L L − r h(r, θ) = h (r) + h (r)P (cos θ), 0 2 2 (2.18) m(r, θ) = m0(r) + m2(r)P2(cos θ), (2.4) m2− = r f −(r)h2−, (2.19) k(r, θ) = k2(r)P2(cos θ), − where C1, C2 and C3 are integration constants. The exte- where equatorial symmetry is assumed and the expan- rior and interior perturbed metric are matched so as to sion is truncated at the l = 2 order, since we are inter- fulfill the junction conditions, as we discuss in the next ested only in quadrupole perturbations. section. In the absence of a perturbing tidal field, we assume Now let us turn our attention to the case of the tidal that the spacetime is asymptotically flat, i.e., the met- deformations of a spherically symmetric gravastar. We ric perturbation must vanish in the limit of r+ ∞. → focus on gravastars whose deformation is induced by a The exterior solutions are therefore given by (see, e.g., static and axially symmetric 1 external tidal field. Thus, Refs. [43, 44]) we have to assume (i) ω = 0 and (ii) the spacetime is 2J not asymptotically flat because of the existence of a tidal ω+ = , (2.5) r3 source on the symmetry axis at a sufficiently large dis- tance. These requirements change the solutions given in J2 m+ = δM , (2.6) Eqs. (2.5)–(2.10) and (2.14)–(2.19) as follows: 0 − r3 δM J2 J = 0, (2.20) h+ = + , (2.7) 0 − r 2M r3(r 2M) r r2 2M − − h+ = D Q 2 1 3D 1 , (2.21) 1 1 r 2 − 1 2 M − − 2 M2 − r + 2 2 h2 = J 3 + 4 + BQ2 1 , (2.8) Mr r M − + 2M 1 r k2 = D1 Q2 1 J2 2M r r(r 2M) M − k+ = B Q 1 1 h+, (2.9) − 2 4 2 2 r − r − r(r 2M) M − − + 6pD 1 h+, (2.22) − 2 M − − 2 r4 dω+ 2 + p + m+ = (r 2M)h+, (2.23) m2 =(r 2M) h2 + , (2.10) 2 2 − − 6 dr ! − − C1 = 0, (2.24) where J and δM are the angular momentum and the spin-induced mass shift, respectively. Here, B is an in- tegration constant, through which the quadrupole mo- 1 In the case of a nonspinning, tidally deformed object, the spheri- ment of the rotating gravastar is defined as cal symmetry of the background configuration guarantees that per- turbations with different azimuthal number m decouple from each J2 8 Q = + BM3 . (2.11) other and are degenerate (i.e., the radial functions are independent M 5 of m). Therefore, we can consider axial symmetry and set m = 0 without loss of generality. After the radial functions are obtained, Note that, for slowly rotating Kerr black holes, regular- it is straightforward to derive the full, nonaxisymmetric, deformed m metric. ity at the horizon imposes B = 0. The function Ql is 5
µ where D1 and D2 are new integration constants. The The unit normal vector to the thin shell, n , is another other perturbation functions remain invariant. The primary quantity characterizing the thin shell. + asymptotic behavior of h (r) as r ∞ reads The first junction condition shows how the spacetimes 2 → induced on the thin shell are matched smoothly in an 8 M 3 M + intrinsic sense. It states that the induced metric hab := h 1 + O D1 µ 2 → − 5 r r g e eν is continuous through the thin shell, namely (2.25) µν a b r 2 M 3 1 + O D2 . − M r [[hab]] = 0, (2.31) Since (1 + gtt)/2 at large distance from an isolated star may− be regarded as the Newtonian gravitational where the double square brackets indicate potential [48], the first and second lines in the right- + hand side of Eq. (2.25) may be interpreted as the in- [[E]] = E E− , (2.32) duced quadrupole component of the gravitational po- − tential and the quadrupole external tidal potential, re- for the generic quantities E± defined on the thin shell. spectively. The tidal Love number, λ, is then defined by From the first junction condition, Eq. (2.31), and from (see, e.g., [2, 48]) the thin-shell location, Eq. (2.28), we get the following 8 D matching conditions (for details, see Ref. [42]); λ = 1 M5. (2.26) 45 D2 + 2 f (R) 2J (A−) = , ω− = C = , (2.33) f (R) 1 R3 A B. Junction conditions for gluing the two spacetimes − − Θ [[ξ0]] = 0, (l−) = 0, (2.34) + To match the two spacetimes in a physically appro- Rξ− Mξ [[h (R)]]+ 0 + 0 = 0, (2.35) priate way at the location of the thin shell, we need to 0 L2 f (R) R2 f +(R) impose the so-called junction conditions [49, 50]. − ξ The location of the thin shell is given by 2 + k (R) = 0, (2.36) R 2 µ µ a (x±) =(x±) (y ), (2.27) + Rξ2− Mξ2 a [[h2(R)]]+ 2 + 2 + = 0. (2.37) where y are the intrinsic coordinate functions of the L f −(R) R f (R) thin shell. Here and henceforth, the roman letter indices (a, b, c, . . . ) are used to indicate the tensor quantities de- The second junction condition determines the stress fined on the three-dimensional hypersurface of the thin energy tensor S on the thin shell, shell. We denote the intrinsic coordinate functions as ab ya =(T, Θ, Φ) and define the thin-shell location as 1 µ 2 2 Θ Sab = ([[Kab]] hab[[K]]) , (2.38) (x±) =(A±T, R + ǫ ξ±(Θ), Θ + ǫ (l±) (Θ), Φ) 8π − + (ǫ3), (2.28) O in terms of the extrinsic curvature Kab, which is defined µ where A± and R are time rescalings and the radius of in terms of the tangent vectors ea and of the unit normal the gravastar for unperturbed spherical states, respec- vector nµ to the thin shell as Θ tively, and ξ± and (l±) stand for functions of Θ related to the displacement of the thin-shell position. The radius α β ab Kab nα;βea eb and K = habK . (2.39) of spherically symmetric gravastars R satisfies L > R ≡ − and 2M < R because the gravastars possess no hori- Here, the semicolon (; ) denotes the covariant derivative zon. By using the degree of freedom of the coordinate associated with the metric function g . The nonzero choice [42], we may assume, without loss of generality, αβ Θ components of the normal vector to the thin shell up to that A+ = 1 and (l+) = 0. To achieve separation of an accuracy of ǫ2 are given by variables, we also decompose the displacement in Leg- endre polynomials, 1 2 m0 + m2P2(cos θ) 4 ξ(Θ) = ξ0 + ξ2P2(cos Θ). (2.29) nr = + ǫ 3 + (ǫ ), (2.40) f r f 2 O A set of three independent tangent vectors to the thin p 2 ξ2 4 µ nθ = ǫ ∂θ P2(cos θ) + (ǫ ). (2.41) shell, ea , is given by − f O µ µ ∂x e = . (2.30) p a ∂ya The nonzero components of the extrinsic curvature up 6 to an accuracy of ǫ are then given by 1. Master equations for spherically symmetric unperturbed solutions T f ′ 2 KT = + (ǫ ), (2.42) − 2 f O In the unperturbed state, the only nonzero component of ua is uT = 1/ f , then we get Θ Φp f 2 KΘ = KΦ = + (ǫ ), (2.43) R + − p O pf − f σ0 = − , (2.50) Φ f (2ω + Rω′) Rω f ′ 3 4πR KT = ǫ − A + (ǫ ), (2.44) p p 2R f O 1 R M 1 2R2/L2 2 p0 = 2 − R − . (2.51) T R ω′ p2 3 8πR f + − f KΦ = ǫ sin Θ + (ǫ ), (2.45) − ! − 2A f O Note that the pressurep of the thin shellp in the unper- 4 f + Rp f ′ 2 turbed state is positive for any positive values of R and K = + (ǫ ). (2.46) + − 2R f O M provided f > 0 and f − > 0.
The ǫ2-order extrinsicp curvature is lengthy and is given in Appendix A. To specify the second junction condi- 2. Master equations for the ǫ-order solutions tions concretely, as argued later, we need to prescribe some equation of state for the thin-shell matter. In the case of a slowly rotating gravastar, we assume the thin shell to be uniformly rotating, i.e., the angular velocity Ω of the shell is constant. Then, the components C. Matter properties of the thin shell: Perfect-fluid thin of ua satisfies shell Φ Φ u d Ω Θ T = = , u = 0. (2.52) The energy density σ of the thin shell is defined by the u dT eigenvalue of the stress-energy tensor S , We further assume that ǫ Ω/Ω 1, where Ω = ab ≡ k ≪ k √M/R3 is the Keplerian frequency of the spherical Saub = σua, (2.47) gravastar. From the T and Φ components of Eq. (2.47), b − we get where ua is the matter velocity tangent to the thin shell a Φ satisfying uau = 1. By using the projection tensor Ω S T 3 q = h + u u , we− may define the projected stress ten- = Φ + (ǫ ). (2.53) ab ab a b − S Φ + σ0 O sor of the thin shell, namely Then the angular momentum is given by γ = Scdq q . (2.48) ab ac bd f +R (R 3M) f J = R2 − − − Ω . (2.54) We consider a perfect-fluid thin shell, whose stress- + k p f − + 2 f p energy tensor reads γab = pqab with p being the isotropic pressure of the fluid thin shell. We also expand Detailed calculationsp can bep found in Ref. [42]. The mo- the energy density and the pressure of the thin shell as ment of inertia of the rotating gravastar, I, is defined by I ǫJ/Ω. Thus, we have follows: ≡ + 2 4 J f R (R 3M) f − σ = σ0 + ǫ δσ + (ǫ ) I = = − − R2 . (2.55) O Ω + 2 4 k p f − + 2 f p σ0 + ǫ (δσ0 + δσ2P2) + (ǫ ) , ≡ O (2.49) Note that this definition of I is different from that used p = p + ǫ2δp + (ǫ4) p p 0 O in Ref. [41], where ω at the surface of the gravastar is p + ǫ2(δp + δp P ) + (ǫ4) , 2 Ω ≡ 0 0 2 2 O used instead of (cf. Eq. (2.33) in the σ0 0 limit). Note that, in the case of a tidally deformed→ gravas- where σ0 and p0 are the energy density and pressure for tar, (ǫ) perturbations vanish because the background the spherical gravastar, respectively, and δσ and δp are is nonspinningO (i.e., ω = 0). Thus, in such case we have the energy density and pressure perturbations about the J = 0 and Ω = 0. spherical gravastar, respectively. The latter perturba- tions are also expanded in Legendre polynomials P2 P (cos θ). By combining the perfect-fluid condition with≡ 2 2 the first junction condition, summarized in Eqs. (2.33)– In the σ0 = 0 limit investigated in Ref. [41], f = f+ and Eq. (2.55) simply yields I = MR2, in contrast with the− value of Ref. [41], (2.37), and by specifying the equation of state for per- I = R3/2. While both definitions agree in the black-hole limit, the turbations δp = δp(δσ)=(dp/dσ) δσ, as argued later, one adopted here has the correct dimensions of a moment of inertia. we can obtain the quadrupolar deformations of slowly This also implies that a quantitative comparison between our results rotating and tidally deformed gravastars. and those presented in Ref. [41] is not possible. 7
2 2 2 3Ω 2 5Ω 3. Master equations for the ǫ -order solutions of quadrupole ǫ sin θ 2(2J R k) JL R k γ− = − − perturbations 2 3 + 2 16πR ( R f − L f − J(3M 2R) + MR3Ω ξp The matter three velocity of the thin shell, ua, satisfies + − k + 3 2 + (ǫ4). the normalization condition R f + ! "" f ##) O p p (2.60) h uaub = 1, (2.56) ab − Since we assume that the shell is composed of a perfect fluid, the projected stress tensor must be proportional to Φ Φ and therefore the non-zero components of ua up to ǫ2- the projection tensor, i.e., γab = pqab. The ( , ) and Θ Θ order read ( , ) components of the projected tensor are Φ Φ γΦ = pqΦ + + Θ 2 Ω + 2 2 Θ Φ 1 ( f )′ξ ( ) + R ( k ω ) sin 2 4 uT = + ǫ2 − − =(p0 + ǫ δp)(1 + u uΦ) + (ǫ ) + 3 O f + 2 Φ 4 2 f = p + ǫ δp + p u uΦ + (ǫ ), (2.61) 0 0 O p h+(R, Θ) p Θ Θ 2 4 + (ǫ4) , (2.57) γΘ = pqΘ = p0 + ǫ δp + (ǫ ). (2.62) − f + O O ! Then, the condition for a perfect-fluid thin shell reduces Φ 1 u = Ω + (ǫ3) . p to f + O Φ 2γ− = p u uΦ. (2.63) − 0 p The energy density of the shell, σ, is calculated through The explicit form of this condition is given by Eq. (2.47). The explicit forms of the second order quan- 2 Ω ξ2 2 4J 2J(M + R) k 2 2 tities δσ0 and δσ2 are given in Appendix B. The non-zero = + MR Ω + R3 − R k components of the projected stress tensor read "" f ## 3 f p p 4πRp JL2 R5Ω 2(2J R3Ω ) 0 Ω − k . T T 4 T T 3 k k 2 3 γ T = p0 q T + (ǫ ), γ Φ = p0 q Φ + (ǫ ), − − 3 − 3L R f − ! Φ Φ O O γ = p q + (ǫ3), (2.64) T 0 T O p Θ + Φ + γΘ = γ + γ− , γΦ = γ γ− , (2.58) The explicit forms of the second order quantities δp0 and − 2 δp2 are given in Appendix B. The ǫ -order solutions for where we have defined quadrupole perturbations of a slowly rotating thin-shell gravastar can be fully computed by determining the + four unknown constants, B, C3, ξ2 , and ξ2−, which are ǫ2 2(Ω R3 2J) JL2 R5Ω γ+ = p k − − k solutions of the set of four coupled linear algebraic equa- 0 2 3 + 2 − 8πR ( 3R f L f − dp tions, Eqs. (2.36), (2.37), (2.64), and δp2 = δσ2 3 dσ J(3M 2R) + MR Ω p ξ− k 2 + + 0 (the explicit form of δp2 and δσ2 is given in Appendix + − R f (h0 )′ 3 R f + ! − − B). Note that, together with the quadrupolar perturba- q f − tions, the junction conditions allow to fully determine 2 p 2 + + (3M 3MR + R )ξ0 + R(R M)m0 p also the spin-induced spherically symmetric deforma- + − 3 − R2 f + ) tions to (ǫ2). We neglect these deformations here, de- tails areO given in Ref. [42] ǫ2P 2(JL2 Rp5Ω )(R3Ω 2J2) + 2 − k k − Let us now consider a spherically symmetric, tidally 2 2 3 + 8πR ( 3L R f − f deformed gravastar. In this case, we have to impose the conditions J = 0 and Ω = 0 in Eq. (2.58). Thus, the con- + R2[[ f (h + k )]] p 2′ 2′ dition for the perfect fluid thin shell, Eq. (2.64), simply 2 2 2 2 (3Rp 2L )ξ2− +(L 2R )m2− reduces to + − 3 − L2 f − ξ2/ f = 0. (2.65) 3Ω 3Ω 2(2J R k)(pJ(3M 2R) + MR k) hh p ii + − − Note that ǫ2 in this case is just a bookkeeping parameter 3R4 f + and can be factored out. To have a unique solution, we (3M2 + 3MR 2Rp2)ξ+ + R(R M)m+ set a value of D [see, Eq. (2.25)] as − 2 − 2 + (ǫ4), 2 − 2 +3 O R f ) M2 (2.59) D2 = . (2.66) p 3R2 8
The solutions with D2 > 0 represent the tidally de- speed of the fluid, vs, is defined by formed stars by the binary companion on the symmetry 3 2 dp 1 p axis As in the spin-induced case, the ǫ -order solutions v2 = 1 + 0 + (ǫ2). (3.2) s ≡ dσ n σ O of quadrupole perturbations for tidally deformed thin- 0 shell gravastars can be found by determining the four The equation of state for the perturbation is also given unknown constants, D , C , ξ+, and ξ , which are so- 1 3 2 2− by the adiabatic relation lutions of the set of coupled four linear algebraic equa- dp 1 p tions, Eqs. (2.36), (2.37), (2.65), and δp2 = δσ2. 2 0 2 dσ δp = vs δσ = 1 + δσ + (ǫ ). (3.3) n σ0 O Tidally-induced monopole terms are not given here but can be computed through the same procedure. 2 Note that the speed of sound vs is necessary to close the 2 system of junction conditions through δp2 = vs δσ2, as discussed in the previous section. III. NUMERICAL RESULTS
When presenting our numerical results, we often em- B. Spherically symmetric unperturbed solutions and their ploy dimensionless quantities in terms of the length stability scale of the de Sitter horizon radius, L. In other words, we use units such that L = 1. In this study, we construct a one-parameter family of equilibrium unperturbed gravastars for fixed equa- tion of state, the parameter being the compactness of the A. Equation of state for the thin-shell matter gravastar, c M/R. In the unperturbed configuration, the energy density≡ and pressure of the thin shell can be To obtain slowly rotating and tidally deformed solu- written in terms of the three physical quantities c, M, tions for gravastars with a perfect-fluid thin shell, we and L, as need to assume some equation of state for the thin-shell matter. In Ref. [42], sequences of equilibrium solutions M2 √ 1 c2 L2 1 2c characterized by a constant fixed value of M/L are con- σ = − − − c, 0 q 4πM sidered and, from these sequences of equilibrium so- 2M2 lutions, the corresponding equation of state is deter- c 1 c 1 c2L2 p0 = − − . (3.4) mined. In practice, for some fixed value of M/L, the 2 8πM √1 2c − 1 M energy density and pressure of the thin shell in an equi- − − c2L2 librium state are, respectively, given by σ = σ (R) and q 0 0 By imposing the equation of state (3.1) and by using the p = p (R), with R being the radius of the thin shell. 0 0 two above relations, we obtain an algebraic relation be- Thus, for each value of M/L, we obtain a relation be- tween M and c. Thus, for fixed values of n, k, and c, we tween σ and p , which is used as the equation of state. 0 0 can compute the mass of the gravastar, M = M(n, k; c). The latter, however, changes with the compactness of The equilibrium solution of the unperturbed spherical the object. In Ref. [41], a thin shell with vanishing en- gravastar is therefore characterized by n, k, c, and M. ergy density is assumed, i.e., σ = 0 is employed as the 0 Giving a set of fixed parameters, n and k, we obtain se- equation of state for the thin shell. quences of equilibrium solutions by changing values of In this study, we consider a more generic and realis- c. The effects of (spin- or tidal-)induced deformations of tic configuration and assume that the thin shell is com- the thin-shell gravastar along several sequences of equi- posed of a polytropic fluid, whose equation of state is librium solutions characterized by values of k and n are given by discussed later. + 1 Let us consider the radial stability of spherical unper- p = k σ1 n , (3.1) turbed gravastars with a polytropic thin shell. Here, we where k and n are positive constants and are kept fixed basically rely on the dynamical stability analysis given along a sequence of equilibrium solutions. The sound by Visser and Wiltshire [33]. In their treatment, the mo- tion of the spherical thin shell is determined by the en- ergy equation 1 3 We recall that, due to the spherical symmetry of the background R˙ 2 + V(R) = 0, (3.5) configuration, nonaxisymmetric perturbations can be easily ob- 2 tained from the axially symmetric ones, the only difference is in the angular decomposition, cf. e.g. [48]. Note also that the Love num- where the shell radius R is extended to be a function ˙ bers do not depend on the value of D2 chosen in Eq. (2.66), but only of the proper time of the thin shell, τ, R denotes the on the ratio D1/D2 (cf. Eq. (2.26)). On the other hand, the ellipticity proper time derivative of R, and V stands for the po- discussed later is proportional to the external tidal field. tential, given in Eq. (40) of Ref. [33]. In terms of V and 9
V′ dV/dR, the energy density and pressure of the dy- ≡ k=3 namical thin shell are given by k=5 0.16 k=7 k=3 R2 2M k=5 1 2V(R) 2 1 2V(R) R k=7 σ = − − L − − − , (3.6) 0 q 4πRq 0.12 M M 1 1 2V(R) RV′(R) R p0 = − − − 0.08 8πR 1 2V(R) 2M − − R q 2R2 1 2V(R) RV′(R) 2 0.04 − − − L , (3.7) − R2 0.2 0.25 0.3 0.35 0.4 0.45 0.5 1 2V(R) 2 − − L c q respectively. As in a standard dynamical analysis of me- 1 k=3 chanical systems, equilibrium states satisfy V(R) = 0 k=5 k=7 and V′(R) = 0 and their dynamical stability is deter- 0.8 2 2 mined by the sign of V′′ d V(R)/dR for the cor- ≡ 0.6 responding equilibrium solutions, namely they are sta- 2 /c > < 2 ble when V′′ 0 and unstable when V′′ 0. Note 0.4 that Eqs. (3.6) and (3.7) reduce to Eqs. (2.50) and (2.51) if WM V(R) = 0 and V′(R) = 0. 0.2 For a polytropic thin shell satisfying Eq. (3.1), we have 0
dp0 1 p0 dσ0 -0.2 = 1 + , (3.8) 0.2 0.25 0.3 0.35 0.4 0.45 0.5 dR n σ0 dR c and the above relation can be used to extract properties of the potential, V(R), from the equation of sate. FIG. 1. Top: The mass M of the spherical unperturbed gravas- For equilibrium solutions, using Eq. (3.8), we obtain tar with the n = 1 polytropic thin shell as a function of the compactness c. Bottom: The radial stability discriminant, W WM2/c2, for the spherical unperturbed gravastar with the V′′ = , (3.9) = H n 1 polytropic thin shell as a function of the compactness c. Each curve corresponds to the sequence of spherical un- where perturbed gravastars characterized by the same values k. The marginally stable solutions are indicated by marks (square, cir- c2 M2 (3c 2)(1 c) f + cle, or cross mark). W − − ≡ M2 c2L2 f − f + − 3 ( − pf − f (3c2 3c + 1) − p results for the spin-induced and tidal deformations for −3 − p f + these six choices of equations of state. In these figures, 2 2 2 the behaviors of M and WM2/c2 are shown as func- 1p (3c 1)(1 c) c L 2M + 1 + − − − tions of c. As shown in the top panels of Figs. 1 and n f + − c2L2 f − 2, the mass displays a local maximum as a function of 2 f + 2M2(1 3c) the compactness along the sequences of equilibrium so- + − , (3.10) 2 2 + lutions. The maximum masses are, respectively, given f − − c L f f − #) p by M 0.11344, 0.15239, and 0.18070 for the sequences + ≈ f p f − p p (n k) = ( ) ( ) ( ) M H + 2. (3.11) with , 1,3 , 1,5 , and 1,7 , and by + 0.12744, 0.25245, and 0.31035 for the sequences with≈ ≡ p f − p f − (n, k)=(3,1), (3,3), and (3,5). The marginally stable To evaluatep the stability,p it is sufficient to check the sign solutions characterized by W = 0 appear at c 0.42843, W, since H > 0 because of the inequality of arithmetic 0.44461 and 0.45372 for the sequences of solutions≈ with and geometric means. Thus, spherical gravastars with a (n, k)= (1,3), (1,5), and (1,7), and at c 0.40183, polytropic thin shell are stable (unstable) against radial 0.45705 and 0.47146 for the sequences of solutions≈ with perturbations when W > 0(W < 0). (n, k)=(3,1), (3,3), and (3,5) (see the bottom panels Figures 1 and 2 show three equilibrium sequences of Figs. 1 and 2). In Figs. 1 and 2, the marginally stable of unperturbed spherical gravastars characterized by solutions are indicated by marks (square, circle, or cross (n, k)=(1,3), (1,5), and (1,7) and by (n, k)=(3,1), mark, depending on the equation of state). Here and (3,3), and (3,5), respectively. Later on, we present the henceforth, the square, circle, and cross marks are used 10
p0 = 0 can be obtained from the intersection points of 0.3 k=1 | | k=3 Mc and M(c) in the top panels in Figs. 1 and 2, and they k=5 are c 0.47813, 0.47459 and 0.46877 for the sequences of 0.25 k=1 ≈ k=3 solutions with (n, k)=(1,3), (1,5), and (1,7), respec- 0.2 k=5 tively, and c 0.47979 and 0.27241 for the sequences of solutions with≈ (n, k)=(3,1) and (3,3), respectively. M 0.15 The solution with (n, k) = (3,5) does not satisfy the 0.1 dominant energy condition in the range of 0.2 c 0.5. Thus, except for the solutions with (n, k)=(≤3,3≤) and 0.05 (n, k)=(3,5), radially stable solutions satisfy the dom- inant energy condition. 0 0.2 0.25 0.3 0.35 0.4 0.45 0.5 c 3 C. Rescaled physical quantities expressing rotational and k=1 k=3 tidal deformations 2 k=5 Yagi and Yunes [10, 11] have found that for neutron
2 1 star models and quark star models, the dimensionless /c 2 quantities
WM 0 I λ QM I¯ , λ¯ , Q¯ , (3.13) -1 ≡ M3 ≡ M5 ≡ J2
-2 obey universal relations within an accuracy of a few per- 0.2 0.25 0.3 0.35 0.4 0.45 0.5 cent for a variety of reasonable equations of state. c As mentioned before, here we are interested in whether similar nearly universal relations among I¯, λ¯ , FIG. 2. Same as Fig. 1, but for the spherical unperturbed ¯ gravastar with the n = 3 polytropic thin shell. and Q exist also for thin-shell gravastars, and we wish to investigate their behaviors in the black-hole limit. For slowly rotating Kerr black holes, the angular momen- Ω to indicate the marginally stable solutions. As shown in tum J and angular velocity at horizon H are, respec- a the top panels of Figs. 1 and 2, these marginally stable tively, given by J = a M and Ω = = BH H 2 2 solutions coincide with the maximum mass solutions for r+ + a the equilibrium sequences characterized by fixed values a 3 2 + (a ), where a, MBH, and r+ are the spin pa- of n and k within the accuracy of the numerical calcula- 4MBH O tion. Thus, the critical point of the radial stability for the rameter, mass, and event-horizon radius of the black spherical thin-shell gravastar seems to coincide with the hole, respectively. Therefore, the moment of inertia maximum mass solutions along the sequences of equi- J 3 2 reads IBH = = 4MBH + (a ). Thus, for slowly librium solutions characterized by a single equation of ΩH O state, or fixed values of n and k, which is the usual situa- rotating Kerr black holes, we have I¯ = 4 + (a2). Fur- tion for generic relativistic polytropic spheres. Although thermore, as argued by Binnington and PoissonO [54], the this result strongly suggests that the same stability crite- tidal Love numbers of a nonrotating black hole are zero. rion that applies to ordinary self-gravitating polytropic Thus, we have λ¯ = 0. Finally, from Eq. (2.11), the di- fluids also works for thin-shell gravastars, it would be mensionless quadrupole moment, Q¯ , is given by interesting to find a generic proof for this behavior. Finally, let us discuss the energy conditions for the 8BM3 Q¯ = 1 + . (3.14) thin shell. We can get the upper limit of gravastar mass 5J2 Mc to satisfy the dominant energy condition by squar- ing σ0 p0 0 twice. The square of the upper limit As mentioned before, regularity at the event horizon of 2 −| | ≥ Mc is given as a slowly rotating Kerr black hole implies B = 0 and, therefore, Q¯ = 1 in the black-hole limit. M2 3 2 √ 4 + 4 2 c = c2 G G G , On the other hand, for standard neutron-star models, L2 4 − 8 − 8 ! (3.12) the centrifugal force makes the star oblate, which cor- 1 c responds to B > 0. Thus, Q¯ > 1 for standard rotating = √1 2c + − . neutron stars. This is not the case for the slowly rotating G − 2√1 2c − thin-shell gravastars as discussed in Refs. [41, 42] and The explicit critical values of the compactness for σ below. 0 − 11
1.5 40 k=3 k=3 k=5 k=5 1 k=7 k=7 k=3 30 k=3 k=5 k=5 0.5 k=7 k=7 2 2 ε ε / /
2 0 2 20 e e
-0.5 10 -1
-1.5 0 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.2 0.25 0.3 0.35 0.4 0.45 0.5 c c
1.5 20 k=1 k=1 k=3 k=3 1 k=5 15 k=5 k=1 k=1 k=3 k=3 k=5 k=5 0.5 10 2 2 ε ε / / 2 2 e e 0 5
-0.5 0
-1 -5 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.2 0.25 0.3 0.35 0.4 0.45 0.5 c c
FIG. 3. The square of ellipticity, e2/ǫ2, of the thin shell as func- FIG. 4. Same as Fig. 3, but for the tidally deformed gravastars. tion of the compactness c for slowly rotating gravastars with n = 1 (top panel) and n = 3 (bottom panel) polytropic thin shells. Each curve corresponds to the sequence of spherical unperturbed gravastars characterized by the same values of k. The marginally stable solutions are indicated by marks (square and circle). For the rotationally deformed thin shell, the square of the ellipticity is negative for less compact models, changes sign as the compactness increases, and finally D. Rotational and tidal deformations of thin-shell gravastars converges to a single positive value in the nearly black- hole limit (see Fig. 3). As shown in Fig. 3, the values and the behavior of e2/ǫ2 in the nearly black-hole limit Let us first examine the deformation of the thin shell strongly depend on the polytrope index n. In particular, on the surface of the gravastar. Such deformation of the the ellipticity is generically discontinuous in the black- spherical thin shell due to centrifugal and tidal forces hole limit, as observed in Ref. [41] for the σ = 0 case. may be well described by the square of the ellipticity, 0 For ordinary rotating perfect-fluid stars, the ellipticity at defined by the stellar surface is always positive, which reflects the fact that the centrifugal force typically makes the star e2 ξ 3 k + 2 + (ǫ2) . (3.15) oblate. As discussed before and in Refs. [41, 42], the sit- ǫ2 ≡ − 2 R O uation for the thin-shell gravastar is different from that for the fluid star, i.e., the centrifugal force makes the Here, we use the standard definition of the ellipticity for thin shell gravastar prolate at small compactness. Al- rotating stars given in, e.g., Ref. [44]. Thus, we have e2 > though this is counterintuitive, a similar situation is also 0(e2 < 0) for oblate (prolate) spheroids. observed in the case of rotating stars with anisotropic In Figs. 3 and 4, the squares of the ellipticity of the pressure [20]. As argued in Ref. [20], we think that thin shell respectively associated with the rotational and a key element is the strongly anisotropic stress of the tidal deformations for the six sequences of unperturbed thin-shell matter (for the case of the thin shell, strongly spherical gravastar solutions previously presented are anisotropic stress means the absence of the radial pres- given as functions of the compactness c. In each figure, sure) though the reason for occurrence of rotating pro- the results for the n = 1 and n = 3 polytropic thin shells late gravastars is not yet clear. A different explanation is are presented in the top and bottom panels, respectively. given in Ref. [41] in terms of the peculiar p = ρ equa- − 12
4 tion of state of the de Sitter interior . 0 The case of tidal deformations is slightly more compli- cated than that of spin-induced deformations (cf. Fig. 4). -100 k=3 k=5
For a polytropic thin shell with n = 1, (cf. top panel of λ k=7 -200 k=3 Fig. 4), the ellipticity is a positive, decreasing monotonic k=5 function of the compactness, and converges to a single k=7 -300 BH value in the nearly black-hole limit. For a polytropic thin shell with n = 3 (cf. bottom panel of Fig. 4), the ellip- 0.2 0.25 0.3 0.35 0.4 0.45 0.5 ticity is positive for less compact models, changes sign c 20 as the compactness increases, and finally converges to a k=3 single negative value in the nearly black-hole limit. For 16 k=5 2 2 k=7 the n = 3 case, behavior of e /ǫ is very different from, k=3 3 12 k=5 and more complicated than, the n = 1 case. I/M k=7 For tidally deformed perfect-fluid stars, the ellipticity 8 BH at the stellar surface is always negative, which means 4 that the tidal force in the present setup makes a perfect 0.2 0.25 0.3 0.35 0.4 0.45 0.5 fluid star prolate (see the arguments given at the end c of Sec. II C). On the other hand, for tidally deformed thin-shell gravastars, the direction of the elongation of 1 the thin shell on the meridional cross section is turned 0.5 k=3 by 90 degrees relative to the standard case. The reason 2 k=5 k=7 0 for occurrence of this unusual effect is likely the same as QM/J k=3 that for the case of the rotationally deformed prolate thin k=5 -0.5 k=7 shell, i.e., two key elements are the strongly anisotropic BH stress of the thin-shell matter (cf. Ref. [20]) and the pe- 0.2 0.25 0.3 0.35 0.4 0.45 0.5 culiar equation of state of the de Sitter fluid making the c gravastar interior [41]. FIG. 5. The dimensionless tidal Love number, λ¯ (top panel), the dimensionless moment of inertia, I¯ (middle panel), and E. I-Love-Q relations for gravastars with polytropic thin the dimensionless quadrupole moment, Q¯ (bottom panel), as shell functions of the compactness c for gravastars with n = 1 poly- tropic thin shell. Each curve corresponds to the sequence of Let us now examine how the I-Love-Q triad depends spherical unperturbed gravastars characterized by the same values of k. The marginally stable solutions are indicated by on the compactness. The quantities λ¯ , I¯, and Q¯ are marks (square and circle). shown as functions of c in Figs. 5 and 6 and Figs. 7 and 8 for the sequences of the spherical gravastars with the n = 1 and n = 3 polytropic thin shells, given in Figs. 1 and 2, respectively. Figures. 6 and 8 are zoomed-in ver- find Q¯ < 1 and λ¯ < 0, which is related to the gravastar sions of the same quantities of Figs. 5 and 7 in the range being prolate for rotational deformations and oblate for 0.45 . c 0.5, in order to better visualize their behavior tidal deformations. Note that a similar behavior is ob- ≤ in the nearly black-hole limit. In each figure, we show served for thin-shell gravastar with σ0 = 0 [41]. This the corresponding quantities for the case of a slowly ro- property is consistent with the fact that less compact tating Kerr black hole with dash-dotted or dotted hori- models of slowly rotating (tidally deformed) thin-shell zontal lines. gravastar have a prolate (oblate) shaped thin shell, as From the middle panels of Figs. 5–8, we see that the previously discussed. From the bottom (top) panels of dimensionless moment of inertia I¯ is always positive, Figs. 5–8, we also see that Q¯ (c)’s (λ¯ (c)’s) obtained in this and monotonically decreases to the black-hole value of study are basically monotonically increasing functions I¯ as the compactness increases. These basic properties which approach the value corresponding to the black- seem to be common also for other slowly rotating ob- hole limit as c 1/2. An exception appears in the jects, e.g. slowly rotating neutron stars. We also observe case of the very compact→ models with (n, k)=(3,1) for that the function I¯(c) is only mildly dependent on the c & 0.45. For this exceptional sequence of unperturbed equation of state for the thin-shell matter. spherical gravastars, Q¯ (λ¯ ) exceeds unity (zero) around The functions Q¯ (c) and λ¯ (c) are respectively given in c 0.46 (c 0.45), but once again becomes less than the bottom and top panels of Figs. 5–8. For n = 1, we unity∼ (zero)∼ around c 0.49 (c 0.495) (cf. bottom (top) panel of Fig. 8). In∼ other words,∼ in this case Q¯ (c) and λ¯ (c) cross the black-hole value twice, and then ap- proach the black-hole limit. 4 We thank Leo Stein for suggesting this explanation. In Fig. 9, we present the I-Love-Q relations for gravas- 13
0 0 -0.1 -50 k=1 -0.2 k=3
λ k=3 λ -100 k=5 -0.3 k=5 k=1 -0.4 k=7 k=3 k=7 -150 k=5 -0.5 BH BH -200 0.45 0.46 0.47 0.48 0.49 0.5 0.2 0.25 0.3 0.35 0.4 0.45 0.5 c c
k=3 20 k=1 4.4 k=5 k=3 k=7 16 k=5 k=7 k=1 3 3 BH 12 k=3 I/M 4.2 I/M k=5 8 BH
4 4 0.45 0.46 0.47 0.48 0.49 0.5 0.2 0.25 0.3 0.35 0.4 0.45 0.5 c c
1 1
k=1 2 0.96 2 k=3 k=3 0.5 k=5
QM/J 0.92 k=5 QM/J k=1 k=7 k=3 k=7 k=5 0.88 BH 0 BH 0.45 0.46 0.47 0.48 0.49 0.5 0.2 0.25 0.3 0.35 0.4 0.45 0.5 c c
FIG. 6. Same as Fig. 5, but for the strongly relativistic regime FIG. 7. Same as Fig. 5, but for the gravastar with the n = 3 (0.45 c 0.5). polytropic thin shell. ≤ ≤ tars with a polytropic thin shell. Top and bottom panels correspond to the I¯– λ¯ and the Q¯ – λ¯ relations, respec- tively. We observe various| | interesting| | | | properties. First, the I-Love-Q relations for thin-shell gravastars are dras- tically different from those of neutron stars and quark havior of matter is well understood. In the case of thin- stars. For example, λ¯ and Q¯ for the gravastar case can shell gravastars, the equation of state changes precisely be negative (note that we use λ¯ and Q¯ instead of λ¯ at the surface and, in this sense, it is reasonable to expect | | and Q¯ to represent the I-Love-Q relations in Fig. 9). Sec- a different behavior. It would be interesting to extend ond, we confirm that the relations approach their cor- our study to the case of finite-thickness gravastars, in responding black-hole limit [41] as the compactness in- order to investigate the region which is (possibly) most creases (or, equivalently, as λ¯ decreases). This property responsible for the emergence of the universality. seems to hold in general, i.e.| | for all values of k and n considered in this study. Finally, we observe that the I-Love-Q relations for In the present study, we allow for a wide range of thin-shell gravastars substantially depend on the equa- polytropic equations of state for the thin shell because tions of state for the thin shell. Although all curves the properties of the latter are basically unconstrained. approach the same black-hole limit, they behave differ- It is therefore natural to expect larger deviations in the ently for different polytropic indices n and even for the universality of the I-Love-Q relations relative to the neu- same polytropic index but different values of k. There- tron star case (whose equation of state, as a qualitative fore, there seems to be no universality in the I-Love-Q reference, can be roughly approximated by a polytrope relations for thin-shell gravastars. However, even in the with 0.5 . n . 1). In fact, we expect even larger de- cases of neutron stars and quark stars the universality viations from universality for polytropic indices n > 3. does not hold for any equation of state, but only for a If the equation of state for the thin shell is moderately subclass of equations of state which nonetheless com- restricted, therefore, it is possible that the I-Love-Q re- prises the majority of models, ranging from soft to stiff lations for thin-shell gravastars will become moderately matter. These models differ at very high densities but universal, the degree of the universality would depend are similar near the neutron-star surface, where the be- on the class of equations of state considered. 14
0.1 0 n=1,k=3 n=1,k=5 -0.1 k=1 n=1,k=7 1 λ k=3 -0.2 k=5 3 10 n=3,k=1 k=3 n=3,k=3
-0.3 I/M k=5 n=3,k=5 -0.4 BH BH 0.45 0.46 0.47 0.48 0.49 0.5 c 100 4.6 k=1 n=1,k=3 k=3 | 2 -1 k=5 10 n=1,k=5 4.4 k=3 3 n=1,k=7 k=5
I/M -2 n=3,k=1 4.2 BH |QM/J 10 n=3,k=3 n=3,k=5 -3 4 10 -2 -1 0 1 2 0.45 0.46 0.47 0.48 0.49 0.5 10 10 10 10 10 c |λ |
1 FIG. 9. The I¯– λ¯ (top panel) and the Q¯ – λ¯ (bottom panel) 0.98 relations for a| gravastar| with a polytropic| | | thin| shell. The 2 0.96 k=1 marginally stable solutions are indicated by marks (triangle k=3 0.94 QM/J k=5 and diamond). 0.92 k=3 k=5 0.9 BH 0.45 0.46 0.47 0.48 0.49 0.5 nal [51]. c Tidal deformations enter the gravitational waveform formally at fifth post-Newtonian order, although they FIG. 8. Same as Fig. 7, but for the strongly relativistic regime are enhanced by terms proportional to c5 [47, 48]. In the (0.44 . c 0.5). ≤ stationary phase approximation, the tidal contribution to the gravitational-wave phase reads [47]
IV. GRAVITATIONAL-WAVE CONSTRAINTS ON THE 5 9 v λ1 TIDAL DEFORMABILITY OF THIN-SHELL δΨ = (11m + m) + 1 2 , (4.1) − 16 µm4 2 m ↔ GRAVASTARS 1 where v is the orbital velocity, m and µ are the total and Since the tidal Love numbers of static [54, 55] (and, reduced mass of the binary, respectively. Clearly, this presumably, also rotating [28, 56]) black holes are iden- term is zero if the two objects are black holes, because tically zero, any measurement of a nonzero tidal Love λ1 = λ2 = 0. number of a compact object implies that the latter is Thus, an interesting question is to what extend a pu- not a black hole. This fact is particularly important for tative upper limit on λ can constrain gravastar mod- gravitational-wave tests of black-hole mimickers and of els. The merger frequency and individual masses of exotic compact objects. In the near future, gravitational- GW150914 suggest a lower limit on the compactness of wave detections of compact-binary inspirals [51, 52] will the two bodies, c & 0.25. Because λ vanishes in the be routine; such observations may be used to put upper black-hole limit, an upper bound on λ can be converted bounds on the tidal Love number of the two bodies, thus into a lower bound on c. From our results in Figs. 5 and constraining exotic alternatives. 7, we observe that λ¯ . 20 100 when c & 0.25, the ex- The first post-Newtonian corrections to the inspiral act value depending| | on the− equation of state of the thin phase of two compact bodies depend only on their shell5. Therefore, an upper bound at the level of masses and spins, but not on their internal structure. For example, the chirp mass 30 M of the event 5 M ≈ ⊙ 18 mi 5 GW150914 places a lower bound on the total mass of λi . (1 7) 10− s , (4.2) m = m + m & M | | − × 30 M the system, : 1 2 70 . On the other hand, ⊙ the merger frequency indicates that⊙ the two bodies have to be very compact [51]. In other words, the objects that produced the gravitational-wave events detected so far are (at least) as compact as a neutron star but much more 5 In the large compactness limit, the Love number is well fitted by massive. This suggests a pair of merging black holes as the expression λ¯ (c) a(c 1/2)b, where a and b are constants that ∼ − the most natural source of the gravitational-wave sig- depends on the equation of state. 15 would provide a lower bound on ci := mi/Ri which state, at the level of a few percent. Although (some de- is more stringent than the approximate one set by the gree of) approximate universality might be restored by merger frequency. restricting to a subclass of equations of state, it is never- A detailed analysis of the possible gravitational-wave theless interesting that thin-shell gravastars provide an constraints on λ goes beyond the scope of this work. example of very compact objects for which the univer- Furthermore, because upper limits on the tidal deforma- sality of the I-Love-Q relations is manifestly broken. bility of the binary systems of GW150914 and GW151226 We also evaluated the stability of the unperturbed have not been published yet, it is difficult to estimate the spherical gravastars with a polytropic perfect-fluid thin constraints on λi for these systems. However, similar shell against small radial perturbations. It is found that analyses for neutron-star binary systems have been per- less compact models are stable and that the stability formed before the recent gravitational-wave detection changes at the maximum mass models of equilibrium (cf., e.g., Ref. [58]). By assuming similar absolute errors, sequences characterized by a single equation of state. one might expect an observational bounds at the level The stability analysis implies that the gravastars become 23 5 of λ . 10− s which, for an object of M 30M , unstable as they approach the black-hole limit. There- | | 4 ≈ ⊙ corresponds to λ¯ . 10− . This would provide a very fore, we inevitably had to use unstable models to inves- stringent lower| bound| on the compactness and can po- tigate the I-Love-Q relations very close to the black-hole tentially rule out several equations of state for which a limit. The instability region shrinks to zero as k increases very compact thin-shell gravastar is radially unstable. for a fixed value of n. The k ∞ limit corresponds → In fact, as clear from the results previously presented, to the case of a thin shell with vanishing energy den- all models of gravastars considered in this paper would sity studied in Ref. [41]. Despite some differences in our 4 be radially unstable when λ¯ . 10− . analysis (we corrected the definition of the moment of | | inertia used in Ref. [41] and kept the properties of the thin-shell matter fixed along different sequences of solu- V. CONCLUSION tions) our results are qualitatively similar to those pre- sented in Ref. [41] and extend the latter to more realistic We have studied the rotational and tidal quadrupole configurations. deformations of a thin-shell gravastar, which are char- In this study, we focus on the case of a perfect-fluid acterized by the tidal Love number and the rotational thin shell. However, a thin shell with anisotropic pres- quadrupole moment, respectively. We worked in a per- sure might be more reasonable. For example, if the turbative regime in which the tidal and the rotational phase transition replacing the would-be horizon [32] is effects are described by small deviations from spheri- due to a scalar field, the latter would naturally give rise cal symmetry. We considered a thin shell made of a to anisotropic stresses. Thus, it would be interesting polytropic fluid and studied in detail the cases when the to investigate in detail the case of gravastars with an polytropic indices are n = 1 and n = 3. We found that anisotropic thin shell. Likewise, it would be interesting the I-Love-Q relations of a thin-shell gravastar are dras- to extend our computation to the case of finite-thickness tically different from those of an ordinary compact star shells, although that would likely require a numerical like a neutron star. The Love number and quadrupole integration of the field equations within the shell. moment are negative for less compact models and the I- A detailed analysis of the observational constraints Love-Q relations continuously approach the black-hole on the tidal deformability of gravastars coming from limit. gravitational-wave measurements goes beyond the The appearance of negative values of the Love num- scope of this work. Nonetheless, our results suggest that ber and quadrupole moment means that the elongation near-future constraints on the tidal Love numbers of direction of the matter distribution on the meridional compact objects from inspiral gravitational waveforms cross section is turned by 90 degrees relative to the case can place very stringent lower limits on the compactness of an ordinary compact star. The reason for the appear- of the two objects, ruling out several gravastar models ance of this counterintuitive deformation is not entirely and giving further support to the fact that events like clear, but seems to be related to the strongly anisotropic GW150914 and GW151226 are coalescences of a pair of stress of the thin shell in the horizontal and vertical di- black holes. rections, and to the peculiar equation of state of the So far, studies of rotating models of gravastars are re- gravastar’s interior. Similar counterintuitive proper- stricted to the case of slow rotation. On the other hand, ties have been found also for others infinitesimally-thin it is likely that a substantial fraction of black holes are shells deformed by rotational or tidal effects [41, 42, 59– rapidly spinning, and this is surely the case for black 61]. holes formed in the coalescences recently detected in We considered a wide range of polytropic equations the gravitational-wave band by the LIGO interferome- of state for the thin-shell matter. Within this range, there ter [51, 52]. Highly-spinning gravastars might be unsta- is no universality in the I-Love-Q relations, unlike the ble against the ergoregion instability [40], and a detailed case of neutron stars and quark stars for which such studied is therefore necessary to assess their viability as relations depend only mildly on the stellar equation of exotic compact objects and black-hole mimickers. These 16 investigations remain as future work. Appendix B: The second order energy density and pressure perturbations of the thin shell
The l = 2 energy density and pressure perturbations are, respectively, given by ACKNOWLEDGMENTS 3Ω 2 1 2ξ2− m2− (2J R k) 2 δσ = − − R [[ fk′ ]] 2 2 3 + 2 N. U. acknowledges financial support from MEXT 4πR ( f − − 3R f f − − p Grant-in-Aid for Scientific Research on Innovative Ar- 3Ω 3 J( 2J + R pk) + R m2 3Mp+ 2R + eas “New Developments in Astrophysics Through + − ξ 3 + − + 2 Multi-Messenger Observations of Gravitational Wave R f R f 3Ω 2 Ω Sources” (Grant Number 24103005). This project has re- (R k 2Jp)(J + R (R 3M)pk) ceived funding from FCT-Portugal through the project + − 3 − , (B1) 3 ) IF/00293/2013. 3R f − p 1 2 δp = R f (h′ + k′ ) 2 8πR2 2 2 n hhp2 ii 2 + Appendix A: The second order coefficients of the extrinsic 1 (2R 3MR + 3M )ξ2 + + 3 − (R M)m2 curvature R f + R − − ! 3Ω Rp k 2J 2J(2R 3M) Ω (0) + − − +(R 3M) k) If we expand the extrinsic curvature as Kab = Kab + +3 R3 − (1) 2 (2) 3R f ǫ Kab + ǫ Kab, the second order coefficients of the 3 2 2 2 2 2 (2) (R Ωp 2J) (L 2R )m− (2L 3R )ξ− extrinsic curvature associated with the thin shell, Kab k − + − 2 − − 2 . is given as follows: − 3R3 f + f 2 3 − L f − ) (B2) p p (2) T 1 2 KT = 3 2m0 f ′ 4R f h0′ The l = 0 energy density and pressure perturbations R f + − h are, respectively, given by 4 f ( f 1) 4R3 p 2 2 + Rξ0 ( f ′) 2− + ωω′ 1 (Ω 2J) 1 − R 3 δσ = k − + J(2J R3Ω ) + m+ 0 5 + + k 0 2 4πR ( 3 f f − f − + 2m2 f ′ 3R f h′ − 2 3Ω 4 Ω (2J R k)(pJ + R (R p3M) k) n 3 + − − 2 4 f ( f 1) 4R + +Rξ ( f ′) − ωω′ P , 3 f 2 − R2 − 3 2 2(σ + p ) (A1) 0 0 ξ , (B3) − R 0 (2) Θ (2) Φ KΘ + KΦ
1 4 2 2 = 6m + 3ξ (2 f R f ′) R ωω′ 1 (3M 3MR + R )ξ + R(R M)m 2 0 0 0 0 3R f − − δp0 = 2 − 3 − h − 8πR ( R2 f + 2 6mp2 + 3ξ2(2 f R f ′ 6) 6R fk2′ 3Ω 2 Ω − − − − (R k 2J)(2J(2R 3M) + R (R 3M) k) + − − p − n4 + 3 + R ωω′ P2 , (A2) 3 f R 3Ω 2 o i 4 (2J R k) 2 p + + ξ0 (2) Θ (2) Φ 6ξ2 + R ωω′ 2 − R f (h )′ . (B4) KΘ KΦ = sin Θ, (A3) 3 + 0 3 2 − 3R f f − − ) − 2R f − q f − (2) 1 p 2 2 p p K = 3 2m0(4 f + R f ′) 4R f h0′ R2 f − h 2 2 + ξ p4 f (1 + f R f ′) + R ( f ′) 0 − 2 2 + 2m (4 f + R f ′) 4R f (h′ + 2k′ ) 2 − 2 2 n 2 2 + ξ 4 f ( 5 + f R f ′) + R ( f ′) P . (A4) 2 − − 2 o i 17
[1] W. G. Laarakkers and E. Poisson, Astrophys. J.512, 282 [33] M. Visser and D. L. Wiltshire, Class. Quantum Grav. 21, (1999). 1135 (2004). [2] E. Poisson and C.M. Will, Gravity: Newtonian, Post- [34] A. E. Broderick and R. Narayan, Astrophys. J. L21 638 Newtonian, Relativistic, Cambridge University Press (2006); Class. Quantum Grav. 24 659 (2007). (2014). [35] C. Chirenti and L. Rezzolla, arXiv:1602.08759 (2016). [3] W. Israel, Phys. Rev. D 164, 1776 (1967). [36] G. F. Giudice, M. McCullough and A. Urbano, [4] B. Carter, J. Math. Phys. 10, 70 (1969). arXiv:1605.01209 [hep-ph]. [5] H. M ¨uller zum Hagen, D. C. Robinson and H. J. Seifert, [37] V. Cardoso, E. Franzin and P. Pani, Phys. Gen. Rel. Grav. 4, 53 (1973). Rev. Lett. 116, no. 17, 171101 (2016) [6] D. C. Robinson, Gen. Rel. Grav. 8, 695 (1977). doi:10.1103/PhysRevLett.116.171101 [arXiv:1602.07309 [7] E. Berti et al., Class. Quant. Grav. 32, 243001 (2015) [gr-qc]]. [8] J. M. Lattimer and M. Prakash, Astrophys. J 550, 426 [38] J. Keir, arXiv:1404.7036, (2014). (2001). [39] V. Cardoso, L. C. B. Crispino, C. F. B. Macedo, H. Okawa, [9] M. Urbanec, J. C. Miller, and Z. Stuchl´ık, Mon. Not. R. and P. Pani, Phys. Rev. D, 90 044069 (2014). astr. Soc. 433, 1903 (2013). [40] V. Cardoso, P. Pani, M. Cadoni and M. Cavaglia, Phys. [10] K. Yagi and N. Yunes, Science 341, 365 (2013). Rev. D 77, 124044 (2008) doi:10.1103/PhysRevD.77.124044 [11] K. Yagi and N. Yunes, Phys. Rev. D 88, 023009 (2013). [arXiv:0709.0532 [gr-qc]]. [12] A. Maselli, V.Cardoso, V.Ferrari, L. Gualtieri, and P. Pani, [41] P. Pani, Phys. Rev. D 92, 124030 (2015). Phys. Rev. D 88, 023007 (2013). [42] N. Uchikata and S. Yoshida, Class. Quantum Grav. 33, [13] D. D. Doneva, S. S. Yazadiev, N. Stergioulas, and K. D. 025005 (2016). kokkotas, Astophys. J 781, L6 (2014). [43] J. B. Hartle, Astro. Phys. J. 150, 1005 (1967); J. B. Hartle [14] G. Pappas and T. A. Apostolatos, Phys. Rev. Lett. 112, and K. S. Thorne, Astro. Phys. J. 153, 807 (1968). 121101 (2014) [44] K. S. Thorne, General relativity and cosmology, Course [15] S. Chakrabarti, T. Delsate, N. G ¨urlebeck and J. Steinhoff, XLVII, Proceedings of the International School of Physics Phys. Rev. Lett. 112, 201102 (2014) ’Enrico Fermi’, ed. R. K. Sachs, Academic Press, New York [16] K. Yagi, K. Kyutoku, G. Pappas, N. Yunes and T. A. Apos- (1971). tolatos, Phys. Rev. D 89, no. 12, 124013 (2014) [45] S. Chandrasekhar and J. C. Miller, Mon. Not. R. astr. Soc. [17] G. Martinon, A. Maselli, L. Gualtieri and V. Ferrari, Phys. 167, 63 (1974). Rev. D 90, no. 6, 064026 (2014) [46] K. S. Thorne, Phys. Rev. D 58, 124031 (1998). [18] K. Yagi and N. Yunes, Phys. Rev. D 91, 103003 (2015). [47] E.´ E.´ Flanagan and T. Hinderer, Phys. Rev. D 77, 021502(R) [19] K. Yagi and N. Yunes, Phys. Rev. D 91, 123008 (2015). (2008). [20] K. Yagi and N. Yunes, Class. Quantum Grav. 33, 095005 [48] T. Hinderer, Astro. Phys. J. 677, 1216 (2008). (2016) (arXiv:1601.02171). [49] W. Israel, Nuovo Cimento 44B, 1 (1966). [21] B. Haskell, R. Ciolfi, F. Pannarale, and L. Rezzolla, Mon. [50] C. Barrab´es and W. Israel, Phys. Rev. D 43, 1129 (1991). Not. R. astr. Soc. 438, L78 (2014). [51] B. P. Abbott et al. [LIGO Scientific and Virgo Collab- [22] Y. H. Sham, L. M. Lin and P. T. Leung, Astrophys. J. 781, 66 orations], Phys. Rev. Lett. 116, no. 6, 061102 (2016) (2014) doi:10.1088/0004-637X/781/2/66 [arXiv:1312.1011 doi:10.1103/PhysRevLett.116.061102 [arXiv:1602.03837 [gr-qc]]. [gr-qc]]. [23] P. Pani and E. Berti, Phys. Rev. D 90, no. 2, 024025 [52] B. P. Abbott et al. [LIGO Scientific and Virgo Collab- (2014) doi:10.1103/PhysRevD.90.024025 [arXiv:1405.4547 orations], Phys. Rev. Lett. 116, no. 24, 241103 (2016) [gr-qc]]. doi:10.1103/PhysRevLett.116.241103 [arXiv:1606.04855 [24] D. D. Doneva, S. S. Yazadjiev, K. V. Staykov and [gr-qc]]. K. D. Kokkotas, Phys. Rev. D 90, no. 10, 104021 [53] A. Buonanno and B. S. Sathyaprakash, arXiv:1410.7832 (2014) doi:10.1103/PhysRevD.90.104021 [arXiv:1408.1641 [gr-qc]. [gr-qc]]. [54] T. Binnington and E. Poisson, Phys. Rev. D80, 084018 [25] Y. H. Sham, T. K. Chan, L. M. Lin and P. T. Leung, (2009). Astrophys. J. 798, no. 2, 121 (2015) doi:10.1088/0004- [55] N. G ¨urlebeck, Phys. Rev. Lett. 114, no. 15, 151102 (2015) 637X/798/2/121 [arXiv:1410.8271 [gr-qc]]. doi:10.1103/PhysRevLett.114.151102 [arXiv:1503.03240 [26] K. Yagi, L. C. Stein, G. Pappas, N. Yunes, and T. A. Apos- [gr-qc]]. tolatos, Phys. Rev. D, 90, 063010 (2014). [56] P. Pani, L. Gualtieri, A. Maselli and V.Ferrari, Phys. Rev.D [27] K. Yagi, Phys. Rev. D, 89, 043011 (2014). 92, no. 2, 024010 (2015) doi:10.1103/PhysRevD.92.024010 [28] P. Pani, L. Gualtieri, and V. Ferrari, Phys. Rev. D 92, [arXiv:1503.07365 [gr-qc]]. 124003 (2015). [57] N. Uchikata and S. Yoshida, Phys. Rev. D 90, 064042 [29] H. A. Buchdahl, Phys. Rev. 116, 1027 (1959). (2014). doi:10.1103/PhysRev.116.1027 [58] W. Del Pozzo, T. G. F. Li, M. Agathos, C. Van [30] J. B. Hartle, Phys. Rep. 46, 201 (1978). Den Broeck and S. Vitale, Phys. Rev. Lett. 111, no. [31] R. L. Bowers and E. P. T. Liang, Astrophys. J. 188, 657 7, 071101 (2013) doi:10.1103/PhysRevLett.111.071101 (1974). [arXiv:1307.8338 [gr-qc]]. [32] P. O. Mazur and E. Mottola, gr-qc/0109035 (2001); Proc. [59] V. de la Cruz and W. Israel, Phys. Rev. 170, 1187 (1968). Natl. Acad. Sci. U.S.A. 101, 9545 (2004). [60] H. Pfister and K. H. Braun, Class. Quantum. Grav 2, 909 (1985). 18
[61] H. Pfister and K. H. Braun, Class. Quantum. Grav 3, 335 (1986).