---
title: Kiselev-Type Metric in Black Hole Physics
url: https://www.emergentmind.com/topics/kiselev-type-metric
type: topic
---

# Kiselev-Type Metric in Black Hole Physics

The Kiselev-type metric is a family of static, spherically symmetric geometries used to model black holes immersed in an anisotropic surrounding medium parameterized by an equation-of-state variable. In the literature represented here, it appears both as a direct generalization of the Schwarzschild or Reissner–Nordström metric and as a broader effective template for spacetimes sourced by quintessence-like matter, dust, radiation, or related anisotropic sectors. The same structure has also been used in analogue gravity, modified-gravity constructions, regular black-hole models, wormholes, and compact-star studies [2506.21639, 1508.04761].

## 1. Canonical form and parameterizations

A standard Kiselev line element is written as
\[
ds^2=-f(r)\,dt^2+f(r)^{-1}dr^2+r^2(d\theta^2+\sin^2\theta\,d\phi^2),
\]
with
\[
f(r)=1-\frac{2M}{r}+\frac{Q^2}{r^2}-\frac{\sigma}{r^{3\omega+1}}.
\]
Here \(M\) is the black-hole mass, \(Q\) is the electric charge, \(\sigma\) is a normalization parameter associated with the surrounding medium, and \(\omega\) is the state parameter entering \(p=\omega\rho\) in the effective description [1508.04761, 2106.03672]. In uncharged cases the \(Q^2/r^2\) term is omitted.

A second notation, used in analogue-gravity constructions, writes
\[
f(r)=1-\frac{D}{r}+\frac{C_\varpi}{r^{3\varpi+1}},
\]
where \(D\) is a draining parameter and \(C_\varpi\) is a fluid-dependent “charge” parameter [2506.21639]. This reflects a different parametrization of the same functional motif: a Schwarzschild-like \(1-r^{-1}\) term supplemented by a fluid term with exponent fixed by the state parameter.

Several special cases recur throughout the literature.

| Surrounding medium | Parameter choice | Metric function \(f(r)\) |
|---|---|---|
| Radiation | \(\omega=\tfrac{1}{3}\) | \(1-\frac{2\mathcal{M}}{r}+\frac{Q^2}{r^2}-\frac{\sigma_r}{r^2}\) |
| Dust | \(\omega=0\) | \(1-\frac{2\mathcal{M}}{r}+\frac{Q^2}{r^2}-\frac{\sigma_d}{r}\) |
| Quintessence | \(\omega=-\tfrac{2}{3}\), \(Q=0\) | \(1-\frac{2\mathcal{M}}{r}-\sigma_q r\) |

These cases are used to interpolate between standard vacuum black holes and “dirty” environments in which the asymptotic and horizon structures are modified by a surrounding field [1508.04761]. In thermodynamic studies, \(\omega=-1\) is also identified with the cosmological-constant case and \(\omega=-4/3\) with phantom energy [2106.03672].

## 2. Stress-energy content and source interpretation

Although the Kiselev family is often introduced through an equation-of-state parameter, the source is generically anisotropic rather than a perfect fluid. A general static, spherically symmetric stress tensor takes the form
\[
T^\mu{}_\nu=\mathrm{diag}(-\rho,\,p_r,\,p_t,\,p_t),
\]
and, in the Kiselev construction, one commonly finds
\[
p_r=-\rho,\qquad p_t=\frac{3w+1}{2}\rho,
\]
with the average pressure satisfying
\[
\bar p=\frac{p_r+2p_t}{3}=w\rho
\]
[2307.09079, 2009.06990]. This distinction between directional pressures is central: the Kiselev geometry is not, in general, a perfect-fluid spacetime.

For the one-component Kiselev black hole, the explicit stresses are
\[
\rho=-p_r=-\frac{3Kw}{8\pi r^{3(1+w)}},\qquad
p_t=-\frac{3Kw(1+3w)}{16\pi r^{3(1+w)}},
\]
so that
\[
\rho+p_r=0,\qquad
\rho+p_t=-\frac{9Kw(1+w)}{16\pi r^{3(1+w)}}.
\]
This directly controls the null energy condition (NEC): it is satisfied for \(Kw(1+w)<0\), marginal for \(Kw(1+w)=0\), and violated for \(Kw(1+w)>0\) [1910.08008].

