Accretion disk around the rotating Damour-Solodukhin wormhole

A new rotating generalization of the Damour-Solodukhin wormhole (RDSWH), called Kerr-like wormhole, has recently been proposed and investigated by Bueno \textit{et al} for echoes in the gravitational wave signal. We show a novel feature of the RDSWH, viz., that the kinematic properties such as the ISCO\ or marginally stable radius $r_{\text{ms}}$, efficiency $\epsilon$ and the disk potential $V_{\text{eff}}$ are \textit{independent} of $\lambda$ (which means they are identical to their KBH counterparts for any given spin). Differences however appear in the emissivity properties for higher values $0.1<\lambda\leq 1$ (say) and for the extreme spin $a_{\star}=0.998$. The kinematic and emissivity are generic properties as variations of the wormhole mass and the rate of accretion within the model preserve these properties. Specifically, the behavior of the luminosity peak is quite opposite to each other for the two objects, which could be useful from the viewpoint of observations. Apart from this, an estimate of the difference $\Delta_{\lambda}$ in the maxima of flux of radiation $F(r)$ shows non-zero values but is too tiny to be observable at present for $\lambda<10^{-3}$ permitted by the strong lensing bound. The broad conclusion is that RDSWH\ are experimentally indistinguishable from KBH by accretion characteristics.


I. INTRODUCTION
There is a renewed interest in wormholes after it has been realized that they can mimic post-merger ringdown initial quasi-normal mode (QNM) spectrum of gravitational waves. To our knowledge, this possibility of "wormhole QNM mode" was first put forward by Damour and Solodukhin [1], where resonances were trapped in a double-hump potential associated with what they termed "black hole foil". This is an example of a horizonless wormhole more commonly known as the Damour-Solodukhin wormhole (DSWH). The DSWH differs from the Schwarzschild black hole (SBH) by a dimensionless real deviation parameter λ and represents a twice-asymptotically flat regular spacetime connected by a throat. The authors showed that if λ is so small that the time scale ∆t = 2GM ln 1 λ 2 is longer than the observational time scale, the signals emitted by a source falling into a wormhole will contain the usual QNM ringing signature of a black hole, in spite of the absence of a true horizon [1]. Völkel and Kokkotas [2] have recently shown, using an inverse method, that the knowledge of the observed QNM spectrum can allow one to also accurately construct the double hump Pöschl-Teller potential approximating that of DSWH. Bueno et al. [3] suitably redefined the static DSWH and further generalized it the into a Kerr-like wormhole, called here the rotating DSWH or RDSWH for brevity. They studied its gravitational wave echo properties. Gravitational deflection of relativistic massive particles by RDSWH has been recently studied by Jusufi et al. [4].
QNM ringing by BHs have been well reviewed, see, e.g. Berti et al. [5] but ringing by WHs still remained an open question. A first step towards answering the question was initiated by Cardoso, Franzin and Pani [6], who * Corresponding author: izmailov.ramil@gmail.com studied the QNM ringing using a wormhole assembled by means of Visser's cut-and-paste surgery of two copies of Schwarzschild black holes (SBH) at a radius close to the horizon [7]. (For future perspectives and new directions of research on ultracompact objects (UCO) including wormholes, see [8].) The authors of [7] showed that, while the time evolution of the early QNMs accurately mimic those from SBH horizon, the differences (if any) would appear only at later times. This work inspired an investigation in [9], where it has been shown that the massless Ellis-Bronnikov wormhole (EBWH) [10,11], made of the minimally coupled exotic scalar field, can also reproduce the black hole QNM spectrum in the eikonal limit (large ). Ringing by massive EBWH was studied in [12]. These developments prompt a natural inquiry as to whether static DSWH can exhibit mimicking of SBH in other phenomenon as well. A recent work [13] shows that it indeed can, e.g., it can very accurately mimic SBH strong field lensing properties for λ ≤ 10 −3 . However, this is only an upper bound, which obviously does not rule out values of λ arbitrarily close to zero. Strong field lensing is an excellent diagnostic for studying the signatures of WHs versus BHs, see, e.g., [14].
There is another important diagnostic, namely, the accretion phenomenon around compact objects that has already been a very active field of research (see, e.g., [15][16][17][18][19][20][21][22][23][24][25]). The first comprehensive study of accretion disks using a Newtonian approach was made in [26]. Later a general relativistic model of thin accretion disk was developed in three seminal papers by Novikov and Thorne [27], Page and Thorne [28] and Thorne [29] under the assumption that the disk is in a steady-state, that is, the mass accretion rateṀ is constant in time and does not depend of the radius of the disk. The disk is further supposed to be in hydrodynamic and thermodynamic equilibrium, which ensure a black body electromagnetic spectrum and properties of emitted radiation. The thin accretion disk model further assumes that individual particles are mov-ing on Keplerian orbits, but for this to be true the central object is assumed to have weak magnetic field, otherwise the orbits in the inner edge of the disk will be deformed.
In this paper, we shall study the kinematic as well as emissivity properties such as the luminosity spectra, flux of radiation, temperature profile, efficiency of a thin accretion disk around a stellar sized RDSWH using the Page-Thorne model [26][27][28][29]. As a numerical example, we shall assume the accretion to take place only in the attractive positive mass mouth with mass 15M and an accretion rateṀ 0 ∼ 10 19 gm.sec −1 . The motive is to analyze how far the accretion profiles are influenced by λ and the dimensionless spin parameter a * .
The paper is organized as follows: In Sec.2, we outline the main formulas relating to the thin accretion disk to be used in the paper. The kinematic accretion formulas for generic spinning spacetime are presented in Sec.3. These are applied to RDSWH in Sec.4. Numerical estimates are made in Sec.5, while Sec.6 concludes the paper. We take units such that G = c = 1 unless restored and metric signature (−, +, +, +).

