HTML conversions sometimes display errors due to content that did not convert correctly from the source. This paper uses the following packages that are not yet supported by the HTML conversion tool. Feedback on these issues are not necessary; they are known and are being worked on.

  • failed: morefloats
  • failed: fonttable
  • failed: cellspace
  • failed: centernot

Authors: achieve the best HTML results from your LaTeX submissions by following these best practices.

License: CC BY 4.0
arXiv:2403.11885v1 [gr-qc] 18 Mar 2024

On dark matter effect on BH accretion disks

D. Pugliese&Z. Stuchlík Research Centre for Theoretical Physics and Astrophysics, Institute of Physics, Silesian University in Opava, Bezručovo náměstí 13, CZ-74601 Opava, Czech Republic
(March 18, 2024)
Abstract

Comparing different dark matter (DM) models, we explore the DM influence on the black hole (BH) accretion disk physics, considering co-rotating and counter-rotating thick accretion tori orbiting a central spinning BH. Our results point out accretion onto a central BH as a good indicator of the DM presence, signaling possible DM tracers in the accretion physics. We analyze accretion around a spinning BH immersed in perfect-fluid dark matter, cold dark matter and scalar field dark matter. Our investigation discusses observational evidence of distinctive DM effects on the toroidal accretion disks and proto-jets configurations, proving that BHs accretion tori, immersed in DM, can present characteristics, such as inter–disks cusp or double tori, which have been usually considered as tracers for super-spinars and naked singularity attractors. Therefore in this context DM influence on the BH geometry could manifest as super-spinars mimickers. DM affects also the central spinning attractor energetics associated with the accretion physics, and its influence on the accretion disks can be searched in a variation of the central BH energetics as an increase of the mass accretion rates.

Dark Matter–Black holes– Accretion disks–Accretion; Hydrodynamics –Galaxies: actives

I Introduction

In this work we study the dark matter (DM) influence on the black hole (BH) accretion disk physics, investigating the accretion disks morphology for co-rotating and counter-rotating geometrically thick accretion tori orbiting a central spinning BH. Comparing three different DM models, our analysis points out possible observational evidence of distinctive DM effects on the accretion disks, which can be traces for the DM presence. Our investigation does not cover all the admissible parametric DM values for the deformed metrics but, taking into account constraints formerly obtained by the study of orbits, the DM BH shadows, or the emission spectra, we perform a comparative analysis of the DM models and an investigation on the constraints imposed on the accretion discs, aimed at further restricting the parameters range and pointing out the possible DM mark in some accretion features. From a methodological viewpoint we take advantage of axial symmetry of DM metrics, studying fully general relativistic models of stationary toroidal orbiting configurations.

The physics of accretion disks around BHs and super massive black hole (SMBH), hosted in quasars and active galactic nuclei (AGN), powers the most energetic processes of our Universe, often accompanied with ejection of matter in jet-like structures with extremely large radiative energy output. We investigate the DM effects on these aspects, focusing on SMBH at the center of galaxies and the accretion discs empowering the emissions. We analyze particularly the limiting situation of static (and spherically symmetric) background (with spin parameter a=0𝑎0a=0italic_a = 0) and the DM deformation on the Kerr extreme BH spacetime (with spin value a=M𝑎𝑀a=Mitalic_a = italic_M), seeing significant qualitative detectable variations with respect to the standard vacuum BH case. Many astrophysical observations lead to SMBH hosted at the galactic center embedded in a DM halo, and in metric models considered here the central BH is surrounded by DM envelope that modifies the geometry around the BH, not interacting directly with the accreting matter or its radiation.

More generally there is a large amount of observational evidence of the presence of some kind of DM component in our Universe, for example from the galactic rotation curves and galaxy cluster dynamics. However, within the variety of the different observations pointing at the DM presence, there is no single metric that encompasses all the DM effects in a single model, and which could also explain the absence of DM observed at different scales111An important issue is then the mass (or space) scale when DM effects became significant. For example DM results missing in some galaxies (ex. AGC 114905 –[1, 2]), and an explanation for this situation is that the galaxy may have been stripped of dark matter from nearby massive galaxies, while the DM presence in the Solar system is still an open problem[3].. In this work we consider a spinning BH immersed in perfect fluid DM (PFDM) [4, 5, 6, 7], cold dark matter (CDM) and scalar field dark matter (SFDM)–[8, 9].

These DM models have been extensively studied in recent literature. PFDM model was considered in [5] and the effects of PFDM on particle motion around a static BH in an external magnetic field were studied in [10]. The shadow of the spinning BH in PFDM was studied in [4], whereas in [6] geodesic motion in PFDM Kerr and Kerr-anti-de Sitter/de Sitter BH were studied–see also [7]. For studies of the geodesic in the Kerr-de Sitter spacetimes see [11, 12, 13, 14, 15]. In [16] superradiance and stability of Kerr DM enclosed by anisotropic fluid matter were studied. Spinning BH solutions with quintessential energy have been discussed in [17]. In [8] BH in DM halo was considered. Rotating black holes with an anisotropic matter field was considered in [18], while in [9] there is a discussion of the BH shadow of Sgr A* in DM halo. Superradiance and instabilities in BHs surrounded by anisotropic fluids were considered in [19]. Galactic dark matter in the phantom field model was considered in [20]. The case of rotating (Kerr) naked singularities was treated in [21, 22, 23, 24]. Here we focus on the influence of DM on BHs governing toroidal accretion structures.

There is an extensive literature exploring the DM effects on the BH and BH accretion physics. Since DM influences different aspects of the singularity, from the characteristics of the horizon (for example it can manifest itself in the BH shadows) to the energetics properties of the surrounding matter, there are various assessments of the DM effects and parameter constraints on the models describing DM presence around BHs. In [26] for example DM clouds of axions around BHs were studied with superradiant instabilities and accretion which could manifest on the gravitational wave signal induced by a small compact object in the field of the central BH–see also [27, 28]. More recently in [29] the formation of SMBHs at high redshifts was studied in connection with ultralight DM, see also [30, 31] for growth of accretion driven scalar DM hair around a Kerr BHs. The DM effect on the quasinormal modes of massless scalar field and electromagnetic field perturbations in a BH spacetime surrounded by PFDM was considered in [32]. An analysis of DM in the M87 core in relation to BH shadows effects was presented in [33], see also [35] and shadows of Sgr A* BH surrounded by superfluid DM halo was studied in [34], and shadow from a charged rotating BH in presence of PFDM is explored in [36].

The DM candidates are various, including string and brane theory effects, boson clouds, hypothetical new particles, primordial BHs and alternative theories of gravity222Dark matter was also explained with diffuse clouds of scalar bosons interacting with gravity and with gravitational waves. However recent results of [25] sets constrains to this hypothesis showing that there are no young scalar boson clouds in our galaxy.. From the observational view-point, dark matter can be detectable from the products of its decay or annihilation in cosmic rays, gamma rays, neutrinos or even gravitons–see also [37], and gravitational–wave and neutrino astronomy can then open different windows in the DM analysis333Concerning possible DM constituents and presence in our Galaxy we mention that the DM Milky Way has been recently studied in [38] and [39]. Whereas, DM and primordial BHs are studied in [40, 41] constraints on DM rule out BHs as constituting only a very small possible fraction of the dark matter. Finally hypothesis suggesting antimatter and DM be linked have been recently studied in [42], posing limits on the interaction of antiprotons with axion-like DM, see also [43].. DM comprehension, particularly focused on sub-galactic DM halos, is also a goal of the Webb Telescope 444https://jwst.nasa.gov..

Nevertheless, despite the variety of DM models, the standard cosmological model is in fact the ΛΛ\Lambdaroman_Λ-CDM, which includes a cosmological constant (ΛΛ\Lambdaroman_Λ) (with negative pressure), encoding Dark energy (DE) in empty space (or vacuum energy) explaining the Universe accelerating expansion. (In this scenario the effects of the cosmological constant are also treated as quintessence555See [44] for a recent analysis constraining the fraction of early dark energy, present during the early ages of the Universe..) Polytropic models of DM halos in ΛΛ\Lambdaroman_Λ-CDM cosmology were individuated in [45]. In this model, the DM velocity is less than the speed of light (in this respect neutrinos component are excluded, being non-baryonic but not necessarily cold) and it is dissipationless as it is not cooled by radiating photons. CDM may be constituted by an hypothetical weakly interacting massive particle (WIMPS), or primordial BHs, or axions.

Although considered a DM standard model, CDM is not exempt from various problems, emerging for example from the observations of galaxies and galaxy clusters and clusterization emerging from the rotation curves and morphological studies (as the cuspy halo problem). There is also a more general problem in describing the effects and presence of DM at large and small scales. The CDM model collides therefore with small-scale structure observations. For all these reasons the search for alternative DM models is still an open issue, and in this respect the SFDM model seems to adapt to both large-scale and small-scale structure observations while the PFDM seems capable to explain the asymptotically flat rotation velocity characterizing the spiral galaxies.

In this analysis we study geometrically thick accretion disk models, Polish Doughnuts (P-D), orbiting the central attractor, whose center coincides with the equatorial plane of the central axisymmetric attractor[62, 47, 48, 49, 50]. These thick accretion tori are characterized by very high (super-Eddington) accretion rates and high optical depth. Tori morphology and stability are essentially governed by the pressure gradients on the equatorial plane [47]. The thin (Keplerian) disks can be considered as P-D limiting configurations regulated by the background geodesic structure. The DM background metric has a characteristic geodetic structure constituting a first major constraint on accretion physics. The tori, described by purely hydrodynamic (barotropic) models, are governed by the equipressure surfaces that can be closed, giving stable equilibrium configurations, and open, giving unstable, jet-like (proto-jets) structures caused by the relativistic instability due to the Paczynski mechanism where the effects of strong gravitational fields are dominant666The time scale of the dynamical processes (regulated by the gravitational and inertial forces) is much lower than the time scale of the thermal ones (heating and cooling processes, radiation) that is lower than the time scale of the viscous processes. The entropy is constant along the flow and, according to the von Zeipel condition, the surfaces of constant angular velocity ΩΩ\Omegaroman_Ω and of constant specific angular momentum \ellroman_ℓ coincide. This implies that the rotation law =(Ω)Ω\ell=\ell(\Omega)roman_ℓ = roman_ℓ ( roman_Ω ) is independent of the equation of state [55, 56]. with respect to the dissipative ones and predominant to determine the unstable phases of the systems [57, 51, 58, 59, 60, 53, 54, 61].

Many features of the tori dynamics and morphology like their thickness, their stretching in the equatorial plane, and the location of the tori are predominantly determined by the geometric properties of spacetime via a fluid effective potential function. Consequently, in models where DM is geometrized as a metric deformation, DM has a clear impact in the tori structure, modifying the fluid effective potential. The gradients of the effective potential on the tori equatorial and symmetry plane regulate the pressure gradient of the fluid in the Euler law governing dynamics of the perfect fluid [62]. The special case of cusped equipotential surfaces is related to the accretion phase onto the central attractor [62, 47, 48, 60, 63]. The outflow of matter through the cusp occurs due to an instability in the balance of the gravitational and inertial forces and the pressure gradients in the fluid, i.e., by the so called Paczynski mechanism of violation of mechanical equilibrium of the tori [48].

DM affects the cusp formation and cusp location with respect to the central singularity, modifying the disk accretion throat, constraining the thickness of the accretionary flow and the maximum amount of matter swallowed by the central BH. Consequently DM will influence the energetic characteristics of the BH in accretion and the disk characteristics, as accretion rates or cusp luminosity [64, 65, 66].

More in details plan of the article is as follows: Thick disks in axially symmetric spacetimes are discussed in Sec. (II). The Kerr metric is introduced in Sec. (II.1). The Polish doughnut tori models are detailed in Sec. (II.2). The fluid effective potential is the subject of Sec. (II.2.1). Extended geodesic structure, constrained the tori modes is explored in Sec. (II.2.2). Dark matter models are discussed in Sec. (III). In Sec. (III.1) the perfect fluid dark matter is considered. Cold and scalar field dark matter models are studied Sec. (III.2). Discussion and conclusions follow in Sec. (IV).

II Thick disks in axially symmetric spacetimes

We study geometrically thick tori in axially symmetric DM–BH spacetimes, considered as a DM-induced deformation of the Kerr geometry. Therefore it is useful here to review the properties of the Kerr metric and the construction of tori in this geometry. More specifically, in Sec. (II.1) the Kerr metric is introduced, while the Polish doughnut tori models are discussed in Sec. (II.2).

II.1 The Kerr metric

The Kerr metric is an axially symmetric, asymptotically flat, vacuum exact solution of the Einstein equation describing the spacetime of central spinning compact object. According to the metric parameter values (dimensionless spin a/M𝑎𝑀a/Mitalic_a / italic_M), the Kerr metric describes naked singularities (NSs) for a>M𝑎𝑀a>Mitalic_a > italic_M and black holes (BHs) for a[0,M]𝑎0𝑀a\in[0,M]italic_a ∈ [ 0 , italic_M ]. The Kerr BH geometry has the limiting static solution of Schwarzschild for a=0𝑎0a=0italic_a = 0 and the extreme Kerr BH spacetime for a=M𝑎𝑀a=Mitalic_a = italic_M.

In the Boyer-Lindquist (BL) coordinates {t,r,θ,ϕ}𝑡𝑟𝜃italic-ϕ\{t,r,\theta,\phi\}{ italic_t , italic_r , italic_θ , italic_ϕ }, the metric tensor reads777We adopt the geometrical units c=1=G𝑐1𝐺c=1=Gitalic_c = 1 = italic_G and the (,+,+,+)(-,+,+,+)( - , + , + , + ) signature, Latin indices run in {0,1,2,3}0123\{0,1,2,3\}{ 0 , 1 , 2 , 3 }. The radius r𝑟ritalic_r has unit of mass [M]delimited-[]𝑀[M][ italic_M ], and the angular momentum units of [M]2superscriptdelimited-[]𝑀2[M]^{2}[ italic_M ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, the velocities [ut]=[ur]=1delimited-[]superscript𝑢𝑡delimited-[]superscript𝑢𝑟1[u^{t}]=[u^{r}]=1[ italic_u start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ] = [ italic_u start_POSTSUPERSCRIPT italic_r end_POSTSUPERSCRIPT ] = 1 and [uϕ]=[uθ]=[M]1delimited-[]superscript𝑢italic-ϕdelimited-[]superscript𝑢𝜃superscriptdelimited-[]𝑀1[u^{\phi}]=[u^{\theta}]=[M]^{-1}[ italic_u start_POSTSUPERSCRIPT italic_ϕ end_POSTSUPERSCRIPT ] = [ italic_u start_POSTSUPERSCRIPT italic_θ end_POSTSUPERSCRIPT ] = [ italic_M ] start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT with [uϕ/ut]=[M]1delimited-[]superscript𝑢italic-ϕsuperscript𝑢𝑡superscriptdelimited-[]𝑀1[u^{\phi}/u^{t}]=[M]^{-1}[ italic_u start_POSTSUPERSCRIPT italic_ϕ end_POSTSUPERSCRIPT / italic_u start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ] = [ italic_M ] start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT and [uϕ/ut]=[M]delimited-[]subscript𝑢italic-ϕsubscript𝑢𝑡delimited-[]𝑀[u_{\phi}/u_{t}]=[M][ italic_u start_POSTSUBSCRIPT italic_ϕ end_POSTSUBSCRIPT / italic_u start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ] = [ italic_M ]. For the seek of convenience, we always consider the dimensionless energy and effective potential [Veff]=1delimited-[]subscript𝑉𝑒𝑓𝑓1[V_{eff}]=1[ italic_V start_POSTSUBSCRIPT italic_e italic_f italic_f end_POSTSUBSCRIPT ] = 1 and an angular momentum per unit of mass [L]/[M]=[M]delimited-[]𝐿delimited-[]𝑀delimited-[]𝑀[L]/[M]=[M][ italic_L ] / [ italic_M ] = [ italic_M ].:

ds2=(12MrΣ)dt2+ΣΔdr2+Σdθ2+[(r2+a2)+2Mra2Σsin2θ]sin2θdϕ24rMaΣsin2θdtdϕ,𝑑superscript𝑠212𝑀𝑟Σ𝑑superscript𝑡2ΣΔ𝑑superscript𝑟2Σ𝑑superscript𝜃2delimited-[]superscript𝑟2superscript𝑎22𝑀𝑟superscript𝑎2Σsuperscript2𝜃superscript2𝜃𝑑superscriptitalic-ϕ24𝑟𝑀𝑎Σsuperscript2𝜃𝑑𝑡𝑑italic-ϕ\displaystyle ds^{2}=-\left(1-\frac{2Mr}{\Sigma}\right)dt^{2}+\frac{\Sigma}{% \Delta}dr^{2}+\Sigma d\theta^{2}+\left[(r^{2}+a^{2})+\frac{2Mra^{2}}{\Sigma}% \sin^{2}\theta\right]\sin^{2}\theta d\phi^{2}-\frac{4rMa}{\Sigma}\sin^{2}% \theta dtd\phi,italic_d italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = - ( 1 - divide start_ARG 2 italic_M italic_r end_ARG start_ARG roman_Σ end_ARG ) italic_d italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + divide start_ARG roman_Σ end_ARG start_ARG roman_Δ end_ARG italic_d italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + roman_Σ italic_d italic_θ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + [ ( italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) + divide start_ARG 2 italic_M italic_r italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG roman_Σ end_ARG roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ ] roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ italic_d italic_ϕ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - divide start_ARG 4 italic_r italic_M italic_a end_ARG start_ARG roman_Σ end_ARG roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ italic_d italic_t italic_d italic_ϕ , (1)

where

Δa2+r22rM;Σa2(1σ)+r2,σsin2θ,formulae-sequenceΔsuperscript𝑎2superscript𝑟22𝑟𝑀formulae-sequenceΣsuperscript𝑎21𝜎superscript𝑟2𝜎superscript2𝜃\displaystyle\Delta\equiv a^{2}+r^{2}-2rM;\quad\Sigma\equiv a^{2}(1-\sigma)+r^% {2},\quad\sigma\equiv\sin^{2}\theta,roman_Δ ≡ italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 2 italic_r italic_M ; roman_Σ ≡ italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 1 - italic_σ ) + italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , italic_σ ≡ roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ , (2)

with G=c=1𝐺𝑐1G=c=1italic_G = italic_c = 1. The horizons r<r+subscript𝑟subscript𝑟r_{-}<r_{+}italic_r start_POSTSUBSCRIPT - end_POSTSUBSCRIPT < italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT are respectively given by

r±M±M2a2,subscript𝑟plus-or-minusplus-or-minus𝑀superscript𝑀2superscript𝑎2\displaystyle r_{\pm}\equiv M\pm\sqrt{M^{2}-a^{2}},italic_r start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ≡ italic_M ± square-root start_ARG italic_M start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , (3)

the horizons can be found888Quantities a±subscript𝑎plus-or-minusa_{\pm}italic_a start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT and aϵ±superscriptsubscript𝑎italic-ϵplus-or-minusa_{\epsilon}^{\pm}italic_a start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT turn to be very useful in the comparison of the dark matter (DM) solutions of Sec. (III) with respect to Kerr solutions in absence of DM. by solving the equation a=a±r(2Mr)𝑎subscript𝑎plus-or-minus𝑟2𝑀𝑟a=a_{\pm}\equiv\sqrt{r(2M-r)}italic_a = italic_a start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ≡ square-root start_ARG italic_r ( 2 italic_M - italic_r ) end_ARG for r[0,2M]𝑟02𝑀r\in[0,2M]italic_r ∈ [ 0 , 2 italic_M ] . The outer and inner stationary limits rϵ±superscriptsubscript𝑟italic-ϵplus-or-minusr_{\epsilon}^{\pm}italic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT (ergosurfaces) (solutions of gtt=0subscript𝑔𝑡𝑡0g_{tt}=0italic_g start_POSTSUBSCRIPT italic_t italic_t end_POSTSUBSCRIPT = 0) are respectively

rϵ±M±M2a2(1σ),superscriptsubscript𝑟italic-ϵplus-or-minusplus-or-minus𝑀superscript𝑀2superscript𝑎21𝜎\displaystyle r_{\epsilon}^{\pm}\equiv M\pm\sqrt{M^{2}-a^{2}(1-\sigma)},italic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ≡ italic_M ± square-root start_ARG italic_M start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 1 - italic_σ ) end_ARG , (4)

