arXiv:1011.4127v2 [gr-qc] 6 Dec 2010 † ∗ r aeyhde eida vn oio,adthere- and horizon, they event occur an singularities behind such if hidden (even safely sense singularity) are weak naked rea- a This a a in have (in or cannot sense horizon? strong we event a space-time in absolute sonable formulated an be can in cosmic conjecture singularities, a one naked exist each of occurrence there clothing the does forbids who conjecture: censor cosmic as known horizons. event by that shrouded states space-times always which flat are asymptotically hypothesis, in counterexamples singularities censorship be curvature cosmic would solutions the singulari- such naked of are found, to there If if lead know that ties. to solutions like collapse also black physical would in One results formation. collapse hole what under gravitational and conditions, whether, the know initial to of like rel- one would general One contemporary is in ativity. gravitation, research of of fields theory active most Einstein frame- the the in of matter, of work distribution regular initially an of lcrncades [email protected] address: Electronic lcrncades [email protected] address: Electronic oe ers 1 a h rtt rps h idea, the propose to first the was [1] Penrose Roger collapse gravitational the of fate final the Investigating a crto ikpoete bevtoal distinguis observationally properties disk accretion Can ASnmes 42.C,0.0 w 47.B,04.80.Cc Bw, 04.70. Dw, 04.20. Cv, 04.20. numbers: PACS Kerr-ty from may singularities signatures naked observational rotating these mecha distinguishing Thus, efficient t more holes. than effici much black higher conversion a than always provide the is singularities that naked singularities shown i.e., naked flu also rotating the is by it magnitude fr tion this, of flux to energy orders addition the several charge, In by scalar exceed nak the can and of conse singularity and holes objects, parameter black compact spin between the of discriminate classes could two that distribut natures these temperature for flux, different (energy t geometry, are disks exterior singularity the the naked in of differences rotating properties the a to Ker Due of and velocity parameter. singularities Einstein- angular naked the the in rotating effect, equations between axially- field difference the around jor of disks solutions accretion as K consider obtained t and we singularities through paper, naked is present rotating holes around symmetr black properties axially from disks and singularities si stationary naked naked majo of a entiating context of is the occurrence co singularities In the a naked astronomy. and forbids proposed holes who has black censor Penrose astrophysical cosmic horizon. a event exists an without gularity ae iglrte r yohtclatohsclobjec astrophysical hypothetical are singularities Naked .INTRODUCTION I. eateto hsc n etrfrTertcladComput and Theoretical for Center and Physics of Department h nvriyo ogKn,PkF a od ogKn,P .C R. P. Kong, Hong Road, Lam Fu Pok Kong, Hong of University The .Kov´acs Z. Dtd etme ,2018) 1, September (Dated: singularities? ∗ n .Harko T. and edadasrn udi 2–6.Teclas fthe of collapse The [24–26]. in fluid radiation string a both a include and to extended field field radiat- been Vaidya Einstein has The metric the [23]. ing of [22], derived solutions large been in- known has a which equations, the and approach of fluids, same II most the Type cludes of including framework solutions, of the class In proposed was approach [21]. inverted an in fluid null a of collapse waves, singular. with remain probed singularities classical if stronger clas- regular while the space-times becomes some singularity For sical [20]. particles in with proposed space- than been in naked rather has waves probing analyzed with of singularities idea been times The has [11]-[19]. papers collapse many gravitational spherical [3]-[10]. solu- found nonstatic been two- exact have tions inhomogeneous the reduction, few this is very by with however, Even case obtained symmetry. gravi- relativity such spherical general construct imposing One of to reduction used dimensional models. are collapse metrics symmetries tational intractable, appears special problem with general the tions, far. so successful Un- been have therein). attempts references such and no to fortunately, [2] attempts (see various conjecture been the have prove there proposal, Since s observers). Penrose’ far-away with communicate cannot fore nodrt banteeeg-oetmtno o the for tensor energy-momentum the obtain to order In the models, physical various of framework the Within equa- Einstein full the of complexity the to due Since, ebakholes. black pe o n qiiru aito spectrum) radiation equilibrium and ion observational day present for challenge r r-yebakhls epciey nthe In respectively. holes, black err-type s hrceie yagaiainlsin- gravitational a by characterized ts, ymti oaignkdsingularities, naked rotating symmetric asessaa edter.Afis ma- first A theory. field scalar massless accretion thin of study comparative he ecneso ffiinyfrbakholes, black for efficiency conversion he rvd h eesr ol rmclearly from tools necessary the provide en nesl rprinlt t spin its to proportional inversely being etemdnmcadelectromagnetic and thermodynamic he imfrcnetn asit radiation into mass converting for nism lc oe si h rm dragging frame the in is holes black r unl iigcerosrainlsig- observational clear giving quently nyo h crtn asit radia- into mass accreting the of ency dsnuaiis o pcfi ausof values specific For singularities. ed clgoere,apsiiiyo differ- of possibility a geometries, ical rmteds faKr lc hole. black Kerr a of disk the from x mteds rudarttn naked rotating a around disk the om glrte.Dsigihn between Distinguishing ngularities. † jcue codn owihthere which to according njecture, lc oe rmnaked from holes black h toa Physics, ational hina 2 quark fluid, described by the bag model equation of state objects, thus giving clear observational signatures, which p = (ρ 4B) /3, with B =constant, has been studied in may allow to discriminate, at least in principle, naked [27], and− the conditions for the formation of a naked sin- singularities from black holes. We would like to point gularity have been obtained. The obtained solution has out that the proposed method for the detection of the been generalized to arbitrary space-time dimensions and naked singularities by studying accretion disks is an in- to a more general linear equation of state in [28, 29]. The direct method, which must be complemented by direct gravitational collapse of a high-density null charged mat- methods of observation of the ”surface” of the consid- ter fluid, satisfying the Hagedorn equation of state, was ered black hole/naked singularity candidates, and/or by considered in the framework of the Vaidya geometry in the study of the lensing properties of the central compact [30]. object. Since theoretical studies alone cannot give an answers about the existence or non-existence of naked singulari- ties in nature, the differences in the observational prop- The mass accretion around rotating black holes was erties of the black holes and naked singularities could studied in general relativity for the first time in [36]. By be used to discriminate between these two classes of ob- using an equatorial approximation to the stationary and jects. Such a distinctive observational feature is repre- axisymmetric spacetime of rotating black holes, steady- sented by the lensing properties, and a research project state thin disk models were constructed, extending the that uses gravitational lensing as a tool to differentiate theory of non-relativistic accretion [37]. In these mod- between naked singularities and black holes was initi- els hydrodynamical equilibrium is maintained by efficient ated in [31]. This pioneering work was extended in [32], cooling mechanisms via radiation transport, and the ac- where the distinctive lensing features of black holes and creting matter has a Keplerian rotation. The radiation naked singularities were studied in detail. The gravi- emitted by the disk surface was also studied under the tational lensing due to a strongly naked singularity is assumption that black body radiation would emerge from qualitatively different from that due to a Schwarzschild the disk in thermodynamical equilibrium. The radiation black hole; a strongly naked singularity gives rise to ei- properties of thin accretion disks were further analyzed in ther two or nil Einstein ring(s), and one radial criti- [38, 39], where the effects of photon capture by the hole cal curve. The gravitational lensing (particularly time on the spin evolution were presented as well. In these delay, magnification centroid, and total magnification) works the efficiency with which black holes convert rest for a Schwarzschild black hole and for a Janis-Newman- mass into outgoing radiation in the accretion process was Winicour naked singularities [33] were studied in [34]. also computed. More recently, the emissivity properties The lensing features are qualitatively similar (though of the accretion disks were investigated for exotic cen- quantitatively different) for Schwarzschild black holes, tral objects, such as wormholes, non-rotating or rotating weakly naked, and marginally strongly naked singular- quark, boson or fermion stars, brane-world black holes, ities. However, the lensing characteristics of strongly f(R) type gravity models, and Horava-Lifshitz gravity naked singularities are qualitatively very different from [40]. In all these cases it was shown that particular those due to Schwarzschild black holes. Gravitational signatures can appear in the electromagnetic spectrum, lensing by rotating naked singularities was considered in thus leading to the possibility of directly testing differ- [35], and it was shown that the shift of the critical curves ent physical models by using astrophysical observations as a function of the lens angular momentum decreases of the emission spectra from accretion disks. slightly for the weakly naked, and vastly for the strongly naked singularities, with the increase of the scalar charge. It is the purpose of the present paper to consider another observational possibility that may distinguish The present paper is organized as follows. The ge- naked singularities from black holes, namely, the study ometry of the considered naked singularity is presented of the properties of the thin accretion disks around ro- in Section II. In Section III we obtain the main physi- tating naked singularities and black holes, respectively. cal parameters (specific energy, the specific angular mo- Thus, in the following we consider a comparative study mentum, and angular velocity) for massive test particles of the physical properties of thin accretion disks around a in stable circular orbits in stationary and axisymmetric particular type of rotating naked singularity, and rotating spacetimes. The motion of test particles around rotat- black holes, described by the Kerr metric, respectively. In ing naked singularities is considered in Section IV. The particular, we consider the basic physical parameters de- frame dragging effect is analyzed in Section V. The prop- scribing the disks, like the emitted energy flux, the tem- erties of standard thin accretion disks are reviewed in perature distribution on the surface of the disk, as well as Section VI. The energy flux, temperature distribution, the spectrum of the emitted equilibrium radiation. Due and radiation spectrum from thin disks around naked sin- to the differences in the exterior geometry, the thermody- gularities and Kerr black holes are discussed in Section namic and electromagnetic properties of the disks (energy VII. Some observational implications of our results are flux, temperature distribution and equilibrium radiation considered in Section VIII. We discuss and conclude our spectrum) are different for these two classes of compact results in Section IX. 3

