Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Kerr-Newman-NUT-Kiselev Black Holes in Rastall Theory of Gravity and Kerr-CFT Correspondence PDF

Download as pdf or txt
Download as pdf or txt
You are on page 1of 25

Kerr-Newman-NUT-Kiselev black holes in Rastall theory of gravity

and Kerr/CFT Correspondence


Muhammad F. A. R. Saktia,∗, Agus Surosoa,b , Freddy P. Zena,b
a
Theoretical Physics Laboratory, THEPI Division, Institut Teknologi Bandung, Jl. Ganesha 10 Bandung,
40132, Indonesia
b
Indonesia Center for Theoretical and Mathematical Physics (ICTMP), Institut Teknologi Bandung, Jl.
Ganesha 10 Bandung, 40132, Indonesia
arXiv:1901.09163v3 [gr-qc] 26 Dec 2019

Abstract
We present a new twisted rotating black hole solution by performing Demiański-Newman-
Janis algorithm to the electrically and dyonically charged black hole with quintessence in
Rastall theory of gravity. Using our black hole solution, we argue that Rastall gravity is
not equivalent with Einstein gravity. For further explanation, the black hole properties such
as the horizon and ergosphere are studied for which there are some different properties for
those theories. Some thermodynamic properties of the black hole solution are also discussed.
At the end, considering the Kerr/CFT correspondence is valid for our black hole solution,
the central charge from the CFT of this extremal solution is derived.
Keywords: Demiański-Newman-Janis algorithm, Rastall gravity, thermodynamics,
Kerr/CFT correspondence

1. Introduction
The Einstein theory of gravity till nowadays has been confirmed by some observational
evidence in our universe. This theory considers the covariant conservation of the energy-
momentum tensor. However, the Einstein theory is also believed as a certain case of some
extended theories of gravity. Since the first time it was formulated, people have been working
for some other alternative theories and developing several modified theories of gravity. One
of the promising general theories of relativity was introduced by Rastall [1, 2]. Within
this theory the usual conservation law of the energy momentum tensor is not obeyed or
on the other hand, the energy-momentum tensor is not conservative (∇µ T µν 6= 0). The
usual conservation law, expressed by the null covariant derivative of the energy-momentum
tensor, is questioned whether it is valid generally in curved space-time or only in certain cases
since this theory is tested only in a weak gravitational field limit or in the flat Minkowski


Corresponding author
Email addresses: m.fitrah@students.itb.ac.id (Muhammad F. A. R. Sakti),
agussuroso@fi.itb.ac.id (Agus Suroso), fpzen@fi.itb.ac.id (Freddy P. Zen)
Preprint submitted to Annals of Physics December 30, 2019
spacetime. It is considered a non-minimal coupling of the matter field to space-time geometry
with the coupling constant κλ or Rastall parameter. This quantity quantifies the deviation
from the Einstein theory of gravity. When Rastall parameter vanishes, it produces the
Einstein theory of gravity. Furthermore, in [3, 4] it is argued that Rastall gravity is a special
case of f (R, T ) gravity.
Rastall theory of gravity might be interpreted as a direct accomplishment of the Mach
principle which suggests that the inertial properties of a mass distribution are determined
by the distribution of mass-energy in the external space-time [5]. Thus, the source of grav-
itation, either mass-density or the elements of energy-momentum tensor, depends on the
gravitational tensor. We can find also some recent research work that use this gravitational
theory to explain the cosmological aspects such as the accelerating expansion of the universe
and the inflationary problems [6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16]. In the smaller scale, some
astrophysical configurations are investigated by performing this theory. For instance, the in-
vestigation of perfect fluid spheres, compact stars, neutron stars, black holes, and wormholes
is done in the context of Rastall gravity [17, 18, 19, 20, 21, 22, 23, 24, 25]. Corresponding
to black hole solutions in Rastall gravity, a fascinating non-commutative inspired black hole
solution is obtained in [26]. In addition, the black hole solution with the source of a Gaussian
matter distribution is obtained in [27]. Several extensions of black hole solution in Rastall
gravity are further investigated in [28, 29, 30, 31, 32]. Besides that, the thermodynamic
properties of black hole solutions in the Rastall gravity are discussed in [33].
Recently, it is claimed by Visser [34] that Rastall theory of gravity is equivalent with
Einstein theory. However, Darabi et al. [35] compare these two gravitational theories and
summarize that Visser’s conclusion is not correct. The argument in [35], indeed, supports
Rastall theory of gravity for which this theory is still an open theory comparing to the
usual general relativity. Henceforth, this theory may face the challenges of cosmological
observations as the general relativity.
In recent years, it is found that our universe experiences accelerated expansion due to
the existence of dark energy that fills more than 70% matters of our universe. A promising
model of dark energy is the quintessence. Quintessence model is a model of a scalar field that
governs the pressure of our universe to be negative then makes the accelerated expansion
[36]. The study involving the scalar field in some frame of modified gravity theories in a
cosmological scale can be discovered in [37, 38, 39, 40, 41, 42, 43, 44, 45]. The existence
of quintessence possibly causes the change of the structure of some locally astrophysical
manifestations of it such as black holes. First investigation of the quintessential effect on
the black hole originally is studied by Kiselev [46]. It is found that the black hole solution is
dependent on the quintessential intensity α and equation of state ω of this type of field. The
equation of state of the quintessence could be varied that depends on the domination of the
matter. The value range of the equation of state ω is −1/3 < ω < 0 which makes the black
hole solutions is asymptotically flat and −1 < ω < −1/3 which represents the accelerating
expansion of universe. But it is still possible to set ω for some special cases. In addition,
black hole solution with the quintessential matter in Lovelock gravity is investigated in [48]
while in Gauss-Bonnet gravity can be found in [49].
In the following, we derive an electrically and dyonically charged black hole solution
2
in Rastall theory of gravity with the existence quintessential matter. We add a magnetic
charge in this solution rather than adding the electric charge only [28] to make it more
general and provide an interesting feature related to the magnetic monopole. In addition,
the twisted rotating solution is more fascinating because there exist spin a and the NUT
charge. Regarding NUT charge as a twisting parameter of the surrounding space-time, one
can find the explanation in [62]. Another definition of this parameter can be discovered in
[63] for which it tells that NUT charge represents a gravo-magnetic monopole parameter
of the central mass. Herein we employ the Demiański-Newman-Janis algorithm [64, 65],
the more general one than the Newman-Janis algorithm, to insert the spin a together with
the NUT charge n. Hence we present a novel solution so-called the Kerr-Newman-NUT-
Kiselev black hole solution in Rastall gravity. The solution that we obtain will be highly
favorable to study the astrophysical aspects of the magnetic monopole in relation with
the quintessential matter. We also give some arguments to say that Rastall gravity is
not equivalent with Einstein gravity. Some properties and thermodynamic quantities of this
solution are also investigated. It will be interesting to make similar studies in other modified
gravity models for example f (R) gravity [50], Hořava-Lifshitz gravity [51] and bumblebee
gravity [52, 53, 54, 55, 56, 57, 58, 59, 60] which have been getting significant attention in
many different aspects.
In a recent investigation, the microscopic entropy of black holes may be derived by
governing the Kerr/CFT correspondence which is originally proposed in [66]. This entropy
merely agrees with the Bekenstein-Hawking entropy that obeys the area law of black holes
[67]-[98]. The famous Cardy entropy formula is applied in this fashion where it is a function
of the central charge and conformal temperature. The central charge comes from the central
term in which it comes from the Dirac bracket of the canonically conserved charges associated
with the non-trivial diffeomorphisms of the near-horizon region. This conserved charge
corresponds with the theory that is implemented in [99, 100] for Einstein-Maxwell system
and for several modified gravitational theories is reviewed in [101] and the references therein.
However, for Rastall theory of gravity, the asymptotic conserved charged has not been
derived. Because of this, we assume that the Kerr/CFT correspondence is true in this
theory to obtain the central charge of the solution that we find from the Bekenstein-Hawking
entropy.
We set up the remaining parts of the paper as follows. In section 2, we derive the
spherically symmetric charged black hole solution in Rastall gravity. After that, to find the
twisted rotating solution, Demiański -Newman-Janis algorithm is implemented. Moreover,
we give some comments regarding the equivalency between Rastall gravity and Einstein
gravity. In section 3, we study the horizon and ergosphere of the black hole solution. In
the next section, several thermodynamic properties are derived. Then the entropy from the
previous section is used to find the central charge of the black hole solution. Finally, the
summary is provided in the next section. The resulting solution can be used to argue that
our solution in Rastall gravity is not equivalent with Einstein gravity.

