The authors construct a model in which a piecewise mass function describes a Schwarzschild black hole surrounded by a finite, electromagnetically inert dark-matter distribution. The dark matter is treated as an effective addition to the black hole’s mass, with separate parameters for the mass ratio ΔM/M, the shell’s inner boundary rs and its thickness Δrs. The authors state that the framework can represent dark matter bound in a galactic halo as well as locally captured objects such as brown dwarfs, stellar-mass black holes, primordial black holes and non-radiating dust.

For the non-rotating case, the analysis identifies a critical shell mass beyond which the metric develops an additional pair of horizons. With rs/M = 2 and Δrs/M = 200, the critical values are r*/M = 152 and ΔM/M = 800/9, or about 88.9. Below that threshold, the original event horizon remains in place, although the photon sphere and shadow radius can change slightly. Above it, the outermost root becomes a new event horizon enclosing the central black hole and part of the dark-matter shell, while the photon-sphere and shadow scales increase substantially.

The authors then generalize the static model to an axisymmetric rotating spacetime using the Newman–Janis Algorithm. The resulting metric returns to the Kerr solution when ΔM = 0 and to the static configuration when the spin parameter is zero. They report that the metric satisfies the Einstein field equations and derive corresponding dark-matter density and pressure profiles. The surface stress-energy at the interfaces of the piecewise mass distribution is found to vanish.

In the rotating model, the critical dark-matter mass depends on spin and shell geometry. For an extremal black hole with a/M = 1 and Δrs/M = 200, the reported threshold is ΔM/M = 88.3 at r*/M = 150.5. Models above this threshold develop an enlarged outer horizon and additional roots in the stationary-limit equation, changing the locations and number of ergosphere layers. The reported changes concern the horizon, ergosphere, critical photon-orbit and shadow scales; the study does not establish a general expansion of every physical scale in the system.

Null geodesics are used to locate the unstable circular photon orbits that define the shadow boundary. In the high-mass regime, the critical photon orbits lie outside the dark-matter shell, where the mass function is constant and the geometry has the Kerr form with a much larger total mass. At a fixed spin, rotation and frame-dragging then have less influence relative to the enlarged shadow scale. The calculated distortion parameter consequently approaches zero, producing a nearly circular shadow even at high spin. This circularization is a model prediction rather than an observed black-hole feature.

The study also estimates the high-frequency absorption cross section from the average shadow radius and uses the surface gravity of the outer horizon to calculate Hawking temperature and thermal energy emission. At fixed dark-matter mass, the model predicts a lower emission rate as spin increases. Dark-matter masses below the critical threshold have a marginal effect, whereas above the threshold the calculated Hawking temperature becomes very small and the emission rate is strongly suppressed.

The authors compare predicted shadow diameters with EHT values of dsh/M = 11 ± 1.5 for M87* and dsh/M = 9.77 ± 0.67 for Sagittarius A*. For M87*, all theoretically allowed dark-matter masses in the plotted compatible region satisfy both the 1σ and 2σ observational intervals. For Sagittarius A*, the full allowed region satisfies the 2σ interval, while a 1σ upper limit emerges and decreases as spin rises. At a/M = 0.9, the reported constraint is ΔM/M ≤ 22.3853.

Because the high-mass branch predicts a shadow diameter beyond the EHT-compatible range, the model disfavors dark-matter configurations large enough to trigger the new outer-horizon branch near the two observed black holes. Within the study’s assumptions, the dark matter must either be absent from the immediate vicinity or have a localized mass below the relevant critical value. These limits are conditional on the piecewise, electromagnetically inert shell model, its chosen shell geometry and the adopted horizon and shadow calculations, rather than being general constraints on all possible dark-matter distributions.

Subject of Research: Theoretical modeling of rotating black holes surrounded by finite dark-matter shells and their observable shadows

Article Title: Shadow of rotating black hole surrounded by dark matter

Article References: Feng, H., Cai, Z., Yang, R.-J., & Zhang, J. (2026). Shadow of rotating black hole surrounded by dark matter. The European Physical Journal C, 86(10), Article 1180. https://doi.org/10.1140/epjc/s10052-026-16437-x

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16437-x

Keywords: black holes, dark matter, black hole shadows, Kerr spacetime, Event Horizon Telescope, null geodesics, Newman–Janis Algorithm, M87*, Sagittarius A*, Hawking radiation, Shadow, rotating

The researchers built their analysis around a family of models in which the interior of a collapsing object is described by a mass function that evolves in time. The key control parameter is a number, called n, which characterises how the mass is distributed near the centre of the collapsing body. When n is larger than two, the configuration is perfectly regular: all curvature invariants remain finite at the centre, and the matter resembles a smooth, quasi-de Sitter core, similar to the interior geometry proposed in many so-called regular black hole models. As the collapse proceeds, however, the null convergence condition of general relativity forces this parameter to decrease monotonically, driving the system through a sequence of increasingly singular configurations.

The crucial insight of the new work is that this descent does not pass through ordinary singularities but through what mathematicians call integrable singularities. In such a region, curvature invariants and the effective energy density do diverge as one approaches the centre, but their volume integrals remain finite. In other words, the total mass contained in any small ball around the centre shrinks to zero with the radius, so the divergence is in a precise sense harmless. The authors showed that this behaviour is exactly what one expects from quantum mechanics: in the hydrodynamic Madelung approximation, where a collective wavefunction describes the collapsing matter, the effective energy density is proportional to the squared modulus of the wavefunction, and the resulting mass function naturally produces integrable singular behaviour for a range of the parameter n between minus one and two.

As the collapse drives n downward, the interior geometry undergoes a dramatic and discontinuous change at n equals zero, an event the authors call Minkowski breaking. Before this transition, the spacetime possesses an inner Cauchy horizon in addition to the outer event horizon, a structure known to be fatally unstable through the process of mass inflation identified by Poisson and Israel in 1989. At the moment of Minkowski breaking, this inner horizon suddenly disappears: the slope of the mass function at the centre jumps from zero to infinity, and the metric function at the origin plunges from a finite value to minus infinity. The collapse continues past this point, with the configuration approaching the exact Schwarzschild geometry as n tends to minus one from above.

A central result of the study concerns whether the Schwarzschild singularity can actually be reached. By analysing the effective energy-momentum tensor of the collapsing interior, the authors found that the proper energy flux measured by a locally inertial observer diverges at the centre unless the parameter n stops changing exactly when it approaches minus one. This means that, within the semiclassical description, the point-like Schwarzschild singularity can only be approached asymptotically in time, never crossed in any finite sense. The total mass of the system, initially distributed throughout the interior, cannot continuously concentrate into a point; the transition to a true point-like source is inherently discontinuous, echoing the well-known fact that general relativity does not admit idealised point-particle sources.

To test whether the collapse might instead settle into a static end state before reaching the singularity, the team turned to the Raychaudhuri equation, the fundamental equation governing the expansion of a congruence of freely falling worldlines. Classically, the analysis confirmed the pessimistic expectation: for all configurations with n greater than minus one, any hypothetical static state is unstable. The tension in the interior matter, which supports the horizon, changes sign across the geometry, causing the congruence of infalling observers to expand in some regions and contract in others. After the Minkowski breaking, the rate of change of the expansion becomes negative everywhere, meaning that the classical dynamics inexorably drives the system toward the singular Schwarzschild geometry. Left to its own classical devices, collapse ends in catastrophe.