II. THIN ACCRETION DISK
We assume geometrically thin accretion disk, which means that the disk height H above the equator is much smaller than the characteristic radius R of the disk, H R. The disk is assumed to be in hydrodynamical equilibrium stabilizing its vertical size, with the pressure and vertical entropy gradient being negligible. An efficient cooling mechanism via heat loss by radiation over the disk surface is assumed to be functioning in the disk, which prevents the disk from collecting the heat generated by stresses and dynamical friction. The thin disk has an inner edge defined by the marginally stable circular radius r ms , while the orbits at higher radii are Keplerian. In the steady-state approximation, the mass accretion rateṀ 0 is assumed to be a constant and the physical quantities describing the accreting matter are averaged over a characteristic time scale ∆t, over the azimuthal angle ∆φ = 2π for a total period of the orbits, and over the height H [27,28].
The orbiting particles with the four-velocity u µ form a disk of an averaged surface density Σ, where the rest mass density ρ 0 , the energy flow vector q µ and the stress tensor t µν are measured in the averaged rest-frame. Then where ρ 0 rest mass density averaged over ∆t and 2π and the torque density with the component t r φ averaged over ∆t and 2π. The time and orbital average of the energy flow vector q µ gives the radiation flux F(r) over the disk surface as For the stress energy tensor T µ ν of the disk, the energy and angular momentum four-vectors are defined by −E µ ≡ T µ ν (∂/∂t) ν and J µ ≡ T µ ν (∂/∂φ) ν respectively. The structure equations of the thin disk can be derived by integrating the conservation laws of the rest mass, of the energy, and of the angular momentum [27,28]. From the rest mass conservation, ∇ µ (ρ 0 u µ ) = 0, it follows that the average rate of the accretion is independent of the disk radius,Ṁ In the steady-state approximation, specific energy E and specific angular momentum L of accreting particles have depend only on the radius of the orbits. Defining black hole rotational velocity Ω = dφ/dt, the energy conservation law ∇ µ E µ = 0 yields the integral This is a balance equation, which states that the energy transported by the rest mass flow,Ṁ 0 E, and the energy transported by the torque in the disk, 2πrΩW r φ , is balanced by the energy radiated away from the surface of the disk, 4πrF E.
The angular momentum conservation law, ∇ µ J µ = 0, states the balance of three forms of angular momentum transport, viz., By eliminating W r φ from Eqs. (5) and (6), and applying the energy-angular momentum relation for circular geodesic orbits in the form E ,r = Ω L ,r , the flux F of the radiant energy, or power, over the disk can be expressed as [27,28], The disk is supposed to be in thermodynamical equilibrium, as explained, so the radiation flux emitted by the disk surface will follow Stefan-Boltzmann law: where σ is the Stefan-Boltzmann constant. The observed luminosity L (ν) has a redshifted black body spectrum [16] where i is the disk inclination angle to the vertical, d is the distance between the observer and the center of the disk, r in and r out are the inner and outer radii of the disc, h is the Planck constant, ν e is the emission frequency, I(ν) is the Planck distribution, and k B is the Boltzmann constant. The observed photons are redshifted and received frequency ν is related to the emitted ones by ν e = (1 + z)ν. The redshift factor (1 + z) has the form: where the light bending effect is neglected [57,58]. Another important characteristic of the thin accretion disk is its efficiency , which quantifies the ability with which the central body converts the accreting mass into radiation. The efficiency is measured at infinity and it is defined as the ratio of two rates: the rate of energy of the photons emitted from the disk surface and the rate with which the mass-energy is transported to the central body. If all photons reach infinity, the Page-Thorne accretion efficiency is given by the specific energy of the accreting particles measured at the marginally stable orbit [28]: As the definition indicates, should be non-negative.

