---
title: Kiselev Black Holes Overview
url: https://www.emergentmind.com/topics/kiselev-black-holes
type: topic
---

# Kiselev Black Holes Overview

Searching arXiv for recent and foundational work on Kiselev black holes to ground the article in the current literature.
Kiselev black holes are static, spherically symmetric black-hole spacetimes sourced by an anisotropic fluid that is often described as “quintessence-like,” but whose stress tensor is not that of a perfect fluid. In its simplest asymptotically flat form, the metric can be written as \(ds^2=-f(r)\,dt^2+f(r)^{-1}dr^2+r^2d\Omega_2^2\) with \(f(r)=1-\frac{2M}{r}-\frac{K}{r^{1+3w}}\), while equivalent conventions write \(f(r)=1-\frac{2m}{r}-\sigma r^{-(3w_q+1)}\) [2512.01265, 1612.07279]. The subject has developed into a broad research program spanning geodesics and lensing, thermodynamics, quasinormal spectra, modified gravity, rotating and regular extensions, and analogue-gravity realizations; at the same time, later work has emphasized that the Kiselev source is intrinsically anisotropic and should not be identified with cosmological quintessence in the strict perfect-fluid sense [2512.01265, 1910.08008].

## 1. Definition and matter source

The defining matter sector of a Kiselev black hole is an anisotropic fluid with nonzero components \(T^t_{\ t}=-\rho(r)\), \(T^r_{\ r}=p_r(r)\), and \(T^\theta_{\ \theta}=T^\varphi_{\ \varphi}=p_t(r)\). In the asymptotically flat form discussed in later analyses, the source satisfies
\[
\rho(r)=-p_r(r)=-\frac{3K\,w}{8\pi\,r^{3(1+w)}},\qquad
p_t(r)=-\frac{3K\,w(1+3w)}{8\pi\,r^{3(1+w)}},
\]
with averaged pressure \(\bar p=(p_r+2p_t)/3=w\rho\) [2512.01265]. This is the origin of the persistent but potentially misleading shorthand “quintessence”: the averaged pressure obeys a barotropic relation, but the local fluid is anisotropic.