Quantum mechanics changes this verdict. Using the Madelung formulation, in which the Schrödinger equation is rewritten as a pair of hydrodynamic equations containing an additional quantum potential, the authors computed the quantum correction to the Raychaudhuri equation for the collapsing matter. The result is remarkable: the sign of this correction flips precisely at the Minkowski breaking. For configurations with n between zero and two, the quantum term diverges negatively near the centre, accelerating the collapse. But for n between minus one and zero, that is, after the inner horizon has vanished, the quantum contribution diverges positively at the centre, producing a strong repulsive effect that opposes further compression. The quantum potential, in effect, begins to fight the very collapse that created it.

The magnitude of this quantum repulsion is not negligible. The authors found that the sign change in the quantum term occurs at a radius of roughly half the horizon radius, for all post-Minkowski-breaking configurations, suggesting that the collapse might terminate in a quantum core of significant size rather than a point. The relative uncertainty in the position of the collapsing matter also grows without bound as n approaches minus one, indicating that the notion of a well-defined central point dissolves in the quantum regime. Meanwhile, near the horizon itself, the quantum term diverges negatively in a manner consistent with the existence of Hawking radiation, which would slowly evaporate the black hole and pull the horizon inward, an effect that becomes significant only at distances comparable to the Planck length from the horizon.

The overall picture that emerges is one of a cascade of discontinuous changes concentrated at the Minkowski breaking. At that single moment, the inner Cauchy horizon disappears, the central slope of the mass function jumps from zero to infinity, the metric function at the origin plunges from one to minus infinity, and the quantum contribution to the expansion of the collapsing fluid flips from attractive to repulsive. The authors emphasise that several of their conclusions are deliberately robust, in that they do not depend on the precise time evolution of the mass function, which remains unspecified in the model. Definitive statements about the final fate of the collapse, however, will require solving the full dynamical equations with the quantum term included, most likely through numerical simulations.

If the preliminary indications hold up, the implications would be profound. Black hole singularities, long regarded as the ultimate breakdown of general relativity, might never actually form; instead, collapse would end in a quantum-deformed core hidden behind the event horizon, with the Schwarzschild limit approached only asymptotically. This would resonate with a growing body of work on regular black holes, Planck stars and quantum gravitational collapse, and it offers a concrete, calculable framework in which the transition from classical to quantum behaviour during collapse can be tracked stage by stage. The Minkowski breaking, a seemingly abstract discontinuity in an interior mass function, may thus mark the exact moment when quantum physics takes over from Einstein’s classical theory and saves the universe from its most notorious infinity.

Subject of Research: Semiclassical and quantum analysis of gravitational collapse, integrable singularities and the possible avoidance of the Schwarzschild singularity

Article Title: On gravitational collapse and integrable singularities

Article References: Casadio, R., Giusti, A., Kamenshchik, A., & Ovalle, J. (2026). On gravitational collapse and integrable singularities. The European Physical Journal C, 86(8), Article 1018. https://doi.org/10.1140/epjc/s10052-026-16258-y

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16258-y

Keywords: black holes, gravitational collapse, singularity, Schwarzschild geometry, quantum potential, Raychaudhuri equation, Madelung approximation, Cauchy horizon, mass inflation, Hawking radiation, general relativity, integrable singularity

The study, published in The European Physical Journal C, examines asymptotically de Sitter black holes in generalized Proca theory, a framework in which gravity is supplemented by a massive vector field with derivative self-interactions. First formulated in its complete form by Lavinia Heisenberg, the theory is carefully constructed so that only the three physical polarizations of the vector propagate, avoiding the ghostlike extra degrees of freedom that plague many modified-gravity models. Crucially, the vector sector can carry what physicists call primary hair: independent charges that make the black hole genuinely different from anything in Einstein’s vacuum solutions. In the branch of solutions studied here, the same Proca sector that supports this hair also generates the effective de Sitter expansion of the universe at large distances, so the cosmological constant is emergent rather than fundamental.

That combination makes these black holes an ideal laboratory for scattering problems. In an asymptotically de Sitter spacetime, the static region is bounded on one side by the event horizon and on the other by the cosmological horizon, the boundary beyond which the accelerated expansion of space sweeps everything away. Waves scattered in this patch must tunnel through a single potential barrier separating the two horizons, and the transmission probability through that barrier is precisely the grey-body factor. The authors, led by Bekir Can Lütfüoğlu of the University of Hradec Králové together with Javlon Rayimbaev, Nuriddin Kurbonov, Sardor Murodov and Faisal Javed, computed these factors for three distinct test fields: a massive scalar, an electromagnetic field, and a massless Dirac field.

The technical heart of the analysis lies in reducing each field’s dynamics to a Schrödinger-type equation with a spin-dependent effective potential. For the massive scalar, the potential contains both a curvature contribution proportional to the derivative of the metric function and an explicit mass term, so the barrier is both raised and broadened as the mass grows. The electromagnetic barrier is purely centrifugal, scaling with the angular momentum quantum number, while the Dirac problem is governed by supersymmetric partner potentials that are isospectral under standard black-hole boundary conditions. In a representative configuration with horizons located at roughly 2.25 and 17.8 in units of the mass parameter, the scalar barrier peaks at about 0.167, more than three times the electromagnetic maximum of roughly 0.055 and nearly five times the Dirac peak of about 0.034. As the scalar mass increases from 0.1 to 1.0, the barrier maximum moves outward from radius 3.23 to 6.93, directly shifting the onset of efficient transmission to higher frequencies.

To obtain the transmission curves, the team employed two complementary methods. The first is a sixth-order WKB approximation, a semiclassical technique that evaluates the local shape of the potential at its maximum and incorporates corrections through sixth order in derivatives of the potential with respect to the tortoise coordinate. The second exploits a recently developed correspondence between grey-body factors and quasinormal modes, the characteristic ringing frequencies of the black hole. In the eikonal regime, the WKB parameter controlling transmission can be reconstructed from the real and imaginary parts of the fundamental quasinormal frequency, with beyond-eikonal corrections refining the estimate. The two approaches mostly agree, with visible differences concentrated near the transition region where the grey-body factor rises from strong reflection to strong transmission, exactly where the barrier-top approximation is most delicate. The authors note an important caveat: the correspondence assumes a single smooth dominant potential maximum, and it fails if the potential develops double-well or double-barrier topology, as can happen in strongly deformed higher-curvature models.

The parameter scans reveal how sensitively the transmission depends on the Proca couplings. In the electromagnetic sector, increasing the coupling alpha hardens the barrier, delaying the rise of the transmission probability to higher frequencies, and the coupling beta produces an even clearer monotonic suppression. Increasing the coupling c1, by contrast, softens the barrier and shifts the transmission window to lower frequencies, reversing the ordering. In the Dirac sector, increasing the primary-hair parameter Q progressively suppresses the grey-body factor up to higher frequencies, with the neutral case being the most transmissive. These trends translate directly into effective absorption cross-sections, which the authors define for the static patch as a measure of how efficiently radiation sent inward from the cosmological horizon is captured by the event horizon. In a striking comparison between a neutral black hole and a near-extreme charged configuration with Q equal to 1.02, the larger hair value visibly reduces both the dominant low-multipole peaks and the high-frequency level of the total cross-section, showing that primary hair does not merely move the horizons but genuinely reduces the hole’s capture efficiency.

