Magnetized Kiselev Black Hole

• Magnetized Kiselev black hole is an exact solution to Einstein's equations that combines a dark energy–like quintessential field with a uniform Melvin magnetic field to modify spacetime geometry.
• The interplay of the magnetic parameter B0 and the quintessence parameter k introduces unique orbital dynamics, photon sphere structures, and lensing effects not seen in standard black holes.
• The metric's rich thermodynamic behavior and altered geodesic structure offer critical insights for testing strong-field gravity and exotic matter fields in both astrophysical and cosmological contexts.

A magnetized Kiselev black hole represents a class of exact solutions to Einstein's equations describing a static or stationary black hole immersed both in a dark energy–like quintessential field and in a uniform (Melvin-type) magnetic field. These geometries generalize the classic Kiselev black holes—which incorporate a stress-energy component mimicking quintessence, parameterized by an equation-of-state variable $w$ and density parameter $k$—by embedding them in an external (test or self-consistent) magnetic field characterized by a parameter $B_0$ (Dariescu et al., 6 Aug 2025, Dariescu et al., 12 Jan 2026, Lungu, 7 Apr 2025, Lungu et al., 2024). Both the magnetic and quintessence modifications introduce qualitative changes in geodesic structure, effective potentials, horizon properties, and observable effects such as lensing, periapsis precession, and the possible existence of multiple photon spheres.

1. Origin, Action, and Metric Structure

The magnetized Kiselev black hole first emerges as a solution to Einstein–Maxwell (or more generally Einstein–power–Maxwell) equations sourced by a nonlinear electromagnetic field in addition to an anisotropic fluid corresponding to quintessence (Dariescu et al., 2022, Dariescu et al., 2023).

The metric for the physically most studied case ($w = -2/3$) is

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

with

$$f(r) = 1 - \frac{2M}{r} - k r, \qquad \Lambda(r, \theta) = 1 + B_{0}^{2} r^{2} \sin^{2}\theta.$$

Here $M$ is the ADM mass, $k$ the quintessence parameter (related to the density and $w$), and $B_{0}$ the asymptotic magnetic field amplitude (Ernst–Melvin parameter).

In the more general case, the quintessential term may take the form $k$0 and the exponent $k$1 (Dariescu et al., 6 Aug 2025). The metric reduces to well-known solutions in limiting cases:

• $k$2: the Schwarzschild–Melvin solution,
• $k$3: the standard Kiselev black hole,
• $k$4: the Melvin magnetic universe.

2. Electromagnetic Field and Stress-Energy Sources

The electromagnetic field is typically implemented via a vector potential with the only nonzero component

$k$5

introducing an asymptotically uniform (Melvin-type) magnetic field aligned with the $k$6-axis (Lungu et al., 2024, Dariescu et al., 6 Aug 2025). The energy-momentum tensor encapsulates both the non-linear magnetic contributions (possibly via a Lagrangian $k$7, with $k$8, $k$9), and the anisotropic fluid with $B_0$0, $B_0$1 (Dariescu et al., 2022, Dariescu et al., 2023).

The quintessential component exhibits an equation of state $B_0$2, with $B_0$3.

3. Horizon Structure and Causal Features

The horizons of the solution are given as real, positive roots of

$B_0$4

which yields

$B_0$5

$B_0$6 corresponds to the event (black-hole) horizon, while $B_0$7 is an outer, cosmological–like horizon generated by quintessence. For $B_0$8, these coalesce in an extremal horizon at $B_0$9; for $w = -2/3$0, the solution describes a naked singularity (Lungu, 7 Apr 2025, Lungu et al., 2024).

The curvature invariants (Ricci, Kretschmann), diverge at $w = -2/3$1, confirming a central singularity. At spatial infinity, the spacetime is neither Ricci-flat nor asymptotically flat; the geometry is dominated by the linear quintessence term and the Melvin field (Lungu et al., 2024).

4. Geodesics and Effective Potentials

4.1 Timelike and Charged Particle Motion

Charged test particles with mass $w = -2/3$2, charge $w = -2/3$3, and specific charge $w = -2/3$4, follow the Lagrangian

$w = -2/3$5

with conserved energy and angular momentum (Dariescu et al., 6 Aug 2025, Dariescu et al., 12 Jan 2026): $w = -2/3$6 The effective potential for equatorial ($w = -2/3$7) motion is

$w = -2/3$8

