Bondi accretion in the spherically symmetric Johannsen–Psaltis spacetime

The Johannsen–Psaltis spacetime explicitly violates the no-hair theorem. It describes rotating black holes with scalar hair in the form of parametric deviations from the Kerr metric. In principle, black hole solutions in any modified theory of gravity could be written in terms of the Johannsen–Psaltis metric. We study the accretion of gas onto a static limit of this spacetime. We utilise a recently proposed pseudo–Newtonian formulation of the dynamics around arbitrary static, spherically symmetric spacetimes. We obtain a potential that generalises the Paczyński–Wiita potential to the static Johannsen–Psaltis metric. We also perform a fully relativistic analysis of the geodesic equations in the static Johannsen–Psaltis spacetime. We find that positive (negative) values of the scalar hair parameter, ϵ3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\epsilon _{3}$$\end{document}, lower (raise) the accretion rate. Similarly, positive (negative) values of ϵ3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\epsilon _{3}$$\end{document} reduce (increase) the gravitational acceleration of radially infalling massive particles.


Introduction
The no-hair theorem [1][2][3] states that black holes in general relativity are uniquely and completely characterised by three parameters viz. their mass, electric charge and angular momentum. An electrically charged black hole will accrete charges of the opposite sign and rapidly neutralise [4]. Consequently the only physically significant parameters of astrophysical black holes are their mass and angular momentum. The Kerr metric is thus the only physical stationary, axisymmetric, asymptotically flat vacuum spacetime in general relativity that possesses an event horizon.
Whilst general relativity has been extremely well constrained on solar system scales [5] the theory has only recently been subjected to tests in strong gravitational fields. The discovery of gravitational waves [6] from the merger of a e-mail: a.john@ru.ac.za b e-mail: c.stevens@ru.ac.za binary black holes provided the first strong field tests of general relativity. Precise determination of the metric of astrophysical black holes allows one to test general relativity in the strong field regime. Alternative theories of gravity can admit black hole spacetimes that do not comply with the no-hair theorem. Observed violations of the no-hair theorem necessarily imply the breakdown of general relativity.
Given the advent of gravitational wave astronomy as well as forthcoming complementary observations of black holes and neutron stars in the electromagnetic spectrum, strong field tests of gravitational theories are becoming feasible [7,8]. A large number of modified theories of gravity have been proposed. Examples of such theories include Modified Gravity (MOG) [9], Tensor-Vector-Scalar gravity (TeVeS) [10] and f (R) gravity [11]. These are motivated by the need to explain the dark matter, dark energy and singularity problems in cosmology as well as the search for a quantum theory of gravity.
Modified theories of gravity can be tested by solving their associated field equations to obtain black hole solutions and comparing their observable signatures with those predicted by general relativity. Sufficiently precise astronomical data can, in principle, discriminate between competing gravitational theories. Alternatively one could devise a model-independent test of general relativity by constructing black hole spacetimes with parametric deviations from the Kerr solution. These black hole spacetimes are generic in the sense that they are devised without appealing to a specific theory of gravity. Null measurements of these extra parameters would mark a successful test of the no-hair theorem. The Event Horizon Telescope recently obtained the first image of a black hole photon sphere or shadow [12]. Black holes in modified theories of gravity are predicted to have shadows with different shapes, sizes and multiple images compared to those occurring in general relativity [13]. The Event Horizon Telescope observation places an upper limit on the fractional quadrupole moment deviation of ΔQ/Q ≤ 4 which is con-sistent with general relativity [12]. This constraint should improve with greater resolution in very long baseline interferometry and provide stringent tests on violations of the no-hair theorem.
The Johannsen-Psaltis (JP) metric [14] facilitates modelindependent tests of gravitational theories. This metric describes the spacetime of a stationary, axisymmetric black hole that violates the no-hair theorem by construction. Significantly the JP metric does not arise as a solution to field equations of any particular theory of gravity. Given any gravitational theory that violates the no-hair theorem it should be possible to parametrise the deviations from the Kerr metric of its black hole solution(s). In principle these deviations can be related to the additional parameters of the JP metric. The JP metric was not the first proposed generic black hole spacetime. Earlier attempts e.g. [15] were plagued by problems like the presence of singularities and closed timelike curves outside the event horizon. The JP metric is a promising alternative for parametrising modified gravity corrections to the Kerr solution as it is free of many of these pathologies.
The JP metric has been applied to a number of problems e.g. strong lensing [16] and the formation of black hole shadows [17]. The measurement of black hole spin via continuum fitting and K α iron line methods has also been analysed in the JP spacetime [18].
The discovery of active galactic nuclei and quasars as well as compact objects like neutron stars and pulsars prompted the realisation that their high energy emission was due to the liberation of the gravitational potential energy of infalling material [19]. Bondi's seminal paper [20] determined the mass accretion rate of a polytrope fluid accelerating towards a star with a Newtonian gravitational potential. This result was later extended to full general relativity by Michel [21] who examined accretion onto a Schwarzschild black hole. Accretion disks form when matter falls onto a rotating body. The structure of accretion disks was modelled both in Newtonian gravity [22] and general relativity [23]. Black hole accretion in various modified theories of gravity has been investigated by a number of authors [24][25][26].
Modelling accreting systems in full general relativity becomes significantly more difficult when one includes phenomena like viscosity, magnetic fields, turbulence and radiative processes. Paczyński and Wiita [27] introduced a gravitational potential that mimics many features of the Schwarzschild black hole in general relativity. The Paczyński-Wiita potential is commonly used to study black holes in a Newtonian manner [28]. The advantage of using this Pseudo-Newtonian potential is that the fluid equations are substantially simpler to analyse than those arising from a full relativistic treatment.
The success of the Paczyński-Wiita potential can be accounted for by its prediction of the same values for the innermost stable circular orbit (ISCO), the marginally bound orbit (r mb ), and the Keplerian angular momentum, L(r ), as those obtained from a fully general relativistic analysis of the Schwarzschild geodesic equations [28]. The associated dynamical systems have the same fixed points [28].
The Paczyński-Wiita potential, Φ PW , is not the only pseudo-Newtonian potential devised, but it is by far the most popular. Other proposals tend to be obtained by fitting formulae. The Φ PW potential has some limitations. It cannot properly account for the dynamics around rotating black holes as it fails to predict the Lense-Thirring effect, nor can it accurately describe self-gravitating systems [28].
Given the successful approximation of static black hole physics in general relativity by the Paczyński-Wiita potential, the question arises as to whether pseudo-Newtonian potentials can be found for modified theories of gravity. A number of authors [29,30] have recently devised methods to generate pseudo-Newtonian potentials that could mimic the behaviour of black holes in any static, spherically symmetric spacetime.
In this paper we determine a pseudo-Newtonian potential for the static limit of the Johannsen-Psaltis metric and solve the accretion problem. Our motivation in doing so is to obtain a model-independent formulation of spherically symmetric accretion in modified theories of gravity. Bondi accretion (where the fluid velocity is purely radial) is an important first approximation to more complex accretion problems. Many of the features of this idealized scenario (e.g. the existence of a transonic solution) survive in the more complicated accretion disk problem. Azimuthal motion of accreting gas onto a Newtonian potential was first considered by treating the angular momentum as a small perturbation about a spherically symmetric Bondi flow [31]. The authors concluded that the accretion rate of a rotating gas is slightly lower than in the original problem, and that this change can be expressed in terms of the boundary conditions at infinity. Generalising the pseudo-Newtonian approach to stationary spacetimes has been explored [32]. The authors determined that within their formalism it was difficult to obtain closed form Lagrangians for spacetimes exhibiting frame-dragging or gravitomagnetic effects. Obtaining pseudo-Newtonian potentials that can successfully describe the physics around rotating black holes is a highly non-trivial problem, which we do not attempt to solve here.
In Sect. 2 we introduce the Johannsen-Psaltis metric and its static limit. In Sect. 3 we obtain a pseudo-Newtonian potential for the JP metric. In Sect. 4 we solve the accretion problem for JP black holes approximated by our potential. In Sect. 5 we analyse the motion of test particles in the JP spacetime. Our results are summarised in Sect. 6. Throughout this paper, unless explicitly stated, we utilise geometric units where c = 1 = G.