The thermodynamic analysis adds another layer of richness. Because de Sitter spacetime lacks an asymptotically flat region in which to normalize time, and because the cosmological horizon carries its own temperature, the team distinguished three temperatures: the bare black-hole temperature from the slope of the metric function at the event horizon, the Bousso-Hawking normalized temperature evaluated at the radius where black-hole attraction and cosmological repulsion balance, and the cosmological-horizon temperature. For the neutral reference configuration, these take the values 0.0337, 0.0455 and 0.00728 respectively. Increasing the hair parameter Q drives the event horizon inward and lowers both black-hole temperatures while leaving the cosmological horizon almost untouched; near the high-charge edge, the bare and normalized temperatures drop to about 0.00845 and 0.0114, explaining why the corresponding emission and absorption channels are so strongly suppressed. The couplings leave their fingerprints here too: increasing lambda or c1 enlarges the effective de Sitter contribution and raises the cosmological-horizon temperature even as the bare event-horizon temperature falls.

Why does this matter beyond the mathematics? Grey-body factors are as important as temperature itself in determining the observable Hawking-radiation spectrum, because the flux reaching a distant observer is the product of thermal occupation and barrier transmission. In de Sitter black holes the problem is considerably more involved than in asymptotically flat spacetimes, since two horizons with different temperatures compete. The present work shows that in generalized Proca theory the entire filtering apparatus is tied to the vector couplings, meaning that any future observational handle on black-hole emission spectra could in principle constrain the same physics that drives cosmic acceleration in these models. The quasinormal spectra of this same geometry had already been shown to be highly sensitive to the hair parameter and to the couplings controlling the effective de Sitter radius, and the new grey-body analysis confirms and extends that picture from the resonance side to the scattering side.

The authors point toward a natural next step: a full direct-scattering computation of the transmission probabilities and Hawking-radiation intensities across the entire parameter space, followed by complete Bose-Einstein and Fermi-Dirac emission sums for all three spin sectors. Such calculations would complete the portrait of these vector-hair black holes as thermal objects and sharpen the interplay between black-hole spectroscopy, modified gravity, and the accelerated expansion of the universe. For now, the message is clear: in a universe whose cosmological constant is written by a massive vector field, even the darkness of a black hole carries the signature of its hair.

Subject of Research: Grey-body factors and thermodynamics of asymptotically de Sitter black holes in generalized Proca theory

Article Title: Grey-body factors and thermodynamics of asymptotically de Sitter black holes in generalized Proca theory

Article References: Lütfüoğlu, B. C., Rayimbaev, J., Kurbonov, N., Murodov, S., & Javed, F. (2026). Grey-body factors and thermodynamics of asymptotically de Sitter black holes in generalized Proca theory. The European Physical Journal C, 86(9), Article 1027. https://doi.org/10.1140/epjc/s10052-026-16289-5

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16289-5

Keywords: black holes, grey-body factors, generalized Proca theory, de Sitter spacetime, Hawking radiation, quasinormal modes, WKB approximation, absorption cross-sections, black hole thermodynamics, primary hair, cosmological horizon, modified gravity

The object at the center of the study is the Canfora–Maeda solution, an exact analytic black hole in a theory where Einstein gravity is coupled to the Skyrme model, a nonlinear field theory of an SU(2)-valued field that supports topologically stable solitons identified with baryons. Exact analytic hairy black holes are exceedingly rare because the Skyrme equations are strongly nonlinear and involve higher-derivative terms, so most known configurations must be found numerically. The Canfora–Maeda solution is characterized by two independent couplings: K, which controls a quadratic sigma-model term and produces a solid-angle-deficit-like deformation of the geometry, and lambda, which governs the quartic Skyrme interaction and enters the metric function with the same positive sign as the charge term of a Reissner–Nordström black hole.

A crucial foundation for the new work is a recently formulated extended first law of black hole thermodynamics, in which the Skyrme couplings K and lambda are promoted to extensive thermodynamic variables with their own conjugate potentials, much as the cosmological constant is treated as pressure in the black hole chemistry paradigm. Building on this framework, the authors constructed the Ruppeiner thermodynamic geometry of the black hole. Ruppeiner geometry builds a Riemannian metric on the space of equilibrium states from the Hessian of the thermodynamic potential, and its scalar curvature is widely interpreted as a measure of the type and strength of interactions among the microscopic constituents of a thermodynamic system. For ordinary fluids, the sign of the curvature distinguishes attractive from repulsive interactions, while vanishing curvature signals an ideal gas.

The result is remarkably compact. The scalar curvature takes the closed form R = 1/(S_ext − S), where S_ext is the extremal entropy at which the Hawking temperature vanishes. Three consequences follow immediately. The curvature diverges in the extremal, zero-temperature limit, where the small-black-hole branch terminates, rather than at the second-order phase transition diagnosed by a divergence of the heat capacity. At that transition entropy, which equals three times the extremal entropy, the curvature remains perfectly finite because the thermodynamic metric itself stays regular there. And throughout the entire physically admissible region, the curvature keeps a single sign, with its magnitude decreasing monotonically as the black hole grows larger.

That absence of a sign change is a finding in its own right, not a technical caveat. In the influential picture developed for charged anti-de Sitter black holes, the Ruppeiner curvature changes sign along the small-black-hole branch, signaling a crossover between repulsive and attractive microscopic interactions that mirrors the Van der Waals-like phase behavior of those systems. The Einstein–Skyrme black hole shows no such crossover: its effective microscopic interactions retain one dominant character across the whole physical domain, strengthening as extremality is approached and weakening toward the large, nearly non-interacting regime. The authors are careful to note that the sign itself is convention dependent, but the location of the divergence, the regularity at the phase transition and the absence of any sign flip are robust, convention-independent statements.

The team also performed a Joule–Thomson-like isenthalpic expansion analysis, treating the quartic coupling lambda as a pressure-like variable through the work term it contributes to the first law. In ordinary fluids, the Joule–Thomson coefficient determines whether a gas cools or heats upon expansion, and an inversion temperature separates the two regimes. For charged AdS black holes, such an inversion curve is a well-known feature. Here, the authors proved analytically that the Joule–Thomson-like coefficient is strictly negative across the entire physical parameter space: the black hole always cools as lambda increases at fixed mass, the equation setting the coefficient to zero has no solution, and no inversion temperature exists. The entire physical region forms a single cooling regime, a distinctive signature of the Skyrme matter sector, since increasing lambda enhances the effective charge-like contribution in the metric and systematically lowers the Hawking temperature.

Turning to dynamics, the researchers derived the effective potential governing massless scalar field perturbations on the Einstein–Skyrme background. The potential takes a Schrödinger-like form in the tortoise coordinate, and the two Skyrme couplings enter through distinct channels. The solid-angle-deficit coupling K dominates: switching it on strongly lowers the potential barrier, with the peak for the l = 2 scalar mode dropping from roughly 0.247 in the Schwarzschild case to about 0.058 at K = 0.015, while the horizon grows outward. By contrast, at fixed K the quartic coupling lambda produces only minor changes, raising the peak by little more than one percent over the range studied, because its contribution enters the metric function repulsively, exactly like an electric charge.

Using leading-order WKB estimates, calibrated against the known Schwarzschild quasinormal frequency, the authors showed that the fundamental oscillation frequencies decrease markedly as K increases, reflecting the lowered and broadened potential barrier, while the dependence on lambda at fixed K is weak. This dominance of K over lambda is itself a nontrivial result: although lambda enters the lapse function exactly like a Reissner–Nordström charge, the Einstein–Skyrme black hole does not behave as a naive charged-black-hole analogue. Its perturbative spectroscopy is governed primarily by the global conical structure of the geometry rather than by the effective charge, a distinction that matters for any attempt to read black hole parameters off gravitational-wave ringdown signals.