the ergosurfaces can be found by solving the equation a=aϵ±a±/1σ𝑎superscriptsubscript𝑎italic-ϵplus-or-minussubscript𝑎plus-or-minus1𝜎a=a_{\epsilon}^{\pm}\equiv{a_{\pm}}/{\sqrt{1-\sigma}}italic_a = italic_a start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ≡ italic_a start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT / square-root start_ARG 1 - italic_σ end_ARG (for σ1𝜎1\sigma\neq 1italic_σ ≠ 1), where r+<rϵ+subscript𝑟superscriptsubscript𝑟italic-ϵr_{+}<r_{\epsilon}^{+}italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT < italic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT on θ0𝜃0\theta\neq 0italic_θ ≠ 0 and rϵ+=2Msuperscriptsubscript𝑟italic-ϵ2𝑀r_{\epsilon}^{+}=2Mitalic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT = 2 italic_M in the equatorial plane θ=π/2𝜃𝜋2\theta=\pi/2italic_θ = italic_π / 2 (σ=1𝜎1\sigma=1italic_σ = 1). Static observers, with four-velocity θ˙=r˙=ϕ˙=0˙𝜃˙𝑟˙italic-ϕ0\dot{\theta}=\dot{r}=\dot{\phi}=0over˙ start_ARG italic_θ end_ARG = over˙ start_ARG italic_r end_ARG = over˙ start_ARG italic_ϕ end_ARG = 0(where q˙˙𝑞\dot{q}over˙ start_ARG italic_q end_ARG indicates the derivative of any quantity q𝑞qitalic_q with respect the proper time (for time-like particles) or a properly defined affine parameter for the light-like orbits) cannot exist inside the (outer) ergoregion999The ergoregion is the range [rϵ,rϵ+]superscriptsubscript𝑟italic-ϵsuperscriptsubscript𝑟italic-ϵ[r_{\epsilon}^{-},r_{\epsilon}^{+}][ italic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT , italic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT ] (where rϵ±superscriptsubscript𝑟italic-ϵplus-or-minusr_{\epsilon}^{\pm}italic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT are functions of the plane σ[0,1]𝜎01\sigma\in[0,1]italic_σ ∈ [ 0 , 1 ]). Here we often intend the outer ergoregion (or simply ergoregion) in the BH spacetimes as the region ]r+,rϵ+]]r_{+},r_{\epsilon}^{+}]] italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT , italic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT ]. Then on the equatorial planes, in the Kerr spacetime there is rϵ=0superscriptsubscript𝑟italic-ϵ0r_{\epsilon}^{-}=0italic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT = 0 and the outer ergosurface is rϵ+=2Msuperscriptsubscript𝑟italic-ϵ2𝑀r_{\epsilon}^{+}=2Mitalic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT = 2 italic_M. , but trajectories r˙0˙𝑟0\dot{r}\geq 0over˙ start_ARG italic_r end_ARG ≥ 0, including particles crossing the stationary limit and escaping outside in the region rrϵ+𝑟superscriptsubscript𝑟italic-ϵr\geq r_{\epsilon}^{+}italic_r ≥ italic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT are possible.

The constants of the geodesic motions are

=(gtϕϕ˙+gttt˙),=gϕϕϕ˙+gtϕt˙,gabuaub=μ2,formulae-sequencesubscript𝑔𝑡italic-ϕ˙italic-ϕsubscript𝑔𝑡𝑡˙𝑡formulae-sequencesubscript𝑔italic-ϕitalic-ϕ˙italic-ϕsubscript𝑔𝑡italic-ϕ˙𝑡subscript𝑔𝑎𝑏superscript𝑢𝑎superscript𝑢𝑏superscript𝜇2\displaystyle\mathcal{E}=-(g_{t\phi}\dot{\phi}+g_{tt}\dot{t}),\quad\mathcal{L}% =g_{\phi\phi}\dot{\phi}+g_{t\phi}\dot{t},\quad g_{ab}u^{a}u^{b}=-\mu^{2},caligraphic_E = - ( italic_g start_POSTSUBSCRIPT italic_t italic_ϕ end_POSTSUBSCRIPT over˙ start_ARG italic_ϕ end_ARG + italic_g start_POSTSUBSCRIPT italic_t italic_t end_POSTSUBSCRIPT over˙ start_ARG italic_t end_ARG ) , caligraphic_L = italic_g start_POSTSUBSCRIPT italic_ϕ italic_ϕ end_POSTSUBSCRIPT over˙ start_ARG italic_ϕ end_ARG + italic_g start_POSTSUBSCRIPT italic_t italic_ϕ end_POSTSUBSCRIPT over˙ start_ARG italic_t end_ARG , italic_g start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT italic_u start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT italic_u start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT = - italic_μ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , (5)