II. SPACETIME GEOMETRY OF THE NAKED geodesic in the equatorial plane connecting the singular- SINGULARITY ity with the future null infinity, i.e., the singularity is naked [42]. As an example of a rotating naked singularity geometry For γ = 1 the line element describes the spacetime we will consider in the following the Kerr-like solution of geometry of the Kerr black hole, where the vanishing the Einstein gravitational field equations with a massless horizon function, ∆(r) = 0, has the roots r± = M(1 2 ± scalar field, R = 8πϕ ϕ , where ϕ is the scalar field, 1 a ). Then r+ gives the horizon radius for a∗ < 1. µν − ,µ ,ν − ∗ | | obtained in [41]. In the coordinate system (t,r,θ,φ), the pIn the case of rotating black holes, the condition gtt =0 line element, adapted to the axial symmetry, is given by does not provide any singular surface, but the ergosphere 2 2 located at re = M(1 + 1 a∗ cos θ). 2 γ 2 − γ ds = f (dt wdφ) 2w(dt wdφ)dφ We note that since forp the metric (1) gtφ = w(1 f ) − 1−γ − −1 2− 2 − 2 2 1−γ 2 2 γ − − +f Σ(∆ dr + dθ + sin θdφ ) , (1) and gφφ = f Σ sin θ + w (2 f ), respectively, the frame dragging frequency of this− rotating solution can be where written as Mr w(1 f γ) f(r, θ) = 1 2 , ω = − . (6) − γΣ f 1−γΣ sin2 θ + w2(2 f γ) − w(θ) = a sin2 θ , For the Kerr black holes (γ = 1) we obtain the familiar 2 2 2 Σ(r, θ) = r + a cos θ , expression ω = 2MarA−1, with A(r, θ) = (r2 + a2) Mr ∆a2 sin2 θ. − ∆(r) = r2 + a2 2 , − γ In the equatorial approximation ( θ π/2 1), the components of the metric (1) reduce| − in the| equatorial≪ and γ is a constant. coordinate system (t,r,z = r cos θ, φ) to the form The energy-momentum tensor of the massless scalar γ field generating the naked singularity is given by gtt = f , (7) − g = a(1 f γ) , (8) 1 tφ − − T ϕ = ϕ ϕ g gαβϕ ϕ . (2) 1−γ 2 2 γ µν ,µ ,ν − 2 µν ,α ,β gφφ = f r + a (2 f ) , (9) 2 − µν 1−γ r The scalar field satisfies the equation g ϕ = grr = f , (10) ∇µ∇ν ∆ 0, which follows from the conservation of the energy- 1−γ momentum tensor, µT ϕ = 0, and is given by gzz = f , (11) ∇ µν where 1 γ2 2Mr ϕ = − ln 1 . (3) 2M p 4  − γ (r2 + a2 cos2 θ)  f(r)=1 . − γr For 0 <γ < 1 this metric describes the spacetime The scalar field is given, in this approximation, by ϕ = geometry of a naked singularity, with a total mass M, 2 2 1 γ ln f/4. In the equatorial plane the frame drag- and an angular momentum J = aM = a∗M . Here − gingp frequency, given by Eq. (6), has the expression a∗ = J/M is the dimensionless spin parameter. It can be shown that the scalar curvature of this geometry diverges −1 −6 −1 ω =2M a∗x A , (12) at a radius where the condition g = f γ = 0 holds, tt − which indicates the existence of a curvature singularity where x = r/M, and at the radius p 1−γ 2 −4 γ f + a∗x (2 f ) −1 2 2 2 A r = γ M(1 + 1 γ a cos θ) . (4) =2 2 γ − , s − ∗ x (1 f ) p − Considering the angular dependence of this expression, respectively, with f(x)=1 2γ−1x−2. For the Kerr we see that the minimal radial coordinate of the singu- solution (γ = 1), A reduces to− the form [36, 38] −1 2 2 larity is at the pole, rs(θ =0)= γ M(1 + 1 γ a∗), − A 2 −4 2 −6 whereas it has a maximal radius in the equatorialp plane, =1+ a∗x +2a∗x . rs(θ = π/2) = 2M/γ. The bigger root of the equation ∆(r) = 0 is III. CIRCULAR GEODESIC MOTION IN −1 2 2 STATIONARY AND AXISYMMETRIC r+ = γ M(1 + 1 γ a ) . (5) − ∗ SPACETIMES p Since for any value of θ, r+ is always less than or equal to rs, the singularity is not hidden behind the event hori- Let us consider, in the coordinate system (t,r,θ,φ), zon. It can also be shown that there is at least one null an arbitrary stationary and axially symmetric geometry, 4 with line element given by Any stationary observer, moving along a world line r = constant and θ = constant with a uniform angular ve- 2 2 2 2 2 ds = gttdt +2gtφdtdφ+grrdr +gθθdθ +gφφdφ . (13) locity Ω, has a four-velocity vector uµ (∂/∂t)µ + Ω(∂/∂φ)µ, which lies inside the surface of the∝ future light The metric (13) is adapted to the symmetries of the cone. Therefore, the condition spacetime, endowed with the time- and space-like Killing vectors (∂/∂t)µ and (∂/∂φ)µ for time translations and spatial rotations, respectively. In the equatorial approxi- mation ( θ π/2 1), and the metric functions g , g , | − | ≪ tt tφ grr, gθθ and gφφ in Eq. (13) depend only on the radial co- ordinate r. Then the geodesic equations in the equatorial µ µ 2 2 [(∂/∂t) + Ω(∂/∂φ) ] = gtt +2Ωgtφ +Ω gφφ 0, (20) plane take the form ≤

dt Eg + Lg dφ Eg + Lg = φφ tφ , = tφ tt , (14) dτ g2 g g dτ −g2 g g etφ − tte φφ etφ − tte φφ and appearing in Eqs. (17) and (18), provides the constraint Ωmin < Ω < Ωmax for the angular velocity of the sta- dr 2 g = V (r), (15) tionary observers, with rr dτ  respectively, where τ is the affine parameter, E and L are the specific energy and specific angular momentum of the particles moving along the time-like geodesics,e ande the potential term V (r) is defined by 2 gtt Ωmin = ω ω , (21) − r − gφφ 2 2 E gφφ +2ELgtφ + L gtt V (r) 1+ . (16) gtt 2 Ω = ω + ω2 , (22) ≡− gtφ gttgφφ max e −e e e r − gφφ where the frame dragging frequency of the spacetime is For circular orbits in the equatorial plane, the condi- defined by ω = gtφ/gφφ. tions V (r) = 0 and V,r(r) = 0, respectively, must hold. − These conditions determine the specific energy E, the specific angular momentum L and the angular velocity e Only circular orbits for which the condition (20) holds Ω of particles moving on circular orbits as e do exist. The limiting case of this condition, gtt +2Ωgtφ + Ω2g = 0, gives r , the innermost boundary of the cir- g + g Ω φφ ph E = tt tφ , (17) cular orbits for particles, called photon orbit. The circu- − 2 gtt 2gtφΩ gφφΩ lar orbits (r > r ), for which the condition E < 1 holds, e − − − ph p gtφ + gφφΩ L = , (18) are bound, and the condition E = 1 gives the radius rmb 2 e gtt 2gtφΩ gφφΩ of the marginally bound circular orbit, i.e., the innermost e − − − e dφp g + (g )2 g g orbits. The marginally stable circular orbits rms around Ω = = − tφ,r tφ,r − tt,r φφ,r . (19) the central object are determined by the condition dt p gφφ,r

E2g +2ELg + L2g (g2 g g ) φφ,rr tφ,rr tt,rr − tφ − tt φφ ,rr V,rr r=rms = 2 =0 , (23) | e e e gtφ egttgtφ − r=rms

where the condition V,rr < 0 holds for all stable cir- IV. EQUATORIAL GEODESIC MOTION cular orbits. The marginally stable orbit is the inner- AROUND ROTATING NAKED SINGULARITIES most boundary of the stable circular orbits of the Keple- rian rotation. By inserting Eqs. (17)-(18) into Eq. (23), Inserting the metric components (7)-(9) into Eqs. (17)- and solving the resulting equation for r, we obtain the (19), we obtain the specific energy, the specific angular marginally stable orbit, once the metric coefficients gtt, momentum, and the angular velocity of the particles or- gtφ and gφφ are explicitly given. biting along circular geodesics in the equatorial plane of 5

10 10 9 9 8 8 7 7 6 6 5 5 4 4

r/M 3 r/M 3

2 2 rms rms rmb γ=1.00 rmb γ=0.99 rph rph r+ rs 1 1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

a* a*

10 10 9 9 8 8 7 7 6 6 5 5 4 4

r/M 3 r/M 3

2 2 rms rms rmb γ=0.70 rmb γ=0.50 rph rph rs rs 1 1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

a* a*

FIG. 1: Radii of circular equatorial orbits around a black hole (γ = 1.00), and a naked singularity (γ = 0.99, 0.7 and 0.5), with a total mass M, as functions of the spin parameter a∗. The dashed curves show the direct orbits, and the dotted curves represent the retrograde ones. The event horizon/singularity is denoted by a solid curve, while the horizontal dotted-dashed line at 2M represents the ergosphere in the equatorial plane. the rotating naked singularity. Hence Ω, E and L can be reduce to the familiar expressions written as e e −1 −3 −1 Ω = M x B , (24) −3 B = 1+ a∗x , E = C −1/2G , (25) −2 −3 C = 1 3x +2a∗x , C −1/2F − −2 −3 Le = Mx , (26) G = 1 2x + a∗x , − −3 2 −4 F = 1 2a∗x + a x , where e − ∗

1−γ −3 B = f h + a∗x , 1−γ −3 −2 obtained for the equatorial approximation of the Kerr C = f [(fh +2a∗x )h x ], −3 − solution [36, 38]. G = fh + a∗x , F 1−γ γ −1 2 −4 = f [1 a∗ (1 f ) x h]+ a∗x , If we compare the locations of the marginally stable, of − − the marginally bound, and of the photon orbits, plotted 2 −1 1/2 with h(x) = [1 (1 γ)(2 γx ) ] . For γ = 1 we in Fig. 1 as functions of the spin parameter a∗, we see have f(x)=1 −2x−2−and h−= 1, and the latter formulae that they have a strong dependence on the parameter γ. −

For static black holes (a∗ = 0 and γ = 1 on the top left hand plot), the orbits with radii rms, rmb and rph, are 6

8 8 a*=0.00 a*=0.00 6 a*=0.40 6 a*=0.40 a*=0.70 a*=0.70 4 a*=0.90 4 a*=0.90 2 a*=0.99 2 a*=0.99 γ γ 2 =1.00 2 =0.99

0 0 V / d (r/M) V / d (r/M) 2 2

d -2 d -2 -2 -2

10 -4 10 -4

-6 -6

-8 -8 2 3 4 5 6 7 8 9 10 2 3 4 5 6 7 8 9 10 r/M r/M

8 8 a*=0.00 a*=0.00 6 a*=0.40 6 a*=0.40 a*=0.70 a*=0.70 4 a*=0.90 4 a*=0.90 2 a*=0.99 2 a*=0.99 γ γ 2 =0.70 2 =0.50

0 0 V / d (r/M) V / d (r/M) 2 2

d -2 d -2 -2 -2

10 -4 10 -4

-6 -6

-8 -8 2 3 4 5 6 7 8 9 10 2 3 4 5 6 7 8 9 10 r/M r/M