Finally, the study established rigorous closed-form lower bounds on the greybody factors, the transmission coefficients that quantify how the curved-spacetime potential barrier filters Hawking radiation. Using Visser’s general bounding method, the tortoise integral of the potential can be evaluated analytically, and the would-be deficit prefactor cancels between the potential and the tortoise measure, so the bound remains finite despite the conical asymptotics. Increasing K raises the bound substantially, because the larger horizon lowers the centrifugal barrier: at a representative frequency, the bound climbs from about 0.01 in the Schwarzschild case to roughly 0.10 at K = 0.015 and 0.40 at K = 0.025. Increasing lambda at fixed K mildly suppresses low-frequency transmission, by only a few percent, while at high frequencies the bound approaches unity for all couplings, as expected for short-wavelength radiation.

Together, these results provide the first comprehensive portrait of the thermodynamic microstructure and perturbative response of one of the very few known exact analytic hairy black holes. They open several avenues for future work, including extending the Ruppeiner analysis to a full non-degenerate thermodynamic space, computing high-precision quasinormal frequencies with higher-order WKB or spectral methods, incorporating electromagnetic and gravitational perturbations, and exploring whether the microscopic interactions identified here connect to the topological properties of the underlying Skyrmion field. As gravitational-wave astronomy matures, understanding how matter fields such as the Skyrme sector reshape the spectra and radiation of black holes may prove essential for distinguishing genuine signatures of hairy black holes from the predictions of classical general relativity.

Subject of Research: Thermodynamic geometry and perturbative spectroscopy of the Einstein–Skyrme black hole

Article Title: Ruppeiner thermodynamic geometry, microstructure, quasinormal modes and greybody factors of the Einstein–Skyrme black hole

Article References: Khasanov, S., Safarov, A., Yang, B., & Tedila, H. M. (2026). Ruppeiner thermodynamic geometry, microstructure, quasinormal modes and greybody factors of the Einstein–Skyrme black hole. The European Physical Journal C, 86(9), Article 1080. https://doi.org/10.1140/epjc/s10052-026-16230-w

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16230-w

Keywords: black hole thermodynamics, Einstein–Skyrme theory, Ruppeiner geometry, quasinormal modes, greybody factors, microstructure, Joule–Thomson expansion, Skyrmions, hairy black holes, phase transitions, Hawking radiation, general relativity

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

The physics behind the proposal traces back to two landmark discoveries. In 1976, William Unruh showed that an observer accelerating through the vacuum of flat spacetime perceives that vacuum as a warm thermal bath, with a temperature proportional to the acceleration. A year earlier, Stephen Hawking had demonstrated that black holes radiate thermally. Shortly afterward, Gibbons and Hawking found that observers in an expanding de Sitter universe—our best mathematical description of a dark-energy-dominated cosmos—also perceive a thermal bath at a temperature set by the cosmological constant. These effects share a common signature: the correlation functions of the quantum field obey the Kubo–Martin–Schwinger (KMS) condition, the mathematical fingerprint of thermal equilibrium. Detecting them directly, however, is extraordinarily difficult because the temperatures involved are vanishingly small for any achievable acceleration.

The new work turns the problem around. Instead of measuring particles directly, the researchers treat the accelerated detector as a quantum battery that can be charged by an external classical driving field while simultaneously interacting with the fluctuating vacuum. A quantum battery stores energy in quantum states, and the ergotropy quantifies how much of that stored energy is genuinely extractable as useful work. For a two-level battery with transition frequency, the ergotropy is determined by the Bloch vector describing the battery’s quantum state: specifically, it equals half of the sum of the vector’s magnitude and its third component. As the battery evolves, vacuum fluctuations in the surrounding spacetime cause dissipation and decoherence, gradually reshaping the state and therefore the extractable work.

Using the open quantum system approach, the team derived a Kossakowski–Lindblad master equation in the weak-coupling limit, valid under the Born–Markov approximation. The coefficients of this equation are fixed by the response function of the detector—the Fourier transform of the field’s Wightman function along the accelerated trajectory. In de Sitter spacetime, this response takes a purely thermal form with a temperature that combines the Unruh contribution from acceleration and the Gibbons–Hawking contribution from curvature: the effective temperature is proportional to the square root of the acceleration squared plus the curvature scale squared. In anti-de Sitter spacetime, by contrast, thermality appears only when the acceleration exceeds the curvature scale—a supercritical condition—because the negative curvature prevents the formation of the relevant horizon below that threshold.

The most striking result concerns the asymptotic behavior. After a long charging time, the battery thermalizes with the field, and its ergotropy settles at a steady value equal to one half of the hyperbolic tangent of the transition frequency divided by twice the effective temperature. Because this expression depends only on the KMS temperature, it is completely independent of the spacetime dimension and, in anti-de Sitter, of the boundary conditions imposed on the quantum field. The hotter the effective bath—whether from higher acceleration or stronger curvature—the smaller the steady ergotropy. Crucially, the same formula unifies the two spacetimes: once the acceleration greatly exceeds the curvature scale, de Sitter and anti-de Sitter batteries converge to identical steady behavior, providing a single witness for Unruh–Hawking thermality across both geometries.

The transient dynamics, however, tell a richer story. In de Sitter spacetime, the ergotropy exhibits pronounced oscillations during its evolution at low accelerations, ringing before settling to its steady value, while at high accelerations it relaxes rapidly and smoothly to equilibrium. In anti-de Sitter spacetime, the picture changes dramatically with the boundary conditions applied to the scalar field at the spacetime boundary. The team examined three cases—Dirichlet, transparent, and Neumann—and found that the oscillatory behavior of the ergotropy is strongest under Dirichlet conditions, weaker for transparent boundaries, and weakest for Neumann boundaries. Notably, choosing the right boundary condition can improve the energy storage of the moving battery, and the boundary-induced corrections become less important as the acceleration grows.

Dimensionality adds another twist. Extending the calculation to six-dimensional anti-de Sitter spacetime, the researchers showed that vacuum fluctuations can modestly amplify the ergotropy during the initial charging stage, because the detector’s response rate increases with the number of spacetime dimensions. A short-time expansion of the ergotropy reveals that it grows quadratically at early times, with a coefficient proportional to the response function, which is larger in six dimensions than in four. Higher dimensions also accelerate thermalization, stabilizing the oscillations more quickly. Yet the asymptotic ergotropy remains the same regardless of dimension, again dictated solely by acceleration and curvature through the KMS condition. The team also notes that in odd-dimensional anti-de Sitter spacetimes, a curious inversion of quantum statistics occurs—bosonic fields behave as if fermionic—but this anomaly does not destroy the underlying thermal nature witnessed by the battery.

The sensitivity of the probe differs between the two geometries in an instructive way. Comparing the derivatives of the steady ergotropy with respect to acceleration, the researchers found that the anti-de Sitter battery responds more sharply to changes in acceleration than its de Sitter counterpart, and this heightened sensitivity persists regardless of boundary conditions. At very short charging times, by contrast, the acceleration has almost no effect, because the thermal response of the vacuum is still too weak and the extractable work is dominated by the external driving field. This separation of timescales suggests a practical strategy: short-time measurements characterize the charging protocol and the geometry-dependent corrections, while long-time measurements isolate the universal thermal fingerprint encoded in the KMS condition.