3
2. Kerr-Newman-NUT-Kiselev black hole in Rastall gravity
Rotational parameter or spin a and twisting parameter or NUT charge n can be tucked
into the spherically symmetric solution using Demiański-Newman-Janis algorithm [65]. This
is the extension of Newman-Janis algorithm. We employ this algorithm to the dyonic
Reissner-Nordström black hole surrounded by the quintessence in Rastall theory of grav-
ity.

2.1. Dyonic Reissner-Nordström black hole with Quintessence in Rastall gravity


Based on Rastall’s theory [1, 2], the energy-momentum tensor is not conserved and
its covariant derivative is proportional to the derivative of Ricci scalar with the coupling
constant λ, namely Rastall parameter that can be written as

∇µ T µν = λ∇ν R. (1)

For the flat Minkowski space-time or specifically in a weak gravitational field limit, this will
reduce to the usual conservation law. So the Einstein field equation will be modified as

Gµν + κλgµν R = κTµν , (2)

where κ = 8πG and G is the Newtonian gravitational coupling constant.


In order to obtain spherically symmetric black hole solutions, we consider the general
space-time metric in the standard Schwarzschild coordinates as

ds2 = −f (r)dt2 + f (r)−1 dr2 + r2 dΩ2 , (3)

where f (r) is a function depending on radial coordinate only and two-dimensional sphere
dΩ2 = dθ2 + sin2 θdφ2 . Using metric (3) to the r.h.s. of Eq. (2), the non-vanishing compo-
nents of gravitational tensor are given by
1 0
G00 + κλR = (f r − 1 + f ) + κλR, G11 = G00 , (4)
r2  
1 0 1 00 2
G2 + κλR = 2 f r + f r + κλR, G33 = G22 ,
2
(5)
r 2

and the Ricci scalar reads as


1 2 00 0

R=− r f + 4rf − 2 + 2f , (6)
r2
in which the prime represents the derivative with respect to the radial coordinate r.
We have two kinds of matter in this case which are the electromagnetic field (Eµν )
and scalar field of quintessence (Qµν ). As derived in [46], the non-zero components of the
energy-momentum tensor of the quintessence are
1
Q00 = Q11 = −ρq (r), Q22 = Q33 = (3ω + 1)ρq (r). (7)
2
4
ω is the parameter of equation of state of the quintessence to determine matter domination
of the solution. Nonetheless, there is an argument from Visser [47] which tells that Kiselev
solution does not represent quintessence with isotropic pressure. However, Kiselev derives
the solution by assuming the isotropic averaging over the angles of the space-like matter
tensor component. The calculations end up with an equation of state like the quintessence
in cosmology.
Then the energy-momentum tensor of the electromagnetic field is given by
 
2 γ 1 γβ
Eµν = Fµγ Fν − gµν F Fγβ , (8)
κ 4
where Fµν = ∂µ Aν − ∂ν Aµ is the trace-free electromagnetic field tensor. Using Bianchi
identity and the variational principle over the electromagnetic potential, we arrive at

∂[σ Fµν] = 0, ∂µ ( −ḡF µν ) = 0. (9)

Regarding the spherical symmetry exists in the space-time metric (3), it is imposed the
only non-vanishing components of the electromagnetic field tensor F µν to be F 01 = −F 10
and F 23 = −F 32 . We impose more non-zero component of the electromagnetic field tensor
because we wish to obtain a more general solution than in [28] which only contains the electric
charge. By performing this assumption to Eq. (9) then, the non-vanishing electromagnetic
potential is given by
e
Aµ dxµ = − dt + g cos θdφ, (10)
r
where e, g are defined as electric and magnetic charges, respectively. The non-zero compo-
nents of the electromagnetic energy-momentum tensor are given by
1 e2 g 2 1 e2 g 2
   
0 1 2 3
E0 = E1 = − + , E2 = E3 = + . (11)
κ r4 r4 κ r4 r4
Making use of Eqs. (5), (7) and (11) to Eq. (2), we may derive two equations that are
given by
1 0 κλ 2 00 0
 e2 g 2
(f r − 1 + f ) − r f + 4rf − 2 + 2f = −ρ q − − , (12)
r2 r2 r4 r4
e2 g 2
 
1 0 1 00 2 κλ 2 00 0
 1
f r + f r − r f + 4rf − 2 + 2f = ρq (3ω + 1) + + . (13)
r2 2 r2 2 r4 r4
Now Eqs. (12) and (13) are able to be solved. The solution to these equations is

2M e2 + g 2 1+3ω−6κλ(1+ω)
f (r) = 1 − + 2
− αr− 1−3κλ(1+ω) , (14)
r r
and the quintessence energy density that reads as
−2Ws α
ρ(r) = , (15)
κrKs
5
where
(1 − 4κλ)[κλ(1 + ω) − ω] 3(1 + ω) − 12κλ(1 + ω)
Ws = − , Ks = .
[1 − 3κλ(1 + ω)]2 1 − 3κλ(1 + ω)

Regarding ρq (r) ≥ 0, hence we have Ws α ≤ 0. On the other hand, for Ws > 0 the quintessen-
tial intensity α < 0 and vice versa. This is the dyonic Reissner-Nordström black hole with
quintessence in Rastall gravity that generalizes the solution obtained in [28] because of the
presence of the magnetic charge g. The dyonic black hole solutions provide special interest
as it allows one to study spherical configuration with magnetic monopole [121, 122].

2.2. Rotating, twisting, charged black hole with quintessence in Rastall theory of gravity
We will derive the more general solution than the electrically charged rotating black
hole solution surrounded by the quintessence in Rastall gravity that is obtained in [31].
Demiański-Newman-Janis algorithm, the extended Newman-Janis algorithm, will be per-
formed to the solution that we have gained in the previous section. It is firstly proposed in
[106] and exposed in detail in [65]. It is calculated in [65] that this algorithm satisfies the
Einstein equation. In addition, Newman-Janis algorithm is valid to derive Kerr-Newman-
AdS black hole solution surrounded by the quintessential matter in Rastall gravity [31]. So
it denotes that this algorithm may be applied in our spherically symmetric solution. In this
manner, the tetrad formalism is not applied but we make use of only the complex coordinate
transformation as shown in [65]. Starting the algorithm, the spherically symmetric metric
(3) is written in more general form as

ds2 = −ft (r)dt2 + fr (r)dr2 + fs (r)dΩ2 , (16)

where here dΩ2 = dθ2 + H 2 (θ)dφ2 and



 sin θ, k = 1,
H(θ) = 1, k = 0, (17)
sinh θ, k = −1.

The constant k is defined as the sign of the surface curvature. The electromagnetic vector
potential is written in the form

Aµ dxµ = fa dt + g cos θdφ. (18)

However because of the vector potential (10) contains no effect from the quintessence and
Rastall parameter, the later potential will remain the same as the usual Kerr-Newman-NUT
black holes [80]. So we stress the derivation only for the metric.
Using this null coordinate transformation
s
fr
dt = du + dr, (19)
ft

6
makes the metric (16) becomes
p
ds2 = −ft du2 − 2 ft fr dudr + fs dΩ2 . (20)

After gaining the above form, we need to make use of the complex coordinate to the coor-
dinate u and r, that is given by

r = r̂ + iF (θ), u = û + iG(θ), (21)

where û, r̂ ∈ R, and F (θ), G(θ) are two arbitrary functions. To find the twisting and rotating
solution, F (θ) and G(θ) have form as

F (θ) = n − a cos θ, G(θ) = a cos θ − 2nln(sin θ), (22)

where a, n are the spin and NUT charge, respectively. Along with the coordinate transfor-
mation, the coordinate r and mass M are required to be complexified by obeying
 
1 2 2 M 1 M M̄ M r̄ + M̄ r̂
r̂ → (r̂ + r̄) = Re(r̂), r̂ → |r̂| , → + = , (23)
2 r̂ 2 r̂ r̄ |r̂|2

where the mass is complexified to M → M + in. Note that r̄ is the conjugate of r̂ and M̄
is the conjugate of M . Next, the differential form of (21) is

dr = dr̂ + iF 0 (θ)dθ, du = dû + iG0 (θ)dθ, (24)

where the prime on the function F (θ), G(θ) denotes the derivative respects to coordinate θ.
Then we have to use Giampieri’s ansatz [107] on the angular coordinate, i.e.

idθ = H(θ)dφ. (25)

So, Eq. (24) becomes

dr = dr̂ + F 0 (θ)H(θ)dφ, du = dû + G0 (θ)H(θ)dφ. (26)