2 2 FIG. 2: The function d V/d(r/M) versus r/M for a rotating black hole, with (γ = 1), and a naked singularity (γ = 0.99, 0.7, and 0.5), with a total mass M. The spin parameter a∗ is set to 0.99, 0.9, 0.7 and 0, respectively. located at 6M,4M and 3M, respectively. As the spin in- γ = 0.99 (see the top right hand plot). After reaching creases from zero they pass through the ergosphere, and this critical value of the rotation, the marginally stable, approach the even horizon r+ at the maximal spin value marginally bound and photon orbits ”jump” to the singu- a∗ = 1. For γ < 1 (see the other three plots in Fig. 1), larity, located at rs =2M/γ. If the rotation exceeds the there is a critical value a∗c of the spin, where for each type value a∗c, corresponding to the marginally stable orbit, of direct orbit the curves have a cut-off. The behavior then all the radii rms, rmb, rph and rs become degener- of the retrograde orbits exhibits only a slight dependence ated. With increasing values of γ, these critical values on γ. The value of a∗c is minimal for the radius rph, and are decreasing and even in the case of slow rotation we maximal for r , e.g., it is a∗ 0.6 and a∗ 0.85 for obtain at the singularity fully degenerated orbits. ms c ∼ c ∼

As an illustration of the dependence of rms on the spin is already approaching the static limit at r = 2M, and a∗ and on the scalar charge parameter γ, we present in for a∗ 1, the marginally stable orbit enters the ergo- Fig. 2 the second order derivative of the effective po- sphere,→ and it approaches the event horizon. The other tential with respect to the radial coordinate, given by panels in Fig. 2, presenting V,rr for the naked singular- Eq. (23), as a function of the radius r. Here the values ity (γ < 1), show that Eq. (23) has no longer solutions of γ are the same as in Fig. (1), and the spin parameter for higher values of the spin parameter, and V,rr remains a∗ takes different values. The top left hand panel shows negative everywhere. As a result, the particles rotat- the case of the black hole (γ = 1). For a static black ing around the singularity have stable circular orbits in hole V,rr vanishes at 6M, and its zero shifts to lower the whole spacetime, up to its boundary, i.e., to the sin- radii with increasing spin parameter. For a∗ 0.9, r gularity itself. With decreasing γ, the critical value of ∼ ms 7 a∗, above which V,rr has no zeros any more, is also de- 8.5 . a*=0.40 creasing. For γ 0.5, only very slowly rotating naked a =0.50 8 * singularities have a marginally stable orbit, which does a*=0.60 a =0.70 not coincide with the singularity. 7.5 * a*=0.80

] a*=0.90 -1 7 a =0.99 s * V. THE FRAME DRAGGING EFFECT 4