Beyond its conceptual appeal, the framework connects several active research frontiers: quantum thermodynamics, relativistic quantum information, and the physics of open quantum systems in curved spacetime. Previous studies had used uncertainty relations, quantum coherence, geometric phases, and Fisher information as probes of the Unruh effect in flat spacetime, and recent work had explored quantum batteries near black holes. The present study extends this program to the two maximally symmetric curved spacetimes of constant positive and negative curvature, offering an operational interpretation of thermality in terms of energy transfer. While an actual laboratory demonstration remains distant—the Unruh temperatures involved are minuscule—the model provides theorists with a concrete, work-based observable that distills the essence of horizon thermodynamics. In the researchers’ view, the relativistic quantum battery is not merely an energy storage device but a functional probe of vacuum fluctuations, opening a new window onto the thermal structure of spacetime itself.

Subject of Research: Quantum thermodynamics of relativistic quantum batteries as probes of Unruh–Hawking thermality in curved spacetimes

Article Title: Quantum thermodynamics of ergotopy for a relativistic battery as a witness to Unruh–Hawking thermality in curved (A)dS spacetimes

Article References: Hao, X., Gan, T.-F., Wu, C.-T., Ren, T.-X., Zhang, W.-W., & Wu, Y.-Z. (2026). Quantum thermodynamics of ergotopy for a relativistic battery as a witness to Unruh–Hawking thermality in curved (A)dS spacetimes. The European Physical Journal C, 86(10), Article 1137. https://doi.org/10.1140/epjc/s10052-026-16400-w

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16400-w

Keywords: quantum battery, ergotropy, Unruh effect, Hawking radiation, de Sitter spacetime, anti-de Sitter spacetime, Unruh–DeWitt detector, KMS condition, quantum thermodynamics, curved spacetime, open quantum systems, relativistic quantum information

The starting point of the analysis is a deceptively simple idea that has haunted general relativity for half a century: the quantum vacuum is not truly empty. Quantum fields fluctuate even in the absence of matter, and those fluctuations carry energy and stress. When a black hole is present, the curved spacetime forces these quantum fields into a strained state, and their averaged energy density, described by the expectation value of the quantum stress-energy tensor, feeds back on the geometry itself. In the semi-classical framework adopted by the authors, Einstein’s field equations are rewritten so that the right-hand side contains this quantum expectation value rather than zero. Solving the modified equations for a spherically symmetric, non-rotating black hole produces a corrected metric that differs from Karl Schwarzschild’s 1916 solution by terms of order the Planck length squared, the scale at which quantum gravity is expected to dominate.

The structure of the quantum stress-energy tensor is central to the result. The authors decompose it into two physically distinct pieces: a trace part governed by the conformal anomaly, a genuine quantum effect in which classically traceless fields acquire a nonzero trace in curved spacetime, and a regular part that itself splits into a vacuum polarization contribution and a component encoding the outgoing flux of Hawking radiation. For massless fields of spin zero, one and two, the relevant numerical coefficients can be computed, yielding a corrected metric function with three small correction constants. Remarkably, the resulting geometry no longer has a single event horizon. The classical horizon splits into an outer horizon, which reduces to the familiar Schwarzschild radius when quantum effects are negligible, and a purely quantum inner horizon that has no classical counterpart.

This two-horizon structure transforms the thermodynamics of the hole. In classical theory, the Hawking temperature of a Schwarzschild black hole is inversely proportional to its mass, so as radiation drains mass away, the temperature climbs without bound and the evaporation accelerates in a runaway process. The quantum-corrected calculation tells a different story. Using the surface gravity formalism, the authors derive a modified temperature that reproduces the classical value for large masses but contains a correction term proportional to the square of the Planck mass divided by the mass cubed. As the hole shrinks, the temperature rises only until the mass reaches a critical value of roughly 0.7 Planck masses, where it attains a finite maximum. Beyond that point the temperature reverses course and falls, dropping to zero at a minimum mass of about 0.4 Planck masses. Evaporation simply stops.

The timing of this final act is extraordinary. Applying the Stefan-Boltzmann law to the corrected temperature shows that in the late stages the evaporation timescale scales as the inverse cube of the distance to the minimum mass, meaning the approach to the remnant becomes ever more sluggish. For a distant observer, the black hole appears to freeze just above the Planck scale, radiating ever more slowly as its temperature sinks. The authors also track the horizons themselves: as mass is lost, the outer and inner horizons converge, squeezing the trapped region between them. When the temperature reaches zero, the two horizons merge into a single extremal surface and the trapped region vanishes entirely. Unlike the extremal remnant of a charged Reissner-Nordstrom black hole, this object carries no electric charge, making it a genuinely stable endpoint of collapse.

The entropy of the black hole acquires an equally significant correction. Integrating the first law of black hole thermodynamics with the corrected temperature yields an entropy that consists of the familiar Bekenstein-Hawking area law, proportional to the horizon area divided by four times Newton’s constant, plus a logarithmic term whose coefficient depends on the quantum fields present. This logarithmic correction is a recurring motif across quantum gravity research, appearing independently in loop quantum gravity, Euclidean quantum gravity, and generalized uncertainty principle models, and its natural emergence here from vacuum backreaction strengthens the case that it reflects a universal feature of quantum spacetime. Intriguingly, the entropy now depends not only on the mass but on the field content of the theory, which the authors interpret as a form of quantum hair. Above the Planck mass the corrected entropy exceeds the classical value, while in the Planckian regime it falls below it, signaling a reduction in the number of available microscopic states.

The heat capacity reveals a dramatic phase transition at the critical mass. Classical Schwarzschild black holes have negative heat capacity: they grow hotter as they lose energy, making them thermodynamically unstable. In the quantum-corrected picture, once the mass drops below the critical value the heat capacity flips sign and becomes positive, meaning the hole can now exchange heat stably with its environment and settle into thermal equilibrium. The black hole thus evolves from an unstable radiator into a stable, cold object, a qualitative transformation with no analogue in the classical theory. The authors emphasize that these features are not model-specific quirks; similar maximal temperatures, slowed evaporation, and remnant formation appear in loop-inspired Planck star models, asymptotically safe gravity, noncommutative geometries, and generalized uncertainty principle frameworks, suggesting the results capture generic consequences of semi-classical backreaction.

To verify the thermodynamic picture independently, the authors turn to the Parikh-Wilczek tunneling method, which treats Hawking radiation as a quantum tunneling process in which a particle emerges from just inside the horizon while energy conservation shrinks the hole by the emitted energy. Working in Painleve-Gullstrand coordinates that are regular at the horizon, they compute the imaginary part of the tunneling action and recover an emission probability whose exponent matches the corrected entropy exactly. In the limit where the emitted energy is small compared with the black hole mass, the radiation reduces to a thermal Boltzmann spectrum at precisely the corrected temperature. But the higher-order terms in the emitted energy, which become important near the Planck scale, render the spectrum distinctly non-thermal, a deviation that could have profound consequences for how information is encoded in the radiation.

That last point touches the deepest open question in the field: the black hole information paradox. Because evaporation in this model terminates at a finite-mass remnant rather than ending in complete disappearance, part of the information that fell in could conceivably remain stored in the remnant configuration, offering a possible escape route from Hawking’s original argument that information is destroyed. The authors are careful not to claim a resolution, noting that a full analysis would require going beyond the semi-classical approximation, but they point out that the modified causal structure could alter the Page time and the entanglement structure of the radiation, opening a connection to recent island and quantum extremal surface constructions. For astrophysical black holes, the corrections remain hopelessly small, with the quantum parameter for a solar-mass hole sitting near ten to the minus seventy-sixth power, so the photon orbits and shadows we observe today are essentially classical. The drama, if these calculations are right, plays out only in the final Planck-scale moments, when a dying black hole cools into an eternal, information-bearing ember rather than vanishing from the universe.

