Kiselev Black Hole Solution

• The Kiselev black hole solution is a static, spherically symmetric model where the Schwarzschild metric is modified by the inclusion of an anisotropic fluid characterized by a free parameter w.
• It generalizes classical black hole models by introducing a power-law matter term in the metric function, linking to nonlinear electrodynamics, modified gravity, and various state equations.
• This framework influences horizon structure, thermodynamic behavior, and gravitational lensing, making it pivotal for understanding complex astrophysical and cosmological phenomena.

The Kiselev black hole solution describes a static, spherically symmetric black hole surrounded by an anisotropic fluid whose equation of state is characterized by a free parameter $w$. It generalizes the Schwarzschild geometry by introducing a power-law matter term in the metric function, providing a flexible framework for modeling black holes embedded in quintessential, radiation, dust, and other fluid environments. Recent developments have connected the Kiselev metric to nonlinear electrodynamics, modified gravity, black hole thermodynamics including Rényi-type entropies, regular black hole models, and extensions to rotation and cosmological backgrounds.

1. Metric Structure and Anisotropic Fluid Source

The canonical Kiselev metric in Schwarzschild-like coordinates $(t,r,\theta,\varphi)$ is

$$ds^2 = -f(r)dt^2 + \frac{dr^2}{f(r)} + r^2(d\theta^2 + \sin^2\theta\, d\varphi^2)$$

with

$f(r) = 1 - \frac{2M}{r} - \frac{c}{r^{3w+1}}$

where $M$ is the ADM mass, $c$ the fluid density normalization, and $w$ is a constant barotropic index $(p_r=w\rho)$ (Dariescu et al., 2022). The stress–energy tensor of the anisotropic fluid is

$$T^\mu{}_\nu = \mathrm{diag}(-\rho, p_r, p_t, p_t),\qquad p_r = w\rho, \quad p_t = \frac{1}{2}(1+3w)\rho$$

(Abbas et al., 2019).

For $w$ within $(t,r,\theta,\varphi)$0 and $(t,r,\theta,\varphi)$1, the energy density is positive and the spacetime admits two horizons—a black hole and a cosmological horizon—analogous to Schwarzschild–de Sitter or quintessence black holes.

2. Horizon Structure, Energy Conditions, and Limiting Cases

Horizons are roots of $(t,r,\theta,\varphi)$2,

$(t,r,\theta,\varphi)$3

with explicit solutions for special $(t,r,\theta,\varphi)$4, e.g., $(t,r,\theta,\varphi)$5 yields a quadratic

$(t,r,\theta,\varphi)$6

with a two-horizon structure for $(t,r,\theta,\varphi)$7 (Boonserm et al., 2019, Abbas et al., 2019).

Energy conditions depend critically on the fluid parameters:

• The null energy condition (NEC) requires $(t,r,\theta,\varphi)$8, and is violated for $(t,r,\theta,\varphi)$9, corresponding to exotic fluids (Boonserm et al., 2019).
• The stress–energy tensor can be decomposed into a perfect fluid plus either an electromagnetic or scalar field for $$ds^2 = -f(r)dt^2 + \frac{dr^2}{f(r)} + r^2(d\theta^2 + \sin^2\theta\, d\varphi^2)$$0-component extensions (Boonserm et al., 2019).

Limiting cases include:

• $$ds^2 = -f(r)dt^2 + \frac{dr^2}{f(r)} + r^2(d\theta^2 + \sin^2\theta\, d\varphi^2)$$1: Cosmological constant, $$ds^2 = -f(r)dt^2 + \frac{dr^2}{f(r)} + r^2(d\theta^2 + \sin^2\theta\, d\varphi^2)$$2 (Schwarzschild–(A)dS).
• $$ds^2 = -f(r)dt^2 + \frac{dr^2}{f(r)} + r^2(d\theta^2 + \sin^2\theta\, d\varphi^2)$$3: Dust, $$ds^2 = -f(r)dt^2 + \frac{dr^2}{f(r)} + r^2(d\theta^2 + \sin^2\theta\, d\varphi^2)$$4.
• $$ds^2 = -f(r)dt^2 + \frac{dr^2}{f(r)} + r^2(d\theta^2 + \sin^2\theta\, d\varphi^2)$$5: Radiation, $$ds^2 = -f(r)dt^2 + \frac{dr^2}{f(r)} + r^2(d\theta^2 + \sin^2\theta\, d\varphi^2)$$6.
• $$ds^2 = -f(r)dt^2 + \frac{dr^2}{f(r)} + r^2(d\theta^2 + \sin^2\theta\, d\varphi^2)$$7: Border of acceleration, resulting in a constant shift (Abbas et al., 2019, Santos, 2024).

3. Kiselev Solution in Nonlinear Electrodynamics and Power–Maxwell Theory

The Kiselev geometry can be derived as an exact solution of Einstein equations coupled to nonlinear electrodynamics, specifically "power–Maxwell" theory: $$ds^2 = -f(r)dt^2 + \frac{dr^2}{f(r)} + r^2(d\theta^2 + \sin^2\theta\, d\varphi^2)$$8 where $$ds^2 = -f(r)dt^2 + \frac{dr^2}{f(r)} + r^2(d\theta^2 + \sin^2\theta\, d\varphi^2)$$9 is a real power. Both electric and magnetic ansatzes can generate the Kiselev form, with distinct relations between the power $f(r) = 1 - \frac{2M}{r} - \frac{c}{r^{3w+1}}$0, the equation-of-state parameter, and the metric exponent (Dariescu et al., 2022). For example, the metric retains

$f(r) = 1 - \frac{2M}{r} - \frac{c}{r^{3w+1}}$1