III. GENERIC ROTATING SPACETIME
The accretion disk is formed by particles moving in circular orbits around a compact object, with the geodesics determined by the space-time geometry around the object, be it a WH, BH or NS. For a rotating geometry the metric is generically given by ds 2 = −g tt dt 2 + 2g tφ dtdφ + g rr dr 2 + g θθ dθ 2 + g φφ dφ 2 .
(12) At and around the equator, i.e., when |θ − π/2| 1,we assume, with Harko et al [41], that the metric functions g tt , g tφ , g rr , g θθ and g φφ depend only on the radial coordinate r. The angular velocity Ω, of the specific energy E, and of the specific angular momentum L of accreting particles in the above geometry are given by g rr dr dτ One may hence define an effective potential term defined as Existence of circular orbits at any arbitrary radius r in the equatorial plane demands that V eff (r) = 0 and dV eff /dr = 0. These conditions allow us to write the kinematic parameters as where, throughout the paper, X ,r ≡ dX/dr. Stability of orbits depends on the signs of d 2 V eff /dr 2 , while the condition d 2 V eff /dr 2 = 0 gives the inflection point or marginally stable (ms) orbit or innermost stable circular orbit (ISCO) at r = r ms . The thin disk is assumed to have an inner edge defined by the marginally stable circular radius r ms , while the orbits at higher radii than r = r ms are Keplerian. Note that marginally stable circular orbits play a crucial role when an accreting BH is surrounded, e.g., by quintessence matter [14,60]. The generic Eqs. (7)(8)(9)(10)(11)(12)(13)(14)(15)(16)(17)(18)(19) are valid for any rotating spacetime that will be explicitly calculated in what follows.

