Kerr-Newman-AdS Black Hole Surrounded by Perfect Fluid Matter in Rastall Gravity

The Rastall gravity is the modiﬁed Einstein general relativity, in which the energy-momentum conservation law is generalized to T µν ; µ = λR ,ν . In this work, we derive the Kerr-Newman-AdS (KN-AdS) black hole solutions surrounded by the perfect ﬂuid matter in the Rastall gravity using the Newman-Janis method and Mathematica package. We then discuss the black hole properties surrounded by two kinds of speciﬁc perfect ﬂuid matter, the dark energy ( ω = − 2 / 3) and the perfect ﬂuid dark matter ( ω = − 1 / 3). Firstly, the Rastall parameter κλ could be constrained by the weak energy condition and strong energy condition. Secondly, by analyzing the number of roots in the horizon equation, we get the range of the perfect ﬂuid matter intensity α , which depends on the black hole mass M and the Rastall parameter κλ . Thirdly, we study the inﬂuence of the perfect ﬂuid dark matter and dark energy on the ergosphere. We ﬁnd that the perfect ﬂuid dark matter has signiﬁcant eﬀects on the ergosphere size, while the dark energy has smaller eﬀects. Finally, we ﬁnd that the perfect ﬂuid matter does not change the singularity of the black hole. Furthermore, we investigate the rotation velocity in the equatorial plane for the KN-AdS black hole with dark energy and perfect ﬂuid dark matter. We propose that the rotation curve diversity in Low Surface Brightness galaxies could be explained in the framework of the Rastall gravity when both the perfect ﬂuid dark matter halo and the baryon disk are taken into account.