The Johannsen-Psaltis spacetime
The Johannsen-Psaltis metric takes the form where and m and a are the black hole mass and spin parameters respectively. Here h(r, θ) is the parametric deviation from the Kerr spacetime and the terms are dimensionless coefficients. It can easily be seen that when h(r, θ) = 0 this metric reduces to the standard Kerr spacetime. We are interested only in the spherically symmetric case, which reduces the form of the metric to where the parametric deviation now takes the form In the absence of any scalar hair i.e. k = 0 we recover the line element for the Schwarzschild black hole.
There are a number of observational constraints on the magnitude of the k parameters [14]. In general relativity stationary and asymptotically flat spacetimes must fall off as 1/r or faster. Spacetimes with slower fall off rates will include gravitational radiation and thus cannot be stationary. Similar arguments hold for spacetimes that do not arise as solutions of the Einstein field equations. Thus the function h(r ) must be of order O(r −n ) where n ≥ 2 and this forces 0 = 1 = 0. Weak field deviations from general relativity can be determined using the parameterized post-Newtonian (PPN) framework. The Lunar Laser Ranging experiment [5] sets an upper bound of | 2 | ≤ 4.6 × 10 −4 . The first unconstrained parameter in h(r ) is thus 3 . Analyses involving the JP spacetime often only consider the leading non-vanishing contribution to h(r ) i.e. the scalar hair is approximated via Throughout this paper we will either use the full expression for h(r ) viz. (6) or the truncated version (7).