with101010The other constant of geodesic motion of the Kerr metric is the Carter constant 𝒬=(cosθ)2[a2(μ22)+(sinθ)2]+(gθθθ˙)2𝒬superscript𝜃2delimited-[]superscript𝑎2superscript𝜇2superscript2superscript𝜃2superscriptsubscript𝑔𝜃𝜃˙𝜃2\mathcal{Q}=(\cos\theta)^{2}\left[a^{2}\left(\mu^{2}-\mathcal{E}^{2}\right)+% \left(\frac{\mathcal{L}}{\sin\theta}\right)^{2}\right]+(g_{\theta\theta}\dot{% \theta})^{2}caligraphic_Q = ( roman_cos italic_θ ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [ italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_μ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - caligraphic_E start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) + ( divide start_ARG caligraphic_L end_ARG start_ARG roman_sin italic_θ end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] + ( italic_g start_POSTSUBSCRIPT italic_θ italic_θ end_POSTSUBSCRIPT over˙ start_ARG italic_θ end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT. In this work, where tori share symmetry plane with the equatorial plane of the central BH, this constant is irrelevant. ua{t˙,r˙,θ˙,ϕ˙}superscript𝑢𝑎˙𝑡˙𝑟˙𝜃˙italic-ϕu^{a}\equiv\{\dot{t},\dot{r},\dot{\theta},\dot{\phi}\}italic_u start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT ≡ { over˙ start_ARG italic_t end_ARG , over˙ start_ARG italic_r end_ARG , over˙ start_ARG italic_θ end_ARG , over˙ start_ARG italic_ϕ end_ARG }. In Eqs (5) quantities \mathcal{E}caligraphic_E and \mathcal{L}caligraphic_L represent the total energy and momentum of the test particle coming from radial infinity, as measured by a static observer at infinity.

The relativistic angular velocity and the specific angular momentum are

Ωuϕut=gϕt+gttgϕϕ+gϕt=gtϕ+gttgϕϕ+gtϕ,=uϕut=gϕϕuϕ+gϕtutgttut+gϕtuϕ=gtϕ+gϕϕΩgtt+gtϕΩ,formulae-sequenceΩsuperscript𝑢italic-ϕsuperscript𝑢𝑡subscript𝑔italic-ϕ𝑡subscript𝑔𝑡𝑡subscript𝑔italic-ϕitalic-ϕsubscript𝑔italic-ϕ𝑡subscript𝑔𝑡italic-ϕsubscript𝑔𝑡𝑡subscript𝑔italic-ϕitalic-ϕsubscript𝑔𝑡italic-ϕsubscript𝑢italic-ϕsubscript𝑢𝑡subscript𝑔italic-ϕitalic-ϕsuperscript𝑢italic-ϕsubscript𝑔italic-ϕ𝑡superscript𝑢𝑡subscript𝑔𝑡𝑡superscript𝑢𝑡subscript𝑔italic-ϕ𝑡superscript𝑢italic-ϕsubscript𝑔𝑡italic-ϕsubscript𝑔italic-ϕitalic-ϕΩsubscript𝑔𝑡𝑡subscript𝑔𝑡italic-ϕΩ\displaystyle\Omega\equiv\frac{u^{\phi}}{u^{t}}=-\frac{\mathcal{E}g_{\phi t}+g% _{tt}\mathcal{L}}{\mathcal{E}g_{\phi\phi}+g_{\phi t}\mathcal{L}}=-\frac{g_{t% \phi}+g_{tt}\ell}{g_{\phi\phi}+g_{t\phi}\ell},\quad\ell\equiv\frac{\mathcal{L}% }{\mathcal{E}}=-\frac{u_{\phi}}{u_{t}}=-\frac{g_{\phi\phi}u^{\phi}+g_{\phi t}u% ^{t}}{g_{tt}u^{t}+g_{\phi t}u^{\phi}}=-\frac{g_{t\phi}+g_{\phi\phi}\Omega}{g_{% tt}+g_{t\phi}\Omega},roman_Ω ≡ divide start_ARG italic_u start_POSTSUPERSCRIPT italic_ϕ end_POSTSUPERSCRIPT end_ARG start_ARG italic_u start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT end_ARG = - divide start_ARG caligraphic_E italic_g start_POSTSUBSCRIPT italic_ϕ italic_t end_POSTSUBSCRIPT + italic_g start_POSTSUBSCRIPT italic_t italic_t end_POSTSUBSCRIPT caligraphic_L end_ARG start_ARG caligraphic_E italic_g start_POSTSUBSCRIPT italic_ϕ italic_ϕ end_POSTSUBSCRIPT + italic_g start_POSTSUBSCRIPT italic_ϕ italic_t end_POSTSUBSCRIPT caligraphic_L end_ARG = - divide start_ARG italic_g start_POSTSUBSCRIPT italic_t italic_ϕ end_POSTSUBSCRIPT + italic_g start_POSTSUBSCRIPT italic_t italic_t end_POSTSUBSCRIPT roman_ℓ end_ARG start_ARG italic_g start_POSTSUBSCRIPT italic_ϕ italic_ϕ end_POSTSUBSCRIPT + italic_g start_POSTSUBSCRIPT italic_t italic_ϕ end_POSTSUBSCRIPT roman_ℓ end_ARG , roman_ℓ ≡ divide start_ARG caligraphic_L end_ARG start_ARG caligraphic_E end_ARG = - divide start_ARG italic_u start_POSTSUBSCRIPT italic_ϕ end_POSTSUBSCRIPT end_ARG start_ARG italic_u start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG = - divide start_ARG italic_g start_POSTSUBSCRIPT italic_ϕ italic_ϕ end_POSTSUBSCRIPT italic_u start_POSTSUPERSCRIPT italic_ϕ end_POSTSUPERSCRIPT + italic_g start_POSTSUBSCRIPT italic_ϕ italic_t end_POSTSUBSCRIPT italic_u start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT end_ARG start_ARG italic_g start_POSTSUBSCRIPT italic_t italic_t end_POSTSUBSCRIPT italic_u start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT + italic_g start_POSTSUBSCRIPT italic_ϕ italic_t end_POSTSUBSCRIPT italic_u start_POSTSUPERSCRIPT italic_ϕ end_POSTSUPERSCRIPT end_ARG = - divide start_ARG italic_g start_POSTSUBSCRIPT italic_t italic_ϕ end_POSTSUBSCRIPT + italic_g start_POSTSUBSCRIPT italic_ϕ italic_ϕ end_POSTSUBSCRIPT roman_Ω end_ARG start_ARG italic_g start_POSTSUBSCRIPT italic_t italic_t end_POSTSUBSCRIPT + italic_g start_POSTSUBSCRIPT italic_t italic_ϕ end_POSTSUBSCRIPT roman_Ω end_ARG , (6)

respectively. The sign of ()\mathcal{L}(\ell)caligraphic_L ( roman_ℓ ) defines the co-rotation/counter-rotation of the particles (fluid). The DM models are axis-symmetric and stationary and we define similarly notion of co-rotating and counter-rotating motions.

II.2 Geometrically thick tori: the Polish doughnut models

We specialize our analysis to the Polish doughnut (P-D) tori, general relativistic hydrodynamic (GRHD) toroidal configurations centered on the central BH equatorial plane, which is coincident with the tori equatorial symmetry plane.

These toroidal models are well known and used in different contexts. They are analytic and general relativistic models defined and integrable in axis symmetric spacetimes, where the results known as the "von Zeipel theorem" hold, ensuring the integration condition on the equations for the fluid. For this reason we here apply these results to the stationary DM metric models111111The toroids are constant pressure surfaces, whose construction in the axis–symmetric spacetimes is based on the application of the von Zeipel theorem, for which the surfaces of constant angular velocity ΩΩ\Omegaroman_Ω and of constant specific angular momentum \ellroman_ℓ coincide and the toroids rotation law =(Ω)Ω\ell=\ell(\Omega)roman_ℓ = roman_ℓ ( roman_Ω ) is independent of the details of the equation of state. More precisely, the von Zeipel theorem reduces to an integrability condition on the Euler equation, in the case of barotropic fluids, where =(Ω)Ω\ell=\ell(\Omega)roman_ℓ = roman_ℓ ( roman_Ω ) and, consequently, in the geometrically thick disks the functional form of the angular momentum and entropy distribution, during the evolution of dynamical processes, depends on the initial conditions of the system and not on the details of the dissipative processes[58]..

Tori are composed by a one particle-species perfect fluid, where

Tab=(ϱ+p)uaub+pgab,subscript𝑇𝑎𝑏italic-ϱ𝑝subscript𝑢𝑎subscript𝑢𝑏𝑝subscript𝑔𝑎𝑏T_{ab}=(\varrho+p)u_{a}u_{b}+pg_{ab},italic_T start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT = ( italic_ϱ + italic_p ) italic_u start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT + italic_p italic_g start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT , (7)

is the fluid energy momentum tensor, ϱitalic-ϱ\varrhoitalic_ϱ and p𝑝pitalic_p are the total energy density and pressure, respectively, as measured by an observer moving with the fluid. The timelike flow vector field uasuperscript𝑢𝑎u^{a}italic_u start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT denotes the fluid four-velocity. The fluid dynamics is described by the continuity equation and the Euler equation respectively:

uaaϱ+(p+ϱ)aua=0,(p+ϱ)uaauc+hbcbp=0,formulae-sequencesuperscript𝑢𝑎subscript𝑎italic-ϱ𝑝italic-ϱsuperscript𝑎subscript𝑢𝑎0𝑝italic-ϱsuperscript𝑢𝑎subscript𝑎superscript𝑢𝑐superscript𝑏𝑐subscript𝑏𝑝0\displaystyle u^{a}\nabla_{a}\varrho+(p+\varrho)\nabla^{a}u_{a}=0,\quad(p+% \varrho)u^{a}\nabla_{a}u^{c}+\ h^{bc}\nabla_{b}p=0,italic_u start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT ∇ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_ϱ + ( italic_p + italic_ϱ ) ∇ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT italic_u start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = 0 , ( italic_p + italic_ϱ ) italic_u start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT ∇ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_u start_POSTSUPERSCRIPT italic_c end_POSTSUPERSCRIPT + italic_h start_POSTSUPERSCRIPT italic_b italic_c end_POSTSUPERSCRIPT ∇ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_p = 0 , (8)

where the projection tensor hab=gab+uaubsubscript𝑎𝑏subscript𝑔𝑎𝑏subscript𝑢𝑎subscript𝑢𝑏h_{ab}=g_{ab}+u_{a}u_{b}italic_h start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT = italic_g start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT + italic_u start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT and agbc=0subscript𝑎subscript𝑔𝑏𝑐0\nabla_{a}g_{bc}=0∇ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_g start_POSTSUBSCRIPT italic_b italic_c end_POSTSUBSCRIPT = 0.

We assume a barotropic equation of state (EoS) p=p(ϱ)𝑝𝑝italic-ϱp=p(\varrho)italic_p = italic_p ( italic_ϱ ) and the stationary and axially symmetric matter distribution moves on circular trajectories. We investigate the case of a fluid toroidal configuration defined by the constraint ur=0superscript𝑢𝑟0u^{r}=0italic_u start_POSTSUPERSCRIPT italic_r end_POSTSUPERSCRIPT = 0, as for the circular test particle motion no motion it is assumed in the θ𝜃\thetaitalic_θ angular direction, which means uθ=0superscript𝑢𝜃0u^{\theta}=0italic_u start_POSTSUPERSCRIPT italic_θ end_POSTSUPERSCRIPT = 0. Because of these symmetries, the continuity equation is identically satisfied and the orbiting configurations are regulated by the Euler equation for the pressure p𝑝pitalic_p only, which can be written as

apϱ+p=aW+Ωa1Ω,withWlnVeff,andVeff=ut,formulae-sequencesubscript𝑎𝑝italic-ϱ𝑝subscript𝑎𝑊Ωsubscript𝑎1Ωwithformulae-sequence𝑊subscript𝑉𝑒𝑓𝑓andsubscript𝑉𝑒𝑓𝑓subscript𝑢𝑡\displaystyle\frac{\partial_{a}p}{\varrho+p}=-{\partial_{a}}W+\frac{\Omega% \partial_{a}\ell}{1-\Omega\ell},\quad\mbox{with}\quad W\equiv\ln V_{eff},\quad% \mbox{and}\quad V_{eff}=u_{t},divide start_ARG ∂ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_p end_ARG start_ARG italic_ϱ + italic_p end_ARG = - ∂ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_W + divide start_ARG roman_Ω ∂ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT roman_ℓ end_ARG start_ARG 1 - roman_Ω roman_ℓ end_ARG , with italic_W ≡ roman_ln italic_V start_POSTSUBSCRIPT italic_e italic_f italic_f end_POSTSUBSCRIPT , and italic_V start_POSTSUBSCRIPT italic_e italic_f italic_f end_POSTSUBSCRIPT = italic_u start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , (9)

where Veffsubscript𝑉𝑒𝑓𝑓V_{eff}italic_V start_POSTSUBSCRIPT italic_e italic_f italic_f end_POSTSUBSCRIPT is the torus effective potential. Tori are regulated by the balance of the hydrostatic and centrifugal factors due to the fluid rotation and by the curvature effects of the background, encoded in the effective potential function Veffsubscript𝑉𝑒𝑓𝑓V_{eff}italic_V start_POSTSUBSCRIPT italic_e italic_f italic_f end_POSTSUBSCRIPT.

Assuming the fluid is characterized by the specific angular momentum \ellroman_ℓ constant (see also discussion [55]), we consider the equation for W:ln(Veff)=c=constant:𝑊subscript𝑉𝑒𝑓𝑓cconstantW:\;\ln(V_{eff})=\rm{c}=\rm{constant}italic_W : roman_ln ( italic_V start_POSTSUBSCRIPT italic_e italic_f italic_f end_POSTSUBSCRIPT ) = roman_c = roman_constant or Veff=K=subscript𝑉𝑒𝑓𝑓𝐾absentV_{eff}=K=italic_V start_POSTSUBSCRIPT italic_e italic_f italic_f end_POSTSUBSCRIPT = italic_K =constant. By setting =absent\ell=roman_ℓ =constant as a torus parameter, the maximum density points in the disk, the pressure gradients (from the Euler equation) are determined by the gradients of the tori effective potential function121212 The procedure adopted here borrows from the Boyer theory on the equipressure surfaces applied to a thick torus [81, 58]. The Boyer surfaces tori are given by the surfaces of constant pressure.. The maximum points of the tori effective potential as function of the radial coordinate provide the minimum points of pressure, where fluid particles are free on unstable circular geodetic orbits.

II.2.1 The fluid effective potential

The fluid effective potential (9) is explicitly [58, 71]

Veff2=(μ)2=gtϕ2gϕϕgttgϕϕ+2gtϕ+gtt2.superscriptsubscript𝑉𝑒𝑓𝑓2superscript𝜇2superscriptsubscript𝑔𝑡italic-ϕ2subscript𝑔italic-ϕitalic-ϕsubscript𝑔𝑡𝑡subscript𝑔italic-ϕitalic-ϕ2subscript𝑔𝑡italic-ϕsubscript𝑔𝑡𝑡superscript2\displaystyle V_{eff}^{2}=\left(\frac{\mathcal{E}}{\mu}\right)^{2}={\frac{g_{t% \phi}^{2}-g_{\phi\phi}g_{tt}}{g_{\phi\phi}+2g_{t\phi}\ell+g_{tt}\ell^{2}}}.italic_V start_POSTSUBSCRIPT italic_e italic_f italic_f end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = ( divide start_ARG caligraphic_E end_ARG start_ARG italic_μ end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = divide start_ARG italic_g start_POSTSUBSCRIPT italic_t italic_ϕ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_g start_POSTSUBSCRIPT italic_ϕ italic_ϕ end_POSTSUBSCRIPT italic_g start_POSTSUBSCRIPT italic_t italic_t end_POSTSUBSCRIPT end_ARG start_ARG italic_g start_POSTSUBSCRIPT italic_ϕ italic_ϕ end_POSTSUBSCRIPT + 2 italic_g start_POSTSUBSCRIPT italic_t italic_ϕ end_POSTSUBSCRIPT roman_ℓ + italic_g start_POSTSUBSCRIPT italic_t italic_t end_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . (10)

The extremes of the pressure are regulated by the angular momentum distributions ±:rVeff=0:superscriptplus-or-minussubscript𝑟subscript𝑉𝑒𝑓𝑓0\ell^{\pm}:\partial_{r}V_{eff}=0roman_ℓ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT : ∂ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT italic_V start_POSTSUBSCRIPT italic_e italic_f italic_f end_POSTSUBSCRIPT = 0, on the equatorial plane θ=π/2𝜃𝜋2\theta=\pi/2italic_θ = italic_π / 2 for co-rotating ()(-)( - ) and counter-rotating (+)(+)( + ) fluids respectively131313Note, in the test particles analysis and accretion tori models, for slowly spinning NSs (a]M,1.29M[a\in]M,1.29M[italic_a ∈ ] italic_M , 1.29 italic_M [), there are circular geodesic orbits with (<0,<0)formulae-sequence00(\mathcal{E}<0,\mathcal{L}<0)( caligraphic_E < 0 , caligraphic_L < 0 ) and (>0,<0)formulae-sequence00(\mathcal{E}>0,\mathcal{L}<0)( caligraphic_E > 0 , caligraphic_L < 0 ) on the equatorial plane of the ergoregion–see Figs (1). These solutions correspond to the relativistic angular velocity (the Keplerian velocity with respect to static observers at infinity Ω=dϕ/dtΩ𝑑italic-ϕ𝑑𝑡\Omega=d\phi/dtroman_Ω = italic_d italic_ϕ / italic_d italic_t) Ω>0Ω0\Omega>0roman_Ω > 0; therefore, in this sense, they are all co-rotating with respect to the static observers at infinity but they can be counter-rotating according to <00\mathcal{L}<0caligraphic_L < 0 and <00\ell<0roman_ℓ < 0 or counter-rotating according to <00\mathcal{L}<0caligraphic_L < 0 but co-rotating according to >00\ell>0roman_ℓ > 0, there can also be orbits with =Ω=0Ω0\ell=\Omega=0roman_ℓ = roman_Ω = 0 –see, for example, [22, 67, 68, 69, 70, 80, 73]. This possibility has not been discussed in the analysis of DM models. .

Torus cusp r×subscript𝑟r_{\times}italic_r start_POSTSUBSCRIPT × end_POSTSUBSCRIPT is the minimum point of pressure and density in the torus corresponding to the maximum point of the fluid effective potential. The torus center rcentersubscript𝑟𝑐𝑒𝑛𝑡𝑒𝑟r_{center}italic_r start_POSTSUBSCRIPT italic_c italic_e italic_n italic_t italic_e italic_r end_POSTSUBSCRIPT is the maximum point of pressure and density in the torus, corresponding to the minimum point of the fluid effective potential. At the cusp (rr×𝑟subscript𝑟r\leq r_{\times}italic_r ≤ italic_r start_POSTSUBSCRIPT × end_POSTSUBSCRIPT) the fluid may be considered pressure-free. Fluid effective potential defines the function K(r)=Veff((r))𝐾𝑟subscript𝑉𝑒𝑓𝑓𝑟K(r)=V_{eff}(\ell(r))italic_K ( italic_r ) = italic_V start_POSTSUBSCRIPT italic_e italic_f italic_f end_POSTSUBSCRIPT ( roman_ℓ ( italic_r ) ). Cusped tori have parameter in the open ranges141414The notation is as follows: the (closed) interval between quantities qlsubscript𝑞𝑙q_{l}italic_q start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT and qrsubscript𝑞𝑟q_{r}italic_q start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT, including qlsubscript𝑞𝑙q_{l}italic_q start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT and qrsubscript𝑞𝑟q_{r}italic_q start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT, is denoted as[ql,qr]={q:qlqqr}subscript𝑞𝑙subscript𝑞𝑟conditional-set𝑞subscript𝑞𝑙𝑞subscript𝑞𝑟[q_{l},q_{r}]=\{q\in\mathbb{R}:q_{l}\leq q\leq q_{r}\}[ italic_q start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT , italic_q start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT ] = { italic_q ∈ blackboard_R : italic_q start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT ≤ italic_q ≤ italic_q start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT }, the notation ]ql,qr[={q:ql<q<qr}]q_{l},q_{r}[=\{q\in\mathbb{R}:q_{l}<q<q_{r}\}] italic_q start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT , italic_q start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT [ = { italic_q ∈ blackboard_R : italic_q start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT < italic_q < italic_q start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT } denotes the open interval. Similarly there is ]ql,qr]={q:ql<qqr}]q_{l},q_{r}]=\{q\in\mathbb{R}:q_{l}<q\leq q_{r}\}] italic_q start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT , italic_q start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT ] = { italic_q ∈ blackboard_R : italic_q start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT < italic_q ≤ italic_q start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT } and [ql,qr[={q:qlq<qr}[q_{l},q_{r}[=\{q\in\mathbb{R}:q_{l}\leq q<q_{r}\}[ italic_q start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT , italic_q start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT [ = { italic_q ∈ blackboard_R : italic_q start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT ≤ italic_q < italic_q start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT }.K=K×K(r×)]Kcenter,1[]Kmso,1[K=K_{\times}\equiv K(r_{\times})\in]K_{center},1[\subset]K_{mso},1[italic_K = italic_K start_POSTSUBSCRIPT × end_POSTSUBSCRIPT ≡ italic_K ( italic_r start_POSTSUBSCRIPT × end_POSTSUBSCRIPT ) ∈ ] italic_K start_POSTSUBSCRIPT italic_c italic_e italic_n italic_t italic_e italic_r end_POSTSUBSCRIPT , 1 [ ⊂ ] italic_K start_POSTSUBSCRIPT italic_m italic_s italic_o end_POSTSUBSCRIPT , 1 [, where KcenterK(rcenter)subscript𝐾𝑐𝑒𝑛𝑡𝑒𝑟𝐾subscript𝑟𝑐𝑒𝑛𝑡𝑒𝑟K_{center}\equiv K(r_{center})italic_K start_POSTSUBSCRIPT italic_c italic_e italic_n italic_t italic_e italic_r end_POSTSUBSCRIPT ≡ italic_K ( italic_r start_POSTSUBSCRIPT italic_c italic_e italic_n italic_t italic_e italic_r end_POSTSUBSCRIPT ). (We adopt the notation qq(r)subscript𝑞𝑞subscript𝑟q_{\bullet}\equiv q(r_{\bullet})italic_q start_POSTSUBSCRIPT ∙ end_POSTSUBSCRIPT ≡ italic_q ( italic_r start_POSTSUBSCRIPT ∙ end_POSTSUBSCRIPT ) for any quantity q𝑞qitalic_q evaluated on a radius rsubscript𝑟r_{\bullet}italic_r start_POSTSUBSCRIPT ∙ end_POSTSUBSCRIPT.)

II.2.2 Extended geodesic structure and notable radii

The geometry equatorial circular geodesic structure constrains the accretion disk physics governing, in the P-D model, the tori cusps and centers locations. In the Kerr geometry the geodesic structure is constituted by the marginally circular orbit for timelike particles rmso±superscriptsubscript𝑟𝑚𝑠𝑜plus-or-minusr_{mso}^{\pm}italic_r start_POSTSUBSCRIPT italic_m italic_s italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT, which is also a photon circular orbit, rmco±rγ±superscriptsubscript𝑟𝑚𝑐𝑜plus-or-minussuperscriptsubscript𝑟𝛾plus-or-minusr_{mco}^{\pm}\equiv r_{\gamma}^{\pm}italic_r start_POSTSUBSCRIPT italic_m italic_c italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ≡ italic_r start_POSTSUBSCRIPT italic_γ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT, the marginally bounded orbit, rmbo±superscriptsubscript𝑟𝑚𝑏𝑜plus-or-minusr_{mbo}^{\pm}italic_r start_POSTSUBSCRIPT italic_m italic_b italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT, and the marginally stable circular orbit, rmso±superscriptsubscript𝑟𝑚𝑠𝑜plus-or-minusr_{mso}^{\pm}italic_r start_POSTSUBSCRIPT italic_m italic_s italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT (see Figs (1))151515In the Kerr spacetime rmco±superscriptsubscript𝑟𝑚𝑐𝑜plus-or-minusr_{mco}^{\pm}italic_r start_POSTSUBSCRIPT italic_m italic_c italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT is the marginal circular orbit (a photon circular orbit) where timelike circular orbits can fill the spacetime region r>rmco±𝑟superscriptsubscript𝑟𝑚𝑐𝑜plus-or-minusr>r_{mco}^{\pm}italic_r > italic_r start_POSTSUBSCRIPT italic_m italic_c italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT. Stable circular orbits are in r>rmso±𝑟superscriptsubscript𝑟msoplus-or-minusr>r_{\mathrm{mso}}^{\pm}italic_r > italic_r start_POSTSUBSCRIPT roman_mso end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT for counter-rotating and co-rotating test particles respectively. The marginal bounded circular orbit is rmbo±superscriptsubscript𝑟mboplus-or-minusr_{\mathrm{mbo}}^{\pm}italic_r start_POSTSUBSCRIPT roman_mbo end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT, where ±(rmbo±)=1superscriptplus-or-minussuperscriptsubscript𝑟mboplus-or-minus1\mathcal{E}^{\pm}(r_{\mathrm{mbo}}^{\pm})=1caligraphic_E start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ( italic_r start_POSTSUBSCRIPT roman_mbo end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ) = 1. More details and the exact forms of these radii can be found for example in [68].. Radii {rmso±,rmbo±,rmco±}superscriptsubscript𝑟𝑚𝑠𝑜plus-or-minussuperscriptsubscript𝑟𝑚𝑏𝑜plus-or-minussuperscriptsubscript𝑟𝑚𝑐𝑜plus-or-minus\{r_{mso}^{\pm},r_{mbo}^{\pm},r_{mco}^{\pm}\}{ italic_r start_POSTSUBSCRIPT italic_m italic_s italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT , italic_r start_POSTSUBSCRIPT italic_m italic_b italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT , italic_r start_POSTSUBSCRIPT italic_m italic_c italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT } constrain the location of the tori cusps (inner edges) with fluid specific angular momentum =±superscriptplus-or-minus\ell=\ell^{\pm}roman_ℓ = roman_ℓ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT respectively

wherermco±<rmbo±<rmso±<r(mbo)±<r(mco)±,wheresuperscriptsubscript𝑟𝑚𝑐𝑜plus-or-minussuperscriptsubscript𝑟mboplus-or-minussuperscriptsubscript𝑟msoplus-or-minussuperscriptsubscript𝑟mboplus-or-minussuperscriptsubscript𝑟𝑚𝑐𝑜plus-or-minus\displaystyle\mbox{where}\quad r_{mco}^{\pm}<r_{\mathrm{mbo}}^{\pm}<r_{\mathrm% {mso}}^{\pm}<r_{\mathrm{(mbo)}}^{\pm}<r_{(mco)}^{\pm},where italic_r start_POSTSUBSCRIPT italic_m italic_c italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT < italic_r start_POSTSUBSCRIPT roman_mbo end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT < italic_r start_POSTSUBSCRIPT roman_mso end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT < italic_r start_POSTSUBSCRIPT ( roman_mbo ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT < italic_r start_POSTSUBSCRIPT ( italic_m italic_c italic_o ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT , (11)
Refer to caption
Refer to caption
Figure 1: Geodesic equatorial circular structure of the Kerr geometry for spin a[0,M]𝑎0𝑀a\in[0,M]italic_a ∈ [ 0 , italic_M ] (left panel) and a>M𝑎𝑀a>Mitalic_a > italic_M (right panel), for co-rotating (()(-)( - )-red-curves) and counter-rotating ((+)(+)( + )- blue curves) orbits. The marginally stable orbits (mso𝑚𝑠𝑜msoitalic_m italic_s italic_o) are the solid curves, the marginally bounded orbits (mbo𝑚𝑏𝑜mboitalic_m italic_b italic_o) are dashed curves, and marginally circular orbits rmco±superscriptsubscript𝑟𝑚𝑐𝑜plus-or-minusr_{mco}^{\pm}italic_r start_POSTSUBSCRIPT italic_m italic_c italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT (which are also photon circular orbits) are the dotted curves. (Note for a>M𝑎𝑀a>Mitalic_a > italic_M there is no last co-rotating circular orbit (i.e. rmco=0superscriptsubscript𝑟𝑚𝑐𝑜0r_{mco}^{-}=0italic_r start_POSTSUBSCRIPT italic_m italic_c italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT = 0.) Radius r+subscript𝑟r_{+}italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT is the outer horizon, rϵ+superscriptsubscript𝑟italic-ϵr_{\epsilon}^{+}italic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT is the outer ergosurface on the equatorial plane. Black region is r<r+𝑟subscript𝑟r<r_{+}italic_r < italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT and gray region is r[r+,rϵ+]𝑟subscript𝑟superscriptsubscript𝑟italic-ϵr\in[r_{+},r_{\epsilon}^{+}]italic_r ∈ [ italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT , italic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT ] for Kerr BHs and r[0,rϵ+]𝑟0superscriptsubscript𝑟italic-ϵr\in[0,r_{\epsilon}^{+}]italic_r ∈ [ 0 , italic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT ] for Kerr NSs. On radii r0±superscriptsubscript𝑟0plus-or-minusr_{0}^{\pm}italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT (black dashed curve-right panel) there is =00\mathcal{L}=0caligraphic_L = 0 (where \mathcal{L}caligraphic_L is the test particle angular momentum) and on rδ±superscriptsubscript𝑟𝛿plus-or-minusr_{\delta}^{\pm}italic_r start_POSTSUBSCRIPT italic_δ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT (black curve-right panel) there is =00\mathcal{E}=0caligraphic_E = 0 (where \mathcal{E}caligraphic_E is the test particle energy).

where we introduced also the radii (r(mbo)±,r(mco)±)superscriptsubscript𝑟mboplus-or-minussuperscriptsubscript𝑟𝑚𝑐𝑜plus-or-minus(r_{\mathrm{(mbo)}}^{\pm},r_{(mco)}^{\pm})( italic_r start_POSTSUBSCRIPT ( roman_mbo ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT , italic_r start_POSTSUBSCRIPT ( italic_m italic_c italic_o ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ) defined as

r(mbo)±:±(rmbo±)=±(r(mbo)±)mbo±,r(mco)±:±(rmco±)=±(r(mco)±)𝐦𝐜𝐨±,\displaystyle r_{\mathrm{(mbo)}}^{\pm}:\;\ell^{\pm}(r_{\mathrm{mbo}}^{\pm})=% \ell^{\pm}(r_{\mathrm{(mbo)}}^{\pm})\equiv\mathbf{\ell_{\mathrm{mbo}}^{\pm}},% \quad r_{(mco)}^{\pm}:\ell^{\pm}(r_{mco}^{\pm})=\ell^{\pm}(r_{(mco)}^{\pm})% \equiv\mathbf{\ell_{mco}^{\pm}},italic_r start_POSTSUBSCRIPT ( roman_mbo ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT : roman_ℓ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ( italic_r start_POSTSUBSCRIPT roman_mbo end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ) = roman_ℓ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ( italic_r start_POSTSUBSCRIPT ( roman_mbo ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ) ≡ roman_ℓ start_POSTSUBSCRIPT roman_mbo end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT , italic_r start_POSTSUBSCRIPT ( italic_m italic_c italic_o ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT : roman_ℓ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ( italic_r start_POSTSUBSCRIPT italic_m italic_c italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ) = roman_ℓ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ( italic_r start_POSTSUBSCRIPT ( italic_m italic_c italic_o ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ) ≡ roman_ℓ start_POSTSUBSCRIPT bold_mco end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ,

significant as governing the location of the tori centers, More precisely ranges (𝐋𝟏,𝐋𝟐,𝐋𝟑)subscript𝐋1subscript𝐋2subscript𝐋3(\mathbf{L_{1},L_{2},L_{3}})( bold_L start_POSTSUBSCRIPT bold_1 end_POSTSUBSCRIPT , bold_L start_POSTSUBSCRIPT bold_2 end_POSTSUBSCRIPT , bold_L start_POSTSUBSCRIPT bold_3 end_POSTSUBSCRIPT ) of fluids specific angular momentum \ellroman_ℓ govern the tori topology, according to the geodesic structure of Eqs (11), as follows:

𝐋𝟏subscript𝐋1\mathbf{L_{1}}bold_L start_POSTSUBSCRIPT bold_1 end_POSTSUBSCRIPT:

for 𝐋𝟏subscript𝐋1\ell\in\mathbf{L_{1}}roman_ℓ ∈ bold_L start_POSTSUBSCRIPT bold_1 end_POSTSUBSCRIPT there are quiescent (i.e. not cusped) and cusped tori–where there is 𝐋𝟏±[mso±,mbo±[\mp\mathbf{L_{1}}^{\pm}\equiv[\mp\ell_{mso}^{\pm},\mp\ell_{mbo}^{\pm}[∓ bold_L start_POSTSUBSCRIPT bold_1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ≡ [ ∓ roman_ℓ start_POSTSUBSCRIPT italic_m italic_s italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT , ∓ roman_ℓ start_POSTSUBSCRIPT italic_m italic_b italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT [. The cusp is r×±]rmbo±,rmso±]r^{\pm}_{\times}\in]r^{\pm}_{mbo},r^{\pm}_{mso}]italic_r start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT start_POSTSUBSCRIPT × end_POSTSUBSCRIPT ∈ ] italic_r start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_m italic_b italic_o end_POSTSUBSCRIPT , italic_r start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_m italic_s italic_o end_POSTSUBSCRIPT ] (with K×±<1superscriptsubscript𝐾plus-or-minus1K_{\times}^{\pm}<1italic_K start_POSTSUBSCRIPT × end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT < 1)) and the center with maximum pressure in rcenter±]rmso±,r(mbo)±]r^{\pm}_{center}\in]r^{\pm}_{mso},r^{\pm}_{(mbo)}]italic_r start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_c italic_e italic_n italic_t italic_e italic_r end_POSTSUBSCRIPT ∈ ] italic_r start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_m italic_s italic_o end_POSTSUBSCRIPT , italic_r start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ( italic_m italic_b italic_o ) end_POSTSUBSCRIPT ].

𝐋𝟐subscript𝐋2\mathbf{L_{2}}bold_L start_POSTSUBSCRIPT bold_2 end_POSTSUBSCRIPT:

for 𝐋𝟐subscript𝐋2\ell\in\mathbf{L_{2}}roman_ℓ ∈ bold_L start_POSTSUBSCRIPT bold_2 end_POSTSUBSCRIPT there are quiescent tori and proto-jets (open-configurations) –where there is 𝐋𝟐±[mbo±,mco±[\mp\mathbf{L_{2}}^{\pm}\equiv[\mp\ell_{mbo}^{\pm},\mp\ell_{mco}^{\pm}[∓ bold_L start_POSTSUBSCRIPT bold_2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ≡ [ ∓ roman_ℓ start_POSTSUBSCRIPT italic_m italic_b italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT , ∓ roman_ℓ start_POSTSUBSCRIPT italic_m italic_c italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT [. The cusp r×±]rmco±,rmbo±]r_{\times}^{\pm}\in]r_{mco}^{\pm},r_{mbo}^{\pm}]italic_r start_POSTSUBSCRIPT × end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ∈ ] italic_r start_POSTSUBSCRIPT italic_m italic_c italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT , italic_r start_POSTSUBSCRIPT italic_m italic_b italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ] is associated to the proto-jets, with K×>1subscript𝐾1K_{\times}>1italic_K start_POSTSUBSCRIPT × end_POSTSUBSCRIPT > 1, and the center with maximum pressure is in rcenter±]r(mbo)±,r(mco)±]r_{center}^{\pm}\in]r_{(mbo)}^{\pm},r_{(mco)}^{\pm}]italic_r start_POSTSUBSCRIPT italic_c italic_e italic_n italic_t italic_e italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ∈ ] italic_r start_POSTSUBSCRIPT ( italic_m italic_b italic_o ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT , italic_r start_POSTSUBSCRIPT ( italic_m italic_c italic_o ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ]. Proto-jets are associated to (not-collimated) open structures, with matter funnels along the BH rotational axis–see [59, 74, 75, 76];

𝐋𝟑subscript𝐋3\mathbf{L_{3}}bold_L start_POSTSUBSCRIPT bold_3 end_POSTSUBSCRIPT:

for 𝐋𝟑subscript𝐋3\ell\in\mathbf{L_{3}}roman_ℓ ∈ bold_L start_POSTSUBSCRIPT bold_3 end_POSTSUBSCRIPT there are only quiescent tori where there is 𝐋𝟑±mco±minus-or-plussuperscriptsubscript𝐋3plus-or-minusminus-or-plusminus-or-plussuperscriptsubscript𝑚𝑐𝑜plus-or-minus\mp\mathbf{L_{3}}^{\pm}\equiv\mp\ell\geq\mp\ell_{mco}^{\pm}∓ bold_L start_POSTSUBSCRIPT bold_3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ≡ ∓ roman_ℓ ≥ ∓ roman_ℓ start_POSTSUBSCRIPT italic_m italic_c italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT and the torus center is at rcenter±>r(mco)±subscriptsuperscript𝑟plus-or-minus𝑐𝑒𝑛𝑡𝑒𝑟superscriptsubscript𝑟𝑚𝑐𝑜plus-or-minusr^{\pm}_{center}>r_{(mco)}^{\pm}italic_r start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_c italic_e italic_n italic_t italic_e italic_r end_POSTSUBSCRIPT > italic_r start_POSTSUBSCRIPT ( italic_m italic_c italic_o ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT.

In the metric models we consider, the DM affects the orbiting fluids modifying the Kerr axially symmetric geometry and the fluid effective potential–see also [78]. Therefore we study the radii limiting the tori construction, defined through the fluid effective potential for the geometries modified by the DM. More precisely we identify the marginally circular orbit, rmco±superscriptsubscript𝑟𝑚𝑐𝑜plus-or-minusr_{mco}^{\pm}italic_r start_POSTSUBSCRIPT italic_m italic_c italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT, as the radius rmco±:K±(r)=:superscriptsubscript𝑟𝑚𝑐𝑜plus-or-minussuperscript𝐾plus-or-minus𝑟r_{mco}^{\pm}:K^{\pm}(r)=\inftyitalic_r start_POSTSUBSCRIPT italic_m italic_c italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT : italic_K start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ( italic_r ) = ∞, the marginally bounded orbit defined by rmbo±:K±(r)=1:superscriptsubscript𝑟𝑚𝑏𝑜plus-or-minussuperscript𝐾plus-or-minus𝑟1r_{mbo}^{\pm}:K^{\pm}(r)=1italic_r start_POSTSUBSCRIPT italic_m italic_b italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT : italic_K start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ( italic_r ) = 1 (asymptotically flat spacetimes) and the marginally stable orbits rmso±:r±=0:superscriptsubscript𝑟𝑚𝑠𝑜plus-or-minussubscript𝑟superscriptplus-or-minus0r_{mso}^{\pm}:\partial_{r}\ell^{\pm}=0italic_r start_POSTSUBSCRIPT italic_m italic_s italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT : ∂ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT = 0.

The orbiting fluid is governed by the geodesic structure of the considered spacetime. The tori are specified by the profile of the distribution of the specific angular momentum of the orbiting matter, (in the equator) and its relation to the radial profile of the specific angular momentum of equatorial circular geodesics. The so called Keplerian distribution of the circular geodesic angular momentum related to a given spacetime is generally given by the relation

=KΦoΦΦ()2Φ,whereΦ()gtϕ2gttgϕϕ,Φ(gtϕ)2gttgϕϕ,formulae-sequencesubscript𝐾minus-or-plussubscriptΦ𝑜ΦsuperscriptsubscriptΦ2subscriptΦwheresubscriptΦminus-or-plusminus-or-plussuperscriptsubscript𝑔𝑡italic-ϕ2subscript𝑔𝑡𝑡subscript𝑔italic-ϕitalic-ϕΦsuperscriptsuperscriptsubscript𝑔𝑡italic-ϕ2superscriptsubscript𝑔𝑡𝑡superscriptsubscript𝑔italic-ϕitalic-ϕ\displaystyle\ell=\ell_{K}\equiv\frac{\Phi_{o}\mp\sqrt{\Phi\Phi_{(-)}^{2}}}{% \Phi_{\star}},\quad\mbox{where}\quad\Phi_{(\mp)}\equiv g_{t\phi}^{2}\mp g_{tt}% g_{\phi\phi},\quad{\Phi}\equiv\left(g_{t\phi}^{\prime}\right)^{2}-g_{tt}^{% \prime}g_{\phi\phi}^{\prime},roman_ℓ = roman_ℓ start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT ≡ divide start_ARG roman_Φ start_POSTSUBSCRIPT italic_o end_POSTSUBSCRIPT ∓ square-root start_ARG roman_Φ roman_Φ start_POSTSUBSCRIPT ( - ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG start_ARG roman_Φ start_POSTSUBSCRIPT ⋆ end_POSTSUBSCRIPT end_ARG , where roman_Φ start_POSTSUBSCRIPT ( ∓ ) end_POSTSUBSCRIPT ≡ italic_g start_POSTSUBSCRIPT italic_t italic_ϕ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ∓ italic_g start_POSTSUBSCRIPT italic_t italic_t end_POSTSUBSCRIPT italic_g start_POSTSUBSCRIPT italic_ϕ italic_ϕ end_POSTSUBSCRIPT , roman_Φ ≡ ( italic_g start_POSTSUBSCRIPT italic_t italic_ϕ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_g start_POSTSUBSCRIPT italic_t italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_ϕ italic_ϕ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , (12)
ΦogtϕΦ(+)gtϕ(gϕϕgtt),Φgtϕ2gtt+gtt2gϕϕgtt(gtϕ2),formulae-sequencesubscriptΦ𝑜superscriptsubscript𝑔𝑡italic-ϕsubscriptΦsubscript𝑔𝑡italic-ϕsuperscriptsubscript𝑔italic-ϕitalic-ϕsubscript𝑔𝑡𝑡subscriptΦsuperscriptsubscript𝑔𝑡italic-ϕ2superscriptsubscript𝑔𝑡𝑡superscriptsubscript𝑔𝑡𝑡2superscriptsubscript𝑔italic-ϕitalic-ϕsubscript𝑔𝑡𝑡superscriptsuperscriptsubscript𝑔𝑡italic-ϕ2\displaystyle\Phi_{o}\equiv g_{t\phi}^{\prime}\Phi_{(+)}-g_{t\phi}\left(g_{% \phi\phi}g_{tt}\right)^{\prime},\quad\Phi_{\star}\equiv g_{t\phi}^{2}g_{tt}^{% \prime}+g_{tt}^{2}g_{\phi\phi}^{\prime}-g_{tt}(g_{t\phi}^{2})^{\prime},roman_Φ start_POSTSUBSCRIPT italic_o end_POSTSUBSCRIPT ≡ italic_g start_POSTSUBSCRIPT italic_t italic_ϕ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT roman_Φ start_POSTSUBSCRIPT ( + ) end_POSTSUBSCRIPT - italic_g start_POSTSUBSCRIPT italic_t italic_ϕ end_POSTSUBSCRIPT ( italic_g start_POSTSUBSCRIPT italic_ϕ italic_ϕ end_POSTSUBSCRIPT italic_g start_POSTSUBSCRIPT italic_t italic_t end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , roman_Φ start_POSTSUBSCRIPT ⋆ end_POSTSUBSCRIPT ≡ italic_g start_POSTSUBSCRIPT italic_t italic_ϕ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_t italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT + italic_g start_POSTSUBSCRIPT italic_t italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_ϕ italic_ϕ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT - italic_g start_POSTSUBSCRIPT italic_t italic_t end_POSTSUBSCRIPT ( italic_g start_POSTSUBSCRIPT italic_t italic_ϕ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ,

for θ=π/2𝜃𝜋2\theta=\pi/2italic_θ = italic_π / 2, where ()(^{\prime})( start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) is for the derivative with respect to r𝑟ritalic_r. The Keplerian profile intersection with the tori profile determines centers and cusps of the tori. In the standard Kerr spacetime it takes the well known form:

K(r,a)=a3r3/2Δa(4M3r)ra2(r2M)2r.subscript𝐾𝑟𝑎superscriptminus-or-plusminus-or-plussuperscript𝑎3superscript𝑟32Δ𝑎4𝑀3𝑟𝑟superscript𝑎2superscript𝑟2𝑀2𝑟\displaystyle\ell_{K}(r,a)=\ell^{\mp}\equiv\frac{a^{3}\mp r^{3/2}\Delta-a(4M-3% r)r}{a^{2}-(r-2M)^{2}r}.roman_ℓ start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT ( italic_r , italic_a ) = roman_ℓ start_POSTSUPERSCRIPT ∓ end_POSTSUPERSCRIPT ≡ divide start_ARG italic_a start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ∓ italic_r start_POSTSUPERSCRIPT 3 / 2 end_POSTSUPERSCRIPT roman_Δ - italic_a ( 4 italic_M - 3 italic_r ) italic_r end_ARG start_ARG italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - ( italic_r - 2 italic_M ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_r end_ARG . (13)

III Dark matter models

We analyze accretion tori orbiting spinning BHs with spacetimes influenced by different DM models. The metrics reduce, for some limiting values of the DM parameters reduces to the Kerr BH geometry. Thus, using Eqs. (10) we consider the BH DM metric components {gtt,gtϕ,gϕϕ}subscript𝑔𝑡𝑡subscript𝑔𝑡italic-ϕsubscript𝑔italic-ϕitalic-ϕ\{g_{tt},g_{t\phi},g_{\phi\phi}\}{ italic_g start_POSTSUBSCRIPT italic_t italic_t end_POSTSUBSCRIPT , italic_g start_POSTSUBSCRIPT italic_t italic_ϕ end_POSTSUBSCRIPT , italic_g start_POSTSUBSCRIPT italic_ϕ italic_ϕ end_POSTSUBSCRIPT } in BL coordinates, and we refer to the literature for details on the metric tensor and the geometry properties. We investigate the equatorial circular geodesic structures for the fluid effective potential, the effective potential function and the tori structure for co-rotating and counter-rotating tori, in three DM models: in Sec. (III.1) we address the perfect fluid DM (PFDM) model of [4]. Cold and scalar field DM models of [9, 8] are discussed Sec. (III.2).

III.1 Perfect fluid dark matter (PFDM)

A rotating BH solution in perfect fluid dark matter (PFDM) has been discussed in [4], with

gtt=[12MrfD(r)Σ],gtϕ=σa[2MrfD(r)]Σ,gϕϕ=σ[a2σ[2MrfD(r)]Σ+(a2+r2)],formulae-sequencesubscript𝑔𝑡𝑡delimited-[]12𝑀𝑟subscript𝑓𝐷𝑟Σformulae-sequencesubscript𝑔𝑡italic-ϕ𝜎𝑎delimited-[]2𝑀𝑟subscript𝑓𝐷𝑟Σsubscript𝑔italic-ϕitalic-ϕ𝜎delimited-[]superscript𝑎2𝜎delimited-[]2𝑀𝑟subscript𝑓𝐷𝑟Σsuperscript𝑎2superscript𝑟2\displaystyle g_{tt}=-\left[1-\frac{2Mr-f_{D}(r)}{\Sigma}\right],\quad g_{t% \phi}=-\frac{\sigma a[2Mr-f_{D}(r)]}{\Sigma},\quad g_{\phi\phi}=\sigma\left[% \frac{a^{2}\sigma[2Mr-f_{D}(r)]}{\Sigma}+\left(a^{2}+r^{2}\right)\right],italic_g start_POSTSUBSCRIPT italic_t italic_t end_POSTSUBSCRIPT = - [ 1 - divide start_ARG 2 italic_M italic_r - italic_f start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT ( italic_r ) end_ARG start_ARG roman_Σ end_ARG ] , italic_g start_POSTSUBSCRIPT italic_t italic_ϕ end_POSTSUBSCRIPT = - divide start_ARG italic_σ italic_a [ 2 italic_M italic_r - italic_f start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT ( italic_r ) ] end_ARG start_ARG roman_Σ end_ARG , italic_g start_POSTSUBSCRIPT italic_ϕ italic_ϕ end_POSTSUBSCRIPT = italic_σ [ divide start_ARG italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_σ [ 2 italic_M italic_r - italic_f start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT ( italic_r ) ] end_ARG start_ARG roman_Σ end_ARG + ( italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) ] , (14)
grrΣΔD,gθθ=Σformulae-sequencesubscript𝑔𝑟𝑟ΣsubscriptΔ𝐷subscript𝑔𝜃𝜃Σ\displaystyle g_{rr}\equiv\frac{\Sigma}{\Delta_{{D}}},\quad g_{\theta\theta}=\Sigmaitalic_g start_POSTSUBSCRIPT italic_r italic_r end_POSTSUBSCRIPT ≡ divide start_ARG roman_Σ end_ARG start_ARG roman_Δ start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT end_ARG , italic_g start_POSTSUBSCRIPT italic_θ italic_θ end_POSTSUBSCRIPT = roman_Σ (15)
whereΔDΔ+fD(r)andfD(r)krlogr|k|,formulae-sequencewheresubscriptΔ𝐷Δsubscript𝑓𝐷𝑟andsubscript𝑓𝐷𝑟𝑘𝑟𝑟𝑘\displaystyle\mbox{where}\quad\Delta_{D}\equiv\Delta+f_{D}(r)\quad\mbox{and}% \quad f_{D}(r)\equiv kr\log\frac{r}{|k|},where roman_Δ start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT ≡ roman_Δ + italic_f start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT ( italic_r ) and italic_f start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT ( italic_r ) ≡ italic_k italic_r roman_log divide start_ARG italic_r end_ARG start_ARG | italic_k | end_ARG ,

see also [5, 6, 7], where k𝑘kitalic_k is the parameter describing the intensity of the PFDM, set in the ranges k]7.18M,2M[k\in]-7.18M,2M[italic_k ∈ ] - 7.18 italic_M , 2 italic_M [. For k=0𝑘0k=0italic_k = 0 the line element reduces to the Kerr metric161616Constrains of k/M𝑘𝑀k/Mitalic_k / italic_M positive were obtained by fitting the rotation curves in spiral galaxies, with values 106107superscript106superscript10710^{-6}-10^{-7}10 start_POSTSUPERSCRIPT - 6 end_POSTSUPERSCRIPT - 10 start_POSTSUPERSCRIPT - 7 end_POSTSUPERSCRIPT [4].. The metric singularities, r±subscript𝑟plus-or-minusr_{\pm}italic_r start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT, defining the DM deformations on the Kerr horizons, can be found by solving the equation a=a±r(2Mr)fD(r),𝑎subscript𝑎plus-or-minus𝑟2𝑀𝑟subscript𝑓𝐷𝑟a=a_{\pm}\equiv\sqrt{r(2M-r)-f_{D}(r)},italic_a = italic_a start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ≡ square-root start_ARG italic_r ( 2 italic_M - italic_r ) - italic_f start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT ( italic_r ) end_ARG , or from the equation k=k±Δ/[rW(Δ/r2)]𝑘subscript𝑘plus-or-minusΔdelimited-[]𝑟𝑊Δsuperscript𝑟2k=k_{\pm}\equiv\Delta/[rW\left({\Delta}/{r^{2}}\right)]italic_k = italic_k start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ≡ roman_Δ / [ italic_r italic_W ( roman_Δ / italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) ], where W𝑊Witalic_W is the Lambert function, such that W(z)𝑊𝑧W(z)italic_W ( italic_z ) gives the principal solution for w𝑤witalic_w in z=wew𝑧𝑤superscript𝑒𝑤z=we^{w}italic_z = italic_w italic_e start_POSTSUPERSCRIPT italic_w end_POSTSUPERSCRIPT.

The deformed ergosurfaces rϵ±superscriptsubscript𝑟italic-ϵplus-or-minusr_{\epsilon}^{\pm}italic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT can be found, for σ1𝜎1\sigma\neq 1italic_σ ≠ 1, as solution of the equation a=aϵ±a±/1σ𝑎superscriptsubscript𝑎italic-ϵplus-or-minussubscript𝑎plus-or-minus1𝜎a=a_{\epsilon}^{\pm}\equiv{a_{\pm}}/{\sqrt{1-\sigma}}italic_a = italic_a start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ≡ italic_a start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT / square-root start_ARG 1 - italic_σ end_ARG while, on the equatorial plane (σ=1)𝜎1(\sigma=1)( italic_σ = 1 ) there is

rϵ±={kW[e2M/ksgn(k)]<2,fork0,kW[1,e2M/k]>2,fork<0,2,fork=0superscriptsubscript𝑟italic-ϵplus-or-minuscasesformulae-sequence𝑘𝑊delimited-[]superscript𝑒2𝑀𝑘sgn𝑘2forless-than-or-greater-than𝑘0𝑜𝑡ℎ𝑒𝑟𝑤𝑖𝑠𝑒formulae-sequence𝑘𝑊1superscript𝑒2𝑀𝑘2for𝑘0𝑜𝑡ℎ𝑒𝑟𝑤𝑖𝑠𝑒2for𝑘0𝑜𝑡ℎ𝑒𝑟𝑤𝑖𝑠𝑒\displaystyle r_{\epsilon}^{\pm}=\begin{cases}kW\left[\frac{e^{2M/k}}{\text{% sgn}(k)}\right]<2,\quad\mbox{for}\quad k\lessgtr 0,\\ kW\left[-1,-e^{2M/k}\right]>2,\quad\mbox{for}\quad k<0,\\ 2,\quad\mbox{for}\quad k=0\end{cases}italic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT = { start_ROW start_CELL italic_k italic_W [ divide start_ARG italic_e start_POSTSUPERSCRIPT 2 italic_M / italic_k end_POSTSUPERSCRIPT end_ARG start_ARG sgn ( italic_k ) end_ARG ] < 2 , for italic_k ≶ 0 , end_CELL start_CELL end_CELL end_ROW start_ROW start_CELL italic_k italic_W [ - 1 , - italic_e start_POSTSUPERSCRIPT 2 italic_M / italic_k end_POSTSUPERSCRIPT ] > 2 , for italic_k < 0 , end_CELL start_CELL end_CELL end_ROW start_ROW start_CELL 2 , for italic_k = 0 end_CELL start_CELL end_CELL end_ROW (16)

(see red curve in Figs (2)–bottom left panel) where the Lambert function W(s,z)𝑊𝑠𝑧W(s,z)italic_W ( italic_s , italic_z ) gives the sthsuperscript𝑠ths^{\text{th}}italic_s start_POSTSUPERSCRIPT th end_POSTSUPERSCRIPT solution for w𝑤witalic_w in z=wew𝑧𝑤superscript𝑒𝑤z=we^{w}italic_z = italic_w italic_e start_POSTSUPERSCRIPT italic_w end_POSTSUPERSCRIPT, and sgn(k)sgn𝑘\text{sgn}(k)sgn ( italic_k ) gives the sign of k𝑘kitalic_k, therefore it is 1,0,1101-1,0,1- 1 , 0 , 1 for k𝑘kitalic_k negative, zero, or positive.

The PFDM horizons are independent of σ𝜎\sigmaitalic_σ (like in the Kerr case). The PFDM ergosurfaces are independent of spin in the equatorial plane (like in the Kerr cases). The horizons r±subscript𝑟plus-or-minusr_{\pm}italic_r start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT (and the ergosurfaces rϵ±superscriptsubscript𝑟italic-ϵplus-or-minusr_{\epsilon}^{\pm}italic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT on the equatorial plane of Eq. (16)) are shown in Figs (2).

The red curve in Figs (2), which represents the BH horizon for a=0𝑎0a=0italic_a = 0 (and the ergosurfaces Eq. (16) on the equatorial plane for a0𝑎0a\neq 0italic_a ≠ 0), bounds the collections of horizons at different a=𝑎absenta=italic_a =constant in the plane (r/M,a/M)𝑟𝑀𝑎𝑀(r/M,a/M)( italic_r / italic_M , italic_a / italic_M ). The PFDM metric describes solutions with 00,1111, and 2222 horizons. For BH solutions, horizons can be shifted outwardly or inwardly with respect to the Kerr BH spacetime depending on the value of k𝑘kitalic_k. There are also spacetime solutions with horizons for a>M𝑎𝑀a>Mitalic_a > italic_M.

Accordingly we select the PFDM metric parameter in the following six cases

𝐤𝟔{ka=2M,kb1.399M,kc=0.1M,kd=0.1M,ke0.82M,kf2M}.subscript𝐤6formulae-sequencesubscript𝑘𝑎2𝑀formulae-sequencesubscript𝑘𝑏1.399𝑀formulae-sequencesubscript𝑘𝑐0.1𝑀formulae-sequencesubscript𝑘𝑑0.1𝑀formulae-sequencesubscript𝑘𝑒0.82𝑀subscript𝑘𝑓2𝑀\displaystyle\mathbf{k_{6}}\equiv\{k_{a}=-2M,k_{b}\approx-1.399M,k_{c}=-0.1M,k% _{d}=0.1M,k_{e}\approx 0.82M,k_{f}\equiv 2M\}.bold_k start_POSTSUBSCRIPT bold_6 end_POSTSUBSCRIPT ≡ { italic_k start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = - 2 italic_M , italic_k start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ≈ - 1.399 italic_M , italic_k start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = - 0.1 italic_M , italic_k start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT = 0.1 italic_M , italic_k start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ≈ 0.82 italic_M , italic_k start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT ≡ 2 italic_M } . (17)
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 2: Horizons and ergosurfaces of the BHs in perfect fluid dark matter (PFDM) of Eq. (14). Different values of the k]7.18M,2M[k\in]-7.18M,2M[italic_k ∈ ] - 7.18 italic_M , 2 italic_M [ parameter, describing the PFDM intensity, are considered. (For k=0𝑘0k=0italic_k = 0 the line element describes the Schwarzschild and Kerr geometries in absence of DM.) Upper left panel: horizons r±subscript𝑟plus-or-minusr_{\pm}italic_r start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT as functions of the k𝑘kitalic_k parameter, for different values of the spin a0𝑎0a\geq 0italic_a ≥ 0. Curves at aMless-than-or-equals-or-greater-than𝑎𝑀a\lesseqgtr Mitalic_a ⋚ italic_M are shown. Red curves are the horizons of the static case (a=0𝑎0a=0italic_a = 0) (and the ergosurfaces rϵ±superscriptsubscript𝑟italic-ϵplus-or-minusr_{\epsilon}^{\pm}italic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT Eqs (16) on the equatorial plane for a>0𝑎0a>0italic_a > 0). Values 𝐤𝟔subscript𝐤6\mathbf{k_{6}}bold_k start_POSTSUBSCRIPT bold_6 end_POSTSUBSCRIPT of Eq. (17) are signed as dotted vertical lines. Yellow vertical line labels the Schwarzschild and Kerr geometries in absence of DM. Upper right panel: the horizons r±/Msubscript𝑟plus-or-minus𝑀r_{\pm}/Mitalic_r start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT / italic_M as functions of the spin a/M𝑎𝑀a/Mitalic_a / italic_M for different values of k>0𝑘0k>0italic_k > 0 (solid curves) and k<0𝑘0k<0italic_k < 0 (dashed curves). There is k=±0.06M𝑘plus-or-minus0.06𝑀k=\pm 0.06Mitalic_k = ± 0.06 italic_M (green), k=±0.006M𝑘plus-or-minus0.006𝑀k=\pm 0.006Mitalic_k = ± 0.006 italic_M (brown), k=±0.6M𝑘plus-or-minus0.6𝑀k=\pm 0.6Mitalic_k = ± 0.6 italic_M (cyan), k=±1.6M𝑘plus-or-minus1.6𝑀k=\pm 1.6Mitalic_k = ± 1.6 italic_M (black), k=1.2M,0.32M𝑘1.2𝑀0.32𝑀k=-1.2M,0.32Mitalic_k = - 1.2 italic_M , 0.32 italic_M (purple), k=1.999M,2M𝑘1.999𝑀2𝑀k=-1.999M,2Mitalic_k = - 1.999 italic_M , 2 italic_M (blue). Curve k=0𝑘0k=0italic_k = 0 (yellow) represents the Schwarzschild and Kerr geometries in absence of DM. Bottom right and left panels: different views of horizons r±subscript𝑟plus-or-minusr_{\pm}italic_r start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT as functions of k/M𝑘𝑀k/Mitalic_k / italic_M and a/M𝑎𝑀a/Mitalic_a / italic_M, mesh-functions are curves with constant r/M𝑟𝑀r/Mitalic_r / italic_M (solid green), k/M𝑘𝑀k/Mitalic_k / italic_M (dotted gray), and a/M𝑎𝑀a/Mitalic_a / italic_M (solid red). Yellow vertical plane labels the case k=0𝑘0k=0italic_k = 0. Function contour colors are according to the a/M𝑎𝑀a/Mitalic_a / italic_M values (left panel) and k/M𝑘𝑀k/Mitalic_k / italic_M values (right panel). We have marked with a black spot the extreme Kerr BH horizon (a=M,r=M)formulae-sequence𝑎𝑀𝑟𝑀(a=M,r=M)( italic_a = italic_M , italic_r = italic_M ) on all panels. To highlight better this point plane a=M𝑎𝑀a=Mitalic_a = italic_M (blue) is included in the bottom panels (the crossing of a=M𝑎𝑀a=Mitalic_a = italic_M and k=0𝑘0k=0italic_k = 0 planes occurs at the point a=M,k=0,r=Mformulae-sequence𝑎𝑀formulae-sequence𝑘0𝑟𝑀a=M,k=0,r=Mitalic_a = italic_M , italic_k = 0 , italic_r = italic_M).
Refer to caption
Refer to caption
Figure 3: Equatorial circular geodesic structures of the Kerr attractor in perfect fluid dark matter (PFDM) of Eq. (14). Different values of the k]7.18M,2M[k\in]-7.18M,2M[italic_k ∈ ] - 7.18 italic_M , 2 italic_M [, describing the PFDM intensity are considered, where for k=0𝑘0k=0italic_k = 0 (vertical yellow line) the line element describes the Schwarzchild (for a=0𝑎0a=0italic_a = 0) and extreme Kerr BH geometry (for a=M𝑎𝑀a=Mitalic_a = italic_M). Values 𝐤𝟔subscript𝐤6\mathbf{k_{6}}bold_k start_POSTSUBSCRIPT bold_6 end_POSTSUBSCRIPT of Eq. (17) are signed with dotted vertical lines. Right panel shows the situation for a=M𝑎𝑀a=Mitalic_a = italic_M, left panel is for a=0𝑎0a=0italic_a = 0. Radii r±subscript𝑟plus-or-minusr_{\pm}italic_r start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT are the horizons, rϵ±superscriptsubscript𝑟italic-ϵplus-or-minusr_{\epsilon}^{\pm}italic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT are the ergosurfaces on the equatorial plane (coincident with the horizons of the static case). mso𝑚𝑠𝑜msoitalic_m italic_s italic_o is for marginally stable orbit, mbo𝑚𝑏𝑜mboitalic_m italic_b italic_o is for the marginally bounded orbit, mco𝑚𝑐𝑜mcoitalic_m italic_c italic_o is for marginally circular orbit, for co-rotating fluids (=superscript\ell=\ell^{-}roman_ℓ = roman_ℓ start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT, solid curves) and counter-rotating fluids (=+superscript\ell=\ell^{+}roman_ℓ = roman_ℓ start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT, dashed curves). Horizontal gray lines show the radii of the geodesic structure for the Kerr and Schwarzschild BH geometry in absence of dark matter (k=0𝑘0k=0italic_k = 0) for a=0𝑎0a=0italic_a = 0 and a=M𝑎𝑀a=Mitalic_a = italic_M respectively. See also Figs (1).
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 4: Fluid specific angular momentum ±superscriptplus-or-minus\ell^{\pm}roman_ℓ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT, energy parameter K±superscript𝐾plus-or-minusK^{\pm}italic_K start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT and (test particles) Keplerian angular momentum ±superscriptplus-or-minus\mathcal{L}^{\pm}caligraphic_L start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT as function of r/M𝑟𝑀r/Mitalic_r / italic_M for co-rotating and counter-rotating fluids, different perfect fluid dark matter (PFDM) parameters of Eq. (14) signed on the panels. Different values of the k]7.18M,2M[k\in]-7.18M,2M[italic_k ∈ ] - 7.18 italic_M , 2 italic_M [ parameter, describing the PFDM intensity, are considered. Columns are a=0𝑎0a=0italic_a = 0 (left) and a=M𝑎𝑀a=Mitalic_a = italic_M (right), and rows are ±superscriptplus-or-minus\ell^{\pm}roman_ℓ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT (top), K±superscript𝐾plus-or-minusK^{\pm}italic_K start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT (middle), ±superscriptplus-or-minus\mathcal{L}^{\pm}caligraphic_L start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT (bottom). (For k=0𝑘0k=0italic_k = 0 the line element describes the Schwarzschild or the Kerr geometry.) Values 𝐤𝟔subscript𝐤6\mathbf{k_{6}}bold_k start_POSTSUBSCRIPT bold_6 end_POSTSUBSCRIPT are in Eq. (17).

The cases (ka,kf)subscript𝑘𝑎subscript𝑘𝑓(k_{a},k_{f})( italic_k start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT , italic_k start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT ) have as horizon r=2M𝑟2𝑀r=2Mitalic_r = 2 italic_M for a=0𝑎0a=0italic_a = 0 (the static case) therefore, in this sense, these solutions can be compared to the Schwarzschild case. Similarly, for k=kb𝑘subscript𝑘𝑏k=k_{b}italic_k = italic_k start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT, the geometry with a=M𝑎𝑀a=Mitalic_a = italic_M has one horizon at r=2M𝑟2𝑀r=2Mitalic_r = 2 italic_M. There is one horizon when k=ke𝑘subscript𝑘𝑒k=k_{e}italic_k = italic_k start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT and a=M𝑎𝑀a=Mitalic_a = italic_M171717From Figs (3) we note that there is one horizon at r1.583M𝑟1.583𝑀r\approx 1.583Mitalic_r ≈ 1.583 italic_M for the special parameters (a1.25776M,k0.595M)formulae-sequencegreater-than-or-approximately-equals𝑎1.25776𝑀𝑘0.595𝑀(a\gtrapprox 1.25776M,k\approx-0.595M)( italic_a ⪆ 1.25776 italic_M , italic_k ≈ - 0.595 italic_M ) and at r0.73M𝑟0.73𝑀r\approx 0.73Mitalic_r ≈ 0.73 italic_M for (a0.855M,k0.27M)formulae-sequence𝑎0.855𝑀𝑘0.27𝑀(a\approx 0.855M,k\approx 0.27M)( italic_a ≈ 0.855 italic_M , italic_k ≈ 0.27 italic_M )..

In Figs (4) we show also the fluid specific angular momentum distribution for the 𝐤𝟔subscript𝐤6\mathbf{k_{6}}bold_k start_POSTSUBSCRIPT bold_6 end_POSTSUBSCRIPT parameters, for the cases a=0𝑎0a=0italic_a = 0 and a=M𝑎𝑀a=Mitalic_a = italic_M, compared to the distribution on the geometry in absence of DM (see also Figs (1)), the associated K𝐾Kitalic_K parameter and the Keplerian (test particle) angular momentum ±±K±superscriptplus-or-minussuperscriptplus-or-minussuperscript𝐾plus-or-minus\mathcal{L}^{\pm}\equiv\ell^{\pm}K^{\pm}caligraphic_L start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ≡ roman_ℓ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT italic_K start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT respectively, which we have defined as related to the thick tori counterparts from definition of K(r)Veff((r))𝐾𝑟subscript𝑉𝑒𝑓𝑓𝑟K(r)\equiv V_{eff}(\ell(r))italic_K ( italic_r ) ≡ italic_V start_POSTSUBSCRIPT italic_e italic_f italic_f end_POSTSUBSCRIPT ( roman_ℓ ( italic_r ) ), showing influence of the DM on the limiting thin Keplerian (geodesic) disk. The thin Keplerian disk is constrained by the geodesic structure of the gravitational background.

In general the geodesic structure is qualitatively similar for any spin–Figs (3). We focus on the two limiting cases a=0𝑎0a=0italic_a = 0 and a=M𝑎𝑀a=Mitalic_a = italic_M, whose equatorial circular geodesic structures for the orbiting configurations are represented in Figs (3). In this model DM couples with the BH rotation, entangling with the frame dragging, evidenced in the deformed rotational law ±superscriptplus-or-minus\ell^{\pm}roman_ℓ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT of the orbiting matter. According to the DM parameter k𝑘kitalic_k and spin a𝑎aitalic_a, relation (11) i.e. rmso±>rmbo±>rmco±superscriptsubscript𝑟𝑚𝑠𝑜plus-or-minussuperscriptsubscript𝑟𝑚𝑏𝑜plus-or-minussuperscriptsubscript𝑟𝑚𝑐𝑜plus-or-minusr_{mso}^{\pm}>r_{mbo}^{\pm}>r_{mco}^{\pm}italic_r start_POSTSUBSCRIPT italic_m italic_s italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT > italic_r start_POSTSUBSCRIPT italic_m italic_b italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT > italic_r start_POSTSUBSCRIPT italic_m italic_c italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT for =±superscriptplus-or-minus\ell=\ell^{\pm}roman_ℓ = roman_ℓ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT cases holds, similarly to the non deformed BH case. According to the discussion of Sec. (II.2.2) for the Kerr background, this relation, for small magnitude of k𝑘kitalic_k, reflects in the relative location of the maximum-minimum points of pressure of the orbiting disks. For values of k𝑘kitalic_k where this relation is not verified, for example for k>0𝑘0k>0italic_k > 0 of Figs (3), the orbiting toroidal structures show large qualitative divergences with respect to the accretion tori formation and dynamics in the Kerr BH spacetime.

–The static attractor (a=0𝑎0a=0italic_a = 0) The geodesic structure is represented in Figs (3)-left panel, compared with the Schwarzschild case. We note that for k>0𝑘0k>0italic_k > 0 the situation changes qualitatively with respect to the case in absence of DM. In Figs (4) is the fluid specific angular momentum ±superscriptplus-or-minus\ell^{\pm}roman_ℓ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT, the tori energy parameter K±superscript𝐾plus-or-minusK^{\pm}italic_K start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT and (test particles) Keplerian angular momentum ±superscriptplus-or-minus\mathcal{L}^{\pm}caligraphic_L start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT as function of r/M𝑟𝑀r/Mitalic_r / italic_M, for different PFDM parameters of Eq. (14), compared with the case k=0𝑘0k=0italic_k = 0 describing the Schwarzschild geometry. Notably, for (kf,ke)subscript𝑘𝑓subscript𝑘𝑒(k_{f},k_{e})( italic_k start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT , italic_k start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ), there is K(r)>1𝐾𝑟1K(r)>1italic_K ( italic_r ) > 1. From the (±,K±)superscriptplus-or-minussuperscript𝐾plus-or-minus(\ell^{\pm},K^{\pm})( roman_ℓ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT , italic_K start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ) analysis it is noted how, for some values of k𝑘kitalic_k, curves are lower with respect to the corresponding curves in absence of DM–Figs (3). From the analysis of the curves L±superscript𝐿plus-or-minusL^{\pm}italic_L start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT, we can note how in some case the test particle angular momentum is qualitatively different from the corresponding in the Schwarzschild case. Tori and effective potentials in this case are represented in Figs (5), constrained by the geodesic structure of Figs (3).

–The spinning attractor (a>0𝑎0a>0italic_a > 0) The geodesic structure for the spinning attractor geometry in PFDM is in Figs (3)–right panel. The fluid specific angular momentum ±superscriptplus-or-minus\ell^{\pm}roman_ℓ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT, energy parameter K±superscript𝐾plus-or-minusK^{\pm}italic_K start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT and (test particles) Keplerian angular momentum ±superscriptplus-or-minus\mathcal{L}^{\pm}caligraphic_L start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT are shown in Figs (4) as functions of r/M𝑟𝑀r/Mitalic_r / italic_M for co-rotating and counter-rotating fluids and PFDM parameters of Eq. (14), compared to the case k=0𝑘0k=0italic_k = 0 describing the Kerr geometry in absence of DM. Tori effective potentials for the counter-rotating (co-rotating) fluids are in Figs (6). Tori are in Figs (7) and Figs (8).

As clear from Figs (3) the geodesic structure in the DM geometry with a=0𝑎0a=0italic_a = 0 is similar to the counter-rotating geodesic structure in the axially symmetric spacetime (a0)𝑎0(a\neq 0)( italic_a ≠ 0 ). The co-rotating case (for a=M𝑎𝑀a=Mitalic_a = italic_M), =superscript\ell=\ell^{-}roman_ℓ = roman_ℓ start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT, especially for k>0𝑘0k>0italic_k > 0, is remarkably different from the geodesic structure in the geometries with a=0𝑎0a=0italic_a = 0 and it is further complicated by the presence of the ergoregion deformed by the PFDM with, however, quantitative discrepancies with respect to the Kerr spacetime for large part of the DM parameter values k<0𝑘0k<0italic_k < 0. For larger k>0𝑘0k>0italic_k > 0 there are also NS solutions (also for a]0,M]a\in]0,M]italic_a ∈ ] 0 , italic_M ]) and, for even large values of k>0𝑘0k>0italic_k > 0, there are BH solutions (also for a>M𝑎𝑀a>Mitalic_a > italic_M).

Refer to caption
Refer to caption
Figure 5: Case a=0𝑎0a=0italic_a = 0. Tori orbiting BHs in perfect fluid dark matter (PFDM) of Eq. (14). Right panel shows the tori for selected values of the k/M𝑘𝑀k/Mitalic_k / italic_M parameter, cusp location r×subscript𝑟r_{\times}italic_r start_POSTSUBSCRIPT × end_POSTSUBSCRIPT and fluid specific angular momentum \ellroman_ℓ, signed on the curves. Left panel shows the associated tori effective potentials. Values 𝐤𝟔subscript𝐤6\mathbf{k_{6}}bold_k start_POSTSUBSCRIPT bold_6 end_POSTSUBSCRIPT of Eq. (17) are considered. (For k=0𝑘0k=0italic_k = 0 the line element describes the Schwarzschild geometry.) There is r=x2+y2𝑟superscript𝑥2superscript𝑦2r=\sqrt{x^{2}+y^{2}}italic_r = square-root start_ARG italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG and σ=y2/(x2+y2)𝜎superscript𝑦2superscript𝑥2superscript𝑦2\sigma=y^{2}/(x^{2}+y^{2})italic_σ = italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / ( italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ), where σsin2θ𝜎superscript2𝜃\sigma\equiv\sin^{2}\thetaitalic_σ ≡ roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ.
Refer to caption
Refer to caption
Figure 6: Case a=M𝑎𝑀a=Mitalic_a = italic_M. Effective potentials for the tori orbiting BHs in perfect fluid dark matter (PFDM) of Eq. (14), for different values of the parameter k]7.18M,2M[k\in]-7.18M,2M[italic_k ∈ ] - 7.18 italic_M , 2 italic_M [, describing the PFDM intensity and fluid angular momentum \ellroman_ℓ, signed on the curves. Values 𝐤𝟔subscript𝐤6\mathbf{k_{6}}bold_k start_POSTSUBSCRIPT bold_6 end_POSTSUBSCRIPT of Eq. (17) are considered. (For k=0𝑘0k=0italic_k = 0 the line element describes the extreme Kerr BH geometry). Left (right) panel shows the effective potentials for the counter-rotating (co-rotating) fluids. There is r=x2+y2𝑟superscript𝑥2superscript𝑦2r=\sqrt{x^{2}+y^{2}}italic_r = square-root start_ARG italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG and σ=y2/(x2+y2)𝜎superscript𝑦2superscript𝑥2superscript𝑦2\sigma=y^{2}/(x^{2}+y^{2})italic_σ = italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / ( italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ), where σsin2θ𝜎superscript2𝜃\sigma\equiv\sin^{2}\thetaitalic_σ ≡ roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ. There are black curves for (k=0,=4.61,=2.14)formulae-sequence𝑘0formulae-sequence4.612.14(k=0,\ell=-4.61,\ell=2.14)( italic_k = 0 , roman_ℓ = - 4.61 , roman_ℓ = 2.14 ), blue curves for (kb,=9.79,=4.84)formulae-sequencesubscript𝑘𝑏9.794.84\left(k_{b},\ell=-9.79,\ell=4.84\right)( italic_k start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT , roman_ℓ = - 9.79 , roman_ℓ = 4.84 ), yellow curves for (ka,=10.8,=3.84)formulae-sequencesubscript𝑘𝑎10.83.84\left(k_{a},\ell=-10.8,\ell=3.84\right)( italic_k start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT , roman_ℓ = - 10.8 , roman_ℓ = 3.84 ), purple curves for (kd,=3.845,=6)formulae-sequencesubscript𝑘𝑑3.8456\left(k_{d},\ell=-3.845,\ell=6\right)( italic_k start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT , roman_ℓ = - 3.845 , roman_ℓ = 6 ), orange curves for (kc,=5.4,=3.05)formulae-sequencesubscript𝑘𝑐5.43.05\left(k_{c},\ell=-5.4,\ell=3.05\right)( italic_k start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT , roman_ℓ = - 5.4 , roman_ℓ = 3.05 ), green curves for (ke,=4,=1.83)formulae-sequencesubscript𝑘𝑒41.83\left(k_{e},\ell=-4,\ell=1.83\right)( italic_k start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT , roman_ℓ = - 4 , roman_ℓ = 1.83 ), red curves for (kf,=3,=2.63)formulae-sequencesubscript𝑘𝑓32.63\left(k_{f},\ell=-3,\ell=2.63\right)( italic_k start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT , roman_ℓ = - 3 , roman_ℓ = 2.63 ).
Refer to caption
Refer to caption
Figure 7: Case a=M𝑎𝑀a=Mitalic_a = italic_M. Effective potentials and tori orbiting BHs in perfect fluid dark matter (PFDM) of Eq. (14), for different values of the parameter k]7.18M,2M[k\in]-7.18M,2M[italic_k ∈ ] - 7.18 italic_M , 2 italic_M [, describing the PFDM intensity and fluid specific angular momentum \ellroman_ℓ signed on the curves. Co-rotating =>0superscript0\ell^{-}=\ell>0roman_ℓ start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT = roman_ℓ > 0 cases are represented. Values 𝐤𝟔subscript𝐤6\mathbf{k_{6}}bold_k start_POSTSUBSCRIPT bold_6 end_POSTSUBSCRIPT of Eq. (17) are considered. For k=0𝑘0k=0italic_k = 0 the line element describes the extreme Kerr BH geometry. There is r=x2+y2𝑟superscript𝑥2superscript𝑦2r=\sqrt{x^{2}+y^{2}}italic_r = square-root start_ARG italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG and σ=y2/(x2+y2)𝜎superscript𝑦2superscript𝑥2superscript𝑦2\sigma=y^{2}/(x^{2}+y^{2})italic_σ = italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / ( italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ), where σsin2θ𝜎superscript2𝜃\sigma\equiv\sin^{2}\thetaitalic_σ ≡ roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ.
Refer to caption
Refer to caption
Refer to caption
Figure 8: Case a=M𝑎𝑀a=Mitalic_a = italic_M. Effective potentials and tori orbiting BHs in perfect fluid dark matter (PFDM) of Eq. (14), for different values of the parameter k]7.18M,2M[k\in]-7.18M,2M[italic_k ∈ ] - 7.18 italic_M , 2 italic_M [, describing the PFDM intensity and fluid specific angular momentum \ellroman_ℓ, signed on the curves. Counter-rotating +=<0superscript0\ell^{+}=\ell<0roman_ℓ start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT = roman_ℓ < 0 cases are represented. Values 𝐤𝟔subscript𝐤6\mathbf{k_{6}}bold_k start_POSTSUBSCRIPT bold_6 end_POSTSUBSCRIPT of Eq. (17) are considered. For k=0𝑘0k=0italic_k = 0 the line element describes the extreme Kerr BH geometry. Inside plot in the right panel is an enlarged view of the inner green torus. There is r=x2+y2𝑟superscript𝑥2superscript𝑦2r=\sqrt{x^{2}+y^{2}}italic_r = square-root start_ARG italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG and σ=y2/(x2+y2)𝜎superscript𝑦2superscript𝑥2superscript𝑦2\sigma=y^{2}/(x^{2}+y^{2})italic_σ = italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / ( italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ), where σsin2θ𝜎superscript2𝜃\sigma\equiv\sin^{2}\thetaitalic_σ ≡ roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ.

In general DM influence manifests also with the existence of extremely large cusped tori located considerably far from the central singularity as clear comparing the geodesic structures in Figs (1) and Figs (3). In some cases, as clear from Figs (7), there are double configurations at equal \ellroman_ℓ (purple and blue curves of Figs (7) and green curve of Figs (8)), with considerable larger tori with respect to the case in absence of dark matter. Furthermore we note the presence of outer cusps (blue curve in Figs (7)) or possibly the emergence of double cusps also in presence of BH solutions, where the co-rotating marginally stable orbit shows some remarkable peculiarities at kke𝑘subscript𝑘𝑒k\geq k_{e}italic_k ≥ italic_k start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT for a=M𝑎𝑀a=Mitalic_a = italic_M and k[kd,ke]𝑘subscript𝑘𝑑subscript𝑘𝑒k\in[k_{d},k_{e}]italic_k ∈ [ italic_k start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT , italic_k start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ] for a=0𝑎0a=0italic_a = 0.

The inter-disk cusp181818Similarly to the outer cusps, the inter–disks cusp is a tori cusp located between two configurations having same (,K)𝐾(\ell,K)( roman_ℓ , italic_K ) parameters, which eventually could be interpreted as an excretion cusp characterizing some cosmological models–see for example the double separated configurations (purple curves) or the blue curv in Figs (7). (see also Figs (7)) can evolve, following the change in one or two of the tori parameters (,K)𝐾(\ell,K)( roman_ℓ , italic_K ) in an inner cusp followed by an inner configuration (as the green curve in Figs (8)) or two separated configurations (as the purple curve in Figs (7)). For k=ke𝑘subscript𝑘𝑒k=k_{e}italic_k = italic_k start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT, in the case a=M𝑎𝑀a=Mitalic_a = italic_M, there are no horizons consequently the geometry, although is not over-spinning, can be considered a naked singularity, and the green curve in Figs (8) can be seen as a typical double configuration characterizing certain NS geometries. Also the presence of excretion cusps could be a DM indicator. Notably these features are usually read as tracers for the possible NSs observations, emerging as consequences of the repulsive gravity effects characterizing NSs solutions, and in this sense the Kerr BH immersed in PFDM could be a "mimicker" of super-spinars solutions.

III.2 Cold and scalar field dark matter

In this section we consider a spinning BH in scalar field dark matter (SFDM), addressed in Sec. (III.2.1) and in cold dark matter (CDM), considered in Sec. (III.2.2)–see for example [8, 9].

III.2.1 Scalar field dark matter (SFDM)

There is

gtt=[12Mr+r2(1ξSFDM)Σ],gtϕ=aσ[2Mr+r2(1ξSFDM)]Σ,gϕϕσ[(a2+r2)2a2σΔSFDM]Σformulae-sequencesubscript𝑔𝑡𝑡delimited-[]12𝑀𝑟superscript𝑟21subscript𝜉𝑆𝐹𝐷𝑀Σformulae-sequencesubscript𝑔𝑡italic-ϕ𝑎𝜎delimited-[]2𝑀𝑟superscript𝑟21subscript𝜉𝑆𝐹𝐷𝑀Σsubscript𝑔italic-ϕitalic-ϕ𝜎delimited-[]superscriptsuperscript𝑎2superscript𝑟22superscript𝑎2𝜎subscriptΔ𝑆𝐹𝐷𝑀Σ\displaystyle g_{tt}=-\left[1-\frac{{2Mr}+r^{2}(1-\xi_{SFDM})}{\Sigma}\right],% \quad g_{t\phi}=-\frac{a\sigma\left[{2Mr}+r^{2}(1-\xi_{SFDM})\right]}{\Sigma},% \quad g_{\phi\phi}\equiv\frac{\sigma\left[\left(a^{2}+r^{2}\right)^{2}-a^{2}% \sigma\Delta_{SFDM}\right]}{\Sigma}italic_g start_POSTSUBSCRIPT italic_t italic_t end_POSTSUBSCRIPT = - [ 1 - divide start_ARG 2 italic_M italic_r + italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 1 - italic_ξ start_POSTSUBSCRIPT italic_S italic_F italic_D italic_M end_POSTSUBSCRIPT ) end_ARG start_ARG roman_Σ end_ARG ] , italic_g start_POSTSUBSCRIPT italic_t italic_ϕ end_POSTSUBSCRIPT = - divide start_ARG italic_a italic_σ [ 2 italic_M italic_r + italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 1 - italic_ξ start_POSTSUBSCRIPT italic_S italic_F italic_D italic_M end_POSTSUBSCRIPT ) ] end_ARG start_ARG roman_Σ end_ARG , italic_g start_POSTSUBSCRIPT italic_ϕ italic_ϕ end_POSTSUBSCRIPT ≡ divide start_ARG italic_σ [ ( italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_σ roman_Δ start_POSTSUBSCRIPT italic_S italic_F italic_D italic_M end_POSTSUBSCRIPT ] end_ARG start_ARG roman_Σ end_ARG
grrΣΔSFDM,gθθ=Σformulae-sequencesubscript𝑔𝑟𝑟ΣsubscriptΔ𝑆𝐹𝐷𝑀subscript𝑔𝜃𝜃Σ\displaystyle g_{rr}\equiv\frac{\Sigma}{\Delta_{SFDM}},\quad g_{\theta\theta}=\Sigmaitalic_g start_POSTSUBSCRIPT italic_r italic_r end_POSTSUBSCRIPT ≡ divide start_ARG roman_Σ end_ARG start_ARG roman_Δ start_POSTSUBSCRIPT italic_S italic_F italic_D italic_M end_POSTSUBSCRIPT end_ARG , italic_g start_POSTSUBSCRIPT italic_θ italic_θ end_POSTSUBSCRIPT = roman_Σ (18)

where

ΔSFDMa22Mr+r2ξSFDMandξSFDMexp[8ρcR2sin(πrR)π(πr)R].formulae-sequencesubscriptΔ𝑆𝐹𝐷𝑀superscript𝑎22𝑀𝑟superscript𝑟2subscript𝜉𝑆𝐹𝐷𝑀andsubscript𝜉𝑆𝐹𝐷𝑀8subscript𝜌𝑐superscript𝑅2𝜋𝑟𝑅𝜋𝜋𝑟𝑅\displaystyle\Delta_{SFDM}\equiv a^{2}-{2Mr}+r^{2}\xi_{SFDM}\quad\mbox{and}% \quad\xi_{SFDM}\equiv\exp\left[-\frac{8\rho_{c}R^{2}\sin\left(\frac{\pi r}{R}% \right)}{\frac{\pi(\pi r)}{R}}\right].roman_Δ start_POSTSUBSCRIPT italic_S italic_F italic_D italic_M end_POSTSUBSCRIPT ≡ italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 2 italic_M italic_r + italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ξ start_POSTSUBSCRIPT italic_S italic_F italic_D italic_M end_POSTSUBSCRIPT and italic_ξ start_POSTSUBSCRIPT italic_S italic_F italic_D italic_M end_POSTSUBSCRIPT ≡ roman_exp [ - divide start_ARG 8 italic_ρ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT italic_R start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_sin ( divide start_ARG italic_π italic_r end_ARG start_ARG italic_R end_ARG ) end_ARG start_ARG divide start_ARG italic_π ( italic_π italic_r ) end_ARG start_ARG italic_R end_ARG end_ARG ] . (19)

The metric components satisfy the asymptotically flatness condition, and the fluid potential is well defined at infinity (r+𝑟r\to+\inftyitalic_r → + ∞), where Veff=1subscript𝑉𝑒𝑓𝑓1V_{eff}=1italic_V start_POSTSUBSCRIPT italic_e italic_f italic_f end_POSTSUBSCRIPT = 1. We introduce the quantity k:ρc=k/R3:𝑘subscript𝜌𝑐𝑘superscript𝑅3k:\rho_{c}={k}/{R^{3}}italic_k : italic_ρ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = italic_k / italic_R start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT. Here ρcsubscript𝜌𝑐\rho_{c}italic_ρ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT is the central density and R𝑅Ritalic_R is the radius at which the pressure and density are zero191919In the SFDM static solution, the Klein-Gordon equation and a quadratic potential for the scalar field have been considered [8, 9]. (where ρc=0subscript𝜌𝑐0\rho_{c}=0italic_ρ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = 0 the metric reduces to the Kerr solution). The Kerr limit occurs for k0𝑘0k\to 0italic_k → 0 or R+𝑅R\to+\inftyitalic_R → + ∞. The zeros of ΔSFDMsubscriptΔ𝑆𝐹𝐷𝑀\Delta_{SFDM}roman_Δ start_POSTSUBSCRIPT italic_S italic_F italic_D italic_M end_POSTSUBSCRIPT distinguish the metric singularities r±subscript𝑟plus-or-minusr_{\pm}italic_r start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT (SFDM deformed horizons), while the zeros of gttsubscript𝑔𝑡𝑡g_{tt}italic_g start_POSTSUBSCRIPT italic_t italic_t end_POSTSUBSCRIPT define a deformation of the Kerr ergosurfaces rϵ±superscriptsubscript𝑟italic-ϵplus-or-minusr_{\epsilon}^{\pm}italic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT , which can be found by solving the equation a=a±𝑎subscript𝑎plus-or-minusa=a_{\pm}italic_a = italic_a start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT and a=aϵ±𝑎superscriptsubscript𝑎italic-ϵplus-or-minusa=a_{\epsilon}^{\pm}italic_a = italic_a start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT respectively202020The accretion tori considered here are geometrically thick and characterized by a pronounced verticality. The tori surfaces can therefore approach the outer ergosurface out of the equatorial plane (i.e. σ<1𝜎1\sigma<1italic_σ < 1). The ergosurface location (on planes σ1𝜎1\sigma\neq 1italic_σ ≠ 1) is an important factor regulating the Lense–Thirring effect on the disks and on the jets flows coming from the disks[77]., where

a±r(2Mrξsfdm1),aϵ±a±1σwithξsfdme8ksin(πrR)π2r.formulae-sequencesubscript𝑎plus-or-minus𝑟2𝑀𝑟superscriptsubscript𝜉𝑠𝑓𝑑𝑚1formulae-sequencesuperscriptsubscript𝑎italic-ϵplus-or-minussubscript𝑎plus-or-minus1𝜎withsubscript𝜉𝑠𝑓𝑑𝑚superscript𝑒8𝑘𝜋𝑟𝑅superscript𝜋2𝑟\displaystyle a_{\pm}\equiv\sqrt{r\left(2M-r\xi_{sfdm}^{-1}\right)},\quad a_{% \epsilon}^{\pm}\equiv\frac{a_{\pm}}{\sqrt{1-\sigma}}\quad\mbox{with}\quad\xi_{% sfdm}\equiv e^{\frac{8k\sin\left(\frac{\pi r}{R}\right)}{\pi^{2}r}}.italic_a start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ≡ square-root start_ARG italic_r ( 2 italic_M - italic_r italic_ξ start_POSTSUBSCRIPT italic_s italic_f italic_d italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) end_ARG , italic_a start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ≡ divide start_ARG italic_a start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT end_ARG start_ARG square-root start_ARG 1 - italic_σ end_ARG end_ARG with italic_ξ start_POSTSUBSCRIPT italic_s italic_f italic_d italic_m end_POSTSUBSCRIPT ≡ italic_e start_POSTSUPERSCRIPT divide start_ARG 8 italic_k roman_sin ( divide start_ARG italic_π italic_r end_ARG start_ARG italic_R end_ARG ) end_ARG start_ARG italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_r end_ARG end_POSTSUPERSCRIPT . (20)

In accordance with the study of the BHs horizons in Figs (9) we explore the following two cases:

p{kp1000M,Rp1200M},g{kg20000M,Rg120000M},formulae-sequencesubscript𝑝formulae-sequencesubscript𝑘𝑝1000𝑀subscript𝑅𝑝1200𝑀subscript𝑔formulae-sequencesubscript𝑘𝑔20000𝑀subscript𝑅𝑔120000𝑀\displaystyle\mathcal{F}_{p}\equiv\{k_{p}\to 1000M,R_{p}\to 1200M\},\quad% \mathcal{F}_{g}\equiv\{k_{g}\to 20000M,R_{g}\to 120000M\},caligraphic_F start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ≡ { italic_k start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT → 1000 italic_M , italic_R start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT → 1200 italic_M } , caligraphic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ≡ { italic_k start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT → 20000 italic_M , italic_R start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT → 120000 italic_M } , (21)

(see Figs (9)).

Refer to caption
Refer to caption
Refer to caption
Figure 9: Horizons and ergosurfaces of BHs geometries in scalar field dark matter (SFDM) of Eq. (19). Parameters psubscript𝑝\mathcal{F}_{p}caligraphic_F start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT and gsubscript𝑔\mathcal{F}_{g}caligraphic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT are defined in Eqs (21). The case k=0𝑘0k=0italic_k = 0 corresponds to Kerr or Schwarzchild geometries. Left panel shows the BH horizons r±subscript𝑟plus-or-minusr_{\pm}italic_r start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT as functions of a/M𝑎𝑀a/Mitalic_a / italic_M for different DM parameters (k,R)𝑘𝑅(k,R)( italic_k , italic_R ). Gray curve is the Kerr BH horizons. Center panel: horizons radii r±subscript𝑟plus-or-minusr_{\pm}italic_r start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT as functions of the DM parameter R𝑅Ritalic_R, for different k𝑘kitalic_k and spin a𝑎aitalic_a, red curve, for a=0𝑎0a=0italic_a = 0, is the static limiting case and a=M𝑎𝑀a=Mitalic_a = italic_M is the blue curve, a=0.5M𝑎0.5𝑀a=0.5Mitalic_a = 0.5 italic_M is the green curve. Dashed curves correspond to kgsubscript𝑘𝑔k_{g}italic_k start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT and solid curves to kpsubscript𝑘𝑝k_{p}italic_k start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT. Right panel shows the ergosurfaces as functions of R/M𝑅𝑀R/Mitalic_R / italic_M for different k/M𝑘𝑀k/Mitalic_k / italic_M, for spin a=0.99M𝑎0.99𝑀a=0.99Mitalic_a = 0.99 italic_M and different planes σsin2θ𝜎superscript2𝜃\sigma\equiv\sin^{2}\thetaitalic_σ ≡ roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ, σ=1𝜎1\sigma=1italic_σ = 1 (red curve) is the equatorial plane, green curve is σ=0.01𝜎0.01\sigma=0.01italic_σ = 0.01, blue curve is σ=0.1𝜎0.1\sigma=0.1italic_σ = 0.1, dashed curves correspond to kgsubscript𝑘𝑔k_{g}italic_k start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT and solid curves to kpsubscript𝑘𝑝k_{p}italic_k start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT defined in Eqs (21).
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 10: Geodesic equatorial structure of scalar field dark matter (SFDM) geometry of Eq. (19). Upper (bottom) panels show the situation for co-rotating (counter-rotating) fluids with superscript\ell^{-}roman_ℓ start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT (+superscript\ell^{+}roman_ℓ start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT) for parameters psubscript𝑝\mathcal{F}_{p}caligraphic_F start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT (purple curves) and gsubscript𝑔\mathcal{F}_{g}caligraphic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT (blue curves), defined in Eqs (21). Radius r+subscript𝑟r_{+}italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT (black curve) is the outer horizon, rϵ+superscriptsubscript𝑟italic-ϵr_{\epsilon}^{+}italic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT (gray curve) is the outer ergosurface on the equatorial plane, mso𝑚𝑠𝑜msoitalic_m italic_s italic_o (solid purple and blue curves) is for marginally stable orbit, mbo𝑚𝑏𝑜mboitalic_m italic_b italic_o (dashed curves) is for the marginally bounded orbits, mco𝑚𝑐𝑜mcoitalic_m italic_c italic_o (dotted curves) is for the marginally circular orbits. The correspondent Kerr geodesic structure is in Figs (1). Center panels are close-up views of the left panels.
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 11: Fluid specific angular momentum ±superscriptplus-or-minus\ell^{\pm}roman_ℓ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT, energy parameter K±superscript𝐾plus-or-minusK^{\pm}italic_K start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT and (test particles) Keplerian angular momentum ±superscriptplus-or-minus\mathcal{L}^{\pm}caligraphic_L start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT as function of r/M𝑟𝑀r/Mitalic_r / italic_M for co-rotating ()(-)( - ) and counter-rotating (+)(+)( + ) fluids and different scalar field dark matter (SFDM) parameters psubscript𝑝\mathcal{F}_{p}caligraphic_F start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT (purple) and gsubscript𝑔\mathcal{F}_{g}caligraphic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT (blue), defined in Eqs (21). The value k=0𝑘0k=0italic_k = 0 corresponds to the Schwarzschild (a=0𝑎0a=0italic_a = 0) extreme Kerr BH case (a=M𝑎𝑀a=Mitalic_a = italic_M). The columns are a=0𝑎0a=0italic_a = 0 (left) and a=M𝑎𝑀a=Mitalic_a = italic_M (right), and rows are for \ellroman_ℓ as function of r𝑟ritalic_r (top), K𝐾Kitalic_K as function of r𝑟ritalic_r (middle), \mathcal{L}caligraphic_L as function of r𝑟ritalic_r (bottom).
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 12: Fluids effective potentials (left panels) and tori (center and right panels) orbiting in scalar field dark matter (SFDM) geometry of Eq. (19), with parameters psubscript𝑝\mathcal{F}_{p}caligraphic_F start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT (purple curves) and gsubscript𝑔\mathcal{F}_{g}caligraphic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT (blue curves), defined in Eqs (21). The geodesic structure is in Figs (10). Rows are a=0𝑎0a=0italic_a = 0 (top), a=0.1M𝑎0.1𝑀a=0.1Mitalic_a = 0.1 italic_M (middle), a=0.9M𝑎0.9𝑀a=0.9Mitalic_a = 0.9 italic_M (bottom). Black curves for k=0𝑘0k=0italic_k = 0 are the configurations for the case of Schwarzschild and Kerr spacetimes in absence of DM. Tori are shown in the correspondent colors association relative to the effective potentials. Fluid specific angular momenta +superscript\ell^{+}roman_ℓ start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT (dashed curves) and superscript\ell^{-}roman_ℓ start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT (solid curves) are signed close to each tori surface (center and right panels). Right panels are an enlarged view of the center panels in the region close to the central attractor. (The integration is in the entire orbital range). There is r=x2+y2𝑟superscript𝑥2superscript𝑦2r=\sqrt{x^{2}+y^{2}}italic_r = square-root start_ARG italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG and σ=y2/(x2+y2)𝜎superscript𝑦2superscript𝑥2superscript𝑦2\sigma=y^{2}/(x^{2}+y^{2})italic_σ = italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / ( italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ), where σsin2θ𝜎superscript2𝜃\sigma\equiv\sin^{2}\thetaitalic_σ ≡ roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ.

However, as clear from Figs (9), BH horizons exist also for a>M𝑎𝑀a>Mitalic_a > italic_M.

The geometry circular geodesic structure is shown in Figs (10) for the parameter gsubscript𝑔\mathcal{F}_{g}caligraphic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT and psubscript𝑝\mathcal{F}_{p}caligraphic_F start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT, emphasizing the differences for the cases (kp,Rp)subscript𝑘𝑝subscript𝑅𝑝(k_{p},R_{p})( italic_k start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT , italic_R start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ) and (kg,Rg)subscript𝑘𝑔subscript𝑅𝑔(k_{g},R_{g})( italic_k start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , italic_R start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ) and the Kerr geometry circular structure of Figs (1). The orbital range locating the proto-jets cusps, bounded by the radii rmbo±superscriptsubscript𝑟𝑚𝑏𝑜plus-or-minusr_{mbo}^{\pm}italic_r start_POSTSUBSCRIPT italic_m italic_b italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT and rmco±superscriptsubscript𝑟𝑚𝑐𝑜plus-or-minusr_{mco}^{\pm}italic_r start_POSTSUBSCRIPT italic_m italic_c italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT, for fluids specific angular momentum =±superscriptplus-or-minus\ell=\ell^{\pm}roman_ℓ = roman_ℓ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT is very narrow, and in general the geodesic structure is shifted considerably outwardly with the respect to the Kerr geometry geodesic radii. Therefore the range for the location of the accreting disks inner edges is considerably larger that the proto-jets cusp range, constituting possibly a constraint on the formation of proto-jets and tori with large angular momentum magnitude. We note also that in the co-rotating case, for a]0,M]a\in]0,M]italic_a ∈ ] 0 , italic_M ], and for psubscript𝑝\mathcal{F}_{p}caligraphic_F start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT differently from the Kerr case in absence of DM, radii are located out of the ergoregion (and partially for gsubscript𝑔\mathcal{F}_{g}caligraphic_F start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT), this could imply a significant difference in the Lense-Thirring effects in presence of DM. (There can be however over-spinning BHs (with a>M𝑎𝑀a>Mitalic_a > italic_M) where these effects could be present). Differently from the PFDM case, qualitatively the geodesic structure is not differentiated with the respect to the case in absence of SFDM for co-rotating or counter-rotating fluids, and for static attractors (a=0𝑎0a=0italic_a = 0) and spinning attractors (a]0,M]a\in]0,M]italic_a ∈ ] 0 , italic_M ]) with SFDM.