[10 6.5 ω An interesting question related to the geodesic motion 6 in the gravitational potential of the rotating naked sin- gularity is the frame dragging effect, i.e., the relation 5.5 between the rotational velocity ω of the spacetime itself, 5 and the azimuthal component of the four-velocity of the 10-8 10-7 10-6 10-5 10-4 γ freely falling (or rotating) massive test particles. The (r/M-rs/M) frame dragging frequency ω at the singularity can be cal- culated from Eq. (12) as FIG. 3: The radial profile of ω near the rotating naked sin- a∗ gularity with γ = 0.9, for different spin parameters. ω(rs) = lim , (27) 3 1−γ 2 x→√2/γ 2M 2/γ f +2Ma∗ p outer boundary of the region between r and r , where where x = rs/M = 2/γ. From a physical point s i the frame dragging effect of the spacetime with the spin of view ω(rs)p can be interpretedp as the angular velocity S parameter a is weaker than the frame dragging effect of the singularity, Ω = ω(rs). We note that the limit 1∗ lim f 1−γ vanishes for γ < 1, but it becomes unity measured for a lower spin value a2∗. Outside this re- x→√2/γ gion, i.e., for all radii x greater than x , the relation for γ = 1. The latter result is not trivial, but if it holds, i ω(x, a ) > ω(x, a ) holds. Fig. 3, where we plotted we reobtain ω(r ) for the Kerr solution, since r r for ∗1 ∗2 e s → e ω as a function of the radial distance for naked singular- γ = 1. For the rotating naked singularity metric given ities with same mass and same scalar charge (γ = 0.9), by Eq. (1), the expression (27) reduces to but different spin parameters, shows this behavior. In 1 Fig. 3, the solid black line, representing the frame drag- ΩS = ω(r )= , s 2Ma ging frequency obtained for the naked singularity with ∗ the spin value of 0.99, has the highest values at radii −4 which is clearly valid only for a rotating singularity, since r/M greater than rs/M + 10 /γ. For lower radii, a S ∼ Ω would diverge in the static case. Even if a∗ = 0 radius ri, where the ω curves for lower spin values inter- holds, the explicit expression for ΩS might be puzzling,6 sect the black solid curve, always exists. The lower the since it does not depend on γ, and shows that the an- spin value is, the closer ri shifts to rs. For radii less than gular velocity of a rotating naked singularity is inversely ri, the frame dragging frequency for a∗ < 0.99 becomes proportional to its spin parameter. We would expect higher than the value obtained for a∗ =0.99. This result that the angular frequency of any rotating object to be holds for any pair of spin values, i.e., there exist radii ri, proportional to the angular momentum of the body - where for any value of a∗1, the curve ω, and the curve and this is indeed true almost over the whole space- for a∗1 >a∗2, intersect each other. In all cases we obtain time, except in a very small domain, close to the naked the relation ω(x, a1∗) <ω(x, a2∗) for x< ri/M. singularity. In a very close vicinity of the naked sin- A similar pathological behavior can bep found for the gularity, at rs, ω becomes bigger for lower values of angular velocity Ω of massive test particles. For the first a∗, as compared to the values corresponding to higher term of B appearing in Eq. (24), as we approach the spins. We can determine the radius ri at which the singularity we obtain the limit proportionality relation ω(x, a∗1) > ω(x, a∗2) for a∗1 > 1/2 a∗2 is inverted by inserting Eq. (12) into the equation 1 γ lim f (1−γ)h = lim f 1−γ 1 − ω(x ,a∗ ) = ω(x ,a∗ ), with x = r /M. The result 2 i 1 i 2 i i x→√2/γ x→√2/γ  − γx f  is a∗1A (xi,a∗2) = a∗2A (xi,a∗1), whichp can be further 1 simplified to the equation = γ(1 γ) lim f 1/2−γ . (29) 2 − x→ 2/γ 4 1−γ γ p √ x f = a∗ a∗ (2 f ) . (28) i 1 2 − For the angular velocity we obtain For γ = 1 this equation reduces to the cubic equation 4 2 x a∗ a∗ x a∗ a∗ = 0 , which has the triple real i − 1 2 i − 1 2 0 , 1/2 <γ< 1 , root xi = 0, showing that, as expected, there is no in- −1  −1 version in the ω a∗ relation for Kerr black holes. For lim Ω= M (2 + a∗) , γ =1/2 , . (30) − x→√2/γ  −1 γ < 1, the biggest real root of Eq. (28) provides the a∗ , 0 <γ< 1/2  8

By supposing that the naked singularity rotates fast the equatorial plane. enough, so that rph = rs, then, if 1/2 <γ < 1, for the orbits located close to the singularity Ω decreases to zero. For γ = 1/2, the angular velocity of the ro- In Fig. 4 we have plotted Ω, Ωmin, Ωmax and ω for tating particles has a non-vanishing value at the sin- the same values of a∗ and γ as we have used in Fig. 2. gularity, and it remains finite even in the static case, The top left hand panels show the frame dragging effect 1/2−γ in the Kerr spacetime. When ω = 0 (the static case), (limx→ 2/γ f = 1 is supposed again). For γ < 1/2, √ for r > rph = 3M the orbiting particles can have a Ke- −1 Ω is proportional to a , and for a∗ 0 it diverges ∗ → plerian rotation, with the frequency Ω. At lower radii to infinity, as ω also does. Since at the singularity there are no circular orbits, and the particles fall freely Ω (r )=2ω(r ) coincides with Ω, the difference Ω ω max s s − onto the event horizon, with a velocity with azimuthal also diverges at rs for a vanishing spin parameter. How- component restricted to values between Ωmin and Ωmax. ever, this divergence does not involve any violation of For the static black hole at the even horizon these values causality, since the velocity of the orbiting particles, mea- tend to ω = 0, i.e., there is no frame dragging. With in- sured in the locally non rotating frame, does not exceed creasing a∗, the frame dragging effect becomes more and the speed of light, i.e., v 1 still holds, where | |≤ more important: as the orbital radius is approaching r+, the particle is dragged with the rotating spacetime. No gφφ v = Ω ω , static observer exists in the ergosphere, since any particle | | 2 | − | g gttgφφ passing through the static limit (g = 0) is forced to ro- tφ − tt q tate with a positive angular velocity. For rapid rotations for circular equatorial orbits. At rs we have gtt = 0, and a∗ 1, rph is also approaching the event horizon, where v = Ω/ω 1 = 1, a relation that provides the upper the→ particles move along circular orbits with the angular | | | − | H limit of the velocity for time-like particles that satisfy the velocity of the horizon, namely, Ω = ω(r+). At r+ we H condition (20) when moving along circular geodesics in obtain Ω = Ω =Ω = ω =Ω , as a∗ 1. min max →

These effects are considerably different for rotating radius has jumped to the singularity (the top right hand naked singularities, presented in the rest of the panels panel in Fig. 1). Particles can therefore move along circu- in Fig. 3. At r = rs, the lapse function for the metric lar orbits in the whole spacetime, although not all these vanishes, gtt = 0, as it does in the static limit for the orbits are stable. However, at the singularity, located at black holes, and Eqs. (21) and (22) give Ω = 0 and 2M, Ω vanishes, even if ω has a finite, non-zero value. min ∼ Ωmax =2ω =1/Ma∗, without fixing the value of Ω at the For γ =1/2, Ω has already a non-zero value at the singu- singularity. The result that at the singularity, depending larity, which is inversely proportional to 2 + a∗, but it is on the values of γ, Ω is either zero, or takes some non- still less than ω, as shown in the middle right hand panel zero values between 0 and 2ω, follows from Eq. (30). For in Fig 4. The panels in the bottom of Fig. 4 show that relatively small scalar charges, γ = 0.99, the static case both Ω and ω increase very rapidly close to the singular- is still similar to the static case of the black holes: the ity - in fact Ω exceeds the frame dragging frequency at marginally stable orbit is located at around 6M (the top r as Ω(r ) 1/Ma∗, and the relation ω(r ) 1/2Ma∗ s s → s → right hand panel in Fig. 2), and for the orbits with radii holds for γ < 1/2. The relation Ω > ω for r rs is greater than r 3M, the Keplerian rotation can still shown in Fig. 5, where we have plotted the ratio→ Ω/ω as ph ∼ be maintained. At lower radii Ω is restricted to values a function of the radius. As mentioned before, for a∗ 0 → between Ωmin and Ωmax, and Ω vanishes at the singu- this ratio tends to 1/2 at the singularity, the particles or- larity (as it does at the event horizon of a Schwarzschild biting at rs have finite velocities, with the speed of the black hole). For rotating singularities with γ =0.99 and light like particles tending towards the speed of light. a∗ = 0.7, rms is still greater than 4M, but the photon

The angular velocity of particles moving along sta- zero from its maximal value at rmax, which is still close ble circular geodesics around black holes is maximal at to the singularity (r r M). By inserting Eq. (24) max − s ≪ the marginally stable orbit, and decreases monotonically, into the condition Ω,r = 0 we obtain the equation with the increasing coordinate radius r. Eq. (30), and Fig. 4, show that this is not the case for the particles rotating in the gravitational potential of rotating naked singularities with γ > 1/2. At the singularity Ω drops to 3γ2x4 6γ(γ + 1)x2 +4γ2 +6γ +2=0, (31) − 9

2 γ a*=0.99 =1.00 1 a*=0.99 γ=0.99 -1 -1

s 1 s 4 4

10 0 10 0

2 1 a*=0.70 a*=0.70 -1 -1

s 1 s 4 4

10 10 0 0 4 Ω 1 3 a*=0.00 Ω a*=0.00 -1 Ω min -1 s 2 s

4 max 4 ω 0

10 1 10 0 1 2 3 4 5 6 7 8 9 10 2 3 4 5 6 7 8 9 10 r/M r/M

1 a*=0.99 γ=0.70 1 a*=0.99 γ=0.50 -1 -1 s s 4 0 4 0 10 10

1 a*=0.70 1 a*=0.70 -1 -1 s s 4 0 4 0 10 10

1 a*=0.00 1 a*=0.00 -1 -1 s s 4 0 4 0 10 10

3 4 5 6 7 8 9 10 4 5 6 7 8 9 10 r/M r/M 1 1 a*=0.99 γ=0.40 a*=0.99 γ=0.20 -1 -1 s s 4 0 4 0 10 10

1 1 a*=0.70 a*=0.70 -1 -1 s s 4 0 4 0 10 10

1 1 a*=0.00 a*=0.00 -1 -1 s s 4 0 4 0 10 10

2 3 4 5 6 7 8 9 10 2 3 4 5 6 7 8 9 10 r/M r/M

FIG. 4: The radial profiles of Ω, Ωmin, Ωmax and ω for black holes with γ = 1 and for the rotating naked singularity with γ = 0.99, 0.9, 0.7, 0.5, 0.4, and 0.2, respectively. The spin parameter a∗ is set to 0.99, 0.7 and 0, respectively. with the solution the derivative of Ω is not zero, indicating a monotonic decrease in the angular velocity with increasing radial 2 2 1+ γ (1 γ )/3 coordinate. For γ = 1 Eq. (32) has two roots, where xmax,± = ± − . (32) 2 6 pγ xmax,+ = rmax/M is the value of the radial coordinate where Ω is maximal. Since Eq. (32) does not depend on For the limiting case of the Kerr black holes (γ = 1), the spin parameter a∗, the dimensionless radius xmax,+ 2 2 Eq. (32) has only one root, xmax,+ = xmax,− = 2, is only determined by the parameter γ of the singular- or rmax = 2M, which is due to multiplication of Ω,r ity. The location of the marginally stable orbit rms still by f(r) in the computation of Eq. (31). At r = 2M, 10

2 2 a*=0.40 a*=0.40 a =0.50 a =0.50 * 1.8 * 1.8 a*=0.60 a*=0.60 a*=0.70 a*=0.70 a =0.80 1.6 a =0.80 1.6 * * a*=0.90 a*=0.90 a*=0.99 1.4 a*=0.99 ω ω / 1.4 / Ω Ω 1.2 1.2 1

1 0.8

0.8 0.6 0 0.2 0.4 0.6 0.8 1 0 1 2 3 4 5 -10γ -3γ 10 (r/M-rs/M) 10 (r/M-rs/M)

FIG. 5: The radial profiles of the ratio Ω/ω for γ = 0.4 (left hand panel) and γ = 0.2 (right hand panel), for different values of the spin parameter.

depends on a∗, as in the case of Kerr black holes, and it size (defined along the z-axis) is negligible, as compared decreases as the singularity rotates faster and faster. At to its horizontal extension (defined along the radial di- orbits with radii higher than rmax, the angular velocity is rection r), i.e, the disk height H, equal to the maximum increasing with decreasing r, but if rmax is greater than half thickness of the disk, is always much smaller than r , then there is an annulus between the inner edge of the characteristic radius R of the disk, H R. The thin ms ≪ the accretion disk and rmax, where for particles moving disk is in hydrodynamical equilibrium, and the pressure along circular orbits Ω is decreasing for small distances gradient and a vertical entropy gradient in the accret- from the singularity. ing matter are negligible. In the steady-state accretion disk models, the mass accretion rate M˙ 0 is supposed to be constant in time, and the physical quantities of the VI. STANDARD ACCRETION DISKS accreting matter are averaged over a characteristic time scale, e.g. ∆t, and over the azimuthal angle ∆φ = 2π, for a total period of the orbits and for the height H. The Accretion discs are flattened astronomical objects plasma moves in Keplerian orbits around the compact made of rapidly rotating gas which slowly spirals onto object, with a rotational velocity Ω, and the plasma par- a central gravitating body, with its gravitational energy ticles have a specific energy E, and specific angular mo- degraded to heat. A fraction of the heat converts into mentum L, which depend only on the radii of the orbits. e radiation, which partially escapes, and cools down the The particles are orbiting with the four-velocity uµ in a accretion disc. The only information that we have about disk havinge an averaged surface density Σ. The accreting accretion disk physics comes from this radiation, when it matter is modeled by an anisotropic fluid source, where reaches radio, optical and X-ray telescopes, allowing as- the rest mass density ρ0 (the specific internal energy is tronomers to analyze its electromagnetic spectrum, and neglected), the energy flow vector qµ and the stress tensor its time variability. The efficient cooling via the radiation tµν are measured in the averaged rest-frame. The energy- over the disk surface prevents the disk from cumulating momentum tensor describing this source takes the form the heat generated by stresses and dynamical friction. µν µ ν (µ ν) µν In turn, this equilibrium causes the disk to stabilize its T = ρ0u u +2u q + t , (33) thin vertical size. The thin disk has an inner edge at µ µν the marginally stable orbit of the compact object poten- where uµq = 0, uµt = 0. The four-vectors of the energy and of the angular momentum flux are defined by tial, and the accreting plasma has a Keplerian motion in µ µ ν µ µ ν higher orbits. E T ν (∂/∂t) and J T ν (∂/∂φ) , respectively. −The four≡ dimensional conservation≡ laws of the rest mass, For the general relativistic case the theory of mass ac- of the energy and of the angular momentum of the plasma cretion around rotating black holes was developed by provide the structure equations of the thin disk. From Novikov and Thorne [36]. They extended the steady- the structure equations the flux of the radiant energy state thin disk models introduced in [37] to the case of over the disk can be expressed as [38, 39] the curved space-times, by adopting the equatorial ap- r proximation for the stationary and axisymmetric geom- M˙ 0 Ω,r F (r)= (E ΩL)L,rdr , (34) etry. For a steady-state thin accretion disk described in 2 −4π√ g (E ΩL) Zrms − the cylindrical coordinate system (r, φ, z = cos θ) most of − − e e e the matter lies close to the radial plane. Hence its vertical where the no-torquee inner boundarye conditions were pre- 11 scribed [38]. This means that the torque vanishes at the accretion disks. The efficiency is defined by the ratio of inner edge of the disk, since the matter at the marginally two rates measured at infinity: the rate of the radiation of stable orbit rms falls freely into the black hole, and can- the energy of the photons escaping from the disk surface not exert considerable torque on the disk. The latter to infinity, and the rate at which mass-energy is trans- assumption is valid as long as strong magnetic fields do ported to the compact object. If all the emitted photons not exist in the plunging region, where matter falls into can escape to infinity, the efficiency depends only on the the hole. specific energy measured at the marginally stable orbit The geometry of the space-time near to the equator, rms, or the metric potential determines the radial dependence of Ω, E and L for the particles moving on circular or- ǫ =1 E . (37) bits around the central object. We can therefore calcu- − r=rms late thee averagede radial distribution of photon emission e For Schwarzschild black holes the efficiency is about 6%, for accretion disks around the rotating singularity in the no matter if we consider the photon capture by the black equatorial approximation, by applying the flux integral hole, or not. Ignoring the capture of radiation by the Eq. (34). Evaluating of the specific energy at the inner black hole, ǫ is found to be 42% for rapidly rotating black edge of the disk, we can also determine the efficiency of conversion of the rest mass into outgoing radiation. holes, whereas the efficiency is 40% with photon capture in the Kerr potential. The accreting matter in the steady-state thin disk model is supposed to be in thermodynamical equilibrium. Then the radiation emitted by the disk surface can be considered as a perfect black body radiation, where the VII. ELECTROMAGNETIC SIGNATURES OF 4 ACCRETION DISKS AROUND ROTATING energy flux is given by F (r) = σSB T (r) (σSB is the Stefan-Boltzmann constant) and the observed luminos- NAKED SINGULARITIES ity L (ν) has a redshifted black body spectrum [43]: After analyzing, in Sections IV and V, the circular 8 rf 2π ν3rdφdr geodesic motion around a rotating naked singularity, we L (ν)=4πd2I (ν)= cos i e . 2 are now ready to discuss the properties of the disk radia- πc Zri Z0 exp(νe/T ) 1 −(35) tion for standard accretion disk models in the spacetime Here d is the distance to the source, I(ν) is the Planck of the naked singularity. In Fig. 6 we present the flux distribution function, i is the disk inclination angle (we profile, calculated from Eq. (34), for the physical param- eters of the configurations already shown in Fig. 2. In set it to zero), and ri and rf indicate the position of the inner and outer edge of the disk, respectively. We take the following we set the total mass to 5M⊙, and the ac- cretion rate to 10−12M /yr. The top left hand panel ri = rms and rf , since we expect the flux over ⊙ the disk surface vanishes→ ∞ at r for any kind of gen- shows F (r) for the static and the rotating black holes. eral relativistic compact object→ geometry. ∞ The emitted The dependence of the flux distribution on the spin has distinct features: the inner edge of the disk is located frequency is given by νe = ν(1 + z), where the redshift factor can be written as at rms = 6M for the static black hole, and shifts to lower radii, approaching M, as the black hole spins up 1+Ωr sin φ sin γ 1+ z = , (36) to a∗ = 1. The radii of the marginally stable orbits are g 2Ωg Ω2g determined by the zeros of V,rr in Fig. 2. With increas- − tt − tφ − φφ p ing spin, the maximal flux is also increasing with at least where we have neglected the light bending [44, 45]. three orders of magnitude, as compared to the cases with The efficiency ǫ with which the central object converts a∗ = 0 and a∗ =0.99. For higher values of the spin, the rest mass into outgoing radiation is the other important locations of the maxima of the spectra also shift to lower physical parameter characterizing the properties of the radii, located closer to the inner edge of the disk.

The rest of the panels in Fig. 6 shows the flux dis- panel in Fig. 2 shows that for a∗ = 0, 0.4 and 0.7 there tribution of the thermal radiation of the accretion disks exists a radius where the second order derivative of the ef- around naked singularities. Depending on the values of fective potential vanishes, whereas for a∗ =0.9 and 0.99 the scalar charge parameter γ, and of the spin parameter V,rr < 0 in the whole spacetime. Hence for γ = 0.99, a∗, naked singularities and black holes could exhibit ei- F (r) has similar characteristics to the flux profiles de- ther similar, or rather different properties. For γ =0.99, rived in the case of the black holes for a∗ = 0, 0.4 and we obtain flux profiles similar to the corresponding cases 0.7., but the flux distribution around the black hole and of the black holes, with the same spin, provided that the the naked singularity is rather different for higher spin equation V,rr = 0 has a real solution. The top right hand values. 12

101 102 a*=0.00 a*=0.00 a*=0.40 a*=0.40 a*=0.70 1 a*=0.70

] ] 10 0 a*=0.90 a*=0.90 -2 10 -2 a*=0.99 a*=0.99