Subject of Research: Semi-classical Hawking radiation and evaporation endpoint of a quantum-corrected Schwarzschild black hole

Article Title: Hawking radiation from a semi-classical Schwarzschild black hole

Article References: Shafiee, M., & Sheykhi, A. (2026). Hawking radiation from a semi-classical Schwarzschild black hole. The European Physical Journal C, 86(9), Article 1096. https://doi.org/10.1140/epjc/s10052-026-16370-z

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16370-z

Keywords: Hawking radiation, black hole thermodynamics, Schwarzschild black hole, quantum corrections, conformal anomaly, vacuum polarization, black hole remnants, logarithmic entropy correction, Parikh-Wilczek tunneling, event horizon, Planck scale, information paradox

For decades, textbooks have located Hawking radiation in an infinitesimally thin layer hugging the event horizon. But that picture has been challenged. Building on arguments by physicist Steven Giddings, the quantum atmosphere proposal holds that the effective source of Hawking radiation extends to a region comparable in size to the horizon radius itself, roughly Δr on the order of the horizon scale r_H. The idea was originally motivated by estimates based on the Stefan-Boltzmann law and the wavelengths of the emitted quanta, and it has since gained support from independent analyses of the stress-energy tensor, the gravitational Schwinger effect, quantum correlations across the horizon, and the thermal character of the radiation. If the atmosphere is real, then the region between roughly 1.4 and 1.5 horizon radii should leave measurable traces in any quantum system probed there.

Zhang and colleagues tested this idea using a tripartite mixed state—a deliberate departure from most earlier work, which focused on pairs of parties or on pure multipartite states. Their starting state blends a generalized GHZ state, the canonical example of genuine three-party entanglement, with a controlled amount of white noise. The state takes the form ρ_ABC = p|Ψ_GHZ⟩⟨Ψ_GHZ| + (1−p)I_8/8, where the mixing parameter p tunes how much genuine quantum structure survives and the state parameter α controls the coherent superposition α|000⟩ + √(1−α²)|111⟩. By varying p, the researchers could systematically examine how noise and gravity jointly reshape multipartite quantum information.

The gravitational stage for this drama is the Schwarzschild spacetime of a non-rotating black hole. The team modeled massless Dirac fields, solving the curved-spacetime Dirac equation near the horizon and constructing the Kruskal modes from the Schwarzschild modes using the Damour-Ruffini analytic continuation method. The resulting Bogoliubov transformation entangles modes outside the horizon with partners trapped inside: the Kruskal vacuum becomes a fermionic two-mode squeezed state with coefficients μ and ν governed by the local temperature. Crucially, rather than using the standard Hawking temperature, the authors substituted the Hartle-Hawking local temperature T_HH, which depends on radial distance, vanishes exactly at the horizon, peaks at a finite radius outside it, and relaxes to the asymptotic Hawking temperature far away.

That local temperature carries a parameter D_HH associated with the stress tensor in the Hartle-Hawking vacuum. Because the Hartle-Hawking boundary conditions alone do not fix this constant, the researchers required the local temperature to remain real everywhere outside the horizon, which imposes D_HH ≥ D_c ≈ 23.03. For such values, the peak of the local temperature sits in the interval 1.43 r_H ≲ r_peak < 1.5 r_H. The whole analysis then hinges on a simple question: do the quantum resources of a tripartite mixed state track this atmospheric peak?

The first quantity they examined, quantum state texture, is a recently proposed structural measure that treats the density matrix as a geometric landscape and quantifies its unevenness relative to a featureless reference state. It is defined as R(ρ) = −ln⟨f_1|ρ|f_1⟩, where |f_1⟩ is the uniform superposition over the computational basis. Because it can be extracted from a single projection probability, it is experimentally meaningful, and it is not equivalent to conventional coherence measures. When the team computed the texture of the physically accessible state—Alice in flat space, with Bob and Charlie hovering near the horizon—they found it rises and then falls with radial distance, exhibiting a clear local extremum. The inaccessible states, built from the modes hidden behind the horizon, show the opposite trend, revealing a genuine redistribution of quantum structure between the two sectors.

The radial position of that extremum is the striking part. For every fixed value of D_HH, the texture extrema of all accessible and inaccessible reduced states occur at the same radius, and as D_HH grows from the critical value 23.03 up to 100, the extremum shifts outward from r/r_h ≈ 1.4322 toward roughly 1.48—always within the same finite near-horizon window in which the local Hawking temperature peaks. In other words, the structural signature of the tripartite mixed state is written exactly where the atmosphere is hottest. The authors interpret this as evidence that quantum state texture serves as a genuine quantum signature of near-horizon Hawking radiation, not merely an artifact of the chosen measure.

Entanglement tells a subtler story. Using the genuine multipartite concurrence, valid for the X-shaped density matrices that arise here, the team showed that the mixing parameter p dominates the survival of genuine three-party entanglement: below moderate p, white noise erases it almost entirely, while as p grows the entanglement re-emerges and spreads. More importantly, the local Hawking effect does not simply destroy this resource. In the physically accessible state, the concurrence dips in a finite region outside the horizon and then recovers, whereas in the inaccessible states it shows a localized peak near the horizon before decaying. Entanglement is thus redistributed between sectors rather than annihilated, and the accessible state remains its main carrier. The team also found that the Hawking effect shifts the optimal value of the state parameter α away from the balanced superposition, so the initial GHZ structure that maximizes entanglement is itself gravity-dependent.

Tripartite nonlocality proved far more fragile. Detected through violation of the Svetlichny inequality, which rules out all bipartite local hidden-variable models, genuine three-party nonlocality survives only when the local Hawking effect is weak. For the accessible state, the Svetlichny parameter stays above the threshold S = 4 for small D_HH, but as D_HH increases at fixed radius it falls rapidly below the threshold and the nonlocality disappears. The inaccessible states never violate the inequality at all, even though they can retain nonzero genuine multipartite entanglement. This hierarchy—texture and entanglement persist while Svetlichny nonlocality is extinguished—establishes a strict ordering in the robustness of tripartite quantum resources under gravitational radiation.

The unifying conclusion is that three conceptually distinct measures—structural texture, entanglement redistribution, and Svetlichny nonlocality—respond most sensitively in the same radial window between roughly 1.432 and 1.5 horizon radii, with their extrema marching outward together as the local Hawking temperature grows. Taken together, the results provide a unified characterization of how the quantum atmosphere restructures density matrices, redistributes multipartite entanglement, and suppresses the strongest forms of nonclassical correlation. They also add independent, quantum-information-based support to the idea that Hawking radiation is not born at the horizon itself but in a finite atmospheric shell. The authors note that their single-mode approximation, while analytically tractable and widely used, leaves a full multi-mode wave-packet treatment for future work—though previous studies suggest the qualitative behavior reported here should survive.

Subject of Research: Tripartite mixed-state quantum correlations in the black hole quantum atmosphere

Article Title: Quantum correlations of tripartite mixed states in the black hole quantum atmosphere