We analyse below the case of the static attractor i.e. a=0𝑎0a=0italic_a = 0, having as limit in the vacuum the Schwarzschild metric, we then studied the influence of the central attractor spin combined with dark matter effects, considering the case a=M𝑎𝑀a=Mitalic_a = italic_M, corresponding to the DM deformation of the vacuum solution (i.e. in absence of DM) of the extreme Kerr BH, and the case of the slowly spinning attractor having a=0.1M𝑎0.1𝑀a=0.1Mitalic_a = 0.1 italic_M. More specifically:

–The static attractor (a=0𝑎0a=0italic_a = 0)

In Figs (11) the fluid specific angular momentum ±superscriptplus-or-minus\ell^{\pm}roman_ℓ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT, tori energy parameter K±superscript𝐾plus-or-minusK^{\pm}italic_K start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT and (test particles) Keplerian angular momentum ±superscriptplus-or-minus\mathcal{L}^{\pm}caligraphic_L start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT are shown as functions of r/M𝑟𝑀r/Mitalic_r / italic_M for co-rotating and counter-rotating fluids, at different spins and SFDM parameters, compared to the case k=0𝑘0k=0italic_k = 0 corresponding to the Schwarzschild geometry. Tori and effective potentials are in Figs (12).

–The spinning attractors (a=M𝑎𝑀a=Mitalic_a = italic_M and a=0.1M𝑎0.1𝑀a=0.1Mitalic_a = 0.1 italic_M)

