Dark Energy Cools Black Holes: Exact Quantum Solution Reveals Quieter Hawking Radiation

Black holes are not perfectly black. In 1975, Stephen Hawking showed that quantum effects near an event horizon force these gravitational monsters to emit a faint thermal glow, a discovery that fused general relativity with quantum theory and gave black holes a measurable temperature. Half a century later, physicists are still asking how that glow changes when the black hole is not sitting in empty space but is embedded in the dark energy that dominates our accelerating universe. A new theoretical study published in The European Physical Journal C by M. D. de Oliveira and Alexandre G. M. Schmidt of the Federal Fluminense University in Brazil delivers one of the most complete answers yet, solving exactly the quantum equation that governs a relativistic particle in the spacetime of a black hole surrounded by a dark-energy-like fluid.

The stage for the calculation is the Kiselev black hole, an exact solution of Einstein’s equations found by V. V. Kiselev in 2003 that describes a static, uncharged black hole wrapped in what is often called a quintessence field. Quintessence is a proposed dynamical form of dark energy, a scalar field carrying negative pressure that could be driving the accelerated expansion of the cosmos. As the authors emphasize, following a definitive analysis by Matt Visser, the Kiselev spacetime is intrinsically anisotropic and cannot be treated as a perfect fluid, so it is best understood as a toy model: a black hole surrounded by a quintessence-like anisotropic fluid rather than a realistic portrait of cosmic dark energy. Even as a toy model, it is a powerful laboratory, because the strength of the dark-energy influence is controlled by a single parameter, alpha, which can be dialed up and down to see how the quantum behavior of the black hole responds.

The Brazilian duo attacked the problem with the Klein–Gordon equation, the relativistic wave equation for spinless particles, written in the curved geometry of the Kiselev metric. The metric carries a correction term that depends on alpha and on an equation-of-state parameter omega-zero, which characterizes how the pressure of the dark-energy fluid relates to its density. The researchers examined two representative choices: omega-zero equal to minus one-third and omega-zero equal to minus two-thirds, both lying in the range allowed for quintessence-like behavior. In each case, the geometry modifies the event horizon. For the first case the horizon sits at twice the black hole mass divided by one minus alpha, while for the second case two horizons emerge, an inner physical horizon and an outer quintessence-like horizon, provided alpha does not exceed one over eight times the mass. If alpha grows beyond that critical value, the horizons vanish entirely and the solution degenerates into a naked singularity, a configuration that general relativity’s cosmic censorship conjecture suggests should be forbidden.

Solving the Klein–Gordon equation in such a background is notoriously difficult, because the radial part of the wave equation becomes a complicated second-order differential equation with position-dependent coefficients. By separating the wave function into a radial piece, the familiar spherical harmonics, and a time-dependent oscillation, the authors reduced the problem to a radial equation whose solutions turn out to be special functions of a rare breed: Heun functions. For the omega-zero equals minus-one-third case, the radial wave function is expressed through the confluent Heun function, a power series whose coefficients obey a three-term recurrence relation. For the omega-zero equals minus-two-thirds case, the solution involves the generalized local Heun function, governed by a more demanding four-term recurrence. These functions are the curved-spacetime cousins of the Legendre and Bessel functions that dominate textbook quantum mechanics, and their appearance here is typical of wave equations in black hole backgrounds.

The real payoff comes when the infinite series is forced to terminate. Demanding that the wave function remain physically acceptable, finite at the horizon and well-behaved at infinity, truncates the Heun series into a polynomial, and that truncation can only happen at special, discrete values of the particle’s energy. These resonant frequencies form what physicists call the quasispectrum, and they are purely imaginary, meaning they describe modes that decay in time rather than oscillate forever. For the first model, the quasispectrum takes the form of the imaginary quantity proportional to n plus one, scaled by factors involving alpha and the black hole mass, and it reduces exactly to the known Schwarzschild result when alpha goes to zero. Intriguingly, the truncation conditions also quantize the orbital angular momentum quantum number l, which in general becomes complex, a hallmark of quasinormal-mode physics around black holes.