with $f(r) = 1 - \frac{2M}{r} - \frac{c}{r^{3w+1}}$2 related to $f(r) = 1 - \frac{2M}{r} - \frac{c}{r^{3w+1}}$3 and the electrodynamic parameters. This demonstrates that the Kiselev metric naturally emerges in the context of nonlinear sources, not merely as a phenomenological fluid.

Thermodynamic properties such as the Hawking temperature,

$f(r) = 1 - \frac{2M}{r} - \frac{c}{r^{3w+1}}$4

and the distinctive heat capacity peak (Schottky anomaly), signal the potential for cyclic heat engine behavior (Dariescu et al., 2022).

4. Extensions: Magnetized Solutions, Rotation, and Modified Gravity

Magnetized Kiselev Black Holes

A Harrison–Ernst transformation applied to the seed Kiselev metric yields a family with asymptotically Melvin–type magnetic fields: $f(r) = 1 - \frac{2M}{r} - \frac{c}{r^{3w+1}}$5 with $f(r) = 1 - \frac{2M}{r} - \frac{c}{r^{3w+1}}$6 and $f(r) = 1 - \frac{2M}{r} - \frac{c}{r^{3w+1}}$7 (Lungu, 7 Apr 2025). The magnetic field parameter $f(r) = 1 - \frac{2M}{r} - \frac{c}{r^{3w+1}}$8 modifies photon spheres and increases lensing effects.

Rotating Kiselev Metrics and Modified Gravity

Rotating reductions generalize the metric via Newman–Janis-type algorithms. In spherical symmetry, the seed metric is extended to

$f(r) = 1 - \frac{2M}{r} - \frac{c}{r^{3w+1}}$9

and rotation introduces additional parameters (spin $M$0), shifting the horizon structure. The effective mass and horizon radii are modified, e.g.,

$M$1

where $M$2 and $M$3 depend on the specific modified gravity theory, such as $M$4 gravity (Ghosh et al., 2023, Santos et al., 2023, Benali et al., 2024, Czinner et al., 30 Aug 2025).

5. Thermodynamic Properties: Rényi Entropy and Cosmic Censorship

By demanding the black hole thermodynamic temperature derived from Rényi entropy matches the surface gravity (Hawking temperature), the metric function is uniquely fixed to the Kiselev form: $M$5 with $M$6 and $M$7 (Czinner et al., 23 Apr 2025). For such solutions,

$M$8

and the Rényi entropy

$M$9

satisfies the first law $c$0. The heat capacity is strictly negative, pointing to local thermodynamic instability. The third law enforces an upper mass bound, $c$1, prohibiting horizonless solutions (cosmic censorship).

6. Regular Black Hole Extensions and Horizon Stability

Allowing the equation-of-state parameter $c$2 to vary with radius provides general regular solutions, avoiding singularities at the origin. For specific choices of $c$3, one obtains metrics

$c$4

with finite curvature invariants at $c$5, satisfying both weak and strong energy conditions over suitable parameter ranges (Santos, 2024).

In asymptotically AdS backgrounds, the Kiselev metric with $c$6 is unstable against brane nucleation, rendering such solutions nonviable for AdS dark energy environments (Xia et al., 1 Dec 2025).

7. Geodesics, Lensing, and Observational Signatures

The null geodesic structure and gravitational lensing signatures for Kiselev black holes, including charge and magnetization, are rich. The parameter $c$7 and the density $c$8 impact photon sphere radii, shadow shapes, and deflection angles. In particular, terms proportional to $c$9 or higher powers (for $w$0 values such as $w$1) introduce nontrivial corrections to the bending of light which may be observationally distinguishable (Lungu, 7 Apr 2025, Abbas et al., 2019, Shchigolev et al., 2016).

8. Summary Table of Metric Variants

Solution Type	Metric Function $w$2	Source/Context
Static Kiselev	$w$3	Original, anisotropic fluid
Magnetized Kiselev	$w$4; $w$5	Harrison–Ernst magnetization
Power–Maxwell Electrodynamics	$w$6	Nonlinear electrodynamics
Rotating Kiselev	See $w$7, $w$8, etc.	Newman–Janis rotation, $w$9 or GR
Rényi–Kiselev	$(p_r=w\rho)$0	Thermodynamics, cosmic censorship
Regular (Generalized $(p_r=w\rho)$1)	$(p_r=w\rho)$2, etc.	Radial-varying $(p_r=w\rho)$3, regular center

• Null geodesics and magnetized Kiselev: (Lungu, 7 Apr 2025)
• Power–Maxwell electrodynamic realization: (Dariescu et al., 2022)
• Stress–energy decomposition and NEC: (Boonserm et al., 2019)
• Strong lensing and shadow structure: (Abbas et al., 2019, Shchigolev et al., 2016, Czinner et al., 30 Aug 2025, Benali et al., 2024)
• Quantum modified gravity and horizon structure: (Hua et al., 2024)
• Regular solution with $(p_r=w\rho)$4: (Santos, 2024)
• Rényi black hole thermodynamics: (Czinner et al., 23 Apr 2025)
• AdS instability: (Xia et al., 1 Dec 2025)
• Bardeen–Kiselev with cosmological constant: (Rodrigues et al., 2022, Wu et al., 2022)
• Rotating extensions: (Ghosh et al., 2023, Santos et al., 2023, Benali et al., 2024)

The Kiselev black hole solution and its generalizations provide a versatile and rigorous framework for analyzing black-hole physics in the presence of general anisotropic fluids, nonlinear electrodynamics, modified gravity, and nonstandard thermodynamics. This continues to have direct implications for lensing, horizon structure, cosmic censorship, and the stability of black holes in complex astrophysical and cosmological settings.