Circular (and bound) orbits are determined by solving $w = -2/3$9, $$ds^{2} = -f(r)\,\Lambda^{2}\,dt^{2} + \frac{\Lambda^{2}}{f(r)}\,dr^{2} + \Lambda^{2} r^{2}\,d\theta^{2} + \frac{r^{2}\sin^{2}\theta}{\Lambda^{2}}\,d\varphi^{2}$$0, with stability set by $$ds^{2} = -f(r)\,\Lambda^{2}\,dt^{2} + \frac{\Lambda^{2}}{f(r)}\,dr^{2} + \Lambda^{2} r^{2}\,d\theta^{2} + \frac{r^{2}\sin^{2}\theta}{\Lambda^{2}}\,d\varphi^{2}$$1 (Lungu et al., 2024, Dariescu et al., 12 Jan 2026). The effective potential exhibits a double-barrier (well) structure between $$ds^{2} = -f(r)\,\Lambda^{2}\,dt^{2} + \frac{\Lambda^{2}}{f(r)}\,dr^{2} + \Lambda^{2} r^{2}\,d\theta^{2} + \frac{r^{2}\sin^{2}\theta}{\Lambda^{2}}\,d\varphi^{2}$$2 and $$ds^{2} = -f(r)\,\Lambda^{2}\,dt^{2} + \frac{\Lambda^{2}}{f(r)}\,dr^{2} + \Lambda^{2} r^{2}\,d\theta^{2} + \frac{r^{2}\sin^{2}\theta}{\Lambda^{2}}\,d\varphi^{2}$$3; bound orbits (including relativistic rosette orbits) and a variety of escape orbits are possible (Dariescu et al., 6 Aug 2025).

4.2 Null Geodesics and Photon Spheres

Null geodesics are governed by an analogous Lagrangian with zero-mass constraint. The equatorial effective potential for photons is

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

Unstable (and, owing to the magnetization, sometimes stable) circular photon orbits (photon spheres) arise as extremal points of $$ds^{2} = -f(r)\,\Lambda^{2}\,dt^{2} + \frac{\Lambda^{2}}{f(r)}\,dr^{2} + \Lambda^{2} r^{2}\,d\theta^{2} + \frac{r^{2}\sin^{2}\theta}{\Lambda^{2}}\,d\varphi^{2}$$5 (Lungu, 7 Apr 2025). The radius of the photon sphere satisfies a quartic equation whose structure and roots depend sensitively on $$ds^{2} = -f(r)\,\Lambda^{2}\,dt^{2} + \frac{\Lambda^{2}}{f(r)}\,dr^{2} + \Lambda^{2} r^{2}\,d\theta^{2} + \frac{r^{2}\sin^{2}\theta}{\Lambda^{2}}\,d\varphi^{2}$$6.

5. Orbital Precession, Epicyclic Frequencies, and Dynamical Phenomena

5.1 Periapsis Shift

The periapsis shift per revolution for charged particles, crucial for probing both magnetic and quintessence effects, admits two complementary derivations:

• Epicyclic-frequency method: The periapsis shift

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

with $$ds^{2} = -f(r)\,\Lambda^{2}\,dt^{2} + \frac{\Lambda^{2}}{f(r)}\,dr^{2} + \Lambda^{2} r^{2}\,d\theta^{2} + \frac{r^{2}\sin^{2}\theta}{\Lambda^{2}}\,d\varphi^{2}$$8, $$ds^{2} = -f(r)\,\Lambda^{2}\,dt^{2} + \frac{\Lambda^{2}}{f(r)}\,dr^{2} + \Lambda^{2} r^{2}\,d\theta^{2} + \frac{r^{2}\sin^{2}\theta}{\Lambda^{2}}\,d\varphi^{2}$$9 the proper-time orbital and radial epicyclic frequencies (Dariescu et al., 6 Aug 2025, Dariescu et al., 12 Jan 2026).

• Direct integral approach: For generic bound orbits,

$$f(r) = 1 - \frac{2M}{r} - k r, \qquad \Lambda(r, \theta) = 1 + B_{0}^{2} r^{2} \sin^{2}\theta.$$0

with the integrand constructed from the previously defined first integrals.