We restrict our analysis to a=M𝑎𝑀a=Mitalic_a = italic_M and a=0.1M𝑎0.1𝑀a=0.1Mitalic_a = 0.1 italic_M, studying the cusped tori limiting the closed configurations, regulated by the effective potential function. In Figs (11), there are the fluid specific angular momentum ±superscriptplus-or-minus\ell^{\pm}roman_ℓ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT, energy parameter K±superscript𝐾plus-or-minusK^{\pm}italic_K start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT and (test particles) Keplerian angular momentum ±superscriptplus-or-minus\mathcal{L}^{\pm}caligraphic_L start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT as function of r/M𝑟𝑀r/Mitalic_r / italic_M for co-rotating and counter-rotating fluids, different spins and SFDM parameters of Eq. (19), with respect to the Kerr vacuum cases. In Figs (12) are the fluids effective potentials and tori compared to the case of Kerr in absence of DM.

From Figs (11) we note that, within this parameter choice, differently from the case in absence of DM, the fluids energy function can not converge to 1111 for large values of r𝑟ritalic_r. In Figs (11) the fluid specific angular momentum distribution, compared to the distribution on the geometry in absence of DM, the associated K𝐾Kitalic_K energy parameter and the Keplerian (test particle) angular momentum ±±K±superscriptplus-or-minussuperscriptplus-or-minussuperscript𝐾plus-or-minus\mathcal{L}^{\pm}\equiv\ell^{\pm}K^{\pm}caligraphic_L start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ≡ roman_ℓ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT italic_K start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT are shown. From Figs (9) it is clear how the horizon curves in the plane aM𝑎𝑀a-Mitalic_a - italic_M are larger and shifted outwardly with respect to the Kerr BH case, constituting a discriminant for the SFDM model and, for some values of the DM parameters, the BH horizons disappear giving rise to a "DM–induced" NS.