cm γ=1.00 cm γ=0.99 -1 -1 100 -1

erg s 10 erg s 20 20 10-1

-2

F(r) [10 10 F(r) [10 10-2

10-3 10-3 1 2 3 4 5 6 8 10 15 20 30 40 60 1 2 3 4 5 6 8 10 15 20 30 40 60 r/M r/M

102 102 a*=0.00 a*=0.00 a*=0.40 a*=0.40 1 a*=0.70 1 a*=0.70

] 10 ] 10 a*=0.90 a*=0.90 -2 -2 a*=0.99 a*=0.99

cm γ=0.70 cm γ=0.50 -1 100 -1 100 erg s erg s 20 10-1 20 10-1

F(r) [10 10-2 F(r) [10 10-2

10-3 10-3 2 3 4 5 6 8 10 15 20 30 40 60 3 4 5 6 8 10 15 20 30 40 60 r/M r/M

104 104 a*=0.00 a*=0.00 a*=0.40 a*=0.40 3 a*=0.70 3 a*=0.70

] 10 ] 10 a*=0.90 a*=0.90 -2 -2 a*=0.99 a*=0.99

cm γ=0.40 cm γ=0.20 -1 102 -1 102 erg s erg s 20 101 20 101

F(r) [10 100 F(r) [10 100

10-1 10-1 4 5 6 7 8 9 10 15 9 10 11 12 13 14 15 r/M r/M

FIG. 6: The energy flux F (r) radiated by the disk around a rotating black hole (γ = 1), and a naked singularity (γ = 0.99, 0.7, 0.5, 0.4 and 0.2), with the total mass 5M⊙, for different values of the spin parameter a∗. The mass accretion rate is set to −12 10 M⊙/yr.

In the slowly rotating case, when the inner edge of the larity. Since stable circular orbits around the naked sin- disk is located at the marginally stable orbit, the radial gularity with a∗ =0.9 and 0.99 exist in the whole equa- profiles of F (r) have the same shape in both the top torial plane, the inner edge of the accretion disk reaches left and the right hand side plots, and the flux maxima the singularity. The radius rs =2M/γ is bigger than rms are somewhat smaller for the black holes as compared to for the rotating black holes with the spin of 0.9 or 0.99, flux maxima for the naked singularity. This difference is and therefore the surface of the accretion disk around enhanced by the fast rotation of the central naked singu- the rotating naked singularity is smaller than the one for 13 its Kerr black hole counterpart. Even this smaller disk Thus, F (r) is much more sensitive for the variation in a∗ surface is not the whole area that can radiate thermal in the case of the naked singularity. The maximal flux is photons in the thin accretion disk model, since the left somewhat higher for a∗ = 0.9 than for a∗ = 0.99, which end of the flux profile is pushed to a radius somewhat indicates that the maximum of F (r) for fast rotation is higher than rs. This is due to the fact that Ω,r changes inversely proportional to the spin parameter. its sign at rmax(> rs), as seen in Fig. 4, and this also changes the sign of the cofactor of the integral in the flux The fact that the flux maximum is higher for the ro- formulae, given by Eq. (34). As a result, we obtain tating naked singularity than for the black hole, even if negative flux values for the radii between rs and rmax, it is integrated over a smaller surface area, is the conse- (since the sign of the integral remains the same at rmax). quence of the considerable difference in the metric deter- The negative flux values involve non-physical states in minant, characterizing the four-volume element in which the framework of the stationary thin disk model. It in- the radiant flux is measured in the vicinity of the equa- dicates that the thermodynamical equilibrium cannot be torial plane. For Kerr black holes √ g = r2 holds in the maintained in this region and other forms of the energy equatorial approximation, but from the− metric (7)-(11) of and momentum transport become dominant over the ra- the rotating naked singularity we obtain √ g = f 1−γr2. − diative cooling such as advection or convection. Then Since the shape factor f vanishes as r rs, the func- 1−γ → the approximation of the steady-state and geometrically tion f has a small value at rmax, which is close to thin disk model breaks down here. Only outside this re- rs. Then the four-volume element is much smaller for gion, which forms a very thin annulus between rs and the naked singularity than for the Kerr black hole, and rmax though, one can apply the standard accretion disk it produces much higher values in the flux integral (34), scheme. We assume the thin disk exists only in the region even if the disk properties determining Ω, E˜ and L˜ are for r < rmax, where the radiation gives the main contri- similar in the two cases, as far as we integrate from rmax bution for energy and momentum transport, and other in the case of the naked singularity (e.g., in the two plots forms of the transport processes becomes negligible. in the top in Fig. 4, the values of Ω have only a moderate Here we consider the radiation properties only for the difference for the naked singularity and the black hole). geometrically thin, relatively cold disk, which is trun- cated at rmax, and we discuss the contribution of the hot We obtain similar trends in the characteristics of the matter in the region r < rmax to the total radiation of radiated flux if we further decrease the values of γ. By the disk-naked singularity system in the next Section. considering the plots for γ = 0.7 or 0.5 in Fig. 2, one Thus we assume that the accretion disk emits ther- can see that the equation V,rr = 0 has a real solution mal photons only at radii higher than rmax, and that for only for a∗ = 0 and a∗ = 0.4 (γ = 0.7), or only for the r > rmax the photon flux can be computed in the thin static case (γ =0.5). For higher values of a∗, the second disk approximation. Since rmax does depend only on γ, order derivative of the effective potential remains nega- but not of a∗, the left edge of F (r) is located always at the tive in the whole equatorial plane, and the inner edge same radii, irrespectively of the rotational speed of the of the disk always jumps to rs. This effect is shown by naked singularity. In spite of this reduction in the effec- the two middle plots in Fig. 6, where the flux profiles tive radiant area of the disk surface, as compared with for the naked singularity can be separated into the two the Kerr black hole case, with the same spin parame- groups (similarly to the case of γ =0.99), formed by the ters (0.9 and 0.99), the flux maxima is much higher for curves similar to those obtained for black holes (V,rr =0 the rotating naked singularity as compared to the black has a real solution), and the curves with higher maxima, hole case. The two plots on the top of Fig. 6 show that located almost at their left edge (V,rr < 0 holds every- the rotating naked singularity with a∗ =0.99 has a disk where). The high flux values obtained for the second ten times more luminous than the disk of a Kerr black group are again the consequence of the rapidly shrinking holes with the same spin. Comparing the rapidly rotat- volume element in the vicinity of the singularity. We also ing naked singularity with its static case we also see a rise see that the flux maxima for the latter group is inversely of 3 orders of magnitude in the flux maximum, whereas proportional to the spin parameter, but the curves fall this rise is only 2 orders of magnitude for the black hole. more rapidly for lower spin.

Since Ω approaches zero at the singularity for 1/2 < the singularity, where the volume element shrinks to zero, γ < 1, the inner boundary of the radiant disk area is and the cofactor of ( g)−1/2 in the flux formula (34) goes − located at rmax, where we obtain finite flux values for to infinity. Their multiplication for γ =1/2 causes F (R) a∗ 0. For γ =1/2, Ω has a non zero value at rs, and to remain finite at the singularity, as seen in Fig. 6. For it is→ a monotonous function in the entire equatorial plane. all non-static cases, plotted in the middle right hand plot Therefore, we cannot use rmax as a truncation radius for in Fig. 6, the peaks at the left edge of the flux profiles the thin disk, and the inner edge of the disk would reach have infinite amplitudes, which cannot represent physical 14

a*=0.00 a*=0.00 a*=0.40 a*=0.40 a*=0.70 a*=0.70 a*=0.90 a*=0.90 a*=0.99 a*=0.99 γ γ 1 =1.00 1 =0.99 K] K] 6 6 T(r) [10 T(r) [10

0.1 0.1 1 2 3 4 5 6 8 10 15 20 30 40 60 1 2 3 4 5 6 8 10 15 20 30 40 60 r/M r/M

a*=0.00 a*=0.00 a*=0.40 a*=0.40 a*=0.70 a*=0.70 a*=0.90 a*=0.90 a*=0.99 a*=0.99 γ γ 1 =0.70 1 =0.50 K] K] 6 6 T(r) [10 T(r) [10

0.1 0.1 2 3 4 5 6 8 10 15 20 30 40 60 3 4 5 6 8 10 15 20 30 40 60 r/M r/M

10 103 a*=0.00 a*=0.00 a*=0.40 a*=0.40 a*=0.70 a*=0.70 a*=0.90 a*=0.90 a*=0.99 a*=0.99 γ=0.40 γ=0.20 K] K] 6 6 T(r) [10 T(r) [10

1 102

4 5 6 7 8 9 10 15 9 10 11 12 13 14 15 r/M r/M

FIG. 7: The temperature distribution on a disk around a rotating black hole (γ = 1), and a naked singularity (γ = 0.99, 0.7, 0.5, 0.4 and 0.2), with the total mass of 5M⊙, for different values of the spin parameter a∗. The mass accretion rate is set to −12 10 M⊙/yr. states. The thin diks model is not a good approximation and angular momentum transport. at the inner edge of the disk, as the steady state of the disk cannot be maintained close to the singularity. We For γ < 1/2, Ω a−2 holds at r , and Ω will diverge as ∼ ∗ s can still state that the accretion disk must be extremely a∗ approaches zero. Its derivative with respect to r also bright in its innermost area, but there must be some other tends to infinity, also giving, together with the vanishing physical mechanisms besides the radiative cooling, such volume element at rs, an infinite flux value in Eq. (34). as advection, which have an important role in the energy This effect is demonstrated by the bottom plots in Fig. 6. Since Ω and its derivative is inversely proportional to a∗, 15 the rate of the rise of the flux profiles as approaching values, but these maxima are still much higher for rotat- rs is also inversely proportional to a∗. Thus, the slowly ing naked singularities as compared with those calculated rotating naked singularities have brighter accretion disks for disks rotating around Kerr black holes. than their fast rotating counterparts do have, but in both In Fig. 7 we present the temperature distribution of cases F (r) tends to infinity at the inner edge of the disk. the accretion disk for the same configurations which we Without any upper boundary in the flux maxima at the have used to study the flux profiles. The disk tempera- singularity, the physical state of the matter in motion ture exhibits a similar dependence on the parameters γ close to the singularity will recede from the thermody- and a∗ as F (r) does. With decreasing γ, and increasing namical equilibrium. The standard accretion disk model a∗, we obtain temperature profiles with much higher and cannot be applied at such small distances from the singu- sharper maxima than those for the Kerr black hole. For larity. At higher radii there must exists a region, where γ 1/2, there still exist configurations with lower spin, the disk matter can attain thermodynamical equilibrium. which≥ give radial temperature profiles similar to those ob- If this radius is still close enough to rs, the flux profiles, tained for the Kerr black holes. Although these tempera- shown in the bottom plots of Fig. 6, still resemble the ture profiles become uncertain in the innermost region of real physical situation, with some uncertainties in the the accretion disk, the disk must still be extremely hot in position of their left edges. Obviously, in a physical situ- this region, as compared to the typical disk temperatures ation, the amplitudes of the peaks must have some finite obtained for Kerr black holes, with the same spin values.

In Fig. 8 we have also plotted, for the same set of values value, and then it starts to decrease again. The rate of of the parameters a∗ and γ, the disk spectra, which were the increase and of the decrease depends on the location calculated from the luminosity formula Eq. (35). These of the inner edge of the disk, as can be seen in Eq. (37). plots show the same trends found in the behavior of the The behavior of the two groups, with the inner disk edge radiated flux distribution F (r), and of the disk tempera- at rms and rs, respectively, is the same in this respect, ture in the black hole and naked singularity geometries, but ǫ varies in different ranges for the two groups. For respectively. For the Kerr black hole, with γ = 1, the cut- the configurations with inner disk edge located at rms, off frequency of the spectra shifts towards higher values, the increase in a∗ causes only a moderate variation in and the maximal amplitudes increase with the increas- ǫ. While the efficiency for the slowly rotating case is ing spin parameter, i.e, the accretion disks of Kerr black somewhat greater for the naked singularity, the rapidly holes become hotter by rotating faster, and they produce rotating black holes have a considerable higher efficiency a blueshifted surface radiation with higher intensity. In than the naked singularities, with the same spin parame- the case of the rotating naked singularity, the two classes ter, do have. Comparing the first two lines of the table we already identified in the previous discussion of the radial see that the static and the slowly rotating configurations distribution of F (r) and of T (r) preserve their different have the same efficiencies, whereas the accreted mass-to- natures. We have seen that there is only a moderate radiation conversion mechanism is about four times more decrease in the disk temperature for the static solution, efficient for extremely rapidly rotating black holes then and for rotating naked singularities with γ = 0.99 and for the fast spinning naked singularities, with γ = 0.99. with spin of a∗ = 0.4 and 0.7, as compared to the Kerr We find relatively small values in the first two columns black hole with the same spin values. As a result, the of the third line, and in the first column of the fourth disk spectra exhibit only negligible differences for black line as well. All the values in the last four lines belong holes and naked singularities. For fast rotation (a∗ =0.9 to the second group of naked singularities, for which the and 0.99), the accretion disk is much hotter in the area inner disk edge is located at the singularity (indicated close to the singularity. Hence, the cut-off frequencies with the same values in the parenthesis in each line as shift blueward, and the maximal amplitudes of the spec- rs, and which does not depend on a∗). In these cases the tra are much higher then for the black hole disk spectra. accretion disk radiates a great amount of thermal energy, The spectral features for the second group are not sensi- and the values of the efficiency are considerably higher tive to the variations in the spin, and the relative shifts than those found for the first group. For the static cases in the cut-off frequencies and the spectral maxima are the efficiency can reach 50% for γ =0.4, but it drops to very small, as shown in the bottom right hand plot in 40%, as we decrease γ to 0.2. With increasing spin, ǫ has Fig. (8). a very mild increase, reaching the values 70% (γ = 0.5) and 61% (γ = 0.4) at a∗ 0.4, and then it starts to In Table I we present the conversion efficiency ǫ of the decrease for faster rotating∼ singularities. For the naked accreted mass into radiation for both rotating naked sin- singularity with the spin of 0.99 and γ = 0.2, the effi- gularities and black holes, for different values of the pa- ciency falls to 34%, which is still higher than the effi- rameters a∗ and γ. This Table demonstrates that with ciency for extremely fast rotating Kerr black holes. We increasing spin parameter, and decreasing γ (increasing can conclude that there is a range of the spin parameter scalar charge), the efficiency is increasing to a maximal 16

1035 1035 γ=1.00 γ=0.99 1034 1034

] 33 ] 33

