The measurements of quasi-periodic oscillations (QPOs) provide a quite powerful tool to test the nature of astrophysical black hole candidates in the strong gravitational field regime. In this paper, we use QPOs within the relativistic precession model to test a recently proposed family of rotating black hole mimickers, which reduce to the Kerr metric in a limiting case, and can represent traversable wormholes or regular black holes with one or two horizons, depending on the values of the parameters. In particular, assuming that the compact object of GRO J1655-40 is described by a rotating black hole mimicker, we perform a χ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\chi $$\end{document}-square analysis to fit the parameters of the mimicker with two sets of observed QPO frequencies from GRO J1655-40. Our results indicate that although the metric around the compact object of GRO J1655-40 is consistent with the Kerr metric, a regular black hole with one horizon is favored by the observation data of GRO J1655-40.
publisher-imprint-nameSpringervolume-issue-count12issue-article-count88issue-toc-levels0issue-pricelist-year2021issue-copyright-holderThe Author(s)issue-copyright-year2021article-contains-esmNoarticle-numbering-styleContentOnlyarticle-registration-date-year2021article-registration-date-month11article-registration-date-day8article-toc-levels0toc-levels0volume-typeRegularjournal-productNonStandardArchiveJournalnumbering-styleContentOnlyarticle-grants-typeOpenChoicemetadata-grantOpenAccessabstract-grantOpenAccessbodypdf-grantOpenAccessbodyhtml-grantOpenAccessbibliography-grantOpenAccessesm-grantOpenAccessonline-firstfalsepdf-file-referenceBodyRef/PDF/10052_2021_Article_9816.pdfpdf-typeTypesettarget-typeOnlinePDFissue-typeRegulararticle-typeOriginalPaperjournal-subject-primaryPhysicsjournal-subject-secondaryElementary Particles, Quantum Field Theoryjournal-subject-secondaryNuclear Physics, Heavy Ions, Hadronsjournal-subject-secondaryQuantum Field Theories, String Theoryjournal-subject-secondaryMeasurement Science and Instrumentationjournal-subject-secondaryAstronomy, Astrophysics and Cosmologyjournal-subject-secondaryNuclear Energyjournal-subject-collectionPhysics and Astronomyopen-accesstrueIntroduction
The first observations of gravitational waves by LIGO [1] and the first image of a black hole in the galaxy M87 [2] have ushered us into a new era of testing general relativity (GR) in the strong gravity regime. On the other hand, quasi-periodic oscillations (QPOs) are observed in the X-ray flux from black hole and neutron star X-ray binary systems, and detected as narrow peaks in the power density spectrum [3]. QPOs are believed to be associated with motion and accretion-related timescales in a region of order the Schwarzschild radius around the compact object, which makes QPOs excellent probes of the strong gravitational field regime [4–8]. In black hole systems, the observed frequencies of QPOs range from mHz to hundreds of Hz. While low frequency QPOs ≲30Hz\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\left( \lesssim 30\,\mathrm {Hz}\right) $$\end{document} are commonly observed from black hole X-ray binaries [9], high frequency QPOs ≳60Hz\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\left( \gtrsim 60\,\mathrm {Hz}\right) $$\end{document} are very rare. In fact, the Rossi X-ray Timing Explorer (RXTE), operational between 1996 and 2012, first detected high frequency QPOs in black hole systems [10, 11]. Interestingly, a pair of simultaneous high-frequency QPOs was first discovered in the X-ray flux from GRO J1655-40 by RXTE [12]. It was noted that the frequencies of the two high-frequency QPOs are in a 3:2 ratio, suggesting a resonance between orbital and epicyclic motion of accreting matter near the innermost stable circular orbit (ISCO) of black holes [13]. Later, three simultaneous QPO frequencies, consisting of two higher frequencies and one lower frequency, were also observed from the X-ray data of GRO J1655-40 [14].
On the theoretical side, various models have been proposed to explain QPOs by relating them to the orbital and epicyclic frequencies of geodesics, such as the relativistic precession model (RPM) [15], the tidal disruption model [16], the parametric resonance model [17], the resonance model [13, 18, 19], the warped disk oscillation model [20] and the non-axisymmetric disk oscillation model [21, 22]. With these models, the observation data of QPO frequencies have been used to constrain the parameters of X-ray binary systems and test the nature of various gravity theories [23–37].
Particularly, it showed that the X-ray data of GRO J1655-40, especially the QPO triplet, fit nicely in the RPM, and the mass and spin of the compact object of GRO J1655-40 can be precisely determined [14]. Remarkably, the inferred mass is in great agreement with the dynamical mass measurement [23]. The RPM was originally proposed to explain QPOs in low-mass X-ray binaries with a neutron star [38], and later extended to systems with stellar-mass BH candidates [15]. In the RPM, QPO frequencies are assumed to be related to fundamental frequencies of a test particle orbiting a central object. The twin higher frequencies are regarded as the azimuthal frequency νϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\phi }$$\end{document} and the periastron precession frequency νper\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\text {per}}$$\end{document} of quasi-circular orbits in the innermost disk region, respectively. The low-frequency QPO is identified as the nodal precession frequency νnod\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\text {nod}}$$\end{document}, which is emitted at the same radius where the twin higher frequencies are generated.
On the other hand, curvature singularities can be formed during a gravitational collapse. It is commonly believed that singularities can be avoided through quantum gravitational effects. Consequently, since Bardeen proposed the first regular black hole [39], constructing and studying classical black holes without singularities have been a topic of considerable interest in GR and astrophysics communities due to their non-singular property [40–47]. Recently, Simpson and Visser proposed a static and spherically symmetric regular spacetime described by the line element (dubbed the SV metric henceforth),ds2=-1-2Mr2+ℓ2dt2+1-2Mr2+ℓ2-1dr2+r2+ℓ2dθ2+sin2θdϕ2,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} ds^{2}= & {} -\left( 1-\frac{2M}{\sqrt{r^{2}+\ell ^{2}}}\right) dt^{2}+\left( 1-\frac{2M}{\sqrt{r^{2}+\ell ^{2}}}\right) ^{-1}dr^{2}\nonumber \\&+\left( r^{2}+\ell ^{2}\right) \left( d\theta ^{2}+\sin ^{2}\theta d\phi ^{2}\right) ~, \end{aligned}$$\end{document}where M≥0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$M\ge 0$$\end{document} represents the ADM mass, and ℓ>0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell >0$$\end{document} is a parameter responsible for regularizing the center singularity [48–50]. The most appealing feature of the SV metric (1) is that it can smoothly interpolate between a regular black hole for ℓ<2M\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell <2M$$\end{document} and a wormhole for ℓ≥2M\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell \ge 2M$$\end{document}. In the limit of ℓ=0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell =0$$\end{document}, the SV metric reduces to the ordinary Schwarzschild spacetime. Subsequently, the properties of the SV metric were discussed, e.g., quasinormal modes [51], precessing and periodic geodesic motions [52], gravitational lensing [53–55] and shadows [56, 57]. To make the SV metric more relevant to realistic situations, the SV metric (1) was generalized to a family of rotating black hole mimickers, dubbed the rotating SV metric, which may represent a rotating wormhole and a rotating regular black hole with one or two horizons [58]. The strong gravitational lensing and the shadows of the rotating SV metric have been investigated in [59, 60], respectively.
In this paper, we test gravity with QPOs frequencies observed from GRO J1655-40 within the RPM for the rotating SV metric. The content of this paper is as follows. After we briefly review the RPM and the rotating SV metric in Sects. 2 and 3, respectively, the epicyclic frequencies of the rotating SV are computed in Sect. 3. In Sect. 4, we use the data of GRO J1655-40 to put constraints on the parameters of the rotating SV metric. Section 5 is devoted to our conclusions. Throughout the paper, we use units in which G=c=1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$G=c=1$$\end{document}.
Epicyclic frequencies
In this section, we consider the timelike geodesic equations in a stationary and axially symmetric spacetime, and then derive the expressions of the epicyclic frequencies. The metric of a stationary and axially symmetric spacetime which satisfies the circularity condition is given by [61]ds2=gttdt2+grrdr2+2gtϕdtdϕ+gθθdθ2+gϕϕdϕ2,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} ds^{2}=g_{tt}dt^{2}+g_{rr}dr^{2}+2g_{t\phi }dtd\phi +g_{\theta \theta }d\theta ^{2}+g_{\phi \phi }d\phi ^{2}, \end{aligned}$$\end{document}where the metric gμν\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$g_{\mu \nu }$$\end{document} is a function of r and θ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\theta $$\end{document}, and we drop the coordinate dependence of the metric functions to simplify the notation. For a massive particle travelling along a time-like world line xμτ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$x^{\mu }\left( \tau \right) $$\end{document} with τ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\tau $$\end{document} being the proper time, the four-velocity Uμ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$U^{\mu }$$\end{document} is defined by Uμ=x˙μ=dxμ/dτ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$U^{\mu }={\dot{x}}^{\mu }=dx^{\mu }/d\tau $$\end{document}, which satisfies UμUμ=-1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$U^{\mu }U_{\mu }=-1$$\end{document}. Due to the stationarity and axisymmetry, the metric (2) admits two Killing vectors,Kμ=(∂t)μ=(1,0,0,0)andRμ=(∂ϕ)μ=(0,0,0,1).\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&K^{\mu }=(\partial _{t})^{\mu }=(1,0,0,0)\quad \text { and } \nonumber \\&R^{\mu }=(\partial _{\phi })^{\mu }=(0,0,0,1). \end{aligned}$$\end{document}The two Killing vectors correspond to two conserved quantities of geodesic motion,E=-Kμdxμdτ=-gttt˙-gtϕϕ˙,Lz=Rμdxμdτ=gtϕt˙+gϕϕϕ˙,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} E&=-K_{\mu }\frac{dx^{\mu }}{d\tau }=-g_{tt}{\dot{t}}-g_{t\phi }\dot{\phi },\nonumber \\ L_{z}&=R_{\mu }\frac{dx^{\mu }}{d\tau }=g_{t\phi }{\dot{t}}+g_{\phi \phi } {\dot{\phi }}, \end{aligned}$$\end{document}which can be interpreted as the energy per unit mass and the angular momentum per unit mass along the axis of symmetry, respectively. In terms of E and Lz\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$L_{z}$$\end{document}, one can express t˙\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{t}}$$\end{document} and ϕ˙\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{\phi }}$$\end{document} ast˙=gϕϕE+gtϕLzgtϕ2-gttgϕϕ,ϕ˙=gtϕE+gttLzgttgϕϕ-gtϕ2.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} {\dot{t}}=\frac{g_{\phi \phi }E+g_{t\phi }L_{z}}{g_{t\phi }^{2}-g_{tt}g_{\phi \phi } }\text {, }{\dot{\phi }}=\frac{g_{t\phi }E+g_{tt}L_{z}}{g_{tt}g_{\phi \phi } -g_{t\phi }^{2}}\text {.} \end{aligned}$$\end{document}Using the above equations, we can rewrite UμUμ=-1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$U^{\mu }U_{\mu }=-1$$\end{document} asgrrr˙2+gθθθ˙2=Veffr,θ,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} g_{rr}{\dot{r}}^{2}+g_{\theta \theta }{\dot{\theta }}^{2}=V_{\text {eff}}\left( r,\theta \right) , \end{aligned}$$\end{document}where the effective potential Veffr,θ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_{\text {eff}}\left( r,\theta \right) $$\end{document} is defined asVeff(r,θ)=E2gϕϕ+2ELzgtϕ+Lz2gttgtϕ2-gttgϕϕ-1.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} V_{\text {eff}}(r,\theta )=\frac{E^{2}g_{\phi \phi }+2EL_{z}g_{t\phi }+L_{z} ^{2}g_{tt}}{g_{t\phi }^{2}-g_{tt}g_{\phi \phi }}-1. \end{aligned}$$\end{document}We consider a circular geodesic at r=r¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$r={\bar{r}}$$\end{document} on the equatorial plane with θ=π/2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\theta =\pi /2$$\end{document}, which means that the effective potential may develop a double root at r=r¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$r={\bar{r}}$$\end{document} on the equatorial plane, i.e., Veff(r¯,π/2)=∂rVeff(r¯,π/2)=0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_{\text {eff}}({\bar{r}},\pi /2)=\partial _{r}V_{\text {eff}}({\bar{r}},\pi /2)=0$$\end{document}. Along the circular orbit, the angular velocity of the particle measured by an observer at infinity is defined byΩϕ=dϕdt=ϕ˙t˙=-gtϕE+gttLzgϕϕE+gtϕLzr=r¯,θ=π/2.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \Omega _{\phi }=\frac{d\phi }{dt}=\frac{{\dot{\phi }}}{{\dot{t}}}=\left. -\frac{g_{t\phi }E+g_{tt}L_{z}}{g_{\phi \phi }E+g_{t\phi }L_{z}}\right| _{r={\bar{r}},\theta =\pi /2}\text {.} \end{aligned}$$\end{document}Solving Veff(r¯,π/2)=0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_{\text {eff}}({\bar{r}},\pi /2)=0$$\end{document} with Eq. (8) gives the specific energy and angular momentum of the particle,E=-gtt+gtϕΩϕ-gtt-2gtϕΩϕ-gϕϕΩϕ2r=r¯,θ=π/2,Lz=gtϕ+gϕϕΩϕ-gtt-2gtϕΩϕ-gϕϕΩϕ2r=r¯,θ=π/2.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} E&=\left. -\frac{g_{tt}+g_{t\phi }\Omega _{\phi }}{\sqrt{-g_{tt}-2g_{t\phi }\Omega _{\phi }-g_{\phi \phi }\Omega _{\phi }^{2}}}\right| _{r={\bar{r}},\theta =\pi /2},\nonumber \\ L_{z}&=\left. \frac{g_{t\phi }+g_{\phi \phi }\Omega _{\phi }}{\sqrt{-g_{tt}-2g_{t\phi }\Omega _{\phi }-g_{\phi \phi }\Omega _{\phi }^{2}}}\right| _{r={\bar{r}},\theta =\pi /2}. \end{aligned}$$\end{document}Solving ∂rVeff(r¯,π/2)=0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\partial _{r}V_{\text {eff}}({\bar{r}},\pi /2)=0$$\end{document} for Ωϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega _{\phi }$$\end{document}, one obtains the angular velocity of the particle,Ωϕ=-∂rgtϕ±∂rgtϕ2-∂rgtt∂rgϕϕ∂rgϕϕr=r¯,θ=π/2,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \Omega _{\phi }=\left. \frac{-\partial _{r}g_{t\phi }\pm \sqrt{\left( \partial _{r}g_{t\phi }\right) ^{2}-\left( \partial _{r}g_{tt}\right) \left( \partial _{r}g_{\phi \phi }\right) }}{\partial _{r}g_{\phi \phi }}\right| _{r={\bar{r}},\theta =\pi /2}, \end{aligned}$$\end{document}where the sign +/-\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$+/-$$\end{document} corresponds to a prograde/retrograde orbit. The stability of the circular orbit is determined by the sign of ∂r2Veff(r¯,π/2)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\partial _{r} ^{2}V_{\text {eff}}({\bar{r}},\pi /2)$$\end{document}, i.e., ∂r2Veff(r¯,π/2)>0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\partial _{r}^{2}V_{\text {eff}} ({\bar{r}},\pi /2)>0$$\end{document}⇔\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Leftrightarrow $$\end{document} unstable and ∂r2Veff(r¯,π/2)<0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\partial _{r}^{2} V_{\text {eff}}({\bar{r}},\pi /2)<0$$\end{document}⇔\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Leftrightarrow $$\end{document} stable. The transition between stable and unstable circular orbits, which is determined by ∂r2Veff(rISCO,π/2)=0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\partial _{r}^{2}V_{\text {eff}}(r_{\text {ISCO}},\pi /2)=0$$\end{document}, is the ISCO, which is located at r=rISCO\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$r=r_{\text {ISCO}}$$\end{document} on the equatorial plane.
To derive the epicyclic frequencies associated with the circular orbit, we consider small perturbations of the orbit in both the radial and the vertical directions,r(t)=r¯+δr(t),θ(t)=π2+δθ(t).\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} r(t)={\bar{r}}+\delta r(t)\text {,}\quad \theta (t)=\frac{\pi }{2}+\delta \theta (t). \end{aligned}$$\end{document}Inserting Eq. (11) into Eq. (6) yields the differential equations for the perturbations δr(t)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\delta r(t)$$\end{document} and δθ(t)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\delta \theta (t)$$\end{document},d2δr(t)dt2+Ωr2δr(t)=0,d2δθ(t)dt2+Ωθ2δθ(t)=0,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \frac{d^{2}\delta r(t)}{dt^{2}}+\Omega _{r}^{2}\delta r(t)=0\text {,}\quad \frac{d^{2}\delta \theta (t)}{dt^{2}}+\Omega _{\theta }^{2}\delta \theta (t)=0, \end{aligned}$$\end{document}where the frequencies of the oscillations areΩr2=-12grrt˙2∂2Veffr,θ∂r2r=r¯,θ=π/2,Ωθ2=-12gθθt˙2∂2Veffr,θ∂θ2r=r¯,θ=π/2.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \Omega _{r}^{2}= & {} \left. -\frac{1}{2g_{rr}{\dot{t}}^{2}}\frac{\partial ^{2}V_{\text {eff}}\left( r,\theta \right) }{\partial r^{2}}\right| _{r={\bar{r}},\theta =\pi /2},\text { }\Omega _{\theta }^{2}\nonumber \\= & {} \left. -\frac{1}{2g_{\theta \theta }{\dot{t}}^{2}}\frac{\partial ^{2}V_{\text {eff}}\left( r,\theta \right) }{\partial \theta ^{2}}\right| _{r={\bar{r}},\theta =\pi /2}. \end{aligned}$$\end{document}We then define νϕ=Ωϕ/2π\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\phi }=\Omega _{\phi }/2\pi $$\end{document}, νr=Ωr/2π\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{r}=\Omega _{r}/2\pi $$\end{document} and νθ=Ωθ/2π\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\theta }=\Omega _{\theta }/2\pi $$\end{document} as the azimuthal, radial and vertical epicyclic frequencies, respectively. The periastron precession frequency νper\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\text {per}}$$\end{document} and the nodal precession frequency νnod\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\text {nod}}$$\end{document} are defined by νper=νϕ-νr\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\text {per }}=\nu _{\phi }-\nu _{r}$$\end{document} and νnod=νϕ-νθ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\text {nod }} =\nu _{\phi }-\nu _{\theta }$$\end{document}, respectively.
Plots of the effective potential Veff\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_{\text {eff}}$$\end{document} of a massive particle travelling on the equatorial plane with θ=π/2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\theta =\pi /2$$\end{document} as a function of radius r in a traversable wormhole with ℓ=12\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell =12$$\end{document} (left panel), a RBH-I black hole with ℓ=1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell =1$$\end{document} (middle panel) and a RBH-II black hole with ℓ=0.08\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell =0.08$$\end{document} (right panel) for Lz=25\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$L_{z}=25$$\end{document} (red lines), Lz=26\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$L_{z}=26$$\end{document} (green lines) and Lz=27\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$L_{z}=27$$\end{document} (blue lines). We have M=6\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$M=6$$\end{document}, a=1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$a=1$$\end{document} and E=1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$E=1$$\end{document}. The horizons of the RBH-I and RBH-II black holes are represented by orange dashed lines. In the three cases, the effective potential always decreases to zero at infinity. For the traversable wormhole, the effective potential has one maximum and one minimum in each universe. For the RBH-I black hole, the effective potential has one maximum and one minimum outside the event horizon. The effective potential of the RBH-II black hole is similar to that of the RBH-I black hole
Rotating Simpson–Visser metric
In [58], a rotating generalization of the static and spherically symmetric metric (1) has been constructed by employing the Newman–Janis procedure [62]. This stationary and axially symmetric metric can describe a rotating traversable wormhole and a rotating regular black hole with one or two horizons. In particular, the rotating SV metric readsds2=-1-2Mr2+ℓ2Σdt2+ΣΔdr2+Σdθ2-4Masin2θr2+ℓ2Σdtdϕ+Asin2θΣdϕ2,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} ds^{2}= & {} -\left( 1-\frac{2M\sqrt{r^{2}+\ell ^{2}}}{\Sigma }\right) dt^{2} \nonumber \\&+\frac{\Sigma }{\Delta }dr^{2}+\Sigma d\theta ^{2}-\frac{4Ma\sin ^{2}\theta \sqrt{r^{2}+\ell ^{2}}}{\Sigma }dtd\phi \nonumber \\&+\frac{A\sin ^{2}\theta }{\Sigma }d\phi ^{2}, \end{aligned}$$\end{document}withΣ=r2+ℓ2+a2cos2θ,Δ=r2+ℓ2+a2-2Mr2+ℓ2,A=r2+ℓ2+a22-Δa2sin2θ,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \Sigma&=r^{2}+\ell ^{2}+a^{2}\cos ^{2}\theta ,\nonumber \\ \Delta&=r^{2}+\ell ^{2}+a^{2}-2M\sqrt{r^{2}+\ell ^{2}},\nonumber \\ A&=\left( r^{2}+\ell ^{2}+a^{2}\right) ^{2}-\Delta a^{2}\sin ^{2} \theta , \end{aligned}$$\end{document}where a is the spin parameter. The rotating SV metric will reduce to the SV metric (1) if a=0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$a=0$$\end{document} and to the Kerr metric if ℓ=0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell =0$$\end{document}. Interestingly, the rotating SV metric is everywhere regular when ℓ>0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell >0$$\end{document} [58].
The horizons of the rotating SV metric are determined by Δ=0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Delta =0$$\end{document}, whose solutions arer±=M±M2-a22-ℓ2.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} r_{\pm }=\sqrt{\left( M\pm \sqrt{M^{2}-a^{2}}\right) ^{2}-\ell ^{2}}. \end{aligned}$$\end{document}As shown in [58], the phases of the rotating SV metric are determined by the existence of r±\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$r_{\pm }$$\end{document}. Specifically, the rotating SV metric represents
a regular black hole with one horizon (RBH-I): M-M2-a2<ℓ<M+M2-a2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$M-\sqrt{M^{2}-a^{2} }<\ell <M+\sqrt{M^{2}-a^{2}}$$\end{document} and M>a\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$M>a$$\end{document};
a regular black hole with two horizons (RBH-II): ℓ<M-M2-a2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell <M-\sqrt{M^{2}-a^{2}}$$\end{document} and M>a\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$M>a$$\end{document};
there limiting cases: a one-way wormhole with a null throat when ℓ=M+M2-a2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell =M+\sqrt{M^{2}-a^{2}}$$\end{document}, a regular black hole with one horizon and a null throat when ℓ=M-M2-a2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell =M-\sqrt{M^{2}-a^{2}}$$\end{document} and an extremal regular black hole when M=a\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$M=a$$\end{document}.
Plots of the azimuthal epicyclic frequency νϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\phi }$$\end{document} (left panel), the periastron precession frequency νper\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\text {per}}$$\end{document} (middle panel) and the nodal precession frequency νnod\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\text {nod}}$$\end{document} (right panel) of prograde (solid lines) and retrograde (dashed lines) orbits at r¯=6\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\bar{r}}=6$$\end{document} as a function of ℓ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell $$\end{document} for a=0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$a=0$$\end{document} (red lines), a=0.4\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$a=0.4$$\end{document} (green lines) and a=0.8\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$a=0.8$$\end{document} (blue lines) in the rotating SV metric with M=1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$M=1$$\end{document}. The a=0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$a=0$$\end{document} case corresponds to the SV metric, which has νnod=0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\text {nod}}=0$$\end{document} due to spherical symmetry. Except νnod\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\text {nod}}$$\end{document} with a=0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$a=0$$\end{document}, the magnitudes of νϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\phi }$$\end{document}, νper\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\text {per}}$$\end{document} and νnod\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\text {nod}}$$\end{document} decrease as ℓ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell $$\end{document} grows with a fixed a. For the retrograde orbits in the rotating SV metric with a>0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$a>0$$\end{document}, the values of νnod\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\text {nod}}$$\end{document} are shown to be negative
Substituting the rotating SV metric (14) into Eq. (7), one can obtain the effective potential Veff\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_{\text {eff}}$$\end{document} for a massive particle in the rotating SV metric. Like a Schwarzschild black, we find that, for a massive particle with a large enough Lz\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\left| \mathbf {L_{z}}\right| $$\end{document}, Veff\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_{\text {eff}}$$\end{document} possesses a minimum and a maximum on the equatorial plane, corresponding to a stable circular orbit and an unstable one, outside the event horizon of a black hole or in the each universe of a wormhole. As Lz\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\left| \mathbf {L_{z}}\right| $$\end{document} increases, the two circular orbits move further away from each other. In Fig. 1, we plot Veff\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$V_{\text {eff}}$$\end{document} as a function of r for various values of Lz\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$L_{z}$$\end{document} in three cases, i.e., a traversable wormhole, a RBH-I black hole and a RBH-II black hole. Here without loss of generality, we have set M=6\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$M=6$$\end{document}, E=1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$E=1$$\end{document}, a=1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$a=1$$\end{document} and θ=π/2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\theta =\pi /2$$\end{document}. Note that the two universes of the wormhole are connected at r=ℓ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$r=\ell $$\end{document}, where the effective potential has a cusp.
For a circular orbit at r=r¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$r={\bar{r}}$$\end{document} on the equatorial plane, substituting the rotating SV metric (14) into Eqs. (8) and (13) gives epicyclic frequencies in the rotating SV metric,νϕ=12πM1/2(r¯2+ℓ2)3/4±aM1/2,νr=r¯νϕ(r¯2+ℓ2)1/21-6M(r¯2+ℓ2)1/2-3a2(r¯2+ℓ2)±8aM1/2(r¯2+ℓ2)3/4,νθ=νϕ1∓4aM1/2(r¯2+ℓ2)3/4+3a2(r¯2+ℓ2),\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \nu _{\phi }&=\frac{1}{2\pi }\frac{M^{1/2}}{({\bar{r}}^{2}+\ell ^{2})^{3/4}\pm aM^{1/2}},\nonumber \\ \nu _{r}&=\frac{{\bar{r}}\nu _{\phi }}{({\bar{r}}^{2}+\ell ^{2})^{1/2}}\nonumber \\&\quad \sqrt{1-\frac{6M}{({\bar{r}}^{2}+\ell ^{2})^{1/2}}-\frac{3a^{2}}{({\bar{r}} ^{2}+\ell ^{2})}\pm \frac{8aM^{1/2}}{({\bar{r}}^{2}+\ell ^{2})^{3/4}}},\nonumber \\ \nu _{\theta }&=\nu _{\phi }\sqrt{1{\mp }\frac{4aM^{1/2}}{({\bar{r}}^{2}+\ell ^{2})^{3/4}}+\frac{3a^{2}}{({\bar{r}}^{2}+\ell ^{2})}}, \end{aligned}$$\end{document}which can be used to constrain the four parameters ℓ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell $$\end{document}, r¯\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\bar{r}}$$\end{document}, M, and a of the rotating SV metric. Here, the top/bottom row of the ± and ∓\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\mp }$$\end{document} signs corresponds to the orbit co-rotating/counter-rotating with the spacetime. To illustrate the dependence of the epicyclic frequencies on ℓ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell $$\end{document}, we plot νϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\phi }$$\end{document}, νper\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\text {per}}$$\end{document} and νnod\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\text {nod}}$$\end{document} as a function of ℓ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell $$\end{document} for various values of a in Fig. 2, where M=1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$M=1$$\end{document} and r¯=6\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\bar{r}}=6$$\end{document}. The prograde and retrograde cases are represented by solid and dashed lines, respectively. The left panel shows that, for a fixed a, the azimuthal epicyclic frequency νϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\phi }$$\end{document} decreases with ℓ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell $$\end{document} increasing in both prograde and retrograde cases, whereas νϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\phi }$$\end{document} of the prograde orbit is smaller than that of the retrograde orbit. The periastron precession frequency νper\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\text {per}}$$\end{document} is displayed in the middle panel, and also decreases as ℓ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell $$\end{document} increases for both prograde and retrograde orbits. Like νϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\phi }$$\end{document}, the retrograde orbits have larger νper\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\text {per}}$$\end{document}. Note that retrograde circular orbits of radius r¯=6\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\bar{r}}=6$$\end{document} do not exist when ℓ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell $$\end{document} is small enough. However, as shown in the right panel, the nodal precession frequency νnod\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\text {nod}}$$\end{document} of the prograde/retrograde orbits decreases/increases as ℓ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell $$\end{document} increases for a given a. More interestingly, when a>0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$a>0$$\end{document}, νnod\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\text {nod}}$$\end{document} of the retrograde orbits is negative while that of the prograde orbits is positive. Note that the nodal precession is defined as the precession of the ascending node, which is the point where the orbit passes through the plane of reference, e.g., the equatorial plane in our paper. If the SV metric is spherically symmetric with a=0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$a=0$$\end{document}, the orbital plane of a particle remains fixed in space, and hence the nodal precession νnod\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\text {nod}}$$\end{document} vanishes, which can be seen from Fig. 2. However for the rotating SV metric with a nonzero a, a particle in a prograde/retrograde orbit has been shown to have a positive/negative nodal precession νnod\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\text {nod}}$$\end{document}, which means that the orbital plane rotates in the same/opposite direction as the direction of the SV metric’s rotation. Interestingly, it is well-known that the nodal precessions of a satellite in prograde and retrograde orbits around Earth have opposite signs [63]. Finally, it is noteworthy that the ISCO radius rISCO\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$r_{\text {ISCO}}$$\end{document} is determined by [58],rISCO2+ℓ2-6MrISCO2+ℓ2±8aMrISCO2+ℓ2=3a2,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} r_{\text {ISCO}}^{2}+\ell ^{2}-6M\sqrt{r_{\text {ISCO}}^{2}+\ell ^{2}}\pm 8a\sqrt{M\sqrt{r_{\text {ISCO}}^{2}+\ell ^{2}}}=3a^{2}, \end{aligned}$$\end{document}where +/-\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$+/-$$\end{document} is associated with the prograde/retrograde ISCO.
Constraining rotating Simpson–Visser metric by quasi-periodic oscillations
In this section, we use the RPM along with the QPO frequencies from GRO J1655-40 to put constraints on the parameters of the rotating SV metric. GRO J1655-40 is an X-ray binary, consisting of a primary star and a compact companion [64]. The measurement of the X-ray spectrum was found to exhibit type-C low-frequency QPOs and simultaneous high-frequency QPOs, which are observed in pairs and therefore dubbed lower and upper high-frequency QPOs [12, 14]. In particular, we consider two sets of QPOs with the observed frequencies based on the RXTE observations [14],ν1U=441Hz,σ1U=2Hz,ν1L=298Hz,σ1L=4Hz,ν1C=17.3Hz,σ1C=0.1Hz\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \nu _{1\mathrm {U}}=441\mathrm {~Hz},\quad&\sigma _{1\mathrm {U}}=2\mathrm {~Hz},\nonumber \\ \nu _{1\mathrm {L}}=298\mathrm {~Hz},\quad&\sigma _{1\mathrm {L}}=4\mathrm {~Hz},\nonumber \\ \nu _{1\mathrm {C}}=17.3\mathrm {~Hz},\quad&\sigma _{1\mathrm {C}} =0.1\mathrm {~Hz} \end{aligned}$$\end{document}andν2U=451Hz,σ2U=5Hz,ν2L=-,ν2C=18.3Hz,σ2C=0.1Hz.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \nu _{2\mathrm {U}}=451\mathrm {~Hz},\quad&\sigma _{2\mathrm {U}}=5\mathrm {~Hz},\nonumber \\ \nu _{2\mathrm {L}}=&~-~,\quad \nonumber \\ \nu _{2\mathrm {C}}=18.3\mathrm {~Hz},\quad&\sigma _{2\mathrm {C}} =0.1\mathrm {~Hz}. \end{aligned}$$\end{document}
Estimate of the mass M/M⊙\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$M/M_{\odot }$$\end{document}, of the spin parameter a/M and of the parameter ℓ/M\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell /M$$\end{document} of the rotating SV metric, which describes the spacetime around the compact object in GRO J1655-40, from interpreting the observations of QPOs in the RPM. Specifically, the best-fit values (red dots) and the contour levels of 1σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1\sigma $$\end{document} (orange lines), 2σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$2\sigma $$\end{document} (green lines) and 3σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$3\sigma $$\end{document} (blue lines) for the parameters are shown in the ℓ/M\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell /M$$\end{document}-a/M (left panel), ℓ/M\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell /M$$\end{document}-M/M⊙\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$M/M_{\odot }$$\end{document} (middle panel) and M/M⊙\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$M/M_{\odot }$$\end{document}-a/M (right panel) planes. In the yellow and blue regions, the rotating SV metric represents regular black hole solutions with one and two horizons, respectively. The observed frequencies of the QPOs in GRO J1655-40 are consistent with a Kerr black hole (the rotating SV metric with ℓ=0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell =0$$\end{document}), but they also allow for large deviations from the Kerr black hole solution. In particular, a regular black hole with one horizon is favored over the other phases of the rotating SV metric, e.g., a regular black hole with two horizons and a traversable wormhole
In the RPM, three simultaneous QPO frequencies are generated at the same radial coordinate in the accretion disk. The upper high-frequency QPOs correspond to the azimuthal epicyclic frequency νϕ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\phi }$$\end{document}, the lower high-frequency QPOs to the periastron precession frequency νper\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\text {per}}$$\end{document}, and the low-frequency QPOs to the nodal precession frequency νnod\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\nu _{\text {nod} }$$\end{document}. Moreover, it is reasonable to assume that the above two sets of QPOs result from two circular orbits of different radii, i.e., r¯1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\bar{r}}_{1}$$\end{document} and r¯2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${{\bar{r}}}_{2}$$\end{document}. In short, we have five free parameters: the mass M, the spin parameter a, the ℓ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell $$\end{document} parameter, and the radii r¯1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\bar{r}_{1}$$\end{document} and r¯2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${{\bar{r}}}_{2}$$\end{document} corresponding to the QPOs with three frequencies and two frequencies, respectively. To obtain the estimate of the five parameters from the observed QPO frequencies, we follow the procedure used in [28, 33, 34] and perform a χ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\chi $$\end{document}-square analysis withχ2(M,a,ℓ,r¯1,r¯2)=(ν1ϕ-ν1U)2σ1U+(ν1per-ν1L)2σ1L+(ν1nod-ν1C)2σ1C+(ν2ϕ-ν2U)2σ2U+(ν2nod-ν2C)2σ2C,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \chi ^{2}(M,a,\ell ,{{\bar{r}}}_{1},\bar{r}_{2})= & {} \frac{(\nu _{1\phi }-\nu _{1\mathrm {U} })^{2}}{\sigma _{1\mathrm {U}}}+\frac{(\nu _{1\text {per}}-\nu _{1\mathrm {L}})^{2} }{\sigma _{1\mathrm {L}}}\nonumber \\&+\frac{(\nu _{1\text {nod}}-\nu _{1\mathrm {C}})^{2} }{\sigma _{1\mathrm {C}}}\nonumber \\&+\frac{(\nu _{2\phi }-\nu _{2\mathrm {U}})^{2}}{\sigma _{2\mathrm {U}}}+\frac{(\nu _{2\text {nod}}-\nu _{2\mathrm {C}})^{2}}{\sigma _{2\mathrm {C}}}~,\nonumber \\ \end{aligned}$$\end{document}the minimum of which, χmin2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\chi _{\text {min}}^{2}$$\end{document}, occurs at the best estimate of M, aℓ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell $$\end{document}, r¯1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${{\bar{r}}}_{1}$$\end{document} and r¯2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${{\bar{r}}}_{2}$$\end{document}. The range of the parameters at a confidence level (C.L.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathrm {C.L.}$$\end{document}) is determined by the interval χmin2+Δχ2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\chi _{\text {min}}^{2}+\Delta \chi ^{2}$$\end{document}. In the case of five degrees of freedom, the intervals with Δχ2=5.89\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Delta \chi ^{2}=5.89$$\end{document}, 11.29 and 17.96 correspond to 68.3%\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$68.3\%$$\end{document}, 95.4%\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$95.4\%$$\end{document} and 99.7%\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$99.7\%$$\end{document}C.L.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathrm {C.L.}$$\end{document}, respectively, which are the probability intervals designated as 1, 2, and 3 standard deviation limits, respectively.
Computing χ2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\chi ^{2}$$\end{document}, we find χmin2=\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\chi _{\text {min}}^{2}=$$\end{document} 0.195 and obtain the best fits of the parameters of the rotating SV metric within 68.3%\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$68.3\%$$\end{document} credibility,M/M⊙=5.305-0.028+0.041,a/M=0.286-0.002+0.003,ℓ/M=0.347-0.347+0.011.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} M/M_{\odot }= & {} 5.305_{-0.028}^{+0.041}\text {, }a/M=0.286_{-0.002}^{+0.003}\text {, }\ell /\nonumber \\ M= & {} 0.347_{-0.347}^{+0.011}. \end{aligned}$$\end{document}Note that our result is consistent with the measurement of the mass M by optical and infrared observations, which give M=5.4±0.3M⊙\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$M=5.4\pm 0.3M_{\odot }$$\end{document} [23]. We present the best estimate and the 1σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1\sigma $$\end{document}, 2σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$2\sigma $$\end{document} and 3σ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$3\sigma $$\end{document} contour levels of M/M⊙\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$M/M_{\odot }$$\end{document}, a/M and ℓ/M\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell /M$$\end{document} in Fig. 3, where the yellow/blue regions represent the RBH-I/RBH-II phases of the rotating SV metric. As shown in the left and middle panels, while the hypothesis that the compact object of GRO J1655-40 is described by a Kerr black hole is consistent with the interpretation of the QPOs’ data in the RPM, significant deviations from the Kerr metric are allowed. In fact, the best-fit values are in the parametric region of the RBH-I phase, and the regions within 1-, 2- and 3-standard deviation limits are almost in the RBH-I region. Therefore, the observation of GRO J1655-40 favors a regular black hole with one horizon if the compact object of GRO J1655-40 is described by the rotating SV metric. The best-fit values of the radii of the circular orbits associated with the two sets of QPOs are found to be r¯1=5.669M=1.1304rISCO\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${{\bar{r}}}_{1} =5.669M=1.1304r_{\text {ISCO}}$$\end{document} and r¯2=5.563M=1.1094rISCO\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\bar{r}_{2}=5.563M=1.1094r_{\text {ISCO}}$$\end{document}, respectively, where rISCO=5.105M\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$r_{\text {ISCO}}=5.105M$$\end{document} is the innermost stable circular orbit evaluated for the rotating SV metric with the best-fit values (22). Consequently, the two circular orbits responsible for generating the two sets of QPOs lie in the close vicinity of the ISCO, and hence are in the strong-field region of the rotating SV metric.
Conclusions
In this paper, we explored potential deviations from the GR predictions of astrophysical black holes using QPOs observed in the power density spectrum of GRO J1655-40. Specially, we modelled the spacetime around the compact object of GRO J1655-40 by the rotating SV metric, and interpreted the observed QPOs within the RPM, which relates the QPO frequencies to epicyclic frequencies of geodesics. The rotating SV metric reduces to a Kerr black hole in the limit of ℓ=0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\ell =0$$\end{document}, and possesses multiple phases, e.g., a regular black hole with one or two horizons and a traversable wormhole. To test the nature of the the compact object of GRO J1655-40, we performed a χ2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\chi ^{2}$$\end{document} analysis by fitting the QPO frequencies computed in the RPM with the observations of two sets of QPOs from GRO J1655-40. Our results show a preference towards a regular black hole with one horizon compared to the Kerr black hole predicted in GR.
Acknowledgements
We thank Guangzhou Guo for his helpful discussions and suggestions. Houwen Wu is supported by the International Visiting Program for Excellent Young Scholars of Sichuan University. This work is supported in part by NSFC (Grant No. 11875196, 11375121, 11947225 and 11005016).
Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors’ comment: Data sharing not applicable to our article as no datasets were generated or analyzed during the current study.]
ReferencesAbbottBPObservation of gravitational waves from a binary black hole merger20161166061102370775810.1103/PhysRevLett.116.061102arXiv:1602.03837AkiyamaKFirst M87 event horizon telescope results. I. The shadow of the supermassive black hole2019875L110.3847/2041-8213/ab0ec7arXiv:1906.11238van der KlisMMillisecond oscillations in x-ray binaries20003871776010.1146/annurev.astro.38.1.717arXiv:astro-ph/0001167BertiETesting general relativity with present and future astrophysical observations20153224300110.1088/0264-9381/32/24/243001arXiv:1501.07274PsaltisDProbes and tests of strong-field gravity with observations in the electromagnetic spectrum200811910.12942/lrr-2008-91166.85302arXiv:0806.1531PsaltisDTwo approaches to testing general relativity in the strong-field regime200918901203310.1088/1742-6596/189/1/012033arXiv:0907.2746JusufiKAzreg-AïnouMJamilMWeiS-WQiangWWangAQuasinormal modes, quasiperiodic oscillations, and the shadow of rotating regular black holes in nonminimally coupled Einstein–Yang–Mills theory20211032024103421242910.1103/PhysRevD.103.024013arXiv:2008.08450Azreg-AïnouMChenZDengBJamilMZhuTQiangWLimY-KOrbital mechanics and quasiperiodic oscillation resonances of black holes in Einstein–Æther theory20201024044028414659210.1103/PhysRevD.102.044028arXiv:2004.02602IngramAMottaSA review of quasi-periodic oscillations from black hole X-ray binaries: observation and theory20198510152410.1016/j.newar.2020.101524arXiv:2001.08758MorganEHRemillardRAGreinerJRXTE observations of QPOs in the black hole candidate GRS 1915+1051997482993101010.1086/304191RemillardRAMorganEHMcClintockJEBailynCDOroszJARxte observations of 0.1–300 Hz quasi-periodic oscillations in the microquasar gro j1655–401999522139710.1086/307606StrohmayerTEDiscovery of a 450 hz quasi-periodic oscillation from the microquasar gro j1655–40 with the rossi x-ray timing explorer20015221L4910.1086/320258AbramowiczMAKluzniakWA Precise determination of angular momentum in the black hole candidate GRO J1655–402001374L1910.1051/0004-6361:20010791arXiv:astro-ph/0105077MottaSEBelloniTMStellaLMuñoz DariasTFenderRPrecise mass and spin measurements for a stellar-mass black hole through X-ray timing: the case of GRO J1655–40201443732554256510.1093/mnras/stt2068arXiv:1309.3652StellaLVietriMMorsinkSCorrelations in the qpo frequencies of low mass x-ray binaries and the relativistic precession model1999524L63L6610.1086/312291arXiv:astro-ph/9907346CadezACalvaniMKosticUOn the tidal evolution of the orbits of low-mass satellites around black holes200848752753210.1051/0004-6361:200809483arXiv:0809.1783RebuscoPTwin peaks kHz QPOs: mathematics of the 3:2 orbital resonance20045655310.1093/pasj/56.3.553arXiv:astro-ph/0403341W. Kluzniak, M.A. Abramowicz. Parametric epicyclic resonance in black hole disks: qpos in micro-quasars (2002). arXiv:astro-ph/0203314M.A. Abramowicz, W. Kluzniak, Z. Stuchlik, G. Torok, The orbital resonance model for twin peak kHz QPOs: measuring the black hole spins in microquasars (2004). arXiv:astro-ph/0401464KatoSFrequency correlations of QPOs based on a disk oscillation model in warped disks20075945110.1093/pasj/59.2.451arXiv:astro-ph/0701085BursaMAbramowiczMAKarasVKluzniakWThe upper kHz QPO: a gravitationally lensed vertical oscillation2004617L45L4810.1086/427167TorokGBakalaPSramkovaEStuchlikZUrbanecMOn mass-constraints implied by the relativistic precession model of twin-peak quasi-periodic oscillations in Circinus X-1201071474875710.1088/0004-637X/714/1/748arXiv:1008.0088BeerMEPodsiadlowskiPThe quiescent light curve and evolutionary state of gro j1655–40200233135110.1046/j.1365-8711.2002.05189.xarXiv:astro-ph/0109136JohannsenTPsaltisDTesting the no-hair theorem with observations in the electromagnetic spectrum. III. Quasi-periodic variability20117261110.1088/0004-637X/726/1/11arXiv:1010.1000KulkarniAKPennaRFShcherbakovRVSteinerJFNarayanRSadowskiAZhuYMcClintockJEDavisSWMcKinneyJCMeasuring black hole spin by the continuum-fitting method: effect of deviations from the Novikov–Thorne disc model2011414118310.1111/j.1365-2966.2011.18446.xarXiv:1102.0010PappasGWhat can quasi-periodic oscillations tell us about the structure of the corresponding compact objects?20124222581258910.1111/j.1365-2966.2012.20817.xarXiv:1201.6071BambiCProbing the space-time geometry around black hole candidates with the resonance models for high-frequency QPOs and comparison with the continuum-fitting method20120901410.1088/1475-7516/2012/09/014arXiv:1205.6348BambiCTesting the nature of the black hole candidate in GRO J1655–40 with the relativistic precession model201575416210.1140/epjc/s10052-015-3396-7arXiv:1312.2228MaselliAGualtieriLPaniPStellaLFerrariVTesting gravity with quasi periodic oscillations from accreting black holes: the case of the Einstein–Dilaton–Gauss–Bonnet theory2015801211510.1088/0004-637X/801/2/115arXiv:1412.3473SuvorovAGMelatosATesting modified gravity and no-hair relations for the Kerr-Newman metric through quasiperiodic oscillations of galactic microquasars201693024004348186310.1103/PhysRevD.93.024004arXiv:1512.02291BambiCNampalliwarSQuasi-periodic oscillations as a tool for testing the Kerr metric: a comparison with gravitational waves and iron line201611633000610.1209/0295-5075/116/30006arXiv:1604.02643ChenSWangMJingJTesting gravity of a regular and slowly rotating phantom black hole by quasi-periodic oscillations20163319195002355200010.1088/0264-9381/33/19/1950021349.83041arXiv:1604.07106A. Allahyari, L. Shao, Testing no-hair theorem by quasi-periodic oscillations: the quadrupole of GRO J1655-40 (2021). arXiv:2102.02232ChenSWangZJingJTesting gravity of a disformal Kerr black hole in quadratic degenerate higher-order scalar-tensor theories by quasi-periodic oscillations202106043432379610.1088/1475-7516/2021/06/043arXiv:2103.11788I. Banerjee, S. Chakraborty, S. SenGupta, Looking for extra dimensions in the observed quasi-periodic oscillations of black holes (2021). arXiv:2105.06636V. De Falco, M. De Laurentis, S. Capozziello, Epicyclic frequencies in static and spherically symmetric wormhole geometries (2021). arXiv:2106.12564K. Rink, I. Caiazzo, J. Heyl, Testing general relativity using quasi-periodic oscillations from X-ray black holes: XTE J1550-564 and GRO J1655-40 (2021). arXiv:2107.06828StellaLVietriMLense-Thirring precession and QPOS in low mass x-ray binaries1998492L5910.1086/311075arXiv:astro-ph/9709085J.M. Bardeen, Non-Singular General-Relativistic Gravitational Collapse, in Proceedings of the International Conference GR5, Tbilisi, USSR (1968), p 174RomanTABergmannPGStellar collapse without singularities?1983281265127771753710.1103/PhysRevD.28.1265HaywardSAFormation and evaporation of regular black holes20069603110310.1103/PhysRevLett.96.031103arXiv:gr-qc/0506126J.M. Bardeen, Black Hole Evaporation without an Event Horizon, vol. 6 (2014). arXiv:1406.4098FrolovVPNotes on nonsingular models of black holes20169410104056374781410.1103/PhysRevD.94.104056arXiv:1609.01758P.A. Cano, S. Chimento, T. Ortín, A. Ruipérez, Regular stringy black holes? Phys. Rev. D 99(4), 046014 (2019). https://doi.org/10.1103/PhysRevD.99.046014. arXiv:1806.08377J.M. Bardeen, Models for the nonsingular transition of an evaporating black hole into a white hole (2018). arXiv:1811.06683Carballo-RubioRDi FilippoFLiberatiSPacilioCVisserMOn the viability of regular black holes201807023386203010.1007/JHEP07(2018)0231395.83041arXiv:1805.02675Carballo-RubioRDi FilippoFLiberatiSVisserMPhenomenological aspects of black holes beyond general relativity20189812124009396970410.1103/PhysRevD.98.124009arXiv:1809.08238SimpsonAVisserMBlack-bounce to traversable wormhole201902042398389510.1088/1475-7516/2019/02/042arXiv:1812.07114SimpsonAMartin-MorunoPVisserMVaidya spacetimes, black-bounces, and traversable wormholes20193614145007398616510.1088/1361-6382/ab28a507416077arXiv:1902.04232A. Simpson, Traversable wormholes, regular black holes, and black-bounces. Master’s thesis, Victoria U., Wellington (2019). arXiv:2104.14055ChurilovaMSStuchlikZRinging of the regular black-hole/wormhole transition202037707514408673310.1088/1361-6382/ab771707418305arXiv:1911.11823ZhouT-YXieYPrecessing and periodic motions around a black-bounce/traversable wormhole20208011107010.1140/epjc/s10052-020-08661-wNascimentoJRPetrovAYPorfirioPJSoaresARGravitational lensing in black-bounce spacetimes20201024044021414620310.1103/PhysRevD.102.044021arXiv:2005.13096ChengX-TXieYProbing a black-bounce, traversable wormhole with weak deflection gravitational lensing20211036064040424752210.1103/PhysRevD.103.064040N. Tsukamoto, Gravitational lensing by two photon spheres in a black-bounce spacetime in strong deflection limits (2021). arXiv:2105.14336BronnikovKAKonoplyaRAPappasTDGeneral parametrization of wormhole spacetimes and its application to shadows and quasinormal modes202110312124062429004610.1103/PhysRevD.103.124062arXiv:2102.10679M. Guerrero, G.J. Olmo, D. Rubiera-Garcia, D. Sáez-Chillón Gómez, Shadows and optical appearance of black bounces illuminated by a thin accretion disk (2021). arXiv:2105.15073MazzaJFranzinELiberatiSA novel family of rotating black hole mimickers202104082429341110.1088/1475-7516/2021/04/082arXiv:2102.01105S. Ul Islam, J. Kumar, S.G. Ghosh, Strong gravitational lensing by rotating Simpson–Visser black holes (2021). arXiv:2104.00696R. Shaikh, K. Pal, K. Pal, T. Sarkar, Constraining alternatives to the Kerr black hole (2021). https://doi.org/10.1093/mnras/stab1779. arXiv:2102.04299S. Chandrasekhar, The mathematical theory of black holes (1985)NewmanETJanisAINote on the Kerr spinning particle metric1965691591717894610.1063/1.17043500142.46305C.D. Brown, Elements of spacecraft design. Aiaa (2002)OroszJABailynCDOptical observations of GRO J1655–40 in quiescence I: a precise mass for the black hole primary199747787610.1086/303741arXiv:astro-ph/9610211