IV. ACCRETION DISK PROPERTIES OF RDSWH
We start with a Kerr-like wormhole spacetime, recently considered by Bueno et al. [3] that generalized the static Damour-Solodukhin wormhole [1], which we call here RDSWH. The spacetime metric in Boyer-Lindquist coordinates is given by where Σ and∆ are expressed by Σ ≡ r 2 + a 2 cos 2 θ,∆ ≡ r 2 − 2M (1 + λ 2 )r + a 2 . It contains a family of parameters where a and M corresponds to the spin and the mass of wormhole, and λ 2 is the deformation parameter 1 . A non vanishing λ 2 differs this metric from Kerr metric, we can recover the Kerr metric when λ 2 = 0. Although the throat of the wormhole can be easily obtained by equating the∆ to zero, which represents a special region that connects two different asymptotically flat regions. For any values of λ = 0 the metric is regular everywhere and with small values of λ 2 ∼ 0, it is practically indistinguishable from a Kerr BH. Special feature of RDSWH metric is that kinematic properties of accretion disk coincide with those of Kerr BH. As an example let us concider the effective potential, given by Eq. (16), which determines the geodesic motion of the test particles in the equatorial plane of RDSWH 1 We clarify that M here is the Keplerian mass as considered by Bueno et al. [3], viz., gtt ∼ 1 − 2M r in their Eq.(36) and explicitly so stated after their RDSWH metric (46). Keplerian mass is the one measured by the orbiting particles in the accretion disk or the one responsible for the echoes [3], both being strong field phenomena. On the other hand, the ADM mass is measured by an asymptotic static observer and is given by [59]: M ADM = M 1 + λ 2 . While the ADM mass is not relevant for studies in accretion or echoes, the existence of two distinct masses here reveal that the RDSWH possesses characteristics of Machian Brans-Dicke solutions, where the Keplerian mass is distinct from the ADM mass. It would thus be rewarding to relate λ to some kind of scalar field and develop a scalar-tensor theory for which RDSWH could be an exact solution. We thank the anonymous reviewer whose query has led us to think of this interesting possibility to be pursued in the future. solution (20)- (21). This is given by (23) Here we see that expression for effective potential of RDSWH solution does not depend on λ and therefore coincides with expression for effective potential of Kerr BH and so the ISCO radius or r ms is the same as that of Kerr BH.
The remaining expressions determining kinematic properties of RDSWH are specific energy, specific angular momentum, angular velocity and radius of the marginally stable orbit which follow from Eqs.(17)- (19). They also coincide with those of Kerr BH: where The reason of coincidence of kinematic properties is that parameter λ appears only in the g rr component of metric which is not used in Eqs.(16)- (19). We shall consider a stellar sized wormhole of mass 15M as a central compact object with an accretion ratė M 0 ∼ 10 19 gm.sec −1 and study the effect of different values of λ in the spacetime described by RDSWH solution (20)- (21). We should ensure that r ms > r + , which is a necessary condition so that the accreting particles do not reach the throat and disappear into the other universe even for extreme limit a = 1. A special feature of the RDSWH is that it does satisfy this condition for λ < 10 −3 , the strong field lensing constraint for mimicking BHs [13] (see Fig.1).
For illustration, we will consider different values of spin parameter a . Choice a = 0 corresponds to static case of DSWH, which in turn goes to Schwarzschild solution when λ = 0. Choice a = 0.75 is expected from the collapse of a maximally rotating polytropic star [61]. Another highly probable value of spin parameter is a = 0.95, which comes from different observational methods [62,63]. Case a = 0.998 is a maximally allowed value coming from the Page-Thorne limit for Kerr BH [28].

V. NUMERICAL ESTIMATES
We present below three tables showing observable characteristics of the thin disk (Table 1), the minimum sta-ble radius and the accretion efficiency (Table2) and the difference of the maxima of flux of radiation between RDSWH and Kerr BH (Table 3).
Figs.2a-2d display the flux of radiation F (r) emitted by the disk between r ms and received at an arbitrary radius r away from the center of the disk [Eq. (7)]. Figs.3a-3d show variation of temperature over the disk from r ms to an arbitrary radius and Figs.4a-4d show observed luminosity variations over different frequency ranges. The peaks values of emissivity properties of accretion disk are given in the Tab.1. In Tab.2, we show the variation of marginally stable orbits r ms and with the spin parameter a in the range used in the previous plots. The conversion efficiency of RDSWH is independent of λ, so is identical with that of Kerr BH.
From Figs 2a-2d we see that the flux emitted from the disk around Kerr BH always bigger than RDSWH depending on a . To show the difference of profiles we will use their maxima given in Table 1. With λ ≤ 0.1 flux of the accretion disk is practically indistinguishable from Kerr BH. Only limiting case a = 0.998 shows the difference of fluxes up to 6 %. With the increase of the spin parameter the difference of fluxes of Kerr BH and RDSWH is also increasing. The most significant difference appears when λ = 1, here we see that maximum value of flux of the wormhole is only 2 % of maximum value of Kerr BH flux. As r → ∞ the fluxes of RDSWH and Kerr BH become indistinguishable. Another interesting feature of RDSWH is the shift of maxima of fluxes with the increase of λ and spin parameter a .
Similar characteristics appear in the disk temperature profiles, depicted in Figs. 3a-3d. For all values of the spin parameter the disks rotating around the Kerr black holes are hotter than those around the wormhole. With the increase of parameter λ the temperature of the disk around rotating WH decreases.
In Figs. 4a-4d, we display the disk spectra for the spinning wormholes compared to Kerr BH. The cut-off frequencies for the Kerr black hole are systematically higher than those for the wormholes. The biggest deviations in radiation emitted by the accretion disk appears in ultraviolet range. Same as flux and temperature profiles the emissivity profile displays the most significant difference between RDSWH and Kerr BH when a = 0.998 and λ = 1.
Recently, by Bueno et al. [3] have considered small values λ = 10 −5 , 10 −10 and a = 0.7, to analyze how the quasinormal modes of the RDSWH differ from those of the Kerr BH. If we consider the accretion disk properties, we can see, that only flux of radiation of wormhole is differ from BH. Others, temperature and emission spectra  indistinguishable then λ = 10 −5 , 10 −10 .