For uncharged particles, the periapsis advance is always prograde ($$f(r) = 1 - \frac{2M}{r} - k r, \qquad \Lambda(r, \theta) = 1 + B_{0}^{2} r^{2} \sin^{2}\theta.$$1), but for charged particles in sufficiently strong magnetic fields, a retrograde ($$f(r) = 1 - \frac{2M}{r} - k r, \qquad \Lambda(r, \theta) = 1 + B_{0}^{2} r^{2} \sin^{2}\theta.$$2) precession regime appears. The boundary between prograde and retrograde is the locus $$f(r) = 1 - \frac{2M}{r} - k r, \qquad \Lambda(r, \theta) = 1 + B_{0}^{2} r^{2} \sin^{2}\theta.$$3. Explicitly, in the weak-field expansion (Dariescu et al., 6 Aug 2025, Dariescu et al., 12 Jan 2026): $$f(r) = 1 - \frac{2M}{r} - k r, \qquad \Lambda(r, \theta) = 1 + B_{0}^{2} r^{2} \sin^{2}\theta.$$4

$$f(r) = 1 - \frac{2M}{r} - k r, \qquad \Lambda(r, \theta) = 1 + B_{0}^{2} r^{2} \sin^{2}\theta.$$5

where $$f(r) = 1 - \frac{2M}{r} - k r, \qquad \Lambda(r, \theta) = 1 + B_{0}^{2} r^{2} \sin^{2}\theta.$$6 is a positive constant dependent on $$f(r) = 1 - \frac{2M}{r} - k r, \qquad \Lambda(r, \theta) = 1 + B_{0}^{2} r^{2} \sin^{2}\theta.$$7.

5.2 Epicyclic and Nodal Precession

Epicyclic frequencies—radial $$f(r) = 1 - \frac{2M}{r} - k r, \qquad \Lambda(r, \theta) = 1 + B_{0}^{2} r^{2} \sin^{2}\theta.$$8, latitudinal $$f(r) = 1 - \frac{2M}{r} - k r, \qquad \Lambda(r, \theta) = 1 + B_{0}^{2} r^{2} \sin^{2}\theta.$$9, azimuthal $M$0—are modified by both $M$1 and $M$2. The nodal (Larmor) precession of slightly tilted orbits is encoded in (Dariescu et al., 12 Jan 2026): $M$3 reflecting a purely magnetic, non–frame-dragging precession.

5.3 Instabilities and 3D Trajectories

The Lorentz force due to the monopolar magnetic field generically confines charged particles to cones of constant half-angle (the "Poincaré cone" structure), with opening determined by $M$4 (Dariescu et al., 2023). Off-equatorial stability bands and resonance structures appear, particularly in the Mathieu-perturbation regime when quintessence dominates (Lungu et al., 2024).

6. Gravitational Lensing, Photon Rings, and Observational Aspects

The presence of quintessence ($M$5) and magnetic field ($M$6) strongly affects null geodesics, and thus the shadow and lensing structure: $M$7 in the weak-field limit. A positive $M$8 enhances deflection at large impact parameter, while magnetic corrections become significant at moderate-to-large $M$9 or $k$0. Stable and unstable photon rings are possible, leading to richer shadow structures compared to the Schwarzschild or Kiselev limits (Lungu, 7 Apr 2025).

The interplay of $k$1 and $k$2 may in principle be observationally probed via multi-wavelength VLBI interferometry, looking for deviations in black hole shadows and ring structure.

7. Thermodynamics and Stability

Thermodynamic analysis reveals that the temperature, entropy, and magnetic potential of the magnetized Kiselev black hole are sensitive to the exponents set by the dimension of the quintessence parameter and the underlying power–Maxwell index q (Dariescu et al., 2022). The solutions are thermally unstable over physically relevant parameter ranges, with negative heat capacity and the typical multi-horizon "Schottky peak." The magnetized Kiselev geometries thus share the instability familiar from de Sitter–like multi-horizon black holes.

---

The magnetized Kiselev black hole represents a theoretically robust and phenomenologically rich class of spacetimes with explicit parameter dependence on both astrophysical (magnetic, mass) and cosmological (quintessence) scales, offering avenues for precision tests of general relativity, exotic matter fields, and electromagnetic interactions in the strong-field regime. The combined influence of $k$3 and $k$4 yields novel features in orbit dynamics, lensing phenomenology, and thermodynamic behaviour, distinguishing this class sharply from both Ernst–Melvin and conventional Kiselev solutions (Dariescu et al., 6 Aug 2025, Dariescu et al., 12 Jan 2026, Lungu, 7 Apr 2025, Lungu et al., 2024, Dariescu et al., 2023, Dariescu et al., 2022).