Large tori orbiting SFDM spinning BHs are shown in Figs (12) as, for example, the purple surface for the case a=0𝑎0a=0italic_a = 0 (from the tori effective potentials we can also note how the tori K𝐾Kitalic_K parameter for tori orbiting in SFDM are generally considerably lower than the K𝐾Kitalic_K parameter in absence of DM).

This feature of the DM model could also be an indication that such extremely huge tori are actually not formed and similarly the back-reaction on the metric is a predominant factor in these configurations, where self-gravity becomes a determinant factor in the tori equilibrium.

III.2.2 Cold Dark Matter (CDM)

The metric components read

gtt=[12MrξCDM+r2Σ],gtϕ=aσ[2MrξCDM+r2]Σ,gϕϕσ[(a2+r2)2a2σΔCDM]Σ;formulae-sequencesubscript𝑔𝑡𝑡delimited-[]12𝑀𝑟subscript𝜉𝐶𝐷𝑀superscript𝑟2Σformulae-sequencesubscript𝑔𝑡italic-ϕ𝑎𝜎delimited-[]2𝑀𝑟subscript𝜉𝐶𝐷𝑀superscript𝑟2Σsubscript𝑔italic-ϕitalic-ϕ𝜎delimited-[]superscriptsuperscript𝑎2superscript𝑟22superscript𝑎2𝜎subscriptΔ𝐶𝐷𝑀Σ\displaystyle g_{tt}=-\left[1-\frac{{2Mr}-\xi_{CDM}+r^{2}}{\Sigma}\right],% \quad g_{t\phi}=-\frac{a\sigma\left[{2Mr}-\xi_{CDM}+r^{2}\right]}{\Sigma},% \quad g_{\phi\phi}\equiv\frac{\sigma\left[\left(a^{2}+r^{2}\right)^{2}-a^{2}% \sigma\Delta_{CDM}\right]}{\Sigma};italic_g start_POSTSUBSCRIPT italic_t italic_t end_POSTSUBSCRIPT = - [ 1 - divide start_ARG 2 italic_M italic_r - italic_ξ start_POSTSUBSCRIPT italic_C italic_D italic_M end_POSTSUBSCRIPT + italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG roman_Σ end_ARG ] , italic_g start_POSTSUBSCRIPT italic_t italic_ϕ end_POSTSUBSCRIPT = - divide start_ARG italic_a italic_σ [ 2 italic_M italic_r - italic_ξ start_POSTSUBSCRIPT italic_C italic_D italic_M end_POSTSUBSCRIPT + italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] end_ARG start_ARG roman_Σ end_ARG , italic_g start_POSTSUBSCRIPT italic_ϕ italic_ϕ end_POSTSUBSCRIPT ≡ divide start_ARG italic_σ [ ( italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_σ roman_Δ start_POSTSUBSCRIPT italic_C italic_D italic_M end_POSTSUBSCRIPT ] end_ARG start_ARG roman_Σ end_ARG ;
grrΣΔCDM,gθθ=Σformulae-sequencesubscript𝑔𝑟𝑟ΣsubscriptΔ𝐶𝐷𝑀subscript𝑔𝜃𝜃Σ\displaystyle g_{rr}\equiv\frac{\Sigma}{\Delta_{CDM}},\quad g_{\theta\theta}=\Sigmaitalic_g start_POSTSUBSCRIPT italic_r italic_r end_POSTSUBSCRIPT ≡ divide start_ARG roman_Σ end_ARG start_ARG roman_Δ start_POSTSUBSCRIPT italic_C italic_D italic_M end_POSTSUBSCRIPT end_ARG , italic_g start_POSTSUBSCRIPT italic_θ italic_θ end_POSTSUBSCRIPT = roman_Σ (22)

where

ΔCDMa22Mr+ξCDMandξCDMr2(rR+1)8πρcR3r.formulae-sequencesubscriptΔ𝐶𝐷𝑀superscript𝑎22𝑀𝑟subscript𝜉𝐶𝐷𝑀andsubscript𝜉𝐶𝐷𝑀superscript𝑟2superscript𝑟𝑅18𝜋subscript𝜌𝑐superscript𝑅3𝑟\displaystyle\Delta_{CDM}\equiv a^{2}-{2Mr}+\xi_{CDM}\quad\mbox{and}\quad\xi_{% CDM}\equiv r^{2}\left(\frac{r}{R}+1\right)^{-\frac{8\pi\rho_{c}R^{3}}{r}}.roman_Δ start_POSTSUBSCRIPT italic_C italic_D italic_M end_POSTSUBSCRIPT ≡ italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 2 italic_M italic_r + italic_ξ start_POSTSUBSCRIPT italic_C italic_D italic_M end_POSTSUBSCRIPT and italic_ξ start_POSTSUBSCRIPT italic_C italic_D italic_M end_POSTSUBSCRIPT ≡ italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG italic_r end_ARG start_ARG italic_R end_ARG + 1 ) start_POSTSUPERSCRIPT - divide start_ARG 8 italic_π italic_ρ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT italic_R start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG start_ARG italic_r end_ARG end_POSTSUPERSCRIPT . (23)