-1 10 -1 10

1032 1032 ) [erg s ) [erg s ν ν

L( 31 L( 31

ν 10 ν 10 a*=0.00 a*=0.00 a =0.40 a =0.40 30 * 30 * 10 a*=0.70 10 a*=0.70 a*=0.90 a*=0.90 a*=0.99 a*=0.99 1029 1029 1015 1016 1017 1018 1015 1016 1017 1018 ν [Hz] ν [Hz]

1035 1035 γ=0.70 γ=0.50 1034 1034

] 33 ] 33

-1 10 -1 10

1032 1032 ) [erg s ) [erg s ν ν

L( 31 L( 31

ν 10 ν 10 a*=0.00 a*=0.00 a =0.40 a =0.40 30 * 30 * 10 a*=0.70 10 a*=0.70 a*=0.90 a*=0.90 a*=0.99 a*=0.99 1029 1029 1015 1016 1017 1018 1015 1016 1017 1018 ν [Hz] ν [Hz]

1035 1035 γ=0.40 γ=0.20 1034 1034

] 33 ] 33

-1 10 -1 10

1032 1032 ) [erg s ) [erg s ν ν

L( 31 L( 31

ν 10 ν 10 a*=0.00 a*=0.00 a =0.40 a =0.40 30 * 30 * 10 a*=0.70 10 a*=0.70 a*=0.90 a*=0.90 a*=0.99 a*=0.99 1029 1029 1015 1016 1017 1018 1015 1016 1017 1018 ν [Hz] ν [Hz]

FIG. 8: The accretion disk spectra for a rotating black hole (γ = 1), and a naked singularity (γ = 0.99, 0.7, 0.5, 0.4 and 0.2), −12 with the total mass of 5M⊙, for different values of the spin parameter a∗. The mass accretion rate is set to 10 M⊙/yr.

a∗, and of the scalar charge γ, where Kerr black holes parameters, where the conversion efficiency for rotating represent more effective engines in the conversion of ac- naked singularities is much higher than for Kerr black creted mass to radiation than the naked singularities do. holes. Nevertheless, we have also found another range of these

Some composite accretion disk models consider a geo- between the marginally stable orbit and the inner edge metrically thick and optically thin hot corona, positioned 17

γ a∗ = 0 0.4 0.7 0.9 0.99 1.00 5.72% (6.00) 7.51% (4.62) 10.4% (3.40) 15.6% (2.32) 26.4% (1.46) 0.99 5.72% (6.00) 7.51% (4.62) 10.4% (3.37) 6.03% (2.02) 5.92% (2.02) 0.70 6.03% (6.15) 8.92% (4.14) 73.4% (2.86) 70.8% (2.86) 69.7% (2.86) 0.50 6.94% (6.00) 70.6% (4.00) 65.1% (4.00) 61.8% (4.00) 60.4% (4.00) 0.40 50.7% (5.01) 61.3% (5.01) 56.4% (5.01) 53.2% (5.01) 51.8% (5.01) 0.20 40.8% (10.0) 41.6% (10.0) 38.0% (10.0) 35.3% (10.0) 34.3% (10.0)

TABLE I: The efficiency, measured in percents, and x2 = r/M (values in parenthesis) at the inner edge of the accretion disk for Kerr black holes, and rotating naked singularities, respectively. The value of γ = 1 corresponds to the Kerr black hole.

of the geometrically thin, and optically thick, accretion Eddington luminosity becomes a coordinate dependent disk, where an inner edge is lying at a few gravitational quantity, LEdd(r) M(r). A similar effect occurs radii ([50]. In other type of composite models, the corona for the case of the∝ naked singularity solution consid- lies above and under the accretion disk, and the soft pho- ered in the present paper. For simplicity in the follow- tons, arriving from the disk, produce a hard emission via ing we will consider only the case of the static naked- their inverse comptonization by the thermally hot elec- singularity. One can associate to the scalar field de- trons in the corona [51]. In both configurations, the disk scribed by this energy-momentum tensor a mass distribu- is truncated at several gravitational radii - reducing the tion M(r) along the equatorial plane of the disk, given by soft photon flux of the disk - and the soft and hard com- M(r)=4π r T ϕ0r2dr =2π r grrϕ ϕ r2dr. By using rs 0 rs ,r ,r ponents of the broad band X-ray spectra of galactic black the explicitR form of the scalarR field we obtain holes were attributed to the thermal radiation of the ac- π M 2 1 γ2 r ∆ cretion disk and the emission mechanism in the corona, M(r)= − f γ−3dr. (38) 2
4 respectively. 2 γ Zrs r Although the metric potential gives stable circular The corresponding coordinate-dependent Eddington lu- geodesic orbits in the entire equatorial plane of the space- minosity can be obtained as time, it cannot guarantee that the matter in motion can maintain configurations with thermodynamical equilib- 1.37 1033 M 2 1 γ2 L (r) = × × − rium in the whole region, which is the basic assumption of Edd γ2
× the geometrically thin steady state accretion disk model. r ∆ A natural limit of the region where the standard accretion f γ−3dr erg/s. (39) Z r4 disk model applies seems to be the radius rmax, where the rs Keplerian rotational frequency becomes inversely propor- The variation of the Eddington luminosity with respect tional to the radial distance from the naked singular- to the coordinate r is represented, for a naked singularity ity. This inversion in the Ω versus r relation indicates with M =5M⊙, in Fig. 9. that the structure equations of the Novikov-Thorne disk It was argued in [47] that any neutron star, composed model have non-physical solutions in the region below by matter described by a more or less general equation rmax. Since this limiting radius depends only on the pa- of state, should experience thermonuclear type I bursts rameter γ and the geometrical mass M via Eq. (32), the at appropriate mass accretion rates. The question asked lowest radial limit of the validity of the thin disk approx- in [47] is whether an ”abnormal” surface may allow such imation is not a free parameter, but it is determined by a behavior. The naked singularities may also have such the charge parameter of the massless scalar field ϕ. a zero velocity, particle trapping, abnormal surface. The presence of a material surface located at the singularity implies that energy can be radiated, once matter col- VIII. OBSERVATIONAL IMPLICATIONS lides with that surface. One can also consider composite accretion disk models as an alternative solution to this An interesting effect, involving the Eddington lumi- problem, and set the truncation radius of the thin disk nosity for the case of a boson star, was pointed out in at rmax. The hot corona of the composite model can lie [43]. The Eddington luminosity, a limiting luminosity in the region between rs and rmax, producing hard X- that can be obtained from the equality of the gravita- ray spectra via the inverse Compton scattering of pho- tional force inwards and of the radiation force outwards, tons radiated from the disk and the electron gas in the 38 is given by L = 4πMm /σ = 1.3 10 (M/M⊙) corona. This hot corona could represent the ”surface” of Edd p T × erg/s, where mp is the proton’s mass, and σT is the the naked singularity. Thus, at least in principle, naked Thompson cross section [46]. Since a boson star- as singularity models, characterized by high mass, normal well as any other transparent object - has a non-constant matter crusts/surfaces and type I thermonuclear bursts mass distribution, with M = M(r), and therefore the can be theoretically constructed. 18

6 12 1020 M = 4 10 M⊙ and radius R = 1.4 10 cm, for γ=0.99 × × γ γ = 1 we obtain a surface redshift of zgrav = 1.55, for 1018 =0.70 γ=0.50 γ = 0.99 we obtain zgrav = 1.6057, while for γ = 0.85, γ 1016 =0.40 zgrav = 9.99. Therefore the radiation coming from the γ=0.20 14 surface of a naked singularity may be highly redshifted Edd 10 (in standard general relativity the redshift obeys the con- 1012 strain z 2).

(r) / L ≤ 10 Type I X -ray bursts show strong spectroscopic ev- Edd 10 L idence for rapid expansion of the radius of the X-ray 108 photosphere. The luminosities of these bursts reach the 6 Eddington critical luminosity at which the outward radi- 10 ation force balances gravity, causing the expansion lay- 104 ers of the star. The touchdown luminosity of radius- 2 3 4 5 6 8 10 15 20 30 40 60 expansion bursts from a given source remain constant r/M between bursts to within a few percent, giving empirical verification to the theoretical expectation that the emerg- FIG. 9: Coordinate dependent Eddington luminosity as a ing luminosity is approximately equal to the Eddington function of the radial coordinate r for the case of a static critical luminosity. The Eddington luminosity at infinity naked singularity with total mass M = 5M⊙, and for differ- of a naked singularity with a thin surface is given by [52] ent values of γ. 2 ∞ 4πmprs √ gtt d LEdd = − √ gtt (40) σT  √grr dr −  r=rs Since naked singularities could be surrounded by a thin shell of matter, the presence of a turning point for matter For the metric given by Eqs. (7)-(11) we obtain (the point where the motion of the infalling matter sud- denly stops) at the surface of the naked singularity may 4πm M ∆ L∞ = p f 3(γ−1)/2 = have important astrophysical and observational implica- Edd σ  r T i r=rs tions. Since the velocity of the matter at the naked sin- M ∆ gularity surface is zero, matter can be captured and de- 1.3 1038 f 3(γ−1)/2 . (41) × M⊙  r  posited on the surface of the naked singularity. Moreover, r=rs because matter is accreted continuously, the increase in For r 2M/γ, L∞ . This shows that the the size and density of the surface will ignite some ther- → Edd → ∞ monuclear reactions [47]. The ignited reactions are usu- luminosity of naked singularities with thin surfaces can ally unstable, causing the accreted layer of gas to burn reach much higher values than in the case of standard explosively within a very short period of time. After the astrophysical objects. Finally, we consider the apparent nuclear fuel is consumed, the naked singularity also re- surface area during burst cooling. Observations of the verts to its accretion phase, until the next thermonuclear cooling tails of multiple type I bursts from a single source instability is triggered. Thus, once a thin material sur- have shown that the apparent surface area of the emit- ∞ 2 4 face is formed, naked singularities may undergo a semi- ting region, defined as S =4πD Fc,∞/σSBTc,∞, where regular series of explosions, called type I thermonuclear Fc,∞ is the measured flux of the source during the cooling bursts, discovered first for X-ray binaries [48, 49]. tail of the burst, Tc,∞ is the measured color temperature of the burst spectrum, D is the distance to the source and The observational signatures indicating the presence of σ is the Stefan-Boltzmann constant, remains approx- X-ray bursts from naked singularities are similar to those SB imately constant during each burst, and between bursts of standard neutron stars, and are the gravitational red- from the same source. The color temperature on the sur- shift of a surface atomic line, the touchdown luminosity of face of the compact object Tc,h is related to the color tem- a radius-expansion burst, and the apparent surface area tt during the cooling phases of the burst [52]. perature measured at infinity by Tc,h = Tc,∞ g [52]. By introducing the color correction factor fc =pTc/Teff , If the thermal radiation with wavelength λe emitted by where Teff is the effective temperature at the surface, we the thin shell of matter at the surface of the naked singu- obtain larity has absorption or emission features characteristic of atomic transitions, these features will be detected at in- r2 S∞ =4π [z(r) + 1]2 . (42) finity at a wavelength λo, gravitationally redshifted with f 4  −γ/2 c r=rs a value zgrav = (λo λe) /λe = (1 2M/γr) 1, where we have assumed,− for simplicity,− that the radia-− Since the radius of the naked singularity as well as its tion emission takes place in the equatorial plane of the redshift may be very large quantities, the apparent area naked singularity. The value of the redshift for a neu- of the emitting region as measured at infinity may be 6 tron star with mass M = 2M⊙ and radius R = 10 cm also very large. Hence all the astrophysical quantities is zNS =0.56. By assuming a naked singularity of mass related to the observable properties of the X-ray bursts, 19 originating at the surface of the naked singularity can be disk surface is neglected in the energy and momentum calculated, and have finite values on the surface of the transport between the different regions of the disk. Once naked singularity and at infinity. any of these conditions is violated, the thin disk model cannot be applied. On the other hand in all the covari- ant general relativistic formulations of disk models the IX. DISCUSSIONS AND FINAL REMARKS physical quantities are obtained upon integrations over the four volume element. In our case, the behavior of the volume element near the singular state gives the domi- In the present paper we have considered the properties nant contribution to the flux, temperature, and spectrum of the accretion disks that could form around naked sin- of the disk, and this contribution is much larger than the gularities, hypothetical singular general relativistic theo- possible effect on the physical parameters from some im- retical objects, characterized by the absence of an event proved disk models. horizon. As for the naked singularity we have considered a rotating solution of the Einstein - massless scalar field It is generally expected that most of the astrophysical equations, which reduces to the Kerr solution when the objects grow substantially in mass via accretion. Re- scalar charge tends to 1. As a first step in our study cent observations suggest that around most of the active we have investigated the motion of the particles in the galactic nuclei (AGN’s) or black hole candidates there gravitational potential of this solution. Depending on exist gas clouds surrounding the central far object, and the values of the mass, scalar charge, and of the spin pa- an associated accretion disk, on a variety of scales from rameter, respectively, there are two types of disks that a tenth of a parsec to a few hundred parsecs [53]. These could exist around naked singularities. For the first type clouds are assumed to form a geometrically and optically there are marginally stable orbits, located outside the thick torus (or warped disk), which absorbs most of the naked singularity, while for the second type the particles ultraviolet radiation and the soft x-rays. The gas exists in can reach, and be in direct contact, with the singularity. either the molecular or the atomic phase. Evidence for Consequently, the properties of the disk radiation are sig- the existence of super massive black holes comes from nificantly different for these two types of disks. While the very long baseline interferometry (VLBI) imaging of the first type shows similarities with the Kerr black hole molecular H2O masers in active galaxies, like NGC 4258 disk, the thermodynamic/electromagnetic properties of [54], and from the astrometric and radial velocity mea- the second type could differ by several orders of mag- surements of the fully unconstrained Keplerian orbits for nitude from the Kerr disks. A very puzzling result is short period stars around the supermassive black hole at represented by the behavior of the angular velocity of the center of our galaxy [55, 56]. The VLBI imaging, the particles at the singularity, which is inversely pro- produced by Doppler shift measurements assuming Kep- portional to the spin parameter of the naked singularity. lerian motion of the masering source, has allowed a quite The frame dragging properties of the naked singularity accurate estimation of the central mass, which has been 7 found to be a 3.6 10 M⊙ super massive dark object, also show significant differences as compared to the Kerr × black hole case. within 0.13 parsecs. In our approach we have used the thin disk model, Hence, important astrophysical information can be ob- which is an idealized physical model. The thin disk model tained from the observation of the motion of the gas is derived under several physical and geometrical assump- streams in the gravitational field of compact objects. tions [38]. The basic assumptions for the geometrically Therefore the study of the accretion processes by com- thin and optically thick model are as follows: a) The pact objects is a powerful indicator of their physical na- space-time geometry is given by the metric of the ro- ture. However, up to now, the observational results have tating naked singularity, and the disk with a negligible confirmed the predictions of general relativity mainly in self-gravity resides in its central plane. b) The disk is a qualitative way. With the present observational preci- geometrically thin. c) There exist a time interval ∆t, sion, one cannot distinguish between the different classes small enough for neglecting any change in the parame- of compact/exotic objects that appear in the theoreti- ters of the naked singularity during ∆t, but large enough cal framework of general relativity [47]. However, im- for measuring the total inward mass flows at any r in the portant technological developments may allow to image disk, as compared with the total disk mass between r and black holes and other compact objects directly [57]. Re- 2r. d) The stress energy tensor of the disk plasma can be cent observations at a wavelength of 1.3 mm have set a algebraically decomposed with respect to the 4-velocity size of microarcseconds on the intrinsic diameter of SgrA* of the plasma. e) The averaged motion of the baryons [58]. This is less than the expected apparent size of the over the azimuthal angle and ∆t can be described as a event horizon of the presumed black hole, thus suggest- circular geodesic motion in the equatorial plane. f) The ing that the bulk of SgrA* emission may not be centered heat flow within the disk is negligible, except in the verti- on the black hole, but arises in the surrounding accre- cal direction. g) The averaged stress-energy propagating tion flow. A model in which Sgr A* is a compact object to the disk surface is carried by thermal photons, and the with a thermally emitting surface was considered in [59]. photons are emitted from the disk surface vertically on Given the very low quiescent luminosity of Sgr A* in the average. h) The role of the photons emitted from the the near-infrared, the existence of a hard surface, even 20 in the limit in which the radius approaches the hori- thermodynamic and electromagnetic properties of the zon, places a severe constraint on the steady mass ac- disks (energy flux, temperature distribution and equi- −12 cretion rate onto the source: M˙ 10 M⊙/yr. This librium radiation spectrum) are different for these two limit is well below the minimum≤ accretion rate needed classes of compact objects, consequently giving clear ob- to power the observed submillimeter luminosity of Sgr servational signatures that could help to identify them −10 A*: M˙ > 10 M⊙/yr. Thus it follows that Sgr A* observationally. More specifically, comparing the en- does not have a surface, i.e., that it must have an event ergy flux emerging from the surface of the thin accretion horizon. This argument could be made more restrictive disk around black holes and naked singularities of similar by an order of magnitude with microarcsecond resolution masses, we have found that for some values of the spin imaging, e.g., with submillimeter very long baseline inter- parameter and of the scalar charge its maximal value is ferometry. Submilliarcsecond astrometry and imaging of much higher for naked singularities. In fact all the ther- the black hole Sgr A* at the Galactic Centre may become modynamical properties of the disks strongly depend on possible in the near future at infrared and submillimetre the values of the spin parameter and on the scalar charge wavelengths [60]. The expected images and light curves, parameter. These effects are confirmed from the analysis including polarization, associated with a compact emis- of the disk temperatures and disk spectra. In addition sion region orbiting the central black hole were computed to this, we have also shown that for a given range of in [61]. From spot images and light curves of the observed the spin parameter and of the scalar charge, the con- flux and polarization it is possible to extract the black version efficiency of the accreting mass into radiation of hole mass and spin. At radio wavelengths, disc opacity naked singularities is generally much larger than the con- produces significant departures from the infrared behav- version efficiency for black holes, i.e., naked singularities ior, but there are still generic signatures of the black hole provide a much more efficient mechanism for convert- properties. Detailed comparison of these results with fu- ing mass into radiation than black holes. Thus, these ture data can be used to test general relativity, and to observational signatures may provide the possibility of improve existing models for the accretion flow in the im- clearly distinguishing rotating naked singularities from mediate vicinity of the black hole. Kerr-type black holes. With the improvement of the imaging observational Acknowledgments techniques, it will also be possible to provide clear obser- vational evidence for the existence of naked singularities, and to differentiate them from other types of compact The work of T. H. was supported by the General Re- general relativistic objects. search Fund grant number HKU 701808P of the govern- Indeed, in the present paper we have shown that the ment of the Hong Kong Special Administrative Region.