That anisotropy is not incidental. One later decomposition shows that the Kiselev stress tensor can be mimicked by a perfect fluid plus either an electromagnetic component or a scalar-field component, depending on the sign of \(K\,w(1+w)\). In the one-component model, \(K\,w(1+w)<0\) corresponds to a perfect fluid plus electromagnetic field, \(K\,w(1+w)=0\) reduces to Schwarzschild or Kottler, and \(K\,w(1+w)>0\) corresponds to a perfect fluid plus scalar field; in the generalized \(N\)-component case the same trichotomy is controlled by the sign of \(\rho'(r)\) [1910.08008]. This decomposition quantifies the precise sense in which the Kiselev geometry fails to be a perfect-fluid spacetime.

A further reinterpretation replaces the anisotropic fluid by nonlinear electrodynamics. In Einstein–power-Maxwell theory, the geometry
\[
f(r)=1-\frac{2M}{r}-k\,r^p,\qquad p=-(3w+1),
\]
is an exact solution sourced either by an electric or by a magnetic power-Maxwell field, with different powers \(q\) for the two realizations. In that formulation, the Kiselev amplitude is fixed by the electric or magnetic charge, while the source still satisfies the characteristic anisotropic relation \(p_r=-\rho\neq p_t\) [2206.12876].

## 2. Metric structure, parameter regimes, and horizons

In the standard spherical form, the Kiselev metric interpolates among several familiar geometries through the equation-of-state parameter. The cases \(w=-1\), \(w=-2/3\), \(w=-1/3\), \(w=0\), and \(w=1/3\) respectively generate Schwarzschild–de Sitter/anti–de Sitter, a linear-in-\(r\) quintessential term, a constant shift, a dust-like \(1/r\) correction, and a radiation-like \(1/r^2\) term [2206.12876, 2307.09079]. In the convention used for null-geodesic calculations,
\[
f(r)=1-\frac{2m}{r}-\sigma r^{-(3w_q+1)},
\]
with \(-1<w_q<-1/3\); in the limit \(w_q\to -1\), the metric reduces to a cosmological-constant term with \(\Lambda=3\sigma/2\) in the conventions of that treatment [1612.07279].

The case \(w_q=-2/3\) is especially important because the metric simplifies to
\[
f(r)=1-\frac{2M}{r}-\sigma r.
\]
The horizons are then
\[
r_\pm=\frac{1\pm\sqrt{1-8M\sigma}}{2\sigma}.
\]
This yields a standard classification into nonextreme black holes for \(0<\sigma<1/(8M)\), an extreme Kiselev black hole at \(\sigma=1/(8M)\) with degenerate horizon \(r_H^e=4M\), and a naked singularity for \(\sigma>1/(8M)\) [1502.01676].

A distinct regime arises in the “reduced Kiselev black hole,” obtained by setting the Schwarzschild mass term to zero and working in \(-1/3<w<0\). Writing \(\alpha=3w+1\in(0,1)\), the metric becomes
\[
f(r)=1-\left(\frac{K}{r}\right)^\alpha,
\]
with a single Killing horizon at \(r=K\), asymptotic flatness, and a spacelike curvature singularity at \(r=0\). Its causal structure is the same as Schwarzschild, but the geometry is supported entirely by the anisotropic fluid sector rather than by an explicit mass term [2307.09079].

The same formalism also admits regularizations. One construction matches a Kiselev exterior to a de Sitter interior across a timelike dust shell, producing a globally regular black hole in which the \(r=0\) singularity is replaced by a de Sitter core. Another allows the equation-of-state parameter to vary radially, \(w\to w(r)\), and derives regular solutions with finite curvature invariants at the origin and de Sitter-like cores [2009.06990, 2411.18804].

## 3. Null geodesics, photon spheres, and lensing

For static spherical Kiselev geometries, null geodesics are conveniently written in terms of \(u=1/r\). In the notation of the homotopy-perturbation treatment, the general geodesic equation becomes
\[
\frac{d^2u}{d\phi^2}+u=3mu^2+\frac{3\sigma(w_q+1)}{2}u^{3w_q+2},
\]
with the quintessence contribution disappearing from this equation in the limit \(w_q\to -1\) because of the prefactor \(w_q+1\) [1612.07279]. For \(w_q=-1/3\), the added term is proportional to \(u\); for \(w_q=-2/3\), it is constant, and this difference drives qualitatively different bending-angle corrections.

The homotopy-perturbation method yields analytic weak-field series for both the orbit and the deflection angle. At lowest order, the Kiselev corrections are
\[
\beta= \frac{4m}{l}+\frac{\pi\sigma}{2}\qquad (w_q=-1/3),
\]
and
\[
\beta=\frac{4m}{l}+\sigma l\qquad (w_q=-2/3),
\]
with higher-order terms involving \(m^2/l^2\), \(\sigma^2\), and mixed \(m\sigma\) contributions. The appearance of the \(\sigma l\) term for \(w_q=-2/3\) is the distinctive feature of that case. The same work also derives a homotopy-based root-finding series for the deflection angle that does not require assuming \(\beta\ll1\) and is extendable to arbitrary order in the weak-field regime [1612.07279].

In the strong-lensing analysis of the \(w_q=-2/3\) metric, the photon effective potential is
\[
V_{\rm eff}=\frac{L^2}{r^2}\left(1-\frac{2M}{r}-\sigma r\right).
\]
Circular null orbits satisfy \(dV_{\rm eff}/dr=0\), giving
\[
r_{c\pm}=\frac{1\pm\sqrt{1-6M\sigma}}{\sigma},
\]
with the unstable photon sphere identified as \(r_{\rm ps}=r_{c-}\). The critical impact parameter is
\[
b_{\rm sc}=\sqrt{\frac{r_{\rm ps}^3}{r_{\rm ps}-2M-\sigma r_{\rm ps}^2}}.
\]
The exact bending angle can then be written in terms of incomplete and complete elliptic integrals, and the comparative ordering of deflections is
\[
\text{naked singularity}>\text{extreme Kiselev black hole}>\text{nonextreme Kiselev black hole}>\text{Schwarzschild black hole}
\]
for the parameter ranges studied [1502.01676].

These geodesic structures persist in rotating and modified-gravity extensions, where they reappear in shadow radii, relativistic images, and photon-sphere topology. A plausible implication is that Kiselev deformations are best viewed not as a single universal lensing correction, but as a family of distinct optical signatures controlled by the exponent \(3w+1\) and by the sign convention chosen for the surrounding medium.

## 4. Thermodynamics, evaporation, and wave dynamics

In the standard Kiselev metric written as
\[
f(r)=1-\frac{2M}{r}-\alpha r^{-(3w_q+1)},
\]
the horizon equation is \(f(r_h)=0\), and the Hawking temperature is
\[
T_H=\frac{1}{4\pi r_h}\left[1+3\alpha w_q\,r_h^{-(3w_q+1)}\right].
\]
Since \(w_q<0\) in the quintessential range, the second term is negative, so the surrounding medium lowers the temperature relative to Schwarzschild at fixed \(r_h\) [2606.19110]. In the evaporation study based on Page’s geometrical-optics approximation, decreasing \(w_q\) lowers the non-final-stage temperature and markedly prolongs the evaporation lifetime. For \(\alpha=10^{-40}\) and \(M_i=5.4\times10^{14}\,\mathrm{g}\), the reported lifetimes are \(T(-0.99)\approx20.77\,\mathrm{Gyr}\), \(T(-0.98)\approx14.79\,\mathrm{Gyr}\), while \(T(w_q\ge -0.90)\approx13.45\,\mathrm{Gyr}\), essentially Schwarzschild-like [2606.19110].

Rainbow-gravity deformations preserve the Kiselev horizon equation but rescale the temperature by the rainbow function. In that framework, no generic remnant appears: for \(n=0,1,2\) the evaporation does not halt, and only a peculiar \(n=4\), \(\eta<0\) configuration with a cosmological horizon yields a \(T_H\to0\) endpoint at the degenerate horizon [2106.03672]. In Einstein–power-Maxwell realizations, the event-horizon temperature is
\[
T_b=\frac{1-(p+1)k r_b^{p}}{4\pi r_b},
\]
and the heat capacity at fixed charge is negative throughout the physical range, but it exhibits a Schottky peak characteristic of multi-horizon systems [2206.12876].

Wave dynamics around Kiselev backgrounds has likewise become a substantial subfield. In quantum fluctuation modified gravity, the metric is written as
\[
Q(r)=1-\frac{2M}{r}+P\,r^{\alpha(\mho,\bar{\bar\omega})},
\]
and scalar and vector perturbations obey Schrödinger-type equations with effective potentials
\[
\hat V_s(r)=Q(r)\left[\frac{l(l+1)}{r^2}+\frac{1}{r}\frac{dQ}{dr}\right],\qquad
\hat V_v(r)=Q(r)\frac{l(l+1)}{r^2}.
\]
Using 6th-order WKB, the real part of the quasinormal frequency decreases as the quantum-fluctuation parameter \(\mho\) rises, while the greybody factors are significantly affected by \(\mho\), \(\bar{\bar\omega}\), and the multipole number \(l\) [2510.02409].

## 5. Modified-gravity, rotating, regular, topological, and analogue extensions

Kiselev black holes have been generalized extensively in modified gravity. In \(f(R,T)=R+2\varkappa T\) gravity, the static spherical solution takes the form
\[
B(r)=1-\frac{2M}{r}+K\,r^{-\gamma(w,\varkappa)},\qquad
\gamma(w,\varkappa)=\frac{8(3\pi w+w\varkappa+\pi)}{3\varkappa+8\pi-w\varkappa}.
\]
Only the \(w=-1\) case reproduces the GR exponent exactly; the other sectors acquire \(\varkappa\)-dependent falloffs, modified energy-condition bounds, and altered Hawking temperatures [2301.02534]. In quantum fluctuation modified gravity, the corresponding exponent becomes
\[
\gamma=\frac{2(1+3\omega-4\omega\alpha)}{2-3\alpha+\omega\alpha},
\]
so that the GR Kiselev exponent is recovered only in the limit \(\alpha\to0\) [2411.15854].

Rotating solutions have also been constructed. In \(f(R,T)\) gravity, the rotating Kiselev metric generated from a spherical seed has
\[
\Delta(r)=r^2+a^2-2Mr+K\,r^{-(d-2)},
\qquad
d=\frac{8[\gamma w+\pi(3w+1)]}{\gamma(3-w)+8\pi},
\]
and includes Kerr for \(K=0\) and Kerr–Newman for \(d=2\), \(K=Q^2\). The horizon structure exhibits one or two critical values of the coupling \(\gamma\), depending on the spin \(a\), and the black-hole domain in \((a,\gamma)\) is correspondingly segmented by extremal branches [2307.11611]. A different rotating construction, applied to the reduced Kiselev geometry, yields
\[
\Delta_{\rm RRK}=r^2+a^2-r_k^\alpha r^{2-\alpha},
\]
with allowed spin values smaller than in the corresponding Kerr family, an approximate shadow-shape degeneracy at fixed \(a/r_k\), and elliptic shadows for certain parameter ranges [2406.00788].

Regularization has proceeded in two complementary ways. One approach lets \(w\) vary radially and solves
\[
-r^2 B''(r)-2r(\bar w(r)+1)B'(r)-2\bar w(r)(B(r)-1)=0,
\]
obtaining regular black holes with finite \(R\), \(R_{\mu\nu}R^{\mu\nu}\), and Kretschmann scalar at the origin [2411.18804]. Another glues a Kiselev exterior to a de Sitter interior across a timelike dust shell, producing a stationary globally regular black hole with a stable shell for specific parameter ranges [2009.06990].

Thermodynamic-topology methods have introduced a different classification layer. For Kiselev–AdS black holes in \(f(R,T)\) gravity, the Duan \(\phi\)-mapping framework gives total topological charge \(Q_{\rm total}=-1\) in the temperature-based method and total topological number \(W=+1\) in the generalized Helmholtz free-energy method, while the coupled parameters \(\omega\) and \(\gamma\) alter the number and arrangement of individual defects and the photon-sphere structure [2410.02262]. A related study of quantum-corrected AdS–Reissner–Nordström black holes in Kiselev spacetime finds that most quintessence-like cases remain in the \(W=+1\) class, while a dust case can give \(W=0\) and a phantom case \(W=-1\) [2408.09782].

The subject has even entered analogue gravity. Starting from Gross–Pitaevskii theory, one acoustic construction reproduces
\[
f(r)=1-\frac{D}{r}+\frac{C_\varpi}{r^{3\varpi+1}},
\]
so that Schwarzschild, dust, radiation-like, quintessence, and \((A)dS\)-type Kiselev geometries all arise from a unified acoustic metric. In that setting, the quasibound-state spectra are analytically accessible, while quasinormal modes are computed by 6th-order WKB [2506.21639].

## 6. AdS viability, controversies, and open directions

The sharpest recent critique concerns anti–de Sitter embeddings. For the Kiselev–AdS metric
\[
f(r)=1+\frac{r^2}{L^2}-\frac{2M}{r}-\frac{K}{r^{1+3w}},
\]
asymptotic AdS behavior requires \(w>-1\), but the Seiberg–Witten brane-nucleation analysis shows that avoiding instability requires
\[
w\ge -\frac13.
\]
This is the opposite of the dark-energy range \(w<-\frac13\). Accordingly, for \(K>0\) and \( -1<w<-\frac13\), the Euclidean brane action becomes negative outside the horizon and is unbounded below, so Kiselev black holes in AdS are declared nonviable in that parameter range [2512.01265]. Spherical horizons are the least unstable; planar and hyperbolic horizons worsen the instability.

This result is deliberately specific. It does not show that no black hole can coexist with an anisotropic dark-energy fluid in AdS; it rules out the particular Kiselev embedding with that stress tensor and metric function. The same paper explicitly leaves open the endpoint of the instability after backreaction from copious brane production [2512.01265].

A broader conceptual controversy concerns nomenclature. Across the later literature, the “Kiselev quintessence” is repeatedly reinterpreted as an anisotropic medium rather than a genuine cosmological quintessence field, whether through perfect-fluid-plus-field decompositions, nonlinear-electrodynamics sourcing, or modified-gravity effective stresses [1910.08008, 2206.12876]. This does not invalidate the geometry as a mathematical or phenomenological model, but it changes the physical reading of lensing, thermodynamic, and observational inferences drawn from it.

Several open directions remain explicit in the current literature: alternative homotopies that remove weak-field dependence in analytic lensing constructions [1612.07279]; alternative anisotropic-fluid black holes consistent with AdS holography [2512.01265]; and further classes generated by tailored radial profiles \(w(r)\), including regular cores and nonstandard asymptotics [2411.18804]. Taken together, these lines of work place Kiselev black holes at the intersection of black-hole optics, anisotropic matter modeling, modified gravity, and analogue spacetime engineering, while keeping the status of the source interpretation and AdS consistency under active scrutiny.

Source: https://www.emergentmind.com/topics/kiselev-black-holes