On dark matter effect on BH accretion disks
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.
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 ) and the DM deformation on the Kerr extreme BH spacetime (with spin value ), 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 -CDM, which includes a cosmological constant () (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 -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 and of constant specific angular momentum coincide. This implies that the rotation law 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 ), the Kerr metric describes naked singularities (NSs) for and black holes (BHs) for . The Kerr BH geometry has the limiting static solution of Schwarzschild for and the extreme Kerr BH spacetime for .
In the Boyer-Lindquist (BL) coordinates , the metric tensor reads777We adopt the geometrical units and the signature, Latin indices run in . The radius has unit of mass , and the angular momentum units of , the velocities and with and . For the seek of convenience, we always consider the dimensionless energy and effective potential and an angular momentum per unit of mass .:
(1) |
where
(2) |
with . The horizons are respectively given by
(3) |
the horizons can be found888Quantities and 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 for . The outer and inner stationary limits (ergosurfaces) (solutions of ) are respectively
(4) |
the ergosurfaces can be found by solving the equation (for ), where on and in the equatorial plane (). Static observers, with four-velocity (where indicates the derivative of any quantity 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 (where are functions of the plane ). Here we often intend the outer ergoregion (or simply ergoregion) in the BH spacetimes as the region . Then on the equatorial planes, in the Kerr spacetime there is and the outer ergosurface is . , but trajectories , including particles crossing the stationary limit and escaping outside in the region are possible.
The constants of the geodesic motions are
(5) |
with101010The other constant of geodesic motion of the Kerr metric is the Carter constant . In this work, where tori share symmetry plane with the equatorial plane of the central BH, this constant is irrelevant. . In Eqs (5) quantities and 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
(6) |
respectively. The sign of 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 and of constant specific angular momentum coincide and the toroids rotation law 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 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
(7) |
is the fluid energy momentum tensor, and are the total energy density and pressure, respectively, as measured by an observer moving with the fluid. The timelike flow vector field denotes the fluid four-velocity. The fluid dynamics is described by the continuity equation and the Euler equation respectively:
(8) |
where the projection tensor and .
We assume a barotropic equation of state (EoS) 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 , as for the circular test particle motion no motion it is assumed in the angular direction, which means . Because of these symmetries, the continuity equation is identically satisfied and the orbiting configurations are regulated by the Euler equation for the pressure only, which can be written as
(9) |
where 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 .
Assuming the fluid is characterized by the specific angular momentum constant (see also discussion [55]), we consider the equation for or constant. By setting 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]
(10) |
The extremes of the pressure are regulated by the angular momentum distributions , on the equatorial plane for co-rotating and counter-rotating fluids respectively131313Note, in the test particles analysis and accretion tori models, for slowly spinning NSs (), there are circular geodesic orbits with and 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 ) ; 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 and or counter-rotating according to but co-rotating according to , there can also be orbits with –see, for example, [22, 67, 68, 69, 70, 80, 73]. This possibility has not been discussed in the analysis of DM models. .
Torus cusp is the minimum point of pressure and density in the torus corresponding to the maximum point of the fluid effective potential. The torus center is the maximum point of pressure and density in the torus, corresponding to the minimum point of the fluid effective potential. At the cusp () the fluid may be considered pressure-free. Fluid effective potential defines the function . Cusped tori have parameter in the open ranges141414The notation is as follows: the (closed) interval between quantities and , including and , is denoted as, the notation denotes the open interval. Similarly there is and ., where . (We adopt the notation for any quantity evaluated on a radius .)
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 , which is also a photon circular orbit, , the marginally bounded orbit, , and the marginally stable circular orbit, (see Figs (1))151515In the Kerr spacetime is the marginal circular orbit (a photon circular orbit) where timelike circular orbits can fill the spacetime region . Stable circular orbits are in for counter-rotating and co-rotating test particles respectively. The marginal bounded circular orbit is , where . More details and the exact forms of these radii can be found for example in [68].. Radii constrain the location of the tori cusps (inner edges) with fluid specific angular momentum respectively
(11) |
where we introduced also the radii defined as
significant as governing the location of the tori centers, More precisely ranges of fluids specific angular momentum govern the tori topology, according to the geodesic structure of Eqs (11), as follows:
- :
-
for there are quiescent (i.e. not cusped) and cusped tori–where there is . The cusp is (with )) and the center with maximum pressure in .
- :
-
for there are quiescent tori and proto-jets (open-configurations) –where there is . The cusp is associated to the proto-jets, with , and the center with maximum pressure is in . Proto-jets are associated to (not-collimated) open structures, with matter funnels along the BH rotational axis–see [59, 74, 75, 76];
- :
-
for there are only quiescent tori where there is and the torus center is at .
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, , as the radius , the marginally bounded orbit defined by (asymptotically flat spacetimes) and the marginally stable orbits .
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
(12) | |||
for , where is for the derivative with respect to . 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:
(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 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
(14) | |||
(15) | |||
see also [5, 6, 7], where is the parameter describing the intensity of the PFDM, set in the ranges . For the line element reduces to the Kerr metric161616Constrains of positive were obtained by fitting the rotation curves in spiral galaxies, with values [4].. The metric singularities, , defining the DM deformations on the Kerr horizons, can be found by solving the equation or from the equation , where is the Lambert function, such that gives the principal solution for in .
The deformed ergosurfaces can be found, for , as solution of the equation while, on the equatorial plane there is
(16) |
(see red curve in Figs (2)–bottom left panel) where the Lambert function gives the solution for in , and gives the sign of , therefore it is for negative, zero, or positive.
The PFDM horizons are independent of (like in the Kerr case). The PFDM ergosurfaces are independent of spin in the equatorial plane (like in the Kerr cases). The horizons (and the ergosurfaces on the equatorial plane of Eq. (16)) are shown in Figs (2).
The red curve in Figs (2), which represents the BH horizon for (and the ergosurfaces Eq. (16) on the equatorial plane for ), bounds the collections of horizons at different constant in the plane . The PFDM metric describes solutions with ,, and horizons. For BH solutions, horizons can be shifted outwardly or inwardly with respect to the Kerr BH spacetime depending on the value of . There are also spacetime solutions with horizons for .
Accordingly we select the PFDM metric parameter in the following six cases
(17) |
The cases have as horizon for (the static case) therefore, in this sense, these solutions can be compared to the Schwarzschild case. Similarly, for , the geometry with has one horizon at . There is one horizon when and 171717From Figs (3) we note that there is one horizon at for the special parameters and at for ..
In Figs (4) we show also the fluid specific angular momentum distribution for the parameters, for the cases and , compared to the distribution on the geometry in absence of DM (see also Figs (1)), the associated parameter and the Keplerian (test particle) angular momentum respectively, which we have defined as related to the thick tori counterparts from definition of , 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 and , 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 of the orbiting matter. According to the DM parameter and spin , relation (11) i.e. for 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 , reflects in the relative location of the maximum-minimum points of pressure of the orbiting disks. For values of where this relation is not verified, for example for 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 () The geodesic structure is represented in Figs (3)-left panel, compared with the Schwarzschild case. We note that for the situation changes qualitatively with respect to the case in absence of DM. In Figs (4) is the fluid specific angular momentum , the tori energy parameter and (test particles) Keplerian angular momentum as function of , for different PFDM parameters of Eq. (14), compared with the case describing the Schwarzschild geometry. Notably, for , there is . From the analysis it is noted how, for some values of , curves are lower with respect to the corresponding curves in absence of DM–Figs (3). From the analysis of the curves , 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 () The geodesic structure for the spinning attractor geometry in PFDM is in Figs (3)–right panel. The fluid specific angular momentum , energy parameter and (test particles) Keplerian angular momentum are shown in Figs (4) as functions of for co-rotating and counter-rotating fluids and PFDM parameters of Eq. (14), compared to the case 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 is similar to the counter-rotating geodesic structure in the axially symmetric spacetime . The co-rotating case (for ), , especially for , is remarkably different from the geodesic structure in the geometries with 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 . For larger there are also NS solutions (also for ) and, for even large values of , there are BH solutions (also for ).
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 (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 for and for .
The inter-disk cusp181818Similarly to the outer cusps, the inter–disks cusp is a tori cusp located between two configurations having same 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 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 , in the case , 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
(18) |
where
(19) |
The metric components satisfy the asymptotically flatness condition, and the fluid potential is well defined at infinity (), where . We introduce the quantity . Here is the central density and 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 the metric reduces to the Kerr solution). The Kerr limit occurs for or . The zeros of distinguish the metric singularities (SFDM deformed horizons), while the zeros of define a deformation of the Kerr ergosurfaces , which can be found by solving the equation and 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. ). The ergosurface location (on planes ) is an important factor regulating the Lense–Thirring effect on the disks and on the jets flows coming from the disks[77]., where
(20) |
In accordance with the study of the BHs horizons in Figs (9) we explore the following two cases:
(21) |
(see Figs (9)).
However, as clear from Figs (9), BH horizons exist also for .
The geometry circular geodesic structure is shown in Figs (10) for the parameter and , emphasizing the differences for the cases and and the Kerr geometry circular structure of Figs (1). The orbital range locating the proto-jets cusps, bounded by the radii and , for fluids specific angular momentum 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 , and for differently from the Kerr case in absence of DM, radii are located out of the ergoregion (and partially for ), this could imply a significant difference in the Lense-Thirring effects in presence of DM. (There can be however over-spinning BHs (with ) 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 () and spinning attractors () with SFDM.
We analyse below the case of the static attractor i.e. , 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 , 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 . More specifically:
- –The static attractor ()
-
In Figs (11) the fluid specific angular momentum , tori energy parameter and (test particles) Keplerian angular momentum are shown as functions of for co-rotating and counter-rotating fluids, at different spins and SFDM parameters, compared to the case corresponding to the Schwarzschild geometry. Tori and effective potentials are in Figs (12).
- –The spinning attractors ( and )
-
We restrict our analysis to and , studying the cusped tori limiting the closed configurations, regulated by the effective potential function. In Figs (11), there are the fluid specific angular momentum , energy parameter and (test particles) Keplerian angular momentum as function of 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 for large values of . In Figs (11) the fluid specific angular momentum distribution, compared to the distribution on the geometry in absence of DM, the associated energy parameter and the Keplerian (test particle) angular momentum are shown. From Figs (9) it is clear how the horizon curves in the plane 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 (from the tori effective potentials we can also note how the tori parameter for tori orbiting in SFDM are generally considerably lower than the 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
(22) |
where
(23) |
We adopt the parametrization , where is the density of the universe at the moment when the DM halo collapsed, is a characteristic radius. The metric is asymptotically flat and we find the Kerr limit in or . We first consider the metric singularities, identifying the space of the parameters used in the tori analysis. The horizons can be written as solutions of the equation , or or where
(24) |
The metric horizons are defined for
(assuming ).
The ergosurfaces can be found for as solutions of the equation and, on the equatorial plane as solutions of or where212121As clear from Figs (13), these relations also represent some portions of the ergosurfaces according to conditions on the DM parameters.
(25) |
see Figs (13). On the equatorial plane, the outer ergosurface, independent on the spin , corresponds to the metric singularity in the static () case–see Eqs (24).
The geodesic structure for this geometry is shown in Figs (14).
Similarly to the SFDM model, the equatorial geodesic structure shows that for the case, the range of the proto-jets cusp location is remarkably narrow. At , for co-rotating fluids, the radii for () do not enter (are partially contained in) the outer ergoregion.
We have analyzed the static and spinning attractors as follows:
- –The static attractor ()
- –The spinning attractor ()
-
Fluid specific angular momentum , energy parameter and (test particles) Keplerian angular momentum as function of are represented in Figs (15) for co-rotating and counter-rotating fluids, different spins and CDM parameters in comparison with the extreme Kerr BH case in absence of DM. Effective potential and tori for CDM parameters and are in Figs (16), for and in Figs (16), for fluid specific angular momentum () 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 for large holds, where . Large tori orbiting CDM spinning BHs are shown in Figs (16) as, for example, the blue surface for the case , dashed -blue curve for and (similarly to the SFDM case, the tori parameter for tori orbiting in CDM are generally considerably lower than the 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 or NSs for . 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 . For example, in all the models presented, the DM affects the disk inner edge which is a tracer of the 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 ) 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 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 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).