Article References: Zhang, A., Zheng, Y., Huang, X., & Zhang, T. (2026). Quantum correlations of tripartite mixed states in the black hole quantum atmosphere. The European Physical Journal C, 86(9), Article 1105. https://doi.org/10.1140/epjc/s10052-026-16377-6

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16377-6

Keywords: black hole quantum atmosphere, Hawking radiation, quantum entanglement, tripartite nonlocality, Svetlichny inequality, quantum state texture, GHZ state, Dirac fields, Schwarzschild spacetime, Bogoliubov transformation, relativistic quantum information, Hartle-Hawking temperature

To appreciate the significance of the work, it helps to understand what quantum discord actually measures. When two quantum systems share correlations, the most famous variety is entanglement, the spooky linkage that Einstein derided as action at a distance. But discord captures something broader: it quantifies the minimum loss of quantum mutual information that occurs when one part of a correlated pair is measured. Crucially, discord can persist even in separable states that carry no entanglement at all, and it is famously robust against environmental decoherence. That resilience makes it an ideal probe for noisy, extreme environments, and few environments are noisier or more extreme than the hot, radiating halo surrounding a black hole. Unlike binary measures such as Bell inequality violations, discord is a continuous quantity, allowing it to register fine-grained, non-monotonic changes in correlation strength that cruder measures would miss entirely.

The theoretical stage for the study is the quantum atmosphere, a concept sharpened by Steven Giddings in 2016. For decades, textbooks depicted Hawking radiation as arising from quantum excitations hugging the event horizon, in a thin shell where the radial distance above the horizon is tiny compared with the horizon radius itself. Giddings’ analysis of total emission rates and stress tensors overturned that picture: the radiation effectively originates from a region extending roughly one horizon radius above the horizon, a thick shell he called the quantum atmosphere. This revised geography matters enormously for quantum information studies, because any observer or quantum system falling toward a black hole passes through this extended region and experiences its local conditions before ever reaching the horizon. Li and Huang set out to chart exactly how quantum correlations evolve as a system traverses this zone.

The authors’ mathematical machinery begins with the Schwarzschild metric, the standard description of a non-rotating, uncharged black hole, and the curved-spacetime Dirac equation governing massless fermions. Solving this equation yields two sets of positive-energy solutions, one appropriate to the region inside the event horizon and one outside. Following the classic Damour-Ruffini method, the researchers connect these solutions across the horizon using Kruskal modes, and the resulting Bogoliubov transformation reveals the hallmark of the Hawking effect: the vacuum state of one observer appears as a thermal mixture to another, with temperature T equal to one over eight pi times the black hole mass. The transformation mixes creation and annihilation operators in a way that entangles modes on either side of the horizon, encoding the thermal radiation in the very structure of the quantum field.

With this framework in place, Li and Huang consider two observers, Alice and Bob, who initially share a Werner state, a standard two-qubit mixed state whose correlation content is tuned by a single parameter p. Alice stays safely in the asymptotically flat region far from the black hole, while Bob falls freely inward. As Bob crosses into the quantum atmosphere, Hawking radiation transforms their shared state into a tripartite one, with Bob’s degrees of freedom split between a mode outside the horizon, which is physically accessible, and a mode inside, which is forever inaccessible. Tracing over the inaccessible interior modes yields the reduced state that governs correlations Alice can actually observe. The authors then apply the geometric measure of discord, introduced by Dakic, Vedral and Brukner, which has the enormous practical advantage of a closed-form analytical expression, avoiding the difficult extremal optimization that plagues the original discord definition in curved spacetime settings.

The central surprise emerges when the discord is plotted as a function of normalized radial distance from the horizon. In the physically accessible region, the geometric quantum discord first decreases and then increases as distance grows, eventually converging to a maximum value of roughly 0.0457. In the physically inaccessible region behind the horizon, the discord does precisely the opposite, rising and then falling. This mirror-image behavior is the signature of redistribution: the Hawking effect drains quantum correlation from the region observers can probe while simultaneously pumping it into the causally disconnected interior. Most strikingly, the extreme values of discord occur exactly where the local Hawking radiation intensity peaks, in the interval where the radial coordinate lies between 1.43 and 1.5 horizon radii. The quantum atmosphere’s hottest zone thus imprints a sharp, measurable fingerprint on the correlation structure of any quantum system passing through it.

The study also dissects how two key parameters control this redistribution process, and finds that they pull in opposite directions. The first is the Hartle-Hawking constant, a parameter appearing in the local temperature profile of the Hartle-Hawking vacuum state, which describes the black hole in thermal equilibrium with its own radiation. The local temperature vanishes exactly at the horizon and asymptotically approaches the standard Hawking temperature far away, but its peak position and height depend on the Hartle-Hawking constant, which must exceed a critical value of about 23.03 for the temperature to remain physically sensible everywhere. The authors show that increasing this constant amplifies the redistribution effect, deepening the discord loss in the accessible region and correspondingly enhancing it in the inaccessible one. Increasing the event horizon radius, by contrast, suppresses the redistribution, smoothing out the variation of discord across the atmosphere.

Extending the analysis to the case where both Alice and Bob fall into the black hole produces an even richer structure. The shared state becomes four-partite, and the researchers derive analytical expressions for the discord of all six possible pairings of the resulting modes. A remarkable trade-off relation emerges: the product of the discords of the two same-region pairs, both outside or both inside the horizon, equals the product of the discords of the two cross-horizon pairs. This conservation of relative ratios demonstrates that Hawking radiation distributes quantum discord unevenly but according to a strict accounting rule between accessible and inaccessible sectors. The cross-horizon discords also obey a firm upper bound of one eighth, regardless of how the initial state parameter is chosen, hinting at fundamental limits on how much quantum correlation can straddle the horizon.

The robustness of strongly correlated initial states offers another practically important insight. As the Werner state parameter p moves away from one half, the initial discord grows, and the analysis shows that such strongly correlated states retain a non-negligible amount of discord in the physically accessible region even where Hawking radiation is at its fiercest. Weakly correlated states, by contrast, are stripped almost bare near the radiation peak. The discord curve is symmetric about p equals one half, reaching its global minimum there, and the accessible-region discord always bottoms out within the same 1.43 to 1.5 horizon-radius window across all values of the Hartle-Hawking constant. For anyone thinking about quantum information protocols in strong gravitational fields, the lesson is clear: correlation depth buys decoherence resistance.

These findings do not resolve the black hole information paradox, but they sharpen the questions at its heart. By demonstrating that the quantum atmosphere is not a passive backdrop but an active arena where quantum correlations are sorted, traded and conserved according to precise rules, the study adds quantitative texture to the information flow around black holes. The redistribution of discord between accessible and inaccessible regions is a concrete, calculable instance of how Hawking radiation mediates the exchange of quantum information across the horizon, and the precise alignment of discord extrema with the radiation peak suggests that local thermodynamics and quantum correlation dynamics are deeply intertwined. As relativistic quantum information theory matures, results like these move the field closer to understanding whether, and how, the quantum information swallowed by black holes might one day be accounted for in full.

Subject of Research: Dynamics of geometric quantum discord in the quantum atmosphere of a Schwarzschild black hole under Hawking radiation

Article Title: Geometric quantum discord in the black hole quantum atmosphere

Article References: Li, S., & Huang, X. (2026). Geometric quantum discord in the black hole quantum atmosphere. The European Physical Journal C, 86(9), Article 1107. https://doi.org/10.1140/epjc/s10052-026-16283-x

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16283-x