[1] R. Penrose, Nuovo Cimento 1, 252 (1969). (1998). [2] P. S. Joshi, Global Aspects in Gravitation and Cosmol- [18] M. K. Parikh and F. Wilczek, Phys. Lett. B449, 24 ogy, Clarendon, Oxford, (1993). (1999). [3] P. C. Vaidya, Proc. Indian Acad. Sci. A33, 264 (1951). [19] B. C. Nolan, Phys. Rev. D60, 024014 (1999). [4] R. C. Tolman, Proc. Natl. Acad. Sci. USA 20, 169 (1934). [20] A. Ishibashi and A. Hosoya, Phys. Rev. D60, 104028 [5] J. P. S. Lemos, Phys. Rev.Lett. 68, 1447 (1992). (1999). [6] W. B. Bonnor, P. C. Vaidya, Gen. Rel. Grav. 1, 127, [21] V. Husain, Phys. Rev. D53, R1759 (1996). (1970); F. de Felice, Astron. Astrophys. 34, 15 (1974); [22] A. Wang, Y. Wu, Gen. Rel. Grav. 31, 107, (1999). F. de Felice, Astron. Astrophys. 45, 65 (1975). [23] J. F. Villas da Rocha, Int. J. Mod. Phys. D11, 113 [7] K. Lake, T. Zannias, Phys. Rev. D43, 1798 (1991). (2002). [8] B. T. Sullivan and W. Israel, Phys. Lett. A79, 371, [24] E. N. Glass and J. P. Krisch, Phys. Rev. D57, R5945 (1980). (1998). [9] E. Poisson and W. Israel, Phys. Rev. D41, 1796 (1990). [25] E. N. Glass and J. P. Krisch, Class. Quantum. Grav. 16, [10] A. Ori, Phys. Rev. D57, 4745 (1998). 1175 (1999). [11] T. P. Singh and P. S. Joshi, Class. Quantum. Grav. 13, [26] K. S. Govinder and M. Govender, Phys. Rev. D68, 559 (1996). 024034 (2003) [12] S. Jhingan, P. S. Joshi and T. P. Singh, Class. Quantum. [27] T. Harko and K. S. Cheng, Phys. Lett. A266, 249, Grav. 13, 3057 (1996). (2000). [13] P. S. Joshi and A. Krolak, Class. Quantum. Grav. 13, [28] S. G. Ghosh and N. Dadhich, Phys. Rev. D65, 127502, 3069 (1996). (2002). [14] J. H. Dwivedi and P. S. Joshi, Class. Quantum. Grav. [29] S. G. Ghosh and N. Dadhich, Gen. Rel. Grav. 35, 359, 14, 1223 (1997). (2003). [15] T. P. Singh, Phys. Rev. D58, 024004 (1998). [30] T. Harko, Phys. Rev. D68, 064005 (2003). [16] S. Barve and T. P. Singh, Mod. Phys. Lett. A12, 2415 [31] K.S. Virbhadra, D. Narasimha, and S. M. Chitre, Astron. (1997). Astrophys. 337, 1 (1998). [17] L. M. Gonzalez-Romero, Class. Quantum. Grav. 15, L59 [32] K. S. Virbhadra and G. F. R. Ellis, Phys. Rev. D65, 21