For the second model, the analysis produced a surprise: the truncation conditions can only be satisfied if the scalar particle is massless. The quasispectrum in that case again collapses to the Schwarzschild value of i times n plus one over four M in the limit of vanishing quintessence. The authors also verified that the radial wave function remains finite far from the black hole, where it behaves as an oscillating sine of a logarithmic phase with a calculable phase shift, a result that opens the door to studying the scattering of both massive and massless scalar fields off the Kiselev geometry using asymptotic methods developed by Abramov and collaborators.

With the exact wave functions in hand, the team then computed the quantity that connects this abstract mathematics to observable physics: the Hawking radiation itself. Following the standard tunneling-style analysis, they examined how outgoing modes behave just inside and just outside the event horizon, performed an analytic continuation around the pole at the horizon, and transformed to Eddington–Finkelstein coordinates, the coordinate system in which light rays travel along simple straight lines. Using Sannan’s formulation, they obtained the decay rate of particles escaping the horizon and, by matching to the Bose–Einstein distribution, extracted the Hawking temperature. For the omega-zero equals minus-one-third case, the temperature is proportional to one minus alpha squared divided by eight pi k-B times M, while for the omega-zero equals minus-two-thirds case it is proportional to alpha times the separation between the two horizons divided by four pi k-B times the inner horizon radius.

Two features of these temperatures stand out. First, both formulas agree exactly with the temperatures derived independently from the surface gravity approach, in which the temperature is proportional to the derivative of the metric function evaluated at the horizon. This cross-check is a strong consistency test, showing that the full quantum wave-function calculation reproduces the thermodynamic result obtained by purely geometric reasoning. Second, when alpha is set to zero, both cases recover the canonical Hawking temperature of a Schwarzschild black hole, one over eight pi k-B M, confirming that the entire dark-energy apparatus smoothly interpolates with known physics. Notably, in the second model the temperature is completely independent of the energy of the escaping particle, a clean separation between the thermodynamics of the hole and the quantum state of its radiation.

The most striking conclusion, however, concerns what dark energy does to the glow. Graphical analysis of the radiation spectrum shows that in both models, the larger the value of alpha, the weaker the Hawking radiation observed outside the event horizon. In the extreme limit of maximal quintessence influence in the first model, the radiation and the temperature both vanish entirely. The authors suggest this suppression may offer a tantalizing theoretical clue: since observations reveal very little matter streaming away from black hole horizons, a universe permeated by dark energy of the kind modeled here would naturally produce quieter black holes, consistent with what astronomers actually see. In other words, the cosmic acceleration that stretches the universe might simultaneously be muffling the quantum whispers of its most extreme objects.

The work, published as an open-access article in The European Physical Journal C and supported by the Brazilian funding agencies CNPq and FAPERJ, adds a precisely solved entry to the growing catalog of quantum fields in exotic spacetimes. Because quasinormal modes, black hole shadows, photon orbits, and thermodynamic stability are all known to shift in the presence of quintessence-like fluids, exact solutions of this kind provide the analytical scaffolding on which future phenomenological predictions can be built. As gravitational-wave detectors and black hole imaging campaigns grow ever more sensitive, distinguishing the fingerprint of dark energy in the behavior of black holes may move from a theorist’s daydream to a testable proposition, and calculations like this one map out exactly where that fingerprint should appear.

Subject of Research: Exact solutions of the Klein–Gordon equation for scalar particles in the spacetime of a Kiselev black hole surrounded by a quintessence-like anisotropic fluid

Article Title: Exact solution of the Klein–Gordon equation in a Kiselev black hole background

Article References: de Oliveira, M. D., & Schmidt, A. G. M. (2026). Exact solution of the Klein–Gordon equation in a Kiselev black hole background. The European Physical Journal C, 86(9), Article 1097. https://doi.org/10.1140/epjc/s10052-026-16325-4

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16325-4

Keywords: Kiselev black hole, Klein–Gordon equation, Hawking radiation, quintessence, dark energy, event horizon, Heun functions, quasinormal modes, general relativity, black hole thermodynamics, Schwarzschild limit, theoretical physics