The Paczyński-Wiita-like potential
In this section we derive the Paczyński-Wiita like potential Φ J P for the Johannsen-Psaltis metric.
The original Paczyński-Wiita potential is pseudo-Newtonian in the sense that it is a potential representing the relativistic Schwarzschild spacetime that can be used in a Newtonian framework. The potential is where m is the mass parameter from the Schwarzschild metric.
Although the form of this potential was essentially guessed by Paczyński, it has been subsequently derived from an algorithm applied to the Schwarzschild spacetime [28]. This method has been extrapolated to general static spherically symmetric spacetimes [29] and we summarise the argument below.
Any static spherically symmetric spacetime written in coordinates {t, r, θ, φ} can be written in the form ds 2 = g 00 (r )dt 2 + g 11 (r )dr 2 + r 2 dΩ 2 , where dΩ 2 = dθ 2 + sin 2 θ dφ 2 . The derivation of the pseudo-Newtonian potential relies on examining the orbits of a test particle in the spacetime which, without loss of generality, can be restricted to the θ = π/2 plane. The effective gravitational potential can be identified in the geodesic equations and takes the form where k ∈ R is an arbitrary constant. When applied to the Schwarzschild spacetime, the choice k = 1 yields the original Paczyński-Wiita potential. Turning our attention to the spherically-symmetric Johannsen-Psaltis metric Eq. (5), we find the corresponding potential when again choosing k = 1 is which clearly reduces to Φ PW when h(r ) = 0. Note that this expression uses the full expansion of h(r ) given by (6).