A persistent misconception is that the Kiselev metric generically represents a perfect fluid. Later analyses show that it can instead be decomposed into a perfect-fluid component plus either an electromagnetic component or a scalar-field component:
\[
T^{\hat a\hat b}_{\text{total}}=
T^{\hat a\hat b}_{\text f}+
T^{\hat a\hat b}_{\text{em}}+
T^{\hat a\hat b}_{\text s}.
\]
The electromagnetic sector is selected when \(p_t-p_r>0\), with
\[
E^2=\max\{p_t-p_r,0\},
\]
while the scalar sector is selected when \(p_r-p_t>0\), with
\[
(\nabla\phi)^2=\max\{p_r-p_t,0\}
\]
[1910.08008]. In the generalized \(N\)-component case, the same criterion is recast in terms of the density gradient \(\rho'(r)\), leading to an “onion-like” radial structure in which the effective support can change with radius [1910.08008].

A further reinterpretation replaces the anisotropic fluid altogether by nonlinear electrodynamics. In power-Maxwell electrodynamics, the Kiselev geometry arises as an exact solution of Einstein’s equations coupled to a power-Maxwell Lagrangian, for both electric and magnetic ansätze, with the metric written as
\[
g(r)=1-\frac{2M}{r}-\frac{k}{r^{3w+1}}.
\]
In that framework the “quintessence” source is recast as a nonlinear electromagnetic sector rather than an ordinary fluid [2206.12876].

## 3. Horizons, causal structure, and geometric variants

Horizons are determined by the roots of \(f(r)=0\), and the number and character of those roots depend on the state parameter and normalization of the surrounding medium. For the widely studied quintessence case \(w=-2/3\),
\[
f(r)=1-\frac{2M}{r}-\sigma r,
\]
and the horizons are
\[
r_\pm=\frac{1\pm\sqrt{1-8M\sigma}}{2\sigma},
\qquad
0<\sigma\le \frac{1}{8M}.
\]
This yields the standard trichotomy: two horizons for \(0<\sigma<1/(8M)\), a degenerate horizon at \(\sigma=1/(8M)\), and no horizon beyond that threshold, producing a naked singularity [1605.02320, 1502.01676].

The horizon analysis extends directly to charged variants. For example, in the radiation case one has
\[
r_\pm=\mathcal{M}\pm\sqrt{\mathcal{M}^2-Q^2+\sigma_r},
\]
while in the dust case
\[
r_\pm=\frac{2\mathcal{M}+\sigma_d\pm\sqrt{(2\mathcal{M}+\sigma_d)^2-4Q^2}}{2}
\]
[1508.04761, 1812.05479]. These formulas make explicit that the ambient field can shift horizon radii and alter extremality conditions.

A less-studied regime, termed the reduced Kiselev black hole, takes \(-1/3<w<0\) and removes the standard mass term. Writing \(\alpha=3w+1\) with \(0<\alpha<1\), the metric function becomes
\[
f(r)=1-\left(\frac{K}{r}\right)^\alpha.
\]
Even with \(r_g=0\), a black-hole-type horizon appears at \(r=K\). The causal structure is Schwarzschild-like, \(r=K\) is a Killing horizon, and the thermodynamic quantities are
\[
\kappa=\frac{\alpha}{2K},\qquad
T=\frac{\alpha}{4\pi K},\qquad
S=\pi K^2
\]
[2307.09079].

Rotating generalizations have been constructed through Newman–Janis-type procedures. In \(f(R,T)\) gravity, the rotating solution is written with
\[
\Sigma=r^2+a^2\cos^2\theta,\qquad
\Delta=r^2+a^2-2rM(r),\qquad
M(r)=M-\frac{K}{2r^{d-1}},
\]
where the exponent
\[
d=\frac{8(\gamma w+\pi(3w+1))}{\gamma(3-w)+8\pi}
\]
depends on both the state parameter \(w\) and the matter-geometry coupling \(\gamma\). The solution reduces to Kerr for \(K=0\) and to Kerr–Newman for \(d=2\), \(K=Q^2\) [2307.11611].

## 4. Accretion, lensing, chaos, and perturbative probes

The Kiselev family has been used extensively as a dynamical background. In spherical accretion without back-reaction, baryon-number conservation gives
\[
r^2nu=C_1,
\qquad
\dot M=4\pi r^2mnu,
\]
while energy-momentum conservation yields a Bernoulli equation,
\[
\left(\frac{\rho+p}{n}\right)^2
\left[1-\frac{2M}{r}-\sigma r+u^2\right]
=
\left(\frac{\rho_\infty+p_\infty}{n_\infty}\right)^2.
\]
The sonic-point conditions are
\[
u_c^2=\frac14\left(\frac{2M}{r_c}-\sigma r_c\right),
\qquad
a_c^2=\frac{u_c^2}{1-\frac{2M}{r_c}-\sigma r_c+u_c^2},
\]
and the quintessence parameter strongly affects the location of the critical point and the mass accretion rate [1605.02320].

For null geodesics, the \(w_q=-2/3\) Kiselev black hole admits a detailed strong-lensing analysis. The effective potential is
\[
V_{\rm eff}(r)=\frac{L^2}{r^2}\left(1-\frac{2M}{r}-\sigma r\right),
\]
the photon sphere is located at
\[
r_{\rm ps}=\frac{1-\sqrt{1-6M\sigma}}{\sigma},
\]
and the exact bending angle can be written in terms of elliptic integrals. In that specific study, the bending angles satisfy the ordering
\[
\text{naked singularity}>\text{extreme Kiselev black hole}>
\text{non-extreme Kiselev black hole}>
\text{Schwarzschild black hole}
\]
for comparable impact parameters [1502.01676].

Charged-particle dynamics exhibits additional structure. For the charged Kiselev black hole,
\[
F(r)=1-\frac{2M}{r}+\frac{Q^2}{r^2}-\frac{\alpha}{r^{3w+1}},
\]
the Lyapunov exponent governing orbital instability is computed from the effective potential and a Jacobian-matrix treatment of perturbations. Near the horizon, the chaos bound is satisfied for fixed charge-to-mass ratio, whereas away from the horizon it can be violated when the black-hole charge \(Q\) and particle angular momentum \(L\) are sufficiently large and the normalization factor \(\alpha\) is small [2204.07983].

Perturbation theory also enters through scalar waves. In analogue Kiselev acoustic black holes, the massless Klein–Gordon equation separates into radial equations of Heun type, quasibound states can be obtained by the Vieira–Bezerra–Kokkotas method, and quasinormal modes can be extracted with sixth-order WKB methods [2506.21639].

## 5. Thermodynamics, holography, and evaporation

Thermodynamic studies of Kiselev black holes emphasize horizon products, heat capacities, and deformations of the Hawking sector. For radiation and dust backgrounds, the products of horizon areas and entropies are mass-independent:
\[
\mathcal A_+\mathcal A_-=16\pi^2(Q^2-\sigma_r)^2,\qquad
\mathcal S_+\mathcal S_-=\pi^2(Q^2-\sigma_r)^2
\]
for radiation, and
\[
\mathcal A_+\mathcal A_-=16\pi^2Q^4,\qquad
\mathcal S_+\mathcal S_-=\pi^2Q^4
\]
for dust. By contrast, for the Schwarzschild black hole surrounded by quintessence, the corresponding products depend on the mass and are therefore not universal [1508.04761].

The same thermodynamic relations have been used in a CFT setting. For Kiselev black holes surrounded by radiation or dust, the universal entropy products support a Kiselev/CFT correspondence in which the left- and right-moving central charges are found to be equal:
\[
c_L=c_R=24M(2M^2-q^2+\sigma_r)
\]
for the radiation case, with an analogous equality in the dust case [1812.05479].

Beyond general relativity, Rainbow-gravity studies use
\[
ds^2=
-\frac{1}{f^2(E/E_p)}
\left(1-\frac{2GM}{r}-\frac{c}{r^{3\omega+1}}\right)dt^2
+\frac{1}{g^2(E/E_p)}
\left(1-\frac{2GM}{r}-\frac{c}{r^{3\omega+1}}\right)^{-1}dr^2
+\frac{r^2}{g^2(E/E_p)}d\Omega^2,
\]
with \(f(E/E_p)=1\) and \(g(E/E_p)=\sqrt{1-\eta(E/E_p)^n}\). In that setting, cosmic-fluid effects impose a maximum allowed horizon radius, the Hawking temperature can vanish at a finite critical horizon where black-hole and cosmological horizons merge, and black-hole remnants do not generically appear except in the peculiar case \(n=4\), \(\eta<0\) [2106.03672].

A distinct thermodynamic construction links the Kiselev metric to Rényi entropy. Requiring the Hawking temperature to equal the Rényi temperature yields
\[
\bar f(r,\lambda)=1-\frac{2M}{r}-2\pi\lambda Mr+f_0(r),
\]
which is of Kiselev form with \(w=-2/3\) and \(K=2\pi\lambda M\). In the Schwarzschild-like case \(f_0=0\), the geometry admits a maximum mass
\[
M_{\max}=\frac{1}{4\sqrt{\pi\lambda}},
\]
at which the temperature vanishes and the horizons coincide; the third law is then interpreted as a cosmic censor preventing dynamical access to the naked-singularity regime [2504.16705].

Evaporation studies introduce a dynamical state parameter \(w_q\) in
\[
f(r)=1-\frac{2M}{r}-\frac{a}{r^{3w_q+1}}.
\]
Lowering \(w_q\) reduces the non-final-stage temperature and markedly prolongs the evaporation lifetime; as \(w_q\to -1\), the metric approaches Schwarzschild–de Sitter behavior and enters an ultra-slowly evaporating regime. That mechanism is explicitly distinguished from the ultra-long lifetimes found in PFDM and Horndeski black holes [2606.19110].

## 6. Analogue, regularized, and matter-coupled extensions

The Kiselev form is sufficiently flexible to reappear in several adjacent research programs. In analogue gravity, Gross–Pitaevskii theory for Bose–Einstein condensates yields an acoustic metric
\[
ds^2=\sqrt{3}\,c_s^2
\left[
-f(r)\,dt^2+\frac{dr^2}{f(r)}+r^2d\vartheta^2+r^2\sin^2\vartheta\,d\phi^2
\right],
\]
with a radial flow engineered so that
\[
v_r\sim\sqrt{\frac{D}{r}-\frac{C_\varpi}{r^{3\varpi+1}}}.
\]
This realizes analogue Kiselev acoustic black holes and provides quasinormal and quasibound spectra for scalar perturbations [2506.21639].

Regularity constructions replace the Kiselev singular core by de Sitter space and match the two regions across a thin shell using Barrabes–Israel junction conditions. The exterior is
\[
f_+(r)=1-\frac{2m}{r}-\frac{c}{r^{3\omega+1}},
\]
the interior is
\[
f_-(r)=1-\frac{\Lambda}{3}r^2,
\]
and stable stationary solutions arise for suitable parameter ranges, yielding nonsingular black-hole spacetimes with a de Sitter core [2009.06990].

The same anisotropic template has been generalized to traversable wormholes. With the Morris–Thorne ansatz
\[
ds^2=-e^{2\Phi(r)}dt^2+\left(1-\frac{b(r)}{r}\right)^{-1}dr^2+r^2d\Omega^2,
\]
and Kiselev-inspired density
\[
\rho(r)=\frac{c}{8\pi r^{3(\omega+1)}},
\]
integration gives the shape function
\[
b(r)=r_0\left(\frac{r_0}{r}\right)^{3\omega}.
\]
The flare-out condition requires \(\omega>-1/3\), while the NEC and WEC are violated in the physically traversable sector, so the configuration remains supported by exotic matter in the usual wormhole sense [2410.17475].

Compact-star applications replace the phenomenological anisotropic fluid by a chameleon scalar field. In that setting the Kiselev-type metric
\[
f(r)=1-\frac{2M}{r}-\frac{a}{r^{3w_q+1}}
\]
is coupled to modified TOV equations, and scalar gradients induce anisotropy through
\[
\Delta p=p_r-p_t=(\partial_r\phi)^2.
\]
The scalar is screened in the high-density core and unscreened outside a critical radius \(r_{\rm crit}\propto m_\phi^{-1}\), producing stable neutron-star configurations consistent with the conservative mass, radius, and tidal-deformability bounds quoted in the study [2508.04744].

Taken together, these developments show that the Kiselev-type metric functions less as a single-source solution than as a reusable geometric ansatz for anisotropic environments. The precise physical interpretation depends on context—quintessence-like matter, nonlinear electrodynamics, effective laboratory media, thin-shell matching, or screened scalars—but the defining signature remains the same: a Schwarzschild-type potential deformed by a state-parameter-dependent power law.

Source: https://www.emergentmind.com/topics/kiselev-type-metric