We adopt the parametrization ρc=k/R3subscript𝜌𝑐𝑘superscript𝑅3\rho_{c}={k}/{R^{3}}italic_ρ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = italic_k / italic_R start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT, where ρcsubscript𝜌𝑐\rho_{c}italic_ρ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT is the density of the universe at the moment when the DM halo collapsed, R𝑅Ritalic_R is a characteristic radius. The metric is asymptotically flat and we find the Kerr limit in k0𝑘0k\to 0italic_k → 0 or R+𝑅R\to+\inftyitalic_R → + ∞. We first consider the metric singularities, identifying the space of the parameters used in the tori analysis. The horizons r±subscript𝑟plus-or-minusr_{\pm}italic_r start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT can be written as solutions of the equation a=a±(CDM)𝑎subscript𝑎plus-or-minus𝐶𝐷𝑀a=a_{\pm}(CDM)italic_a = italic_a start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ( italic_C italic_D italic_M ), or R=R±(CDM)𝑅subscript𝑅plus-or-minus𝐶𝐷𝑀R=R_{\pm}(CDM)italic_R = italic_R start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ( italic_C italic_D italic_M ) or k=k±(CDM)𝑘subscript𝑘plus-or-minus𝐶𝐷𝑀k=k_{\pm}(CDM)italic_k = italic_k start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ( italic_C italic_D italic_M ) where

a±(CDM)r[2Mr8πkr],R±(CDM)r[11ςr8πk1],k±(CDM)rlogς18πlog,formulae-sequencesubscript𝑎plus-or-minus𝐶𝐷𝑀𝑟delimited-[]2𝑀𝑟superscriptWeierstrass-p8𝜋𝑘𝑟formulae-sequencesubscript𝑅plus-or-minus𝐶𝐷𝑀𝑟delimited-[]11superscript𝜍𝑟8𝜋𝑘1subscript𝑘plus-or-minus𝐶𝐷𝑀𝑟superscript𝜍18𝜋Weierstrass-p\displaystyle a_{\pm}(CDM)\equiv\sqrt{r\left[2M-r\wp^{-\frac{8\pi k}{r}}\right% ]},\quad R_{\pm}(CDM)\equiv r\left[\frac{1}{1-\varsigma^{\frac{r}{8\pi k}}}-1% \right],\quad k_{\pm}(CDM)\equiv\frac{r\log\varsigma^{-1}}{8\pi\log\wp},italic_a start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ( italic_C italic_D italic_M ) ≡ square-root start_ARG italic_r [ 2 italic_M - italic_r ℘ start_POSTSUPERSCRIPT - divide start_ARG 8 italic_π italic_k end_ARG start_ARG italic_r end_ARG end_POSTSUPERSCRIPT ] end_ARG , italic_R start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ( italic_C italic_D italic_M ) ≡ italic_r [ divide start_ARG 1 end_ARG start_ARG 1 - italic_ς start_POSTSUPERSCRIPT divide start_ARG italic_r end_ARG start_ARG 8 italic_π italic_k end_ARG end_POSTSUPERSCRIPT end_ARG - 1 ] , italic_k start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ( italic_C italic_D italic_M ) ≡ divide start_ARG italic_r roman_log italic_ς start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT end_ARG start_ARG 8 italic_π roman_log ℘ end_ARG ,
withrR+1,andς2rMa2r2.formulae-sequencewithWeierstrass-p𝑟𝑅1and𝜍2𝑟𝑀superscript𝑎2superscript𝑟2\displaystyle\mbox{with}\quad\wp\equiv\frac{r}{R}+1,\quad\mbox{and}\quad% \varsigma\equiv\frac{2rM-a^{2}}{r^{2}}.with ℘ ≡ divide start_ARG italic_r end_ARG start_ARG italic_R end_ARG + 1 , and italic_ς ≡ divide start_ARG 2 italic_r italic_M - italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . (24)

The metric horizons are defined for

a]0,M]:r]r,r[;a[0,M]:r>r+;a>M:r>r,wherera22M,\displaystyle a\in]0,M]:\quad r\in\left]r_{\circledast},{r_{-}}\right[;\quad a% \in[0,M]:r>{r_{+}};\quad a>M:r>r_{\circledast},\quad\mbox{where}\quad r_{% \circledast}\equiv\frac{a^{2}}{2M},italic_a ∈ ] 0 , italic_M ] : italic_r ∈ ] italic_r start_POSTSUBSCRIPT ⊛ end_POSTSUBSCRIPT , italic_r start_POSTSUBSCRIPT - end_POSTSUBSCRIPT [ ; italic_a ∈ [ 0 , italic_M ] : italic_r > italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ; italic_a > italic_M : italic_r > italic_r start_POSTSUBSCRIPT ⊛ end_POSTSUBSCRIPT , where italic_r start_POSTSUBSCRIPT ⊛ end_POSTSUBSCRIPT ≡ divide start_ARG italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_M end_ARG ,

(assuming a>0,R>0,k>0formulae-sequence𝑎0formulae-sequence𝑅0𝑘0a>0,R>0,k>0italic_a > 0 , italic_R > 0 , italic_k > 0).

The ergosurfaces rϵ±superscriptsubscript𝑟italic-ϵplus-or-minusr_{\epsilon}^{\pm}italic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT can be found for σ1𝜎1\sigma\neq 1italic_σ ≠ 1 as solutions of the equation a=aϵ±(CDM)a±(CDM)/1σ𝑎superscriptsubscript𝑎italic-ϵplus-or-minus𝐶𝐷𝑀subscript𝑎plus-or-minus𝐶𝐷𝑀1𝜎a=a_{\epsilon}^{\pm}(CDM)\equiv{a_{\pm}(CDM)}/{\sqrt{1-\sigma}}italic_a = italic_a start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ( italic_C italic_D italic_M ) ≡ italic_a start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ( italic_C italic_D italic_M ) / square-root start_ARG 1 - italic_σ end_ARG and, on the equatorial plane (σ=1)𝜎1(\sigma=1)( italic_σ = 1 ) as solutions of R=Rϵ±𝑅superscriptsubscript𝑅italic-ϵplus-or-minusR=R_{\epsilon}^{\pm}italic_R = italic_R start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT or k=kϵ±𝑘superscriptsubscript𝑘italic-ϵplus-or-minusk=k_{\epsilon}^{\pm}italic_k = italic_k start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT where212121As clear from Figs (13), these relations also represent some portions of the ergosurfaces according to conditions on the DM parameters.

kϵ±rlogr2M8πlogandRϵ±r[11(2Mr)r8πk1],formulae-sequencesuperscriptsubscript𝑘italic-ϵplus-or-minus𝑟𝑟2𝑀8𝜋Weierstrass-pandsuperscriptsubscript𝑅italic-ϵplus-or-minus𝑟delimited-[]11superscript2𝑀𝑟𝑟8𝜋𝑘1\displaystyle k_{\epsilon}^{\pm}\equiv\frac{r\log\frac{r}{2M}}{8\pi\log\wp}% \quad\mbox{and}\quad R_{\epsilon}^{\pm}\equiv r\left[\frac{1}{1-\left(\frac{2M% }{r}\right)^{\frac{r}{8\pi k}}}-1\right],italic_k start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ≡ divide start_ARG italic_r roman_log divide start_ARG italic_r end_ARG start_ARG 2 italic_M end_ARG end_ARG start_ARG 8 italic_π roman_log ℘ end_ARG and italic_R start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ≡ italic_r [ divide start_ARG 1 end_ARG start_ARG 1 - ( divide start_ARG 2 italic_M end_ARG start_ARG italic_r end_ARG ) start_POSTSUPERSCRIPT divide start_ARG italic_r end_ARG start_ARG 8 italic_π italic_k end_ARG end_POSTSUPERSCRIPT end_ARG - 1 ] , (25)

see Figs (13). On the equatorial plane, the outer ergosurface, independent on the spin a𝑎aitalic_a, corresponds to the metric singularity in the static (a=0𝑎0a=0italic_a = 0) case–see Eqs (24).

We consider the following sets of parameters

1{k1000M,R120000M},2{k100M,R1200M},andk(p)100Mformulae-sequencesubscript1formulae-sequence𝑘1000𝑀𝑅120000𝑀formulae-sequencesubscript2formulae-sequence𝑘100𝑀𝑅1200𝑀andsubscript𝑘𝑝100𝑀\displaystyle\mathcal{F}_{1}\equiv\{k\to 1000M,R\to 120000M\},\quad\mathcal{F}% _{2}\equiv\{k\to 100M,R\to 1200M\},\quad\mbox{and}\quad k_{(p)}\equiv 100Mcaligraphic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ≡ { italic_k → 1000 italic_M , italic_R → 120000 italic_M } , caligraphic_F start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ≡ { italic_k → 100 italic_M , italic_R → 1200 italic_M } , and italic_k start_POSTSUBSCRIPT ( italic_p ) end_POSTSUBSCRIPT ≡ 100 italic_M (26)

–Figs (13).

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 13: Horizons r±subscript𝑟plus-or-minusr_{\pm}italic_r start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT and ergosufaces rϵ±superscriptsubscript𝑟italic-ϵplus-or-minusr_{\epsilon}^{\pm}italic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT of the cold dark matter CDM geometry of Eqs (22). Dark matter parameters 1subscript1\mathcal{F}_{1}caligraphic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, 2subscript2\mathcal{F}_{2}caligraphic_F start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT and k(p)subscript𝑘𝑝k_{(p)}italic_k start_POSTSUBSCRIPT ( italic_p ) end_POSTSUBSCRIPT are in Eqs (26), parameters (kp,kg,Rg,Rp)subscript𝑘𝑝subscript𝑘𝑔subscript𝑅𝑔subscript𝑅𝑝(k_{p},k_{g},R_{g},R_{p})( italic_k start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT , italic_k start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , italic_R start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , italic_R start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ) are in Eqs (21) . Black curves, k=0𝑘0k=0italic_k = 0, correspond to the Kerr horizons. There is σsin2θ𝜎superscript2𝜃\sigma\equiv\sin^{2}\thetaitalic_σ ≡ roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ, where σ=1𝜎1\sigma=1italic_σ = 1 is the equatorial plane. In the bottom left panel blue curves are a=M𝑎𝑀a=Mitalic_a = italic_M, red curves are a=0𝑎0a=0italic_a = 0, green curves are a=0.5M𝑎0.5𝑀a=0.5Mitalic_a = 0.5 italic_M, dotted curves are for k=k(p)𝑘subscript𝑘𝑝k=k_{(p)}italic_k = italic_k start_POSTSUBSCRIPT ( italic_p ) end_POSTSUBSCRIPT, dashed curves are for k=kg𝑘subscript𝑘𝑔k=k_{g}italic_k = italic_k start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT, solid curves are for k=kp𝑘subscript𝑘𝑝k=k_{p}italic_k = italic_k start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT, black horizontal lines correspond to k=0𝑘0k=0italic_k = 0 for the Kerr spacetimes. Bottom right panel shows the ergosurfaces rϵ±superscriptsubscript𝑟italic-ϵplus-or-minusr_{\epsilon}^{\pm}italic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT as function of the dark matter parameter R/M𝑅𝑀R/Mitalic_R / italic_M for spin a=0.99M𝑎0.99𝑀a=0.99Mitalic_a = 0.99 italic_M and for different planes σ[0,1]𝜎01\sigma\in[0,1]italic_σ ∈ [ 0 , 1 ] and dark matter parameter k𝑘kitalic_k, where red curves correspond to σ=1𝜎1\sigma=1italic_σ = 1, blue curves correspond to σ=0.1𝜎0.1\sigma=0.1italic_σ = 0.1, brown curves to σ=0.01𝜎0.01\sigma=0.01italic_σ = 0.01, dashed curves are for k=kg𝑘subscript𝑘𝑔k=k_{g}italic_k = italic_k start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT, solid curves for k=kp𝑘subscript𝑘𝑝k=k_{p}italic_k = italic_k start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT and dotted curves are for k=k(p)𝑘subscript𝑘𝑝k=k_{(p)}italic_k = italic_k start_POSTSUBSCRIPT ( italic_p ) end_POSTSUBSCRIPT.

The geodesic structure for this geometry is shown in Figs (14).

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 14: Geodesic equatorial structure of the cold dark matter CDM geometry of Eqs (22). Dark matter parameters 1subscript1\mathcal{F}_{1}caligraphic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT (green curves) and 2subscript2\mathcal{F}_{2}caligraphic_F start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT (red curves) are in Eqs (26). Upper (below) panels show the structure for fluid specific angular momentum =>0superscript0\ell=\ell^{-}>0roman_ℓ = roman_ℓ start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT > 0 (+<0superscript0\ell^{+}<0roman_ℓ start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT < 0) for co-rotating (counter-rotating) fluids. Center panels are an enlarged view of the right panels, showing the marginally stable orbits (mso𝑚𝑠𝑜msoitalic_m italic_s italic_o solid curves) marginally bounded orbits (mbo𝑚𝑏𝑜mboitalic_m italic_b italic_o dashed curves) and marginally circular orbits rmco±superscriptsubscript𝑟𝑚𝑐𝑜plus-or-minusr_{mco}^{\pm}italic_r start_POSTSUBSCRIPT italic_m italic_c italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT (dotted curves), where r+subscript𝑟r_{+}italic_r start_POSTSUBSCRIPT + end_POSTSUBSCRIPT (black curves) is the outer horizon, rϵ+superscriptsubscript𝑟italic-ϵr_{\epsilon}^{+}italic_r start_POSTSUBSCRIPT italic_ϵ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT (gray curves) is the outer ergosurface on the equatorial plane. The correspondent Kerr geodesic structure is in Figs (1).

Similarly to the SFDM model, the equatorial geodesic structure shows that for the 2subscript2\mathcal{F}_{2}caligraphic_F start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT case, the range [rmco±,rmbo±]superscriptsubscript𝑟𝑚𝑐𝑜plus-or-minussuperscriptsubscript𝑟𝑚𝑏𝑜plus-or-minus[r_{mco}^{\pm},r_{mbo}^{\pm}][ italic_r start_POSTSUBSCRIPT italic_m italic_c italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT , italic_r start_POSTSUBSCRIPT italic_m italic_b italic_o end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ] of the proto-jets cusp location is remarkably narrow. At a]0,M]a\in]0,M]italic_a ∈ ] 0 , italic_M ], for co-rotating fluids, the radii for 1subscript1\mathcal{F}_{1}caligraphic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT (2subscript2\mathcal{F}_{2}caligraphic_F start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT) do not enter (are partially contained in) the outer ergoregion.

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 15: Fluid specific angular momentum ±superscriptplus-or-minus\ell^{\pm}roman_ℓ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT (upper row), energy parameter K±superscript𝐾plus-or-minusK^{\pm}italic_K start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT (middle row) and (test particles) Keplerian angular momentum ±superscriptplus-or-minus\mathcal{L}^{\pm}caligraphic_L start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT (bottom row) as function of r/M𝑟𝑀r/Mitalic_r / italic_M for co-rotating and counter-rotating fluids, different cold dark matter CDM parameters of Eqs (22). Dark matter parameters 1subscript1\mathcal{F}_{1}caligraphic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT (purple curves) and 2subscript2\mathcal{F}_{2}caligraphic_F start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT (blue curves) are in Eqs (26). The columns are a=0𝑎0a=0italic_a = 0 (left) and a=M𝑎𝑀a=Mitalic_a = italic_M (right). Kerr and Schwarzschild case are the black curves for k=0𝑘0k=0italic_k = 0.
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 16: Effective potential and tori of the cold dark matter CDM geometry. Dark matter parameters 1subscript1\mathcal{F}_{1}caligraphic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT (purple curves) and 2subscript2\mathcal{F}_{2}caligraphic_F start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT (blue curves) are in Eqs (26). Rows are a=0𝑎0a=0italic_a = 0 (top line), a=0.7M𝑎0.7𝑀a=0.7Mitalic_a = 0.7 italic_M (top second and third lines), a=M𝑎𝑀a=Mitalic_a = italic_M (top fourth and fifth lines). Black curves for k=0𝑘0k=0italic_k = 0 are the configurations for the case of Schwarzschild and Kerr spacetimes in absence of DM. Tori are shown in the correspondent colors association relative to the effective potentials. Fluid specific angular momenta +superscript\ell^{+}roman_ℓ start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT (dashed curves) and superscript\ell^{-}roman_ℓ start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT (solid curves) are signed close to each tori surface. Regions close to the central attractor are also shown in enlarged views panels. There is r=x2+y2𝑟superscript𝑥2superscript𝑦2r=\sqrt{x^{2}+y^{2}}italic_r = square-root start_ARG italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG and σ=y2/(x2+y2)𝜎superscript𝑦2superscript𝑥2superscript𝑦2\sigma=y^{2}/(x^{2}+y^{2})italic_σ = italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / ( italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ), where σsin2θ𝜎superscript2𝜃\sigma\equiv\sin^{2}\thetaitalic_σ ≡ roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ.