Bondi accretion
We now compute the mass accretion rate of a gas falling into a black hole described by the line element (5). In the spirit of Paczyński and Wiita we model this problem using classical hydrodynamics and the effective gravitational potential, Φ J P . We model the gas as a perfect fluid with a polytrope equation of state, p = Kρ γ , where the adiabatic index, γ , is given by the ratio of specific heats i.e. γ = c p /c v .
The dynamical motion of the fluid is determined by the conservation of mass and momentum viz.
where u is the fluid 4-velocity, p is the isotropic pressure and ρ is the mass density. The effective pseudo-Newtonian potential Φ J P is given by (11). Imposing spherical symmetry, and considering only radially infalling steady-state flows gives us the simpler system u du dr where u(r ) is the radial velocity of the fluid. For convenience we rewrite the equation of state using the polytropic index, p = Kρ 1+ 1 n . Now we can integrate Eqs. (14) and (15) to find the conserved mass flux and specific energy respectively aṡ where we have denoted a = dp dρ as the adiabatic sound speed. Our boundary conditions are determined at a point, r ∞ , which is sufficiently far away from the black hole that the fluid velocity (u(r ∞ ) =: u ∞ ) and the gravitational potential both vanish. The fluid sound speed "at infinity" is denoted by a ∞ . Further, we rewrite (16) to obtain the accretion ratė which is related to the mass flux viaṀ =Ṁγ n K n [33]. We now combine Eqs. (17) and (18) to obtain To ensure that the accretion flow is smooth, the denominator of (19) must vanish at the same point as the numerator. This critical point, r c is determined by the conditions: The first condition tells us that the critical point is also a sonic point. If the fluid accelerates from rest far away from the black hole (u ∞ = 0) it must reach its local sound speed at the critical point i.e. u c = a c . Solving (20) for the critical points is equivalent to finding solutions to the equation This equation is valid for the full form of h(r ) given by (6). This expression is quite complicated in general so we only consider the leading contribution to h(r ) given by (7), which is commonly done in studies of the JP spacetime [14,34]. The specialisation reduces Eq. (21) to a 8th order polynomial equation in α := r c /m. To simplify this further for the allowance of analytic solutions, we linearize the expression in h(r ) (equiv. 3 ) which reduces (21) to a quartic polynomial equation given by where we must compute the solutions to relate r c to a c . It is not surprising that even in the simplified case of h(r ) considered here, the solutions of this quartic equation are lengthy expressions. We observed that for a wide range of values for To obtain a more tractable analytical expression for this root, we continue to linearize about 3 for the rest of this section.
The physically relevant root is given by with X := √ 1 + 16a 2 . This is further simplified by linearising our equations in a 2 c as done in [35]. This gives the simple expression which reduces to the well known result for the Schwarzschild spacetime found in the literature [21,35]. Using Eq. (24) and Eq. (17) and linearising about 3 , a 2 c and a 2 we find the relationship between the critical sound speed and the boundary condition which is the same as in the Schwarzschild case.
Using the above together with Eq. (18) we obtain the Bondi accretion rate, which is kept in this form (i.e. not being completely linear in a 2 ∞ ) for simplicity.

Comparison with the Schwarzschild case
To compare pseudo-Newtonian accretion in the Johannsen-Psaltis and Schwarzschild spacetimes we investigate a range of different system parameters and tabulate how 3 affects the accretion rate. These parameters are m, γ ⇔ n and a ∞ . For the below we fix the boundary condition a ∞ = 3 × 10 −6 which follows from the temperature of ionised interstellar gas [24]. We consider two fiducial values for the black hole mass viz. a solar mass black hole, m = 1M , and a supermassive black hole, m = 10 6 M . We restrict γ to have values 1 < γ < 5/3. We writeṀ = A + B 3 in order to tabulate the results. One sees from the Tables 1 and 2 that 3 affects the accretion rate A predicted by general relativity by adding the number B which, if A has order a, then B has order roughly where e is the order of 3 . Note that we took 3 to be small to arrive at Eq. (4), and thus e should mimic this restriction. When 3 > 0 the accretion rate is lowered, while when 3 < 0 the accretion rate is raised. This behaviour is consistent with a back-of-the-envelope calculation of the gravitational force of our system. For our potential, Φ J P , we have F ∼ −Φ J P ∼ − 1 r 2 + 3 1 r 4 . For positive values of 3 the gravitational force becomes less negative and hence less attractive. Similarly negative values of 3 result in a more negative and thus larger attractive force. Further, when γ increases, so does the accretion rate. This is consistent with other results in the literature [36].