VI. CONCLUSIONS
Accretion around spinning wormholes is an active area of current research. Along this line, gravitational wave echoes of RDSWH have been recently studied by Bueno et al. [3]. We studied in this work the thin accretion disk profiles of the same wormhole using the steady-state Page-Thorne model. Our aim was to examine how far the attractive positive mass mouth of RDSWH mimicks profiles of the Kerr BH for small values of λ inspired by strong field lensing, which in general is an excellent diagnostic for comparison between two categories of objects (BH and WH) with differing topology, see. e.g., [13,14].
We chose for illustration a rotating stellar sized wormhole as the accreting object and obtained certain non-trivial generic characteristics of the disk. Specifically, we show that RDSWH parameter λ has no influence on the kinematic profiles of the disc. That is, the accretion efficiency , radius of the marginally stable orbit r ms (Table 2), the specific angular momentum L, the specific energy E, angular velocity Ω for a test particle in a circular equatorial orbit around RDSWH are independent of λ.
However, differences appear in the emissivity properties for higher values 0.1 < λ ≤ 1 (say) and for the extreme spin a = 0.998. We numerically computed the time averaged flux F (r), temperature distribution T (r) and emission spectra νL(ν) of the accretion disk for different λ for the same dimensionless spin parameter a . The relevant results are displayed in Figs. 2a-2d, 3a-3d, and 4a-4d for an accreting wormhole with assumed mass 15M and accretion rateṀ 0 ∼ 10 19 gm.sec −1 for λ = 0 (Kerr BH), λ = 0.1, λ = 0.5 and λ = 1 (RDSWH) respectively. Specifically, we observe certain characteristics of the emissivity profile, viz., the maxima (peak) of the F (r), temperature distribution T (r) and emission spectra νL(ν) shift in the following manner as can be read off from Table 1: For a fixed a , an increment of λ leads to an increasing shift in r = r crit at which the peaks of F max (r), T max (r) decrease, while there is a decreasing shift in the critical frequency ν crit at which the peak of νL(ν) max decreases. On the other hand, for a fixed λ, an increment of a leads to a decreasing shift in r = r crit at which the peaks of F max (r), T max (r) increase, while there is an increasing shift in the critical frequency ν crit at which the peak of νL(ν) max increases. These are generic features as variations of the wormhole mass and the rate of accretion within the RDSWH model preserve these characteristics. Thus the behavior of the luminosity peak is generically quite opposite to each other in the indicated variations, which in principle could be useful from the viewpoint of observations. Apart from luminosity, it follows that an estimate of the difference ∆ λ in the maxima of flux of radiation F (r) between RDSWH and Kerr BH shows non-zero values ( Table 3) but is unobservable at present for small values λ < 10 −3 permitted by the strong lensing bound. Given the fact that the accretion profiles are not known with great precision, the overall conclusion is that RDSWH are experimentally indistinguishable from KBH by accretion characteristics as of now. Nevertheless, microlensing perturbations to the flux ratios of gravitationally lensed quasar images can constrain the temperature profiles of their accretion disks (connected by Stefan-Boltzmann law), which could be a more accurate test for the accretion process [64] and may effectively lead to a tighter constraint on λ. Work is underway.