We have analyzed the static and spinning attractors as follows:

–The static attractor (a=0𝑎0a=0italic_a = 0)

The fluid specific angular momentum ±superscriptplus-or-minus\ell^{\pm}roman_ℓ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT, energy parameter K±superscript𝐾plus-or-minusK^{\pm}italic_K start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT and (test particles) Keplerian angular momentum ±superscriptplus-or-minus\mathcal{L}^{\pm}caligraphic_L start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT as function of r/M𝑟𝑀r/Mitalic_r / italic_M, and CDM parameters 1subscript1\mathcal{F}_{1}caligraphic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and 2subscript2\mathcal{F}_{2}caligraphic_F start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT are in Figs (15), compared to the Schwarzschild case. The effective potential and tori are shown in Figs (16) compared with the Schwarzschild case.

–The spinning attractor (a0𝑎0a\neq 0italic_a ≠ 0)

Fluid specific angular momentum ±superscriptplus-or-minus\ell^{\pm}roman_ℓ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT, energy parameter K±superscript𝐾plus-or-minusK^{\pm}italic_K start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT and (test particles) Keplerian angular momentum ±superscriptplus-or-minus\mathcal{L}^{\pm}caligraphic_L start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT as function of r/M𝑟𝑀r/Mitalic_r / italic_M are represented in Figs (15) for co-rotating and counter-rotating fluids, different spins and CDM parameters 1subscript1\mathcal{F}_{1}caligraphic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT in comparison with the extreme Kerr BH case in absence of DM. Effective potential and tori for CDM parameters 1subscript1\mathcal{F}_{1}caligraphic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and 2subscript2\mathcal{F}_{2}caligraphic_F start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT are in Figs (16), for a=M𝑎𝑀a=Mitalic_a = italic_M and a=0.7M𝑎0.7𝑀a=0.7Mitalic_a = 0.7 italic_M in Figs (16), for fluid specific angular momentum =>0superscript0\ell=\ell^{-}>0roman_ℓ = roman_ℓ start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT > 0 (+<0superscript0\ell^{+}<0roman_ℓ start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT < 0) for co-rotating (counter-rotating) fluids, compared with the case in the vacuum Kerr geometry

It can be proved that, in all the cases considered, the limit K1𝐾1K\to 1italic_K → 1 for large r𝑟ritalic_r holds, where ±±K±superscriptplus-or-minussuperscriptplus-or-minussuperscript𝐾plus-or-minus\mathcal{L}^{\pm}\equiv\ell^{\pm}K^{\pm}caligraphic_L start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ≡ roman_ℓ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT italic_K start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT. Large tori orbiting CDM spinning BHs are shown in Figs (16) as, for example, the blue surface for the case a=0𝑎0a=0italic_a = 0, dashed -blue curve for a=0.7M𝑎0.7𝑀a=0.7Mitalic_a = 0.7 italic_M and a=M𝑎𝑀a=Mitalic_a = italic_M (similarly to the SFDM case, the tori K𝐾Kitalic_K parameter for tori orbiting in CDM are generally considerably lower than the K𝐾Kitalic_K parameter in absence of DM).

From Figs (16) we note, as for SFDM, the presence of larger cusped tori located far from the central spinning attractor, distinguishing the DM deformed geometry from the Kerr case.

IV Discussion and Final Remarks

In DM models considered here there are NS solutions, solutions with one horizon and two horizons, according to the DM parameters. There are also BH spacetime solutions with horizons at a>M𝑎𝑀a>Mitalic_a > italic_M or NSs for a<M𝑎𝑀a<Mitalic_a < italic_M. The geodesic structure regulating the accretion physics and the tori location around the central spinning attractor can be shifted considerably outwardly with the respect to the Kerr geometry. DM effects mimic Kerr attactors with altered spin to mass ratio a/M𝑎𝑀a/Mitalic_a / italic_M. For example, in all the models presented, the DM affects the disk inner edge which is a tracer of the a/M𝑎𝑀a/Mitalic_a / italic_M in the Kerr geometry. The presence of an excretion cusp, double cusps, or double tori, which are also typical of Kerr NSs solutions, could indicate the presence of DM. Consequently DM affects BH horizon physics, considering DM (models) as NSs mimickers, or vice versa DM (models) as BHs mimickers for super-spinars (cosmological) solutions. DM could also affect the jet emission. The orbital range locating the proto-jets cusps can be also very small, as discussed in Sec. (III.2.1) for SFDM and in Sec. (III.2.2) for CDM. The open cusped solutions (constraining also the jet emission) are very different from their counterparts in the Kerr spacetime in absence of DM. In general DM manifests also with the existence of extremely large cusped tori orbiting very far from the central singularity. From Figs (12) we see the large dimensions of the cusped tori orbiting SFDM spinning BHs. The equilibrium of these tori may be hugely affected by the their self-gravity.

In all DM models considered here, however, DM affects the geometric and causality properties, while there is no coupling with ordinary matter (nor an hypothetical accretion disk consisting of dark matter in orbit), considering gravity modified by the effects of the dark matter–see also [78]. We addressed three models, drawing qualitative and comparative considerations, ruling out some solutions and tracing some common patterns. We have taken as a selection criterion in the space of the metric parameters the observation that there are two expected regimes, where there is a fully modified geometry, qualitatively divergent with respect to the general relativistic onset, as a strongly different horizons structures compared to the reference Kerr solution, and the second scenarios consisting in an appreciable quantitative deformation of the orbiting structures, but not a qualitatively significant change of the background geometry defined by the spinning BH. The current methods of measuring and identifying BHs are based also on the physics of accretion, being related to the accretion disk inner edge, which we prove to be distorted by the DM treated, in the metric models considered here, as a background deformation.

Spherically symmetric black hole solutions in PFDM have been considered to be adapted to the observed asymptotically flat rotation velocity in spiral galaxies and a possible interaction between the DM halo and central BH has been differently theorized. However it has been supposed that SMBHs could enhance the DM density significantly222222Producing a so-called spike” phenomenon [79].. The results of our analysis prove accretion to be a good indicator of the divergences induced by the DM presence and the study of the accretion disks in DM models to represent a valid DM models discriminant. The tori dimensions provide an indication of the possible effects of DM for the energetics associated with the physics of accretion around BHs. In the P-D models for example, the thickness of the accretion throat (opening of the cusp for tori with specific fluid angular momentum in L1 with K]K×,1[K\in]K_{\times},1[italic_K ∈ ] italic_K start_POSTSUBSCRIPT × end_POSTSUBSCRIPT , 1 [) determines (in the assumptions of vanishing pressure at the inner edge), many characteristics of tori energetics such as mass accretion rates and cusp luminosity, the rate of the thermal-energy carried at the cusp, the mass flow rate through the cusp (i.e., mass loss accretion rate), the fraction of energy produced inside the flow and not radiated through the surface but swallowed by central BH, mass-flux, the enthalpy-flux (related to the temperature parameter), depending also on the EoS, as the polytropic index and constant232323Configurations considered here have been often adopted as the initial conditions in the set up for simulations of the general relativistic magnetohydrodynamic (GRMHD) accretion structures [51, 52, 53, 54]. The geometrically thick axial symmetric hydrodynamical models are widely adopted in many contexts showing a remarkably good fitting with the more complex dynamical models as discussed for example in [55]. In the current analysis of dynamical systems of both general relativistic hydrodynamic (GRHD) and GRMHD set-up, these tori are commonly adopted as initial configurations for the numerical analysis–[53, 54, 82] constituting also a comparative model in many numerical analysis of complex situations sharing the same symmetries. Indeed the general relativistic thick tori morphological features, related to the equilibrium (quiescent) and accretion phases as the cusp emergence, are predominantly determined by the centrifugal and gravitational components of the force balance in the disks rather then the dissipative ones.. It has been shown in [73, 80] that the maximum of flow thickness and the maximum amount of matter swallowed by the central BH is determined by the attractor spin–mass ratio only, being defined by the location of the marginally stable circular orbit, and therefore DM influences on the tori dimension and the marginally stable orbits in these models can be searched in a variation of the central BH energetics242424 A further relevant issue in the analysis of the DM effects on BHs accretion, is if the features shown as track of the DM presence may be used to distinguish the DM models. An answer to this question is immediate from the comparison of the horizons structures for the PFDM model in Figs (2), SFDM model in Figs (9), CDM model in Figs (13), and of the geodesic structure (constraining tori morphology and formation) for the PFDM model in Figs (3), SFDM model in Figs (10), CDM model in Figs (14). Thus, it is immediate to determine that in general the main differences are between PFDM model on the one side and CDM and SFDM models on the other, also for small values of DM parameters. In these DM models we have pointed out ”DM-induced” NSs (slower spinning attractors without BH horizons) and in all cases also ”DM-induced” BHs, (over spinning solutions with one or two horizons). Here, we take into account DM models differentiation by means of the tori characteristics determining possiply the DM presence, as the inter-disk cusps, double accretion tori, presence of extremely large and far tori, limited proto-jets ranges, DM differentiation according to fluid rotation orientation, Lense–Thirring effects in DM presence. While the SFDM and CDM models show qualitatively similar characteristics, not distinguishing substantially DM effects for fluid rotation orientation and slowly spinning from faster spinning attractors, the situation for PFDM is clearly different. PFDM model shows remarkable differences also for a small variation of the DM k𝑘kitalic_k parameter, distinguishing DM effects on co-rotating and counter-rotating fluids and between slowing spinning attractors and faster spinning attractors. The CDM and SFDM models show, for the considered parameters ranges (a[0,M]𝑎0𝑀a\in[0,M]italic_a ∈ [ 0 , italic_M ]) a qualitatively similar geodesics structure compared with the BH case in absence of DM. There are very large tori located far from the attractor, and proto-jets cusps constrained in a narrow orbital range around the attractor, constraining proto-jets emission and the formation of tori with large angular momentum in magnitude. We have also noticed indications of a possible alteration of the Lense–Thirring effects on the disks and flows with respect to the Kerr case without DM. In the considered parameters ranges, major differences of the PFDM models with respect to the case in absence of DM appear in the formation of the inter-cusps, double configurations and possible excretion tori. It must be stressed however that, while we have drawn here a DM models comparative analysis, an in-depth exploration of more extensive DM parameters ranges in all models, would further narrow the DM parameters with the DM effects on the BHs accretion. [64, 65, 80].

It should be emphasized then that since DM-BHs can exhibit features associated with Kerr NSs, it has implications on cosmic censorship, in the fact that observing a compact object with such tracers (excretion cusp, double cusps, double tori) would not require the breaking of cosmic censorship (viewing a Kerr NSs), but instead could mean one is observing a BH surrounded by DM. Finally in this work we developed a comparative analysis of accretion disks in different dark matter models while we reserve the in-depth explorations of different DM parametric values for future analysis.

References

  • [1] P. E. Mancera Piña, F. Fraternali, T. Oosterloo, et al., MNRAS, 512, 3230 (2022).
  • [2] P. E. Mancera Pina et al., ApJL, 883, L33 (2019).
  • [3] E. Belbruno, J. Green, MNRAS, 510, 4, 5154-5163 (2022).
  • [4] X. Hou, Z. Xu and J. Wang, JCAP, 12, 04 (2018).
  • [5] F. Rahaman, K.K. Nandi, A. Bhadra, et al, Physics Letters B, 694, 10–15 (2010).
  • [6] Z. Xu, X. Hou and J. Wang, Class. Quantum Grav. 35, 115003 (2018).
  • [7] A. Das et al. Class. Quantum Grav., 38, 065015 (2021).
  • [8] Z. Xu, X. Hou, X. Gonga, J. Wanga, JCAP,09, 038 (2018).
  • [9] X. Hou, Z. Xu, M. Zhou, and J. Wang, JCAP,07, 015 (2018).
  • [10] S. Shaymatov, D. Malafarina, B. Ahmedov, Physics of the Dark Universe, 34, 100891 (2021).
  • [11] Z. Stuchlik, BAIC 34, 129, (1983).
  • [12] Z. Stuchlík & S. Hledík, Physical Review D, 60, 044006 (1999).
  • [13] Z. Stuchlík, M. Kološ, J. Kovář, Slaný et al., Univ, 6, 26 (2020).
  • [14] Z. Stuchlík, S. Hledík, CQGra, 17, 4541 (2000).
  • [15] Z. Stuchlík, MPLA, 20, 561 (2005).
  • [16] M. Khodadi, R. Pourkhodabakhshi, Physics Letters B, 823, 136775 (2021).
  • [17] B. Toshmatov, Z. Stuchlik, B. Ahmedov, Eur. Phys. J. Plus, 132, 98 (2017).
  • [18] H. Kim, B. Lee, W. Lee, and Y. Lee, Physcal Review D, 101, 064067, (2020).
  • [19] B. Cuadros-Melgar, R. D. B. Fontana, J. de Oliveira, Physical Review D, 104, 104039 (2021).
  • [20] M. Li and K. Yang, Physical Review D 86, 123015, (2012).
  • [21] Z. Stuchlík, J. Schee, CQGra, 30, 075012 (2013).
  • [22] Z. Stuchlik, BAICz, 31, 129 (1980).
  • [23] M. Blaschke, Z. Stuchlík, PhRvD, 94, 086006 (2016).
  • [24] Z. Stuchlík, S. Hledík, K. Truparová, CQGra, 28, 155017 (2011).
  • [25] R. Abbott et al, arXiv:2111.15507v1 [astro-ph.HE].
  • [26] D. Traykova et al. Phys. Rev. D 104, 103014 (2021)
  • [27] K. Clough, P. G. Ferreira, and M. Lagos, Phys. Rev. D 100, 063014, (2019).
  • [28] J. Bamber, K. Clough, P. G. Ferreira, L. Hui, and M. Lagos, Phys. Rev. D 103, 044059 (2021).
  • [29] H. Davoudias, P. B. Denton and J. Gehrlein Physical Review Letters 128, 081101 (2022).
  • [30] L. E. Padilla, T. Rindler-Daller, P. R. Shapiro, et al.Phys. Rev. D 103, 063012 (2021).
  • [31] J. Bamber, K. Clough, P. G. Ferreira, L. Hui, and M. Lagos, Phys. Rev. D 103, 044059, (2021).
  • [32] K. Jusufi Phys. Rev. D ,101, 084055 (2020).
  • [33] T. Lacroix, M. Karami, A. E. Broderick, J. Silk, C. Boehm Physical Review D 96, 063008 (2017).
  • [34] K. Jusufi, M. Jamil, T. Zhu, Eur. Phys. J. C, 80, 354 (2020).
  • [35] K. Jusufi, M. Jamil, P. Salucci, T. Zhu, and S. Haroon, Physical Review D 100, 044012 (2019).
  • [36] F. Atamurotov et al. Class. Quantum Grav. 39, 025014 (2022).
  • [37] A. Das et al., Phys. Rev. Lett. 128, 021101 (2022).
  • [38] T. S. Li, A. P. Ji, A. B. Pace, et al., Astrophys. J. , 928, 30 (2022).
  • [39] R. P. Naidu, C. Conroy, A. Bonaca, et al., Astrophys. J. , 923, 92 (2021).
  • [40] N. Cappelluti, G. Hasinger, P. Natarajan, ApJ, 926, 205 (2022).
  • [41] S. Basak et al,ApJ Lett., 926, L28 (2022).
  • [42] C. Smorra, Y. V. Stadnik, P. E. Blessing,et al. Nature, 575, 7782, 310-314 (2019).
  • [43] S. Afach, B.C. Buchler, D. Budker, et al. Nat. Phys. 17, 1396–1401 (2021).
  • [44] A. Gomez-Valent, Z. Zheng, L. Amendola, et al. Phys. Rev. D, 104, 083536 (2021).
  • [45] Z. Stuchlí, S. Hledík, J. Novotný, PhRvD, 94, 103513 (2016).
  • [46] M. Kozłowski, M. Jaroszyński, M. A. Abramowicz Astron. Astrophys., 63, 209 (1998).
  • [47] M. A. Abramowicz,M. Jaroszyński, M. Sikora Astron. Astrophys, 63, 221 (1978).
  • [48] M. Jaroszynski, M. A.Abramowicz, B. Paczynski, Acta Astron., 30, 1 (1980).
  • [49] D. Pugliese, G. Montani & M. G. Bernardini, Mon. Not. R. Astron. Soc., 428 (2), 952 (2013).
  • [50] D. Pugliese& G. Montani, Europhys. Lett., 101, 19001 (2013).
  • [51] I. V. Igumenshchev, M. A. Abramowicz, Astrophys. J. Suppl.,130, 463 (2000).
  • [52] R. Shafee, J. C McKinney, R. Narayan, et al. Astrophys. J. , 687, L25 (2008).
  • [53] P. C. Fragile, O. M. Blaes, P. Anninois, J. D. Salmonson, Astrophys. J. , 668, 417-429 (2007).
  • [54] J-P. De Villiers,& J. F. Hawley, Astrophys. J. , 577, 866 (2002).
  • [55] Q. Lei, M. A. Abramowicz, P. C. Fragile, et al., A&A., 498, 471 (2008).
  • [56] M. A. Abramowicz, arXiv:astro-ph/0812.3924 (2008).
  • [57] J. A. Font& F. Daigne, Astrophys. J. , 581, L23–L26 (2002).
  • [58] M. A. Abramowicz& P.C. Fragile, Living Rev. Relativity, 16, 1 (2013).
  • [59] D. Pugliese and G. Montani, Phys. Rev. D 91, 8, 083011 (2015).
  • [60] B. Paczyński, Acta Astron., 30, 4 (1980).
  • [61] J. A. Font, Living Rev. Relat., 6, 4 (2003).
  • [62] M. Kozłowski, M. Jaroszyński, M. A. Abramowicz, Astron. Astrophys., 63, 209 (1998).
  • [63] M. A. Abramowicz, M. Calvani, L. Nobili, Astrophys. J., 242, 772 (1980).
  • [64] M. A. Abramowicz, Astronomical Society of Japan, 37, 4, 727-734 (1985).
  • [65] D. Pugliese&Z. Stuchlik , Eur. Phys. J. C 79 4, 288, (2019).
  • [66] D. Pugliese & Z. Stuchlík,Class. Quant. Grav.  35, 18, 185008 (2018).
  • [67] Z. Stuchlik, BAICz, 32, 68 (1981).
  • [68] D. Pugliese, H. Quevedo and R. Ruffini, Phys. Rev. D, 84, 044030 (2011).
  • [69] K. Adamek, Z. Stuchlik, CQGra, 30, 205007, (2013).
  • [70] P. Slaný, Z. Stuchlík, CQGra, 22, 3623, (2005).
  • [71] D. Pugliese&Z. Stuchlík, Astrophys. J.s, 221, 2, 25 (2015).
  • [72] D. Pugliese, G. Montani, Gen. Rel. Grav, 53,5,51, (2021).
  • [73] D. Pugliese, Z. Stuchlik, submitted 2022.
  • [74] D. Pugliese&Z.Stuchlík, Astrophys. J.s, 223, 2, 27 (2016).
  • [75] D. Pugliese & Z. Stuchlík, Class. Quant. Grav.  35,10, 105005 (2018).
  • [76] D. Pugliese, Z. Stuchlik, PASJ,73, 5, 1333-1366, (2021).
  • [77] D. Pugliese, Z. Stuchlik, PASJ- 73, 6, 1497-1539 (2021)
  • [78] E. Kurmanov, K. Boshkayev, R. Giambo, et al.Astrophys. J. , 925, 210, (2022)
  • [79] L. Sadeghian, F. Ferrer and C.M. Will, Phys. Rev. D 88, 063522 (2013).
  • [80] D. Pugliese, Z. Stuchlik, MNRAS, 512, 4, 5895–5926, (2022).
  • [81] R. H. Boyer, Proc. R. Soc. London A, 311, 245 (1969).
  • [82] O. Porth, H. Olivares, Y. Mizuno, Z. Younsi, et al.arXiv:1611.09720 [gr-qc] (2017).