Hence we may replace ft , fr , fs by f¯t , f¯r , f¯s where it shows that the metric functions
are transformed by complex coordinate now. Finally we obtain the twisting and rotating
solution in Eddington-Finkelstein coordinates as follows
q q
ds = −2 ft fr (dûdr̂ + G Hdr̂dφ) − 2(ft G H + f¯t f¯r F 0 H)dûdφ
2 ¯ ¯ 0 ¯ 0
q
− ft dû + fs dθ − (ft G H + 2 f¯t f¯r F 0 G0 H 2 − f¯s H 2 )dφ2 .
¯ 2 ¯ 2 ¯ 02 2
(27)

In order to find the solution in Boyer-Linquist coordinates, we need to transform metric (27)
with the following coordinates transformation

dû = dt̂ − g(r̂)dr̂, dφ = dϕ − h(r̂)dr̂, (28)


7
where
1
(f¯t f¯r )− 2 f¯s − F 0 G0 F0
g(r̂) = f¯s
, h(r̂) = ¯ . (29)
f¯r
+ F 02 H fs
+ F 02
f¯r

It is worth mentioning that the functions g(r̂), h(r̂) cannot be dependent of coordinate θ
because it can make the transformations (28) unintegrable [65]. At the end, after performing
(28), the general twisted rotating metric in Boyer-Linquist coordinates reads as

f¯s
ds2 = −f¯t (dt + χHdϕ)2 + dr2 + f¯s (dθ2 + σ 2 H 2 dϕ2 ), (30)

where
s
f¯r 0 2 f¯r 02 f¯s
χ = G0 + F ,σ = 1 + F ,∆ = + F 02 ,
f¯t f¯s f¯r

and the hat (ˆ) on t, r has been omitted. This metric is invariant under transformation
(F, G, ϕ) → −(F, G, ϕ).
Regarding the dyonic Reissner-Nordström solution with the quintessence in Rastall grav-
ity with f (r) that is given in Eq. (14), we obtain the Kerr-Newman-NUT-Kiselev black hole
solution in Rastall gravity that is given by

∆ 2 ρ2 2 sin2 θ  2
2
ds = − 2 dt − (a sin θ + 2n cos θ)dϕ + dr + ρ dθ + 2 adt − (r2 + a2 + n2 )dϕ ,
2 2 2
ρ ∆ ρ
(31)
where
1 − 3ω
∆ = r2 − 2M r + a2 + e2 + g 2 − n2 − αrυ , υ = , ρ2 = r2 + (n − a cos θ)2 .
1 − 3κλ(1 + ω)

The more general twisting and rotating solution can be gained by governing the different
function of F (θ) and G(θ) but it is not implemented in [65]. By performing
 
θ
F (θ) = −n − a cos θ, G(θ) = a cos θ + 2nln(sin θ) − 2nln tan , (32)
2

we find g(r̂) = [r2 + (a + n)2 ]∆−1 and h(r̂) = a∆−1 , so our coordinate transformation (29)
is integrable. Then the resulting metric reads as

∆ 2 ρ 2 2
ds2 = − dt − {a sin 2
θ + 2n(1 − cos θ)}dϕ + dr
ρ2 ∆
2 
sin θ 2
+ρ2 dθ2 + 2
adt − {r2 + (a + n)2 }dϕ , (33)
ρ

8
where ρ2 = r2 + (n + a cos θ)2 . The electromagnetic potential related to the metric (33) is
given by [111]

−er adt − a2 sin2 θ + 2an(1 − cos θ) dϕ


  
µ
Aµ dx =
aρ2
g(n + a cos θ) [adt − {r2 + (a + n)2 } dϕ]
− . (34)
aρ2
For g, n = 0, the space-time metric reduces to the case explained in [31]. One can check that
this is the solution of Rastall (modified Einstein) field equation using Mathematica package
RGTC. The Einstein tensor is shown in Appendix A. This metric is more general than the
previous metric (31) and than the metrics that are shown in several articles [28, 31, 108]
which also contains magnetic charge and NUT charge (for vanishing cosmological constant).
The existence of magnetic charge provides an opportunity to study the axially symmetric
black hole solutions with magnetic monopole while the NUT charge is related to the gravo-
magnetic monopole [63]. More special about NUT charge is that its presence constructs
asymptotically non-flat solution, so it will be interesting for further study to compute the
mass and angular momentum using Komar integral whether it will be well-defined or not.
We also construct the energy-momentum tensor from the gravitational tensor as
1
Tµν = (Gµν + κλgµν R) . (35)
κ
The non-zero components of this tensor are shown in Appendix B. Obviously we can see
that all non-zero components are dependent of parameters a, M, e, g, n, α, ω and κλ. As we
know, NUT charge is also known as a twisting parameter [62]. With vanishing spin, the
momentum density component of energy-momentum tensor does not vanish. However, the
magnetic charge has similar effect as electric charge in this matter tensor but it dyonically
charges the solution. Regarding the discussion of gravo-magnetic monopole [63] of NUT
parameter, we may expect that our solution contains two magnetic charges. These charges
come from the electromagnetic theory and gravitational theory which give another feature
of magnetic monopole. It is believed also that in the early universe, the magnetic monopole
existed. The distribution of the quintessential matter around black hole solution (33) is also
determined by all parameters especially NUT charge and Rastall parameter that give more
terms.

2.3. Rastall gravity is not equivalent with Einstein gravity


One of the analysis of Visser [34] to say that Rastall gravity and Einstein gravity are
equivalent is that the term containing Rastall parameter is just a redefinition of cosmological
constant term. For vacuum equation he shows that 14 λRgαβ = 0 for λ 6= 1 and 14 λRgαβ =
−Λgαβ for λ = 1. When matter is added, he shows that 14 λRgαβ = − 4λ−4 λ
TR gαβ for λ 6= 1
1
and 4 λRgαβ = −Λgαβ for λ = 1 where TR is the trace of Rastall’s energy momentum tensor
used by Visser and in his paper, he uses normalization κ = 1/4. If this argument is right,
the solution of those conditions should represent a metric tensor containing cosmological
9
constant or, at least, terms that can be redefined as a cosmological constant. Nonetheless,
for example, the following Kerr-NUT-dS(AdS) black hole
2
ρ2 2 ρ2 2

2 ∆r 2
ds = − 2 2 dt − {a sin θ − 2n(1 − cos θ)}dφ + dr + dθ
Ξρ ∆r ∆θ
2
∆θ sin2 θ

2 2
adt − {r + (a + n) }dφ , (36)
Ξ2 ρ2
where
r2 (r2 + 6n2 + a2 ) 3n2 (a2 − n2 )
∆r = r2 − 2M r + a2 − n2 − − ,
l2 l2
a cos θ(4n + a cos θ) a2
∆θ = 1 + , ρ2 = r2 + (n + a cos θ)2 , Ξ = 1 + 2 ,
l2 l
could be used to analyze whether those two gravitational theories are equivalent or not
because the most important feature is the metric tensor. The reason we choose to com-
pare with the Kerr-NUT-dS(AdS) metric is because this metric contains three important
parameters that we need which are spin a, NUT charge n and cosmological constant 1/l2
(or negative for AdS space).
If Rastall gravity is just a redefinition of matter or cosmological constant, it should be
obvious that our metric could be described as Kerr-NUT-dS(AdS) metric for arbitrary value
of ω, κλ. From this metric, we can see that the cosmological constant is coupled with spin
and NUT charge in the functions ∆r , ∆θ . Let us take a look at the ∆r first. We may write
∆ as !
X 1
∆ = r2 − 2M r + a2 − n2 − αn rυ , υ = P , (37)
n
1 − m κλm
where we set e, g, ω = 0 in this case to simplify the explanation. The constants αn , λm can
be identified as any constants to reproduce the function ∆r , for example
3n2 (a2 − n2 ) r2 (6n2 + a2 ) r4
α0 r υ = , α1 rυ = 0, α2 rυ = , α3 rυ = 0, α4 rυ = . (38)
l2 l2 l2
However, in our metric the function ∆θ does not appear. It cannot be produced using the
redefinition of α, λ. Another argument is that R ∼ T is not a constant but it is dependent
on the coordinate. Ricci scalar of our solution can be seen in Eq. (A.6). So from this first
point of view, Rastall gravity is not equivalent to Einstein gravity. Note that for non-linear
form of scalar Ricci, as the Gauss-Bonnet gravity or f (R) gravity, on the non-conservative
matter tensor equation (39) might produce the equivalence between these two theories.
Nevertheless, the non-linear term is beyond the Ratall theory that can be proposed for
further study as
∇µ T µν = λ∇ν K(R), (39)
where K(R) is the function containing non-linear terms of R [61]. Related to the energy-
momentum tensor (35), the role of the Rastall theory also can be seen. If the statement that
10
Einstein gravity is equivalent to Rastall gravity is right, κλ only describes the rearrangement
quintessential matter, however it does not. Moreover, Rastall gravity might correspond to
Mach principle, if this theory is not equivalent to Einstein gravity. Recently, in [3, 4] it
is argued that Rastall gravity is a special case of f (R, T ) gravity. If we also consider the
argument in [3, 4], so once again, Rastall gravity is not equivalent to Einstein gravity.