103004 (2002). Gravity 26, 055009 (2009). [33] K. S. Virbhadra, Int. J. Mod. Phys. A12, 4831 (1997). [47] Y. F. Yuan, R. Narayan and M. J. Rees Astrophys. J. [34] K. S. Virbhadra and C. R. Keeton, Phys. Rev. D77, 606, 1112 (2004). 124014 (2008). [48] J. Grindlay, H. Gursky, H. Schnopper, D. R. Parsignault, [35] G. N. Gyulchev and S. S. Yazadjiev, Phys. Rev. D78, J. Heise, A. C. Brinkman, and J. Schrijver, Astrophys. 083004 (2008). J. 205, L127 (1976). [36] I. D. Novikov and K. S. Thorne, in Black Holes, ed. C. [49] D. Tournear, E. Raffauf , E. D. Bloom, W. Focke, B. DeWitt and B. DeWitt, New York: Gordon and Breach Giebels, G. Godfrey, P. M. Saz Parkinson, K. T. Reilly, (1973). K. S. Wood, P. S. Ray, M. T. Wolff, R. M. Bandyopad- [37] N. I. Shakura and R. A. Sunyaev, Astron. Astrophys. 24, hyay, M. N. Lovellette, and J. D. Scargle, Astrophys. J. 33 (1973). 595, 1058 (2003). [38] D. N. Page and K. S. Thorne, Astrophys. J. 191, 499 [50] K. S. Thorne and R. H. Price, Astrophys. J. 195, L101 (1974). (1975). [39] K. S. Thorne, Astrophys. J. 191, 507 (1974). [51] E. T. P. Liang and R. H. Price, Astrophys. J. 218, 247 [40] S. Bhattacharyya, A. V. Thampan and I. Bombaci, As- (1977). tron. Astrophys. 372, 925 (2001); F. S. Guzman, Phys. [52] D. Psaltis, Phys. Rev. D77, 064006 (2008). Rev. D 73, 021501 (2006); T. Harko, Z. Kovacs and [53] C. M. Urry and P. Padovani, Publ. Astron. Soc. of the F. S. N. Lobo, Phys. Rev. D D78, 084005 (2008); Pacific 107, 803 (1995). T. Harko, Z. Kov´acs and F. S. N. Lobo, Phys. Rev. D79, [54] M. Miyoshi, J. Moran, J. Herrnstein, L. Greenhill, N. 064001 (2009); C. S. J. Pun, Z. Kov´acs and T. Harko, Nakai, P. Diamond and M. Inoue, Nature 373, 127 Phys. Rev.D78, 084015 (2008); Z. Kov´acs, K. S. Cheng (1995). and T. Harko, Astron. Astrophys., 500, 621 (2009); [55] A. M. Ghez, S. Salim, N. N. Weinberg, J. R. Lu, T. Do, C. S. J. Pun, Z. Kov´acs and T. Harko, Phys. Rev. D78, J. K. Dunn, K. Matthews, M. Morris, S. Yelda, E. E. 024043 (2008); Z. Kovacs, K. S. Cheng, and T. Harko, Becklin, T. Kremenek, M. Milosavljevic, and J. Naiman, MNRAS 400, 1632 (2009); T. Harko, Z. Kovacs, and F. Astrophys. J. 689, 1044 (2008). S. N. Lobo, Phys. Rev. D80, 044021 (2009); C. Bambi, [56] S. Gillessen, F. Eisenhauer, S. Trippe, T. Alexander, R. K. Freese, T. Harada, R. Takahashi, N. Yoshida, Phys. Genzel, F. Martins, and T. Ott, Astrophys. J. 692, 1075 Rev. D80, 104023, (2009); T. Harko, Z. Kovacs, and F. (2009). S. N. Lobo, Class. Quant. Grav. 27, 105010 (2010); T. [57] H. Falcke, F. Melia, and E. Agol, Astrophys. J. 528, L13 Harko, Z. Kovacs, and F. S.N. Lobo, arXiv:1009.1958 (2000). (2010); C. Bambi, T. Harada, R. Takahashi, N. Yoshida, [58] S. Doeleman, J. Weintroub, A. E.E. Rogers, R. Plam- Phys. Rev. D81, 104004 (2010); P. Pani, E. Barausse, E. beck, R. Freund, R. P.J. Tilanus, P. Friberg, L. M. Ziurys, Berti, and V. Cardoso, Phys. Rev D82, 044009 (2010); J. M. Moran, B. Corey, K. H. Young, D. L. Smythe, M. M. Heydari-Fard, Class. Quantum Grav. 27, 235004 Titus, D. P. Marrone, R. J. Cappallo, D. C.J. Bock, G. C. (2010). Bower, R. Chamberlin, G. R. Davis, T. P. Krichbaum, J. [41] K. D. Krori and D. R. Bhattacharjee, J. Math. Phys. 23, Lamb, H. Maness, A. E. Niell, A. Roy, P. Strittmatter, 637 (1982). D. Werthimer, A. R. Whitney, and D. Woody, Nature [42] K. S. Virbhadra, S. Jhingan and P. S. Joshi, Int. J. Mod. 455, 78 (2008). Phys. D6, 357 (1997); K. S. Virbhadra, Phys. Rev. D60, [59] A. E. Broderick and R. Narayan, Astrophys. JL. 636, 104041 (1999). L109 (2006). [43] D. Torres, Nucl. Phys. B 626, 377 (2002). [60] A. E. Broderick and A. Loeb, Mon. Not. Roy. Astron. [44] J. P. Luminet, Astron. Astrophys. 75, 228 (1979). Soc. 367, 905 (2006). [45] S. Bhattacharyya, R. Misra and A. V. Thampan, Astro- [61] A. E. Broderick and A. Loeb, Mon. Not. Roy. Astron. phys. J. 550, 841 (2001). Soc. 363, 353 (2005). [46] D. Bini, R. T. Jantzen, and L. Stella, Class. Quantum