INTRODUCTION
The dark matter and dark energy are two unsolved problems in cosmology and particle physics. Today many observations, including Type Ia supernovae, cosmic microwave background (CMB), baryon acoustic oscillations (BAO), weak lensing, rotation curve, larger scale structure etc., have revealed that the dark energy accounts for 73%, the dark matter for 23% and ordinary baryonic matter only for 4% of the total mass-energy of the universe (e.g., Riess et al. 2004;Komatsu et al. 2011). Following these observations, many theoretical models have been proposed (e.g., Peebles & Ratra 2003;Copeland et al. 2006;Wang et al. 2017;Li et al. 2011). For example, the phenomenological dark energy models include quintessence, phantom, quintom and other models. The dark matter includes Cold Dark Matter (CDM) (e.g., Dubinski & Carlberg 1991;Navarro et al. 1997Navarro et al. , 1996, Warm Dark Matter (WDM) and Scalar Field Dark Matter (SFDM) (e.g., Viel et al. 2005;Amendola & Barbieri 2006;Goodman 2000;Hu et al. 2000;Hui et al. 2017;Marsh 2016;Peebles 2000;Press et al. 1990;Schive et al. 2014;Sin 1994;Turner 1983). The perfect fluid dark matter (PFDM) model has also been proposed recently (e.g., Guzman et al. 2003;Rahaman et al. 2010Rahaman et al. , 2011. For spiral galaxies, under the assumption of the spherical symmetric mass distribution and perfect fluid matter, one can obtain the solution of the Einstein equation for spherical symmetric gravitational field. When the equation of state for spiral galaxies is around −1/3 (Guzman et al. 2003), the PFDM model can explain the asymptotically flat rotation curves.
The covariant conservation of the energy-momentum tensor plays an important role in Einstein's general relativity (GR). Taking this conservation as starting point, one can derive the conservation of some global defined physical quantities through Noether symmetry theorem. These conserved quantities take the form of integrals of the components of the energy-momentum tensor over appropriate space-like surfaces, which allow at least one of the Killing vectors of the background space-time as their normal. Thus, the total rest energy and rest mass of a physical system are conserved in the context of GR ). On the other hand, some modified Einstein GR have relaxed conservation condition of the energy-momentum tensor. One of the modified gravity theory was introduced by Rastall (1972Rastall ( , 1976, in which the conservation of the energy-momentum tensor changes to be T µν ;µ = λR ,ν . When the space-time is flat, the conservation law becomes the usual formalism. This theory could be understood as Mach principle which represents the mass distribution depending on the energy and mass in the space-time (Majernik & Richterek 2006), and has many interesting applications in cosmology and strong gravity space-time (e.g., Batista et al. 2012;Fabris et al. 2012Fabris et al. , 2015Heydarzade et al. 2016;Lobo et al. 2017). However, the physical relevance of the Rastall gravity is under dispute. Recently, Visser (2017) pointed out that the Rastall gravity is equivalent to the Einstein GR, while  suggested that two gravitational theories are not equivalent, although both work admit that the parameter κλ always represents the rearrangement of the matter sector of the Einstein GR. In this context, it is still worth studying the Rastall gravity in both observational and theoretical aspects.
The physical effects of the dark matter and dark energy on the black hole properties have been investigated recently and attracted increasing attention. The cosmological constant in Schwarzschild black hole space-time has been discussed in Kottler (1918), and later generalized to the rotational black hole by Carter (1973). The dynamics of the quintessence dark energy in Schwarzschild space-time have been obtained by , in which the equation of state ω for perfect fluid matter is constant. That work has been generalized to the Kerr black hole and Kerr-Newman-AdS (KN-AdS) black hole (Toshmatov et al. 2015;. In the case of spherical symmetry, the Reissner-Nordstrom black hole space-time surrounded by perfect fluid matter in the Rastall gravity has been obtained by . In this work, we further generalize the Reissner-Nordstrom black hole space-time to the rotational black hole (KN-AdS black hole) space-time in the Rastall gravity and investigate accordingly the KN-AdS black hole properties.
The outline of the paper is as follows. In Section 2, we introduce the Reissner-Nordstrom black hole solution with perfect fluid matter in the Rastall gravity. In Section 3, we derive the KN-AdS black hole solution using Newman-Janis method and Mathematica package. In Section 4, we analyze the properties of the KN-AdS black hole surrounded by perfect fluid matter in the Rastall gravity, including the energy condition, horizon structure, ergosphere and singularity. In Section 5, we calculate the rotation velocity in the equatorial plane and discuss its applications. The summary is given in Section 6.

REISSNER-NORDSTROM BLACK HOLE IN PERFECT FLUID MATTER
In this section, we introduce the Reissner-Nordstrom black hole solution surrounded by perfect fluid matter in the Rastall gravity. The Rastall gravity is based on the Rastall hypothesis which generalizes the energy-momentum conservation law to the following formalism (Rastall 1972(Rastall , 1976 T where T µν is the energy-momentum tensor, λ is the Rastall parameter which represents the level of the energy-momentum conservation law in gravity theory and whether the Mach principle is satisfied. In this theory, the generalized Einstein field equation is given by where κ = 8πG N is the gravitational constant of the Rastall gravity and G N is the gravitational constant of the Newton gravity. When λ → 0, the field equation reduces to the Einstein field equation in GR.
In , the Reissner-Nordstrom black hole solution in perfect fluid matter is given by where where α is the perfect fluid matter intensity around the black hole, M is the mass of the black hole, Q is the charge of the black hole and ω describes the equation of state defined by ω = p/ρ with p and ρ the pressure and density of the perfect fluid matter, respectively.
The electromagnetic field is encoded by the gauge potential A which is Although the real dark matter or dark energy could not satisfy the perfect fluid condition, the perfect fluid is a good approximation to describe the behaviour of the dark matter and dark energy. When λ → 0, the black hole solution Eq. (3) reduces to the Kiselev formalism . For −1 < ω < −1/3, the solution represents the black hole surrounded by the dark energy Toshmatov et al. 2015;. For ω = −1/3, the solution represents the black hole surrounded by the perfect fluid dark matter (Guzman et al. 2003;Rahaman et al. 2010Rahaman et al. , 2011.

Newman-Janis Method and Kerr-Newman Solution in Perfect Fluid Matter
Using the Newman-Janis method (Azreg-Aïnou 2014; Capozziello et al. 2010;Erbin 2015;Newman & Janis 1965), we generalize the Reissner-Nordstrom black hole to the Kerr-Newman black hole in perfect fluid matter in the Rastall gravity. First, we perform trans-formation for the spherically symmetric black hole space-time metric (Eq. 3) from Boyer-Lindquist (BL) coordinates (t, r, θ, φ) to Eddington-Finkelstein (EF) coordinates (u, r, θ, φ) through and the space-time metric (Eq. 3) becomes For the electromagnetic field part, the gauge potential A in the (u, r) plane becomes ).
The non-zero components of the inverse space-time metric (Eq. 7) are given by In the null (EF) frame, the metric matrix can be written as where the corresponding components are For any point in the black hole space-time, the null vectors of the null tetrad satisfy the relations of l µ l µ = n µ n µ = m µ m µ = l µ m µ = n µ m µ = 0 and l µ n µ = −m µ m µ = 1. In the plane of (u, r), the coordinate transformations are We then perform the transformations of f (r) → F (r, a, θ), g(r) → G(r, a, θ), and Σ 2 = r 2 + a 2 cos 2 θ. In the (u, r) space, the null vectors become From the definition of the null tetrad, the metric tensor g µν in the EF coordinates are given by Through calculation, the covariant metric tensor in the EF coordinates are Using Eq.12 we can obtain the electromagnetic field gauge potential A in the EF coordinates as A = Qr Σ 2 (du − asin 2 θdφ).
Finally we perform the coordinate transformation from the EF to BL coordinates as where In the BL coordinates, the Kerr-Newman black hole solution in perfect fluid matter is therefore given by where and the gauge potential A is accordingly From the Newman-Janis method, we know that this Kerr-Newman black hole space-time metric in the Rastall gravity (Eq. 19) satisfies the Einstein-Maxwell field equation. If there is no perfect fluid matter around the black hole (α = 0), this solution will reduce to the usual Kerr-Newman black hole without dark matter. If κλ = 0, it will reduce to the Kerr-Newman black hole in perfect fluid matter in the Einstein GR .

KN-AdS Solution in Perfect Fluid Matter
Now we extend the Kerr-Newman black hole to the KN-AdS black hole surrounded by the perfect fluid matter. Since the Newman-Janis method does not include the cosmological constant Λ, we employ other method to obtain the solution with the presence of the cosmological constant, as presented in . First, we rewrite the above Kerr-Newman black hole metric (Eq. 19) as We then guess the KN-AdS black hole solution as where The solution should satisfy the Einstein-Maxwell field equations which are Through calculations using the Mathematica package ( inserting Eq. (23) into G µν (Eq. 25)), we obtain the Einstein tensors in the Rastall gravity as where Generically the Kerr black hole in GR is a vacuum solution and hence the energymomentum tensor is zero. A non-zero cosmological constant Λ in GR incorporates a static geometry which is known to only shift its vacuum energy by a constant, while has no effect on the energy-momentum tensor. Analogically, in the Rastall gravity, the non-zero Λ should not appear in the explicit expression of the Einstein tensor (or energy-momentum tensor) components either, in the case of the KN-AdS black hole space-time surrounded by the perfect fluid matter. This is exactly what we obtained (Eq. 26). Therefore, our guessed solution (Eq. 23) satisfies the Einstein-Maxwell field equation (Eq. 25) with the perfect fluid matter in the Rastall gravity. Similarly, the gauge potential A in the KN-AdS black hole does not involve Λ and has same expression as in the case of Kerr-Newman black hole (Eq. 21), therefore it also satisfies the Einstein-Maxwell equation (Eq. 25).
In order to compute the energy-momentum tensor, we give the following tetrad The non-zero components of the energy-momentum tensor are given by Rg rr + κλg rr R + Λg rr ), We then find that It is clear that the distribution of the perfect fluid matter around the black hole (Eq. 29) is determined by the parameters of a, Λ, Q, α, ω and κλ. The space-time metric (Eq. 23) describes the perfect fluid matter around the KN-AdS black hole in the Rastall gravity. The equation of state ω determines the material form. As described in Section 2, if −1 < ω < −1/3, the perfect fluid matter represents the dark energy. If ω = −1/3, the perfect fluid matter stands for the perfect fluid dark matter. Using this space-time metric, we can study the interaction between the perfect fluid dark matter (or dark energy) and black hole. Such topic is very interesting in regard of the dark matter (dark energy)-black hole systems. From Eq. (29), we can see that the Rastall parameter κλ plays an important role in the distribution of perfect fluid matter around the black hole. Specifically, the Rastall gravity makes the perfect fluid matter to be rearranged around the black hole. If the Rastall gravity is equivalent to the Einstein GR, κλ only describes the re-arrangement of the perfect fluid matter. If the Rastall gravity is not equivalent to the Einstein GR, the meaning of κλ is the same, but the physical origin of κλ is related to the Mach principle.

Energy condition
In this part, we study the energy condition of the KN-AdS black hole in perfect fluid matter in the Rastall gravity. First, we consider the diagonal energy-momentum tensor T µν in the standard locally non-rotating frame (LNRF, Bardeen et al. 1972). The weak energy condition (WEC) can be written as T µν u µ u ν ≥ 0, where u ν is the time-like vector. The physical meaning of the WEC is that the measured total energy density of all matter fields for any observer traversing a time-like curve is never negative. In the LNRF, the WEC corresponds to ρ ≥ 0 and ρ + P i ≥ 0, where i = r, θ, φ. From Eq. 29 the WEC requires that which is independent of the black hole spin a. Thus, the WEC for the KN-AdS black hole in perfect fluid matter in the Rastall gravity is the same as that in the case of the spherically symmetric black hole ). Since α is always positive, the condition (3κλ(1 + ω) − 3ω)(1 − 4κλ) ≥ 0 will always be satisfied. If considering dark energy (ω = −2/3), then this condition will result in −2 ≤ κλ ≤ 1/4. If considering dark matter (ω = −1/3), then −1/2 ≤ κλ ≤ 1/4. The strong energy condition (SEC) requires that the Raychaudhuri equation is (T µν − 1 2 T g µν )u µ u ν ≥ 0. In the LNRF, the Raychaudhuri equation becomes Y = ρ + P r + P θ + P φ .
From (Eq. 29) the SEC implies that This equation describes the condition that α and κλ should satisfy when SEC is satisfied. From Eq . 31 we find that ρ + P r + P θ + P φ = 0 has one zero point, then Y = ρ + P r + P θ + P φ is in the range [r 1 , ∞) in the KN-AdS space-time with perfect fluid matter in the Rastall gravity (where r 1 is the root of Eq . 31), but at the same time, κλ must satisfy the condition (3κλ(1 + ω) − 3ω)(1 − 4κλ) ≥ 0 which is similar to the WEC condition (Eq . 30). Therefore, the range of κλ under the SEC is consistent with the case of WEC.

Horizon structure
Black hole properties are determined by the horizon structure of the black hole. For stationary axisymmetric black hole, the horizon definition is where the property of horizon depends on α, a, Q, Λ, κλ and ω. Especially the parameters of κλ and ω change the horizon significantly. The number of horizon is determined by ω. As examples, we discuss the cases of the dark energy (ω = −2/3) and dark matter (ω = −1/3). We can set the cosmological constant Λ = 0 given that it is small.
Case I: for the dark energy, the horizon equation reduces to Three horizons exist, including inner horizon (r − ), event horizon (r + ) and cosmological horizons (r q ), where r q is determined by the dark energy. Then the above equation should have two extreme value points. Through calculation, we find that α satisfies the condition of When M = 1 and κλ = 0, the above condition reduces to the Kiselev situation. In other cases, the maximum value of α decreases with the increasing κλ when 0 < κλ and κλ = 1, and increases with the increasing of | κλ | when −2 < κλ < 0.
Case II: for the perfect fluid dark matter, the horizon equation reduces to Dark matter does not produce new horizon and only two normal horizons exist, i.e., inner horizon r − and event horizon r + . Then the horizon equation should have one extreme value point. Through calculation, we find that α satisfies the condition of When M = 1 and κλ = 0, the above condition reduces to the Kiselev situation. Similarly, the maximum value of α decreases with the increasing κλ when 0 < κλ < 1 and κλ = 1/2, and increases with the increasing | κλ | when κλ < 0.
These results reveal that κλ would constrain α in the Rastall gravity. As mentioned in Section 1 and 3, the physical origin could have two possibilities. If the Rastall gravity is equivalent to the Einstein GR, κλ only describes the re-arrangement of the perfect fluid matter. If the Rastall gravity is not equivalent to the Einstein GR, the meaning of κλ is the same, but its physical origin is related to the Mach principle. Furthermore, the Mach principle would determine the value of α.

Stationary limit surfaces
The stationary limit surface of the KN-AdS black hole is important because the rotational energy of black hole relates to the ergosphere size. From the black hole metric (Eq. 23), the ergosphere is defined by Through calculation, Eq. (37) becomes Λ 3 r 2 (r 2 + a 2 ) − a 2 cos 2 θ + Λ 3 a 4 sin 2 θcos 2 θ − r 2 + 2Mr − Q 2 + αr 1 − 3ω 1 − 3κλ(1 + ω) = 0. (38) This equation has two roots. There is an ergosphere between the event horizon and the static limit surface. The parameters a, ω, Q, α, Λ, θ and κλ could affect the shape of the ergosphere. In Fig. 1 and Fig. 2, we show the shapes of the ergosphere for dark energy (ω = −2/3) and perfect fluid dark matter (ω = −1/3), respectively. We find following properties of the ergosphere: (1) Its size decreases with the increasing α, indicating that dark energy and perfect fluid dark matter reduce the rotational energy of the black hole; (2) Its size increases with the increasing κλ in the effective range. These results imply that the Mach principle or the re-arrangement of matter makes the rotational energy of the black hole larger.

Singularity
According to the Einstein GR, the stationary axisymmetric black hole space-time has the singularity. We consider whether or not the singularity changes with the presence of the where H is the polynomial function. The singularity occurs at r = 0 and θ = π/2 when Σ 2 = r 2 + a 2 cos 2 θ = 0. In the Boyer-Lindquist coordinates, the singularity is a ring with a radius of a in the equatorial plane. For different ω and α, there are different polynomial functions, but Σ 2 is the same as in the usual Kerr black hole, so the perfect fluid dark matter and dark energy do not change the singularity of the black hole.

ROTATION VELOCITIES IN THE EQUATORIAL PLANE AND ROTATION CURVE DIVERSITY
Basing on the KN-AdS solution in perfect fluid matter in the Rastall gravity (Section 3), we now discuss the rotation curves in the equatorial plane (θ = π 2 ). We first derive the rotation velocity expression using ordinary method, then we calculate the rotation velocities with the presence of the dark energy and perfect fluid dark matter. Especially for the case of perfect fluid dark matter, we can explain the diversity of rotation curve if considering both the perfect fluid dark matter and a baryon disk which is expressed as where Σ 0 is the central surface density and r d is the scale radius.
In the four-dimensional space-time, the four-velocity satisfies the normalized condition for zero angular momentum observer (ZAMO) as g µν u µ u ν = −1. In the axisymmetric spacetime, the space-time symmetry leads to two conserved quantities L and E. In our KN-AdS black hole space-time, we write the normalized condition from the expressions of u µ and u ν as We then obtain Combining the equations of (40) and (41), we get In the black hole space-time, the general conditions for a stable circular orbit are dr dτ = 0 and dV 2 dr = 0, we then get the expressions of L and E as , Ω φ = −g tφ,r + (g tφ,r ) 2 − g tt,r g φφ,r g φφ,r .
From the definition of rotation velocity, we obtain Fig. 3 and Fig. 4 show the curves of the rotation velocity for the KN-AdS black hole surrounded by the dark energy and perfect fluid dark matter, respectively, in the Rastall gravity. We find that: (1) In the case of dark energy, the rotation velocity always decreases with the increasing distance from the black hole and becomes flatter with the decreasing κλ.
(2) In the case of perfect fluid dark matter, the influence of κλ on the rotation velocity is minor at small distances and significant at larger distances. Note that the effect at larger distances is not shown in the figures, but can be inferred from the explicit expression of the rotation velocity (Eq. 46).
As an example of the application of the Rastall gravity, we try to explain the rotation curve diversity of Low Surface Brightness (LSB) galaxies. In the Rastall gravity, the parameter κλ represents the re-arrangement of the perfect fluid matter around the black hole. From Section 3, for the perfect fluid dark matter halo (ω = −1/3), we can calculate the energy density ρ through the Einstein equation. Because the motion velocity of the dark matter particle is much smaller than the speed of light, the energy density of the perfect fluid dark matter then approximates to the mass density. On the other hand, the rotation curve is usually observed at large distances from the black hole, so that the black hole spin a then can be approximated to zero. From Kamada et al. (2017), the baryon matter can be considered as an index disk, i.e., ρ b = Σ 0 exp[−r/r d ]δ(z), in galaxies. Using the perfect fluid dark matter halo and the baryon disk, we can derive the rotation curve expression. In the equatorial plane, the general Mass function is M(r) = 4π r 0 r 2 ρ DM dr + 2π r 0 rρ b dr. We then obtain the Mass function in our case as The rotation velocity is determined by α, Σ 0 , r d and κλ. Varying κλ would result in different expressions of v, which may explain the rotation velocity diversity in LSB galaxies. We stress two points on the Rastall gravity. If the Rastall gravity is equivalent to the Einstein GR, the parameter κλ only describes the re-arrangement of the perfect fluid dark matter. If the Rastall gravity is not equivalent to the Einstein GR, the meaning of κλ is the same, but the physical origin of κλ is related to the Mach principle.

SUMMARY
We have obtained the KN-AdS black hole solutions surrounded by perfect fluid matter in the Rastall gravity using the Newman-Janis method and complex calculations. In this work, we focus on the cases of the dark energy (ω = −2/3) and the perfect fluid dark matter (ω = −1/3). We discuss different black hole properties in the Rastall gravity with the presence of the perfect fluid matter. We find that if the condition (3κλ(1 + ω) − 3ω)(1 − 4κλ) ≥ 0 is satisfied, the WEC and SEC can always be hold everywhere in the KN-AdS black hole.
By analyzing the horizon equation, we can constrain the values of the parameter α, which is determined by the equation of state ω and the Rastall parameter κλ. In the case of the dark energy, the maximum value of α decreases with the increasing κλ when 0 < κλ and κλ = 1, and increases with the increasing | κλ | when −2 < κλ < 0. In the case of the perfect fluid dark matter, the maximum value of α decreases with the increasing κλ when 0 < κλ < 1 and κλ = 1/2, and increases with the increasing | κλ | when κλ < 0. The perfect fluid dark matter and dark energy reduce the ergosphere size, but they don't change the singularity of the black hole.
We also derive the rotation velocity for the KN-AdS black hole surrounded by the perfect fluid matter in the Rastall gravity. In the case of dark energy, the rotation velocity decreases with the increasing distance from the black hole and becomes flatter with the decreasing κλ. In the case of perfect fluid dark matter, the effect at small distances is minor and become significant at galactic scale. Considering both the perfect fluid dark matter halo and the baryon disk, we have the possibility of explaining the rotation curve diversity in LSB galaxies.
We acknowledge the anonymous referee for a constructive report that has significantly improved this paper. We acknowledge the financial support from the National Natural Science Foundation of China through grants 11503078, 11573060 and 11661161010.