3. Some black holes structures


We have explained that our solution is a new solution in Rastall theory. For further
understanding, in this section, we discuss the black holes structures of the black hole solution
obtained in the previous section.

3.1. Black holes horizons


Kerr-Newman-NUT-Kiselev black hole solution in Rastall gravity has the number of the
horizon which is dependent on the value of ω and κλ. Coordinate singularity determines
the horizon which corresponds with the space-time. The horizon is a null hypersurface of
constant r. Thus, the horizons are the roots of the following equation
−1
g11 = r2 − 2M r + a2 + q 2 − n2 − αrυ = 0, (40)

where q 2 = e2 + g 2 . The horizon structure for g, n, κλ = 0 is shown in [108] but here we


provide the more general solution. Performing similar fashion in [110, 31], it is found that
for general value of ω, the quintessential intensity satisfies

2 − 6κλ(1 + ω) 3ω−3κλ(1+ω)
0<α≤ 2 1−3κλ(1+ω) , (41)
1 − 3ω
where we have put M = 1. For κλ = 0, it reduces to the result obtained in [108]. Note
that even the NUT charge is present, it does not affect the constraint of the quintessential
intensity as the spin [108]. Because of the existence of quintessence and besides the inner
and outer horizons, our black hole solutions might possess a so-called cosmological horizon
but it differs from the cosmological horizon when the cosmological constant is non-vanishing.
However, it might exist only one horizon or naked singularity. Each horizon radius of the
black holes will depend on a, q, n, α, ω, κλ and M . In Fig. 1, we show the possibilities of the
number of the horizon.
The top three pictures shows the variation of α and κλ. Negative or positive value of
α, κλ is related to the energy condition of the solution, for instance, the explanation for
non-rotating solution is given in [28]. We compare the horizons for negative and positive
κλ and also with the solution in Einstein general relativity (κλ = 0). For the given value of
each parameter, these pictures show four cases of black hole solutions. Firstly, the red solid
line describes a black hole solution with naked singularity because no horizon is shown. The
yellow line denotes the existence of extremal black hole solution where the inner and outer
horizons coincide each other. The purple line clearly tells us that the black holes only have
two horizons, yet it is generic Kerr-Newman-NUT black hole because α = 0. The green line
11
1.6 1.6 1.6
= -0.2 =0 = 0.2

0.8 0.8 0.8

0.0 0.0 0.0

= -0.5 = -0.5 = -0.5


-0.8 = -0.1 -0.8 = -0.1 -0.8 = -0.1
= 0.5 = 0.5 = 0.5
= 0.1 = 0.1 = 0.1
=0 =0 =0
-1.6 -1.6 -1.6
0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0
r r r
1.6 1.6 1.6
= -0.2 =0 =0 =0 = 0.2 =0
=0 = -1/2 =0 = -1/2 =0 = -1/2
= -2/3 = -2/3 = -2/3
1.2 = -5/6
1.2 = -5/6
1.2 = -5/6

0.8 0.8 0.8

0.4 0.4 0.4

0.0 0.0 0.0

0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0
r r r
1.6 1.6 1.6

0.8 0.8 0.8

0.0 0.0 0.0

-0.8 n = 0.1 -0.8 n = 0.1 -0.8 n = 0.1


n = 0.35 n = 0.35 n = 0.35
n = 0.55 = -0.2 n = 0.55 =0 n = 0.55 = 0.2
n = 0.9 =0 n = 0.9 =0 n = 0.9 =0
-1.6 -1.6 -1.6
0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0
r r r
1.6 1.6
a = 0.2 = 0.01
a = 0.5 = 0.1 0.8
a = 0.77 = 0.4
0.8 a = 0.9 0.8 = 0.8
0.4

0.0 0.0
0.0

-0.8 -0.8 = -1
-0.4 =1
=3
=5
-1.6 -1.6 -0.8
0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0
r r r

Figure 1: The plot showing the behavior of ∆ for M = 1. For the top pictures, we use a = 1, q = 0.64, n =
0.7, ω = −1/2 and vary α, κλ. The second line pictures, we vary ω, κλ and use a = 0.8, q = 0.8, n = 0.4, α =
0.1, except the purple line α = 0. The third line pictures, a = 1, q = 0.65, α = 0.1, ω = −1/2 is taken,
except the purple line α = 0 and we vary n, κλ. The bottom pictures are intended to show a possibility of
having three horizons for q = 0.7, n = 0.65, ω = −1/2. For left bottom picture, we take α = 0.1, κλ = 3 and
vary a. For middle bottom picture, we use a = 0.9, κλ = 3 and vary α. The last is right bottom picture for
which we take a = 0.9, α = 0.1 and vary κλ.

12
Table 1: Analytical inner and outer horizon for some values of ω, κλ.
ω, κλ Horizon (r±q
)
2
M + α2 ± M + α2 + n2 − a2 − q 2

0,0
√ 2 2 2 2
M M −(a +q −n )(1−α)
-1/3,0 1−α
±√ 1−α
M M 2 −(a2 +q 2 −n2 )(1−α)
0, 1/6 1−α
± 1−α
α
 q 2
-1/3, -1/2 M+2 ± M + α2 + n2 − a2 − q 2
p
1/3,0 M ± M 2 + n 2 + α − a2 − q 2

in the left and middle pictures has similar meaning as the purple but not for the right one.
The green line in the right picture, together with the blue line in every top picture shows
that the black hole solution only has a single horizon.
The second and third line pictures are intended for the similar reasons such as the
top pictures. Yet we vary the different parameters which are ω and n. We also compare
with the Kerr-Newman-NUT black holes (purple line). If we carefully compare the middle
pictures with the right pictures for both variation, we can see that the existence of Rastall
coupling constant makes the number of the horizon decreases (see the yellow, blue, green
lines). Note that ω = 0 is for the dust domination and ω = −1/2, −2/3, −5/6 is for
quintessential field domination. The bottom pictures explain that we can have the black
holes with three horizons. On the other hand, there might exist the cosmological horizon
(rq ) with r− ≤ r+ ≤ rq . If there exist three horizons, so the probabilities of the extremal
case are r− = r+ , r− = rq , r+ = rq , and r− = r+ = rq . It is clear from the left and middle
bottom pictures, the green line shows that it fulfills an extremal case r− = r+ . Though
Rastall coupling constant could make the number of horizon decreases, we may see from the
bottom right picture when κλ = −1, the horizons are only two. However, when κλ = 1 or
κλ = 3, the cosmological horizon then appears.
For some values of ω, κλ, the horizon coincides to have only inner and outer parts or
has no cosmological horizon. In this case, some thermodynamic properties are easier to be
investigated. That is why in [112], they consider only two roots of the horizon. For more
than two roots of the horizon, it will bring to more complicated computation for example to
study the entropy product. We show some circumstances that will give rise to two analytic
roots in Table 1. For cases ω, κλ = 0, 0; −1/3, −1/2; 1/3, 0, the existence of the quintessence
increase the horizon radius because of the positive quintessential intensity α and vice versa.
For cases −1/3, 0; 0, 1/6, α plays a role as the denominator, so α should not be equal to 1
or on the other hand, the horizon radius will be infinite.

3.2. Ergoregion of the black holes


For rotating black hole, there exists a surface besides the event horizon, i.e. the static
limit surface. It corresponds to the time-like hypersurface where the time translational
Killing vector become null, so it means g00 = 0. For the Kerr-Newman-NUT with the

13
a=0.2,q=0.4,n=0.4, =0.1, =-1/3, =1/6 a=0.6,q=0.4,n=0.4, =0.1, =-1/3, =1/6 a=0.9,q=0.4,n=0.4, =0.1, =-1/3, =1/6
3 3 3

2 2 2

1 1 1
z/M

z/M

z/M
0 0 0

-1 -1 -1

-2 -2 -2

-3 -3 -3
-3 -2 -1 0 1 2 3 -3 -2 -1 0 1 2 3 -3 -2 -1 0 1 2 3
x/M x/M x/M
a=0.4,q=0.1,n=0.4, =0.1, =-1/3, =1/6 a=0.4,q=0.6,n=0.4, =0.1, =-1/3, =1/6 a=0.4,q=0.9,n=0.4, =0.1, =-1/3, =1/6
3 3 3

2 2 2

1 1 1
z/M

z/M

z/M
0 0 0

-1 -1 -1

-2 -2 -2