Keywords: quantum discord, black holes, Hawking radiation, quantum atmosphere, Schwarzschild black hole, quantum information, event horizon, Werner states, relativistic quantum information, Hartle-Hawking temperature, information paradox, quantum correlations

The signature in question is a subtle shift in what cosmologists call the effective number of relativistic species, denoted N_eff. Any black hole with a mass below roughly 10^15 grams would, by quantum physics, have completely evaporated by today through Hawking radiation, the process by which black holes slowly leak particles and energy. Tiny primordial black holes formed in the early universe would have evaporated almost immediately, injecting entropy and new particles into the primordial plasma of photons, electrons, and neutrinos. This injection of energy subtly raises the radiation content of the universe, an effect that would be recorded today as an excess in N_eff above its standard value of 3.044. Precision measurements of the cosmic microwave background, the relic glow of the Big Bang, can in principle detect such an excess, making evaporating primordial black holes accessible not through telescopes or detectors, but through cosmological bookkeeping.

The difficulty, as the authors emphasize, is that most theoretical studies have modeled primordial black hole populations as monochromatic, meaning every black hole has exactly the same mass. This is a convenient simplification, but it erases the rich variety of mass functions that different formation mechanisms actually predict. A brief ultra-slow-roll phase during inflation, when the universe expanded at a nearly frozen rate, produces a log-normal mass function peaked at a characteristic mass. A scale-invariant spectrum of primordial fluctuations collapsing during a radiation-dominated era produces a power-law distribution. Critical gravitational collapse, in which density perturbations hover just above the threshold for forming a black hole, yields a function with a low-mass power-law tail and an exponential cutoff. And metric preheating, a process in which the oscillations of the inflaton field at the end of inflation resonantly amplify perturbations, produces a numerically determined, sharply peaked distribution.

The new work goes further still by combining all four mechanisms into a single multimodal population, in which several formation channels operate simultaneously at different cosmic epochs and mass scales. Crucially, the authors do not treat the relative weights of the four sub-populations as free parameters to be fitted by hand. Instead, each channel’s weight is derived from the primordial collapse probability, which itself is fixed by the amplitude of the primordial curvature power spectrum at the scale associated with that channel. This ties the resulting mass function directly to inflationary model building: for any specific multi-feature model of the early universe, the population fractions, and therefore the predicted imprint on N_eff, are in principle calculable rather than assumed. The framework was implemented in a new public code called FRISHBEE, an extension of the existing FRISBHEE package, which solves the coupled Friedmann-Boltzmann equations governing the evaporation of the black hole population and the heating of the cosmic plasma.

The numerical results deliver a striking message: the monochromatic approximation systematically underestimates the cosmological imprint of evaporating primordial black holes. Working with initial black hole masses of 10^5, 10^7, and 10^8 grams, all safely within the window where evaporation completes before Big Bang nucleosynthesis and well before neutrino decoupling, the team computed the excess in the effective number of relativistic species for five distributions, three spin configurations, and three weighting schemes. Extended mass functions enhance the signal over the monochromatic benchmark by factors ranging from about 1.03 for the critical collapse case to roughly 1.84 for the log-normal case, with the power-law, metric preheating, and multimodal mixtures falling in between. The hierarchy among the distributions is preserved across the entire mass window, confirming that the shape of the mass function, rather than its characteristic mass alone, is the primary determinant of the evaporation signal.

The physical reason for the enhancement is intuitive. Broad distributions contain a population of lighter black holes that evaporate earlier, injecting their energy into the plasma when the universe was hotter and the number of available particle degrees of freedom was larger. Each unit of deposited energy therefore produces a bigger effect on the radiation content. A monochromatic population, in contrast, dumps all of its energy at one characteristic epoch and misses this compounding advantage. The multimodal mixture, dominated in the mean-mass-weighted scenario by its log-normal component at nearly seventy percent of the total weight, produces an enhancement of about 77 percent over the monochromatic case, and the result proves robust across all three physically motivated weighting scenarios the authors tested, never dropping below roughly 43 percent.

Detectability is where the shape effect becomes potentially decisive. The forthcoming CMB-S4 experiment and the Simons Observatory are expected to measure the excess in the effective number of relativistic species with sensitivities of about 0.06 and 0.05 respectively. For a population of primordial black holes with initial mass 10^7 grams and a scalar dark radiation species emitted by evaporation, the monochromatic prediction of approximately 0.057 falls below the detection threshold and is effectively invisible. But the log-normal, power-law, and multimodal distributions yield values between roughly 0.096 and 0.105, crossing the threshold at the 1.6 to 1.7 sigma level. In other words, the same underlying black hole population can flip from undetectable to marginally observable simply because of the shape of its mass function. For a spin-2, graviton-like dark radiation species, however, the absolute signal drops by more than an order of magnitude, and none of the distributions considered would be detectable at this mass scale.

Spin adds a further twist, and one of the study’s most surprising results. Rotating black holes emit Hawking radiation more efficiently through a process called superradiance, in which co-rotating wave modes are amplified rather than absorbed, extracting both energy and angular momentum from the hole. The amplification grows steeply with the spin of the emitted quantum: negligible for scalars, a few percent for photons, and more than one hundred percent for gravitons, while fermions are protected by Pauli blocking. The authors find that near-extremally spinning black holes, modeled with a Gaussian spin distribution centered at a spin parameter of 0.99, can boost the evaporation signal by factors of six to twelve for spin-2 dark radiation. Yet here the logic reverses: in broad mass distributions, lighter black holes spin down and shed their angular momentum long before they finish evaporating, averaging away the superradiant advantage. For near-extremal spin and spin-2 dark radiation, this effect can outweigh the mass-broadening enhancement entirely, causing the monochromatic approximation to actually overestimate the signal, an inversion the authors quantify in detail.

Taken together, these results elevate the effective number of relativistic species from a mere bound on the existence of primordial black holes to a diagnostic probe of their formation history. A future measurement at the sensitivity of CMB-S4 or the Simons Observatory would discriminate between formation scenarios: a signal near the extended-distribution predictions would disfavor a monochromatic or critical-collapse-dominated population at the 10^7 gram scale, while a non-detection would constrain the broader scenarios. Because the multimodal framework connects the population weights directly to the primordial power spectrum through the collapse probability, precision cosmology could in principle probe the number and relative amplitude of features in the inflaton potential, the physics that governed the universe’s earliest moments. The authors caution that their analysis is confined to the pre-Big Bang nucleosynthesis mass window, and that extended distributions with heavy tails approaching the boundary may require a treatment of competing dilution effects, which they defer to future work. But the central conclusion stands: in the hunt for primordial black holes through their Hawking afterglow, the shape of the population is everything.

Subject of Research: Cosmological imprint of evaporating primordial black holes with multimodal mass and extended spin distributions on the effective number of relativistic species

Article Title: Evaporation of primordial black holes with multimodal mass and extended spin distributions: cosmological imprints on the effective number of relativistic species

Article References: Toghrai, T., Daassou, A., Ouchhaine, Y., Laassiri, H., & Benbrik, R. (2026). Evaporation of primordial black holes with multimodal mass and extended spin distributions: cosmological imprints on the effective number of relativistic species. The European Physical Journal C, 86(9), Article 1087. https://doi.org/10.1140/epjc/s10052-026-16372-x

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16372-x

Keywords: primordial black holes, Hawking radiation, effective number of relativistic species, cosmic microwave background, inflation, mass distribution, black hole spin, superradiance, dark radiation, CMB-S4, Simons Observatory, early universe