Test particle orbits
To further investigate the accretion properties of the previous section, we perform a fully relativistic calculation of test particle orbits in the spherically-symmetric Johannsen-Psaltis spacetime and examine how they are effected by a non-zero 3 . Note that we do not linearise about 3 in this section.
First, we obtain the constants of motion of a test particle travelling in the θ = π/2 plane with 4-velocity u a = {t (λ), r (λ), 0, φ (λ)} parametrised by λ and satisfying u a u a = −κ 2 : Here E and L are the conserved energy and angular momentum, V (λ) is the potential and κ = 1, κ = 0 correspond to timelike and null orbits respectively.
To get an idea of how timelike test particles behave outside the event horizon under different choices of 3 , we can numerically solve Equation (29) with potential given by Equation The results are plotted in Fig. 1 for different values of the 3 parameter. 2 One can see that 3 > 0 has the effect of slowing down the particles' approach toward the central object. This is in direct agreement with the slower accretion rate calculated in Sect. 4. For the choice 3 < 0 the opposite is true and is again consistent with a faster accretion rate.
It is useful to present the linearised (in 3 ) unstable circular photon orbit which is found by solving dV dr = 0 and is in agreement with [14] once linearised and appropriate limits have been taken. This simple expression shows that with 3 > 0 and | 3 | 1 the photon orbit is slightly closer to the black hole. This again supports the paradigm of a slower accretion rate. We note that the Innermost Stable Circular Orbit (ISCO) does not change from r = 6m, even without linearising 3 . The fully relativistic analysis of geodesic motion supports the results of our pseudo-Newtonian calculation of the accretion rate, M as well as the back-of-the-envelope calculation in the previous section.

Summary
We examined the problem of accretion onto a static black hole that generically violates the no-hair theorem. The Johannsen-Psaltis spacetime describes a rotating black hole with scalar hair and does not arise from a specific theory of gravity. We restricted our focus to the static limit of the JP metric. This limit contains parametric deviations from the Schwarzschild solution of general relativity.
We studied a Newtonian formulation of the accretion problem using the method of Faraoni et al [29]. This technique generalises the Paczyński-Wiita potential to general static spherically symmetric black hole spacetimes. We obtained a pseudo-Newtonian potential that approximates many of the gravitational features of the static JP metric.
We modelled the matter accreting onto a static JP black hole as a polytropic fluid. The fluid is at rest far from the black hole and accelerates radially inwards at subsonic speeds. The fluid's speed reaches its local sound speed at a critical point. Thereafter the fluid continues to accelerate towards the event horizon supersonically.
In order to determine the critical point we solved a fourth order polynomial equation numerically. Only one of the four roots was found to be real, positive and outside the event horizon. This root was identified as the critical point of the transonic flow. We linearised the quartic equation about 3 and obtained an approximate analytical expression for the critical point's position. We then obtained an analytical expression for the critical velocity of the fluid in terms of its sound speed at "infinity". The critical velocity is independent of 3 .
We obtained an analytical expression for the accretion rate,Ṁ. Our accretion rate is proportional to the square of the black hole mass i.e.Ṁ ∼ m 2 . This is a common feature in Bondi accretion. The accretion rate depends on the black hole mass, m, the fluid's polytrope index, n, the fluid sound speed at infinity and the JP parameter, 3 . Positive (negative) values of 3 were found to reduce (increase) the accretion rate and the accretion rate increased with γ . These results were tabulated.
We obtained and solved the geodesic equations for massive and massless particles orbiting a static JP black hole. Positive values of 3 were found to slow down massive particles in radially infalling trajectories around a static JP black hole. Similarly, negative values of 3 accelerated massive particles. Our fully relativistic analysis of geodesic trajectories corroborates our findings on the effect of 3 on the accretion rate viz. that small, positive values of 3 lowered the accretion rate and vice versa.
There is no a priori reason to expect an overall positive or negative sign to scalar hair contributions to a black hole metric. Black holes with scalar hair arising in modified theories of gravity will thus either increase or decrease the efficiency of the conversion of gravitational energy into radiant energy. Our results can be used to vindicate the validity of the pseudo-Newtonian approach to black holes. The problem we investigated was an idealised study of transonic accretion onto a JP black hole. This framework can used to incorporate various important physical phenomena like radiative processes, viscous dissipation, magnetic fields and accretion disks. Further, our results are independent of the theory of gravity, as opposed to comparable results in the literature, such as [37] for example, which assumes the Einstein vacuum equations hold. Perhaps one disadvantage in taking this theory-independent philosophy is that one cannot compute multipole moments as done in [37], and thus cannot deduce anything about the quadrupole moment. This is because the calculation of multipole moments rests fundamentally on a chosen theory of gravity. However, as stated above, we still obtain a framework to analyse a variety of important physical processes, supplementing existing literature with different results from a different perspective.