-3 -3 -3
-3 -2 -1 0 1 2 3 -3 -2 -1 0 1 2 3 -3 -2 -1 0 1 2 3
x/M x/M x/M
a=0.4,q=0.4,n=0.1, =0.1, =-1/3, =1/6 a=0.4,q=0.4,n=0.6, =0.1, =-1/3, =1/6 a=0.4,q=0.4,n=0.9, =0.1, =-1/3, =1/6
3 3 3

2 2 2

1 1 1
z/M

z/M

z/M

0 0 0

-1 -1 -1

-2 -2 -2

-3 -3 -3
-3 -2 -1 0 1 2 3 -3 -2 -1 0 1 2 3 -3 -2 -1 0 1 2 3
x/M x/M x/M
a=0.4,q=0.4,n=0.4, =0 a=0.4,q=0.4,n=0.4, =0.01, =-1/3, =1/6 a=0.4,q=0.4,n=0.4, =0.1, =-1/3, =1/6
3 3 3

2 2 2

1 1 1
z/M

z/M

z/M

0 0 0

-1 -1 -1

-2 -2 -2

-3 -3 -3
-3 -2 -1 0 1 2 3 -3 -2 -1 0 1 2 3 -3 -2 -1 0 1 2 3
x/M x/M x/M
a=0.4,q=0.4,n=0.4, =0.1, =-2/3, =-2 a=0.4,q=0.4,n=0.4, =0.1, =-2/3, =0 a=0.4,q=0.4,n=0.4, =0.1, =-2/3, =2
3
3 3
2 2 2
1 1 1
z/M

z/M

z/M

0 0 0

-1 -1 -1

-2 -2 -2

-3 -3
-3
-3 -2 -1 0 1 2 3 -3 -2 -1 0 1 2 3 -3 -2 -1 0 1 2 3
x/M x/M x/M
a=0.6,q=0.6,n=0.6, =0.1, =1/3, =-1/2 a=0.6,q=0.6,n=0.6, =0.1, =-1/3, =-1/2 a=0.6,q=0.6,n=0.6, =0.1, =-2/3, =-1/2

2 2 2

1 1 1
z/M

z/M

z/M

0 0 0

-1 -1 -1
14
-2 -2 -2

-2 -1 0 1 2 -2 -1 0 1 2 -2 -1 0 1 2
x/M x/M x/M

Figure 2: The ergoregion of Kerr-Newman-NUT-Kiselev black holes in Rastall gravity with the variation of
a, q, n, α, ω, κλ. The chosen value of each parameter is given on the picture.
quintessence in Rastall gravity, the static limit surface comes from the equation that is
given by
g00 = r2 − 2M r + a2 cos2 θ + q 2 − n2 − αrυ = 0. (42)
These surfaces meet at the poles. The region between them gives rise to ergoregion or
ergosphere admitting negative energy orbits where the asymptotic time-like Killing vector
sls
becomes space-like. The region of ergosphere is between r+ < r < r+ . It turns out that the
shape of the ergosphere also depends on such parameters as the horizon. As an example for
the perfect fluid dark matter dominated solution with κλ = −1/2, we have
r
sls
 α α 2
r+ = M + ± M+ + n2 − a2 cos2 θ − q 2 . (43)
2 2
The plot of ergosphere is illustrated in Fig. 2 with the variation of all parameters. Note
that when the ergoshepere becomes smaller, it means the rotational energy of the black
holes is decreasing and vice versa. From Fig. 2, we can see that the bigger a, the wider
ergosphere. But the behavior of the ergosphere due to the presence of q is different because
the bigger charge, the smaller its size. Then if we enlarge the value of NUT charge n and
quintessential intensity α, it makes the ergosphere larger. Remember that g and n are related
to the magnetic monopole. However, from the ergosphere size, both parameters possesses
opposite effect. Then when we vary κλ, it is obvious that in Einstein theory the ergosphere
is bigger than in Rastall theory because the denominator on the power of r that contains
Rastall coupling constant is equal to 1. However, for general κλ 6= 0, the increasing of its
value causes the increasing of ergosphere’s size. Then the bottom pictures show the effect
of ω for which the decreasing of it makes the ergosphere bigger (clearer when ω = −1/3 to
ω = −2/3).

4. Some Thermodynamic Properties


In this section, some thermodynamic quantities are investigated. For the black hole
solution that is indicated by Eq. (33), the number of horizon depends on ω and κλ as the
power of r. To study the thermodynamics of this solution, we need to know the event horizon
but it is very complicated to find it for non-arbitrary ω, κλ. Hence, here we will define the
thermodynamics in term of event horizon r+ . Firstly, the mass is determined from Eq. (40)
as
1 2 υ
+ a2 + q 2 − n 2 .

M= r+ − αr+ (44)
2r+
To calculate the Hawking temperature, the tunneling method can be considered as in
[114, 115, 116]. In this fashion, the metric is allowed to be diagonal so, we put dθ = dϕ = 0
and take the angular coordinate to be θ = 0. Hence, the metric and Hawking temperature
is expressed respectively as

2 2 1 2 ∂r f (r)
ds = −f (r)dt + dr , TH = . (45)
f (r) 4π r=r+
15
To make use of the tunneling method, we apply the above assumption such that the metric
(33) becomes
∆ r2 + (a + n)2 2
ds2 = − dt2
+ dr . (46)
r2 + (a + n)2 ∆
Then for metric (46), the Hawking temperature is computed and straightforwardly yields
υ−1
2(r+ − M ) − αυr+
TH = 2
, (47)
4π[r+ + (a + n)2 ]
where the mass is given in Eq. (44). The angular momentum of solution (33) is then

g03 a
ΩH = − = 2 . (48)
g33 r=r+ [r+ + (a + n)2 ]

The Coulomb electromagnetic potential is given by [80]


er+ + g(a + n) g(a + n)
ΦH = −K µ Aµ r=r+ = Φe + Φg = 2 2
− 2 , (49)
r+ + (a + n) r+ + (a + n)2
where K = ∂t + ΩH ∂ϕ . The resulting Hawking temperature is general for every value of
ω, κλ. So for specific matter domination as in [28], we just have to put the value of the
equation of state. For α = 0, the temperature reduces to Kerr-Newman-NUT black hole
solution [80]. The form of angular momentum and Coulomb potential is also similar to
Kerr-Newman-NUT one but it differs on the definition of r+ . Then it is obvious that for
α, n = 0, it reduces to Kerr-Newman black hole properties. It is worth pointing out that the
existence of the quintessence causes another thermodynamic variable to appear [117, 118].
Hence, we can derive the new thermodynamic variable as

∂M 1 υ−1
ΘH = = − r+ , (50)
∂α r=r+
2
when S, J, Qe , Qg are held constant. We need to add this generalized force in order to satisfy
the first law of black hole’s thermodynamics [117, 118].
The area of the black hole is given by

Z
2
ABH = g22 g33 dθdϕ = 4π[r+ + (a + n)2 ]. (51)

Using the area law of the black hole, the Bekenstein-Hawking entropy can be found as
ABH 2
SBH = = π[r+ + (a + n)2 ]. (52)
4
For constant value of J, Qe , Qg , we can also obtain the heat capacity using the relation that
is given by [48, 49, 119]

dM dM/dr+
C= = , (53)
dTH r=r+ dTH /dr+ r=r+
16
Hence from Eq. (53), we find that
2
2π[r+ + (a + n)2 ]2 2
r+ − α(υ − 1)2 r+υ
+ a2 + q 2 − n2
C= 2 2
, where δ = 2 υ
. (54)
[r+ + (a + n)2 ]δ − 2r+ r+ − α(υ − 1)r+ − a2 − q 2 + n 2

The dependence of α, ω and κλ on area and entropy is implicitly contained on r+ while


on heat capacity, it is obviously seen in Eq. (54). The thermodynamic stability of any
system is related to the sign of the heat capacity [109]. For black holes, when the heat
capacity is positive (C > 0), it will be stable. Whereas when the heat capacity is negative
(C < 0), the black holes are said to be unstable. On the Kerr-Newman-NUT-Kiselev black
hole solution in Rastall gravity, the heat capacity is dependent on many parameters. So,
to find the thermodynamic stability, those parameters need to be determined. Note that
numerator of the heat capacity is always positive so the stability is determined by the sign of
the denominator. It can be obviously seen that the Hawking temperature, entropy, angular
momentum, Coulomb potential, and heat capacity are explicitly dependent of NUT charge.
We give examples of the entropy for some arbitrary values of ω, κλ. The values are chosen
from the Table 1 where the event horizons are provided. Then the extremal Bekenstein-
Hawking entropies are given by
 
(1) α 2 2
SBH = π M + + (a + n) , (55)
2
" 2 #
(2) M
SBH = π + (a + n)2 , (56)
1−α

for the black hole solution with perfect fluid dark matter domination ω, κλ = −1/3, −1/2
and dust domination ω, κλ = 0, 1/6, respectively. The investigation of every specific value
of ω, κλ is studied in [28] for non-rotating case related to the violation of the strong energy
condition. We can note that in Eq. (56), the quintessential intensity is restricted on α = 1 to
gain finite entropy for the given value of ω, κλ. We will also have the same value of entropy
for ω, κλ = 0, 0 and ω, κλ = −1/3, 0 for dust domination and perfect fluid dark matter
domination in Einstein theory, respectively. When the quintessence is absent, it reduces to
the extremal Kerr-Newman-NUT black hole’s entropy.

5. Central charge from CFT


The Kerr/CFT correspondence, firstly proposed in [66], has been successful to compute
the microscopic origin of extremal black hole entropy. To find the microscopic entropy, the
famous Cardy formula is applied that is given by

π2
SCF T = (cL TL + cR TR ), (57)
3
where cL , cR are the left-moving and right-moving central charges while TL , TR are the left-
moving and right-moving conformal temperatures. Using this correspondence, it is proved
17
that entropy of extremal black holes from CFT (57) agrees with the Bekenstein-Hawking one
[67]-[80]. In the previous section, the Bekenstein-Hawking entropy for Kerr-Newman-NUT-
Kiselev black hole in Rastall gravity has been calculated. However, the microscopic origin of
this entropy (52) is not calculated yet using the Kerr/CFT correspondence. One reason is
because the canonically conserved charge for Rastall theory of gravity is not calculated yet.
But, for the other modified gravity theory can be found in [101]. The asymptotic charge that
is related to the asymptotic symmetry is used to find the central charge. The central charge
shows up from the central extension that comes from the Dirac bracket of asymptotic charges
[99, 100]. In this section, we are going to find the central charge considering the Kerr/CFT
correspondence that relates entropy from CFT and the Bekenstein-Hawking entropy. We
argue that this correspondence may be implemented because when κλ vanishes, it becomes
the usual Einstein-Maxwell system with the scalar field. In [102, 103, 104, 105], this relation
is also performed for deformed black hole solutions because the central charge is not cal-
culated yet. For the solutions consisting of the electromagnetic and scalar field besides the
gravity, this correspondence is able to prove the relation between the Bekenstein-Hawking
entropy and CFTs entropy such as for Kerr-Sen black hole solution [68] and black hole with
dilaton-axion field [70]. Generally, this relation is valid not only for extremal black holes,
but also for non-extremal black holes [81]-[97].
Before we derive the central charge, it is common to show the isometry of the near-horizon
extremal black hole solutions. In addition, the isometry should be U (1)×SL(2, R) to employ
the asymptotic symmetry group. So we will show that the near-horizon extremal Kerr-
Newman-NUT-Kiselev black hole solutions in Rastall gravity have that isometry. Firstly,
we need to define the following coordinates transformation based on [66, 79] which is given
by
r0 r0
r = r+ + r0 y, t = τ, ϕ = φ + ΩH τ, (58)
 
where r02 = r+2
+ (a + n)2 . After taking  → 0, the near-horizon extremal metric of (33) reads
2
ρ2+ dy 2 r04 sin2 θ
  
2 2 2 2 2ar+
ds = −y dτ + 2 + V dθ + dφ + ydτ , (59)
V y ρ2+ V r02
where ρ2+ = r+
2
+ (n + a cos θ)2 . Straightforwardly, we could perceive that the near-horizon
metric (59) has AdS2 × S 2 topology. It possesses an isometry U (1) × SL(2, R) which is
generated by the following vector fields
τ2
 
1 2ar+
ζ0 = ∂φ , X1 = ∂τ , X2 = τ ∂τ − y∂y , X3 = + ∂τ − τ y∂y − ∂φ . (60)
2y 2 2 V r02 y
That isometry exhibits that Kerr/CFT should be valid for our extremal solutions. Fur-
thermore, to obtain the central charge, we need to compute the the conformal temperatures.
In computing this temperatures, we apply the fashion shown in [78, 79, 80] for which it is
given by

TH r0 ∂TH /∂r+
TR = , TL = − , (61)
 r=r+ ∂ΩH /∂r+ r=r+
18
where r0 is a constant to factor out the overall scale of the near-horizon geometry and 
is an infinitesimal constant for extremal black holes. It is clear that for extremal black
holes, the right-moving temperature will vanish because it is proportional to the Hawking
temperature. Hence, the remaining temperature that is required to be calculated is the
left-moving temperature based on Eq. (57). Using the formula of left-moving temperature
in Eq. (61), we finally derive
υ−2 2
[2 − αυ(υ − 1)r+ ][r+ + (a + n)2 ]
TL = . (62)
8πar+
The conformal temperatures come from the construction of the Frolov-Thorne vacuum for
generic rotating black holes. When α = 0, it reduces to the left-moving temperature of Kerr-
Newman-NUT black hole. After getting the left-moving temperature, the Cardy formula
(57) is applied. Since we consider SCF T = SBH , finally the left-moving central charge is
given as follows
24ar+
cL = υ−2 . (63)
2 − αυ(υ − 1)r+

In general, the central charge and left-moving temperature show up as a function of the
quintessential intensity, equation of state, and Rastall parameter besides the other parame-
ters. It is interesting because any specific of matter domination as the quintessential model,
along with the Rastall coupling constant affects the value of these quantities. It is obvious
υ−2
that αυ(υ −1)r+ may not be equal to 2 to produce finite left-moving central charge. So the
value of ω, κλ are required to be constrained but we will not explain in more detail within
this paper.
As the entropy, we give examples of the central charge for some arbitrary values of ω, κλ
given in the Table 1. We find the following central charges
(1)
 α
cL = 12a M + , (64)
2
(2) 12aM
cL = , (65)
(1 − α)2

for the black hole solution with perfect fluid dark matter domination ω, κλ = −1/3, −1/2
and dust domination ω, κλ = 0, 1/6, respectively. We will also have the same value of central
charge for ω, κλ = 0, 0 and ω, κλ = −1/3, 0, respectively. When the quintessential intensity
vanishes, the central charge coincides to the central charge of extremal Kerr-Newman-NUT
black hole [80]. In addition, when q vanishes, it reduces to the result obtained in [73, 120].

6. Summary
By performing the Demiański-Newman-Janis algorithm, we obtained Kerr-Newman-
NUT-Kiselev black hole solutions in Rastall theory of gravity which is the extension of
Kerr-Newman-Kiselev black holes in Rastall theory. This new solution contains magnetic
19
charge and NUT parameter which are related to the magnetic monopole. It is believed that
in the early universe, the magnetic monopole existed. The horizon analysis and ergosphere
were also studied to see the possibilities of having several horizons and the dependence of
each parameter. For some arbitrary values of ω, κλ, this solution coincides to possess inner
and outer horizons. The size of the ergosphere is dependent of each parameter related to
the rotational energy of the black hole. We also could map several values of ω in Rastall
gravity to the other values of ω in Einstein general relativity. But we did not say that these
two theories are equivalent because Rastall theory seems more Machian and needs to face
the challenges of the observational evidence. From the argument that we pointed out, we
believed that Rastall gravity is not equivalent to Einstein gravity.
Several thermodynamic properties were investigated. We defined the mass from the
horizon equation. It is useful to define the Hawking temperature that we computed by
applying tunneling method. The angular momentum, Coulomb electromagnetic potential
and heat capacity were also calculated. Because of the quintessential field, we needed to
add the generalized force in order to satisfy the first law of black hole’s thermodynamics.
Then the Bekenstein-Hawking entropy was also derived which is equal with the area law
of the black hole. It could be seen that the dependence of the quintessential intensity α
implicitly lies on the horizon r+ but for the temperature, its dependence is obvious. The
Hawking temperature, entropy, angular momentum, Coulomb potential, and heat capacity
are explicitly dependent of NUT charge also, besides the spin. It will be interesting to study
the thermodynamic stability of this solution from the heat capacity, but we did not explain
in this article.
Instead of using the asymptotic symmetry group in Kerr/CFT correspondence, we cal-
culated the central charge of the extremal black hole solutions by equating the entropy from
the CFT with the Bekenstein-Hawking one which was already proved. Before finding the
central charge, we showed the isometry of the near-horizon extremal form of our solution.
Furthermore, the conformal temperatures were computed first where the temperatures that
we needed are right- and left-moving temperatures. But from the fact that the right-moving
part was proportional to the Hawking temperature, it remained zero in extremal case. Fi-
nally, the remaining central charge was obtained. The left-moving central charge depends
υ−2
on the quintessential intensity explicitly. Generally, αυ(υ − 1)r+ may not be equal to 2
to produce finite left-moving central charge. The central charges for ω, κλ = −1/3, −1/2
and ω, κλ = 0, 1/6 were also shown as the examples. In the future, we want to study this
solution using the famous Kerr/CFT correspondence to confirm the result that we obtained
in this article.

Acknowledgments
We gratefully acknowledge the support from Ministry of Research, Technology, and
Higher Education of the Republic of Indonesia and PMDSU Research Grant. We also want
to express many thanks all members of Theoretical Physics Laboratory, Institut Teknologi
Bandung for the valuable support.

20
Appendix A. Einstein tensor of Kerr-Newman-NUT black holes with quintessence
in Rastall gravity
We derive the Einstein tensor of Kerr-Newman-NUT black holes with the quintessence
in Rastall gravity using Mathematica package RGTC as
2 (r4 − 2r3 + a2 r2 − a4 x2 (1 − x2 )) W 0 ra2 sin2 θW 00
Gtt = −
ρ2 ρ4
2n [n {r2 + a2 (1 − x2 )} + 2a3 x(1 − x2 )]
− , (A.1)
ρ6

4nr2 (1 − x)(a2 − ∆)W 0 a(1 − x2 )(r2 + n2 )2 X


Gtφ = +
ρ2 2r2 ρ6
2ar (1 − x )(2r(M − r) − (q − αrυ ))W 0
2 2 2
+
ρ6
an(1 − x) [4r (x − 2)W 0 + (1 − x)(n2 + r2 )X]
4
+
r2 ρ6
a (1 − x ) [(1 + x )X − 8r2 W 0 ] a3 n2 (1 − x2 )(1 + 4x + x2 )X
3 2 2
+ +
2ρ6 2r2 ρ6
(A.2)
2r2 W 0
Grr = − , (A.3)
ρ2 ∆

2a2 x2 W 0 2(n2 + 2anx)W 0 rW 00


Gθθ = − − − , (A.4)
(1 − x2 )ρ2 (1 − x2 )ρ2 (1 − x2 )

a(1 − x2 )(r2 + a2 )(a2 + (2r2 + a2 )(2x2 − 1)W 0


Gφφ = −
ρ2
2r3 a(1 − x2 )2 W W 0 r2 a2 (1 − x2 )2 (r2 + a2 − n2 − ∆)W
− +
ρ2 ρ6
a2 (r2 + a2 )(1 − x2 ) [2a2 x2 + r2 (4x2 − 2)] W 0
+
ρ6
∆(x − 1)2 (2n + a + ax)2 (2r2 W 0 + ρ2 X)Y a3 x(2n + 2nx + ax)X
+ +
2r2 ρ8 2r2 ρ8
∆(x − 1)2 (2n + a + ax)2 (r2 + n2 )(2r5 W 00 + n2 X)

2r2 ρ8
a[4r2 n(x − 2)W 0 + (1 + x)X] a2 [r2 (1 + x2 ) + n2 (1 + 4x + x2 )]
+ +
ρ8 2r2 ρ8
(r2 + n2 )(1 − x2 )(2r5 W 00 + n2 x)Y 2 2r2 a(x2 − 2)W 0
+ + . (A.5)
2r2 ρ8 ρ8
21
2
where x = cos θ, W = M − 2rq 2 + α2 rυ−1 , q 2 = e2 + g 2 , X = αυ(υ − 1)rυ and Y = r2 + (a + n)2 .
The Ricci scalar is
X
R = 2 2. (A.6)
r ρ
Note that the total gravitational tensor is

Gµν + κλgµν R. (A.7)

Appendix B. Energy-momentum components of Kerr-Newman-NUT black holes


with quintessence in Rastall gravity
Using the gravitational tensor, we may construct the energy-momentum tensor of the
Kerr-Newman-NUT black holes solution with the quintessence in Rastall gravity. We can
get the following equation
1
Tνα = (Gαν + κλδνα R) , (B.1)
κ
by rising one index from Eq. (35). Finally, we can find

2 (r4 + r2 (2a2 + n2 ) + 4an − r2 a(a + 2n)x) a(x − 1)(a + 2n + ax)X λX


E = Ttt = − 6
− 2 4
+ 2 2,
κρ 2r κρ r ρ
(B.2)
2 0
2r W λX
Pr = Trr = − + 2 2, (B.3)
κρ4 r ρ
2r2 W 0 X λX
Pθ = Tθθ = 4
− 2 2 + 2 2, (B.4)
κρ 2r κρ r ρ

r3 (r2 + n2 )W 00 2ar2 [a(x2 − 2) + 2n(x − 2)]W 0 a3 x(2n + ax + 2nx)X


Pφ = Tφφ = − − −
κρ6 κρ6 2r2 κρ6
(r2 + n2 )[n2 + 2an(1 + x)]X a2 [r2 (1 + x2 ) + n2 (1 + 4x + x2 )] λX
− 2 6
− 2 6
+ 2 2 (, B.5)
2r κρ 2r κρ r ρ

4ar2 W 0 aX
S = Ttφ = − 6
+ 2 4. (B.6)
κρ 2r κρ

References
[1] Rastall P 1971 Phys. Rev. D 6 3357
[2] Rastall P 1976 Can. J. Phys. 54 66
[3] Fisher S B and Carlson E D 2019 Phys. Rev. D 100 064059
[4] De Moraes W A G and Santos A F 2019 Gen. Relativ. Grav. 51 167
[5] Majernı́k V and Richterek L Rastalls gravity equations and Machs Principle arXiv:gr-qc/0610070
[6] Campos J P, Fabris J C, Perez R, Piattella O F and Velten H 2013 Eur. Phys. J. C 73 2357
[7] Fabris J C, Daouda M H and Piattella O F 2012 Phys. Lett. B 711 232
[8] Moradpour H 2016 Phys. Lett. B 757 187
[9] Al-Rawaf A S and Taha M O 1996 Gen. Relativ. Gravit. 28 935
22
[10] Batista C E M, Daouda M H, Fabris J C, Piattella O F and Rodrigues D C 2012 Phys. Rev. D 85
084008
[11] Caramês T R P, Daouda M H, Fabris J C, Oliveira A M, Piattella O F and Strokov V 2014 Eur. Phys.
J. C 74 3145
[12] Salako I G, Houndjo M J S and Jawad A 2016 Int. J. Mod. Phys. D 25 1650076
[13] Smalley L L 1974 Phys. Rev. D 9 1635
[14] Smalley L L 1975 Phys. Rev. D 12 376
[15] Wolf C 1986 Phys. Scr. 34 193
[16] Moradpour H, Heydarzade Y, Darabi F, Salako I G 2017 Eur. Phys. J. C 77 259
[17] Hansraj S, Banerjee A and Channuie P 2019 Ann. Phys. 400 320
[18] Abbas G and Shahzad M R 2018 Eur. Phys. J. A 54 211
[19] Abbas G and Shahzad M R 2019 Astrophys. Space Sci. 364 50
[20] Oliveira A M, Velten H E S, Fabris J C and Casarini L 2015 Phys. Rev. D 92 044020
[21] Moradpour H and Sadeghnezhad N 2017 Can. J. Phys. 95 1257
[22] Moradpour H and Salako I G 2016 Adv. High Energy Phys. 2016 3492796
[23] Oliveira A M, Velten H E S and Fabris J C 2016 Phys. Rev. D 93 124020
[24] Bronnikov K A, Fabris J C, Piattella O F and Santos E C 2016 Gen. Relativ. Gravit. 48 162
[25] Prihadi H L, Sakti M F A R, Hikmawan G and Zen F P Dynamics of Charged and Rotating NUT
Black Holes in Rastall Gravity arXiv:1908.09629 [gr-qc]
[26] Ma M-S and Zhao R 2017 Eur. Phys. J. C 77 629
[27] Spallucci E and Smailagic A 2018 Int. J. of Mod. Phys. D 27 1850003
[28] Heydarzade Y and Darabi F 2017 Phys. Lett. B 771 365-373
[29] Heydarzade Y, Moradpour H and Darabi F 2017 Can. J. Phys. 95 1253-1256
[30] Kumar R and Ghosh S G 2018 Eur. Phys. J. C 78 750
[31] Xu Z, Hou X, Gong X and Wang J 2018 Eur. Phys. J. C 78 513
[32] Kumar R, Singh B P, Ali M S and Ghosh S G Rotating black hole shadow in Rastall theory
arXiv:1712.09793 [gr-qc]
[33] Lobo I P, Moradpour H, Graca J P M and Salako I G 2018 Int. J. Mod. Phys. D 27 1850069
[34] Visser M 2018 Phys. Lett. B 782 8386
[35] Darabi F, Moradpour H, Licata I, Heydarzade Y and Corda C 2018 Eur. Phys. J. C 78 25
[36] Copeland E J, Sami M and Tsujikawa S 2006 Int. J. Mod. Phys. D 15 1753
[37] Arianto, Zen F P, Gunara B E, Triyanta and Supardi 2007 J. High Energy Phys. 09 048
[38] Arianto, Zen F P, Triyanta and Gunara B E 2008 Phys. Rev. D 77 123517
[39] Zen F P, Arianto, Gunara B E, Triyanta and Purwanto A 2009 Eur. Phys. J. C 63 477
[40] Arianto, Zen F P and Gunara B E 2010 Gen. Relativ. Gravit. 42 909
[41] Feranie S, Arianto and Zen F P 2010 Phys. Rev. D 81 084058
[42] Arianto, Zen F P, Feranie S, Widyatmika I P and Gunara B E, Phys. Rev. D 84 044008
[43] Suroso A and Zen F P 2013 Gen. Relativ. Gravit. 45 799
[44] Suroso A and Zen F P 2015 Adv. Stud. Theor. Phys. 9 423
[45] Hikmawan G, Soda J, Suroso A and Zen F P 2016 Phys. Rev. D 93 068301
[46] Kiselev V V 2003 Class. Quantum Grav. 20 1187
[47] Visser M The Kiselev black hole is neither perfect fluid, nor is it quintessence arXiv:1908.11058 [gr-qc]
[48] Ghosh S G, Amir M and Maharaj S D 2017 Eur. Phys. J. C 77 530.
[49] Ghosh S G, Maharaj S D, Baboolal D and Lee T 2018 Eur. Phys. J. C 78 90
[50] Buchdahl H A 1970 MNRAS 150 1-8
[51] Hořava P 2009 Phys. Rev. D 79 084008
[52] Engels T, Kolomenskiy D, Schneider K, Lehmann F-O and Sesterhenn J 2016 Phys. Rev. Lett. 116
028103
[53] Maluf R V, Silva J E G and Almeida C A S 2015 Phys. Lett. B 749 304-308
[54] Kanzi S and Sakalli İ 2019 Nuc. Phys. B 946 114703
[55] Jesus W D R and Santos A F 2019 Mod. Phys. Lett. A 34 1950171

23
[56] Oliveira R, Dantas D M, Santos V and Almeida C A S 2019 Class. Quantum Grav. 36 105013
[57] Övgün A, Jusufi K and Sakalli İ 2019 Phys. Rev. D 99 024042
[58] Capelo D and Páramos J 2015 Phys. Rev. D 91 104007
[59] Guiomar G and Páramos J Phys. Rev. D 90 082002
[60] Casana R, Cavalcante A, Poulis F P and Santos E B 2018 Phys. Rev. D 97 104001
[61] Lin K and Qian W-L 2019 Chinese Phys. C 43 083106
[62] Al-Badawi A and Halilsoy M 2006 Gen. Relativ. Gravit. 38 1729
[63] Lynden-Bell D and Nouri-Zonoz M 1998 Rev. Mod. Phys. 70 427
[64] Erbin H 2015 Gen. Relativ. Gravit. 47 19
[65] Erbin H 2016 Gen. Relativ. Gravit. 48 56
[66] Guica M, Hartman T, Song W and Strominger A 2008 Phys. Rev. D 80 124008
[67] Hartman T, Murata K, Nishioka T and Strominger A 2009 J. High Energy Phys. 04 019
[68] Ghezelbash A M 2009 J. High Energy Phys. 08 045
[69] Lü H, Mei J and Pope C N 2009 J. High Energy Phys. 04 054
[70] Li R and Ren JR 2010 J. High Energy Phys. 09 039 (2010)
[71] Anninos D and Hartman T 2010 J. High Energy Phys. 03 096
[72] Ghodsi A and Garousi M R 2010 Phys. Lett. B 687 79 2010
[73] Ghezelbash A M 2012 Mod. Phys. Lett. A 27 1250046
[74] Astorino M 2015 J. High Energy Phys. 10 016
[75] Astorino M 2015 Phys. Lett. B 751 96
[76] Siahaan H M 2016 Class. Quantum Grav. 33 155013
[77] Astorino M 2016 Phys. Lett. B 760 393
[78] Sinamuli M and Mann R B 2016 J. High Energy Phys. 08 148
[79] Compére G 2017 Living Rev. Relativ. 20
[80] Sakti M F A R, Suroso A and Zen F P 2018 Int. J. Mod. Phys. D 27 1850109
[81] Castro A, Maloney A and Strominger A 2010 Phys. Rev. D 82 024008
[82] Chen B, Long J and Zhang J-J 2010 Phys. Rev. D 82 104017
[83] Chen C-M and Sun J-R 2010 J. High Energy Phys. 08 034
[84] Chen D, Wang H, Wu H and Yang H 2010 J. High Energy Phys. 11 002
[85] Wang Y-Q and Liu Y-X 2010 J. High Energy Phys. 08 087
[86] Chen B and Long J 2010 J. High Energy Phys. 08 065
[87] Setare M R and Kamali V 2010 J. High Energy Phys. 10 074
[88] Ghezelbash A M , Kamali V and Setare M R 2010 Phys. Rev. D 82 124051
[89] Chen C-M, Huang Y-M, Sun J-R, Wu M-F and Zou S-J 2010 Phys. Rev. D 82 066004
[90] Chen B, Chen C-M and Ning B 2011 Nuc. Phys. B 853 196209
[91] Chen B, Ghezelbash A M, Kamali V and Setare M R 2011 Nuc. Phys. B 848 108120
[92] Chen B and Zhang J-J 2011 J. High Energy Phys. 08 114
[93] Huang Y-C and Yuan F-F 2011 J. High Energy Phys. 03 029
[94] Shao K-N and Zhang Z 2011 Phys. Rev. D 83 106008
[95] Chen D, Wang P, Wu H and Yang H 2011 Gen. Relativ. Gravit. 43 181190
[96] Ghezelbash A M and Siahaan H M 2013 Class. Quantum Grav. 30 135005
[97] Siahaan H M 2018 Class. Quantum Grav. 35 155002
[98] Sakti M F A R, Suroso A and Zen F P 2019 Eur. Phys. J. Plus 134 580
[99] Barnich G and Brandt F 2002 Nucl. Phys. B 633 3
[100] Barnich G and Compere G 2008 J. Math. Phys. 49 042901
[101] Adami H, Setare M R and Sisman T C Conserved Charges in Extended Theories of Gravity
arXiv:1710.07252 [hep-th]
[102] Lowe D A and Skanata A 2012 J. Phys. A: Math. Theor. 45 475401
[103] Ghezelbash A M and Siahaan H M 2013 Class. Quantum Grav. 30 135005
[104] Ghezelbash A M and Siahaan H M 2014 Gen. Relativ. Gravit. 46 1783
[105] Sakti M F A R, Ghezelbash A M, Suroso A and Zen F P 2019 Gen. Relativ. Gravit. 51 151

24
[106] Demianski M 1972 Phys. Lett. A 42 157
[107] Giampieri G 1990 Introducing Angular Momentum into A Black Hole using Complex Variables Gravity
Research Foundation http://gravityresearchfoundation.org/pdf/awarded/1990/giampieri
[108] Xu Z and Wang J 2017 Phys. Rev. D 95 064015
[109] Sakti M F A R, Prihadi H L, Suroso A and Zen F P Rotating and Twisting Charged Black Holes with
Cloud of Strings and Quintessence arXiv:1911.07569 [gr-qc]
[110] Toshmatov B, Stuchlı́k Z and Ahmedov B 2017 Eur. Phys. J. Plus 132 98
[111] Podolský J and Kadlecová H 2006 Class. Quantum Grav. 26 105007
[112] Xu Z and Wang J Thermodynamics and Phase Transition in Rotational Kiselev Black Hole
arXiv:1610.00376 [gr-qc]
[113] Ghosh S G 2016 Eur. Phys. J. C 76 222
[114] Parikh M K and Wilczek F 2000 Phys. Rev. Lett. 85 5042
[115] Banerjee R and Majhi B R 2008 J. High Energy Phys. 06 095
[116] Ma Z Z 2008 Phys. Lett. B 666 376
[117] Chen S, Wang B and Su R 2008 Phys. Rev. D 77 124011
[118] Sekiwa Y 2006 Phys. Rev. D 73 084009
[119] Cai R G 2004 Phys. Lett. B 582 237
[120] Sakti M F A R, Suroso A and Zen F P 2019 J. Phys.: Conf. Ser. 1204 012009
[121] Gibbons G W 1977 Phys. Rev. D 15 3530
[122] Lee K, Nair V P and Weinberg E J 1992 Phys.Rev. D 45 2751-2761

25

You might also like