Brans-Dicke theory, first proposed in the 1960s, is one of the oldest and most celebrated extensions of Einstein’s general relativity. In it, the strength of gravity itself is not fixed but governed by a scalar field that varies across space and time. The original theory contains a single adjustable parameter, the Brans-Dicke parameter, which controls how strongly this scalar field interacts with the curvature of spacetime. Solar system experiments, notably radio tracking of the Cassini spacecraft, constrain this parameter to be enormous, roughly greater than ten thousand, which pushes ordinary Brans-Dicke theory very close to general relativity and limits its usefulness as a cosmological alternative. The new work revives the framework by generalising it in two directions at once: the scalar field is non-minimally coupled to gravity through a function of the field, and its kinetic term is non-canonical, meaning the field’s motion contributes to the cosmic energy budget in a field-dependent way rather than in the standard quadratic form.

Technically, the authors work with an action in which the Ricci scalar is multiplied by a function f of the scalar field, taken to be one plus a coupling constant xi times the field, while the scalar sector takes the form of a kinetic term divided by that same function, minus a potential. This particular structure is inspired by Noether symmetry arguments, a mathematical technique for selecting physically consistent Lagrangians. The construction is deliberately broad: setting the coupling xi to zero and choosing a constant potential recovers standard LambdaCDM, while other choices reproduce quintessence and the original Brans-Dicke theory. The non-canonical kinetic term, in which the usual kinetic energy is scaled by the inverse of the gravitational coupling function, is what distinguishes this model from most scalar-tensor studies and, as the authors show, noticeably changes the resulting cosmic dynamics even when the dynamical variables are defined in ways similar to those used in simpler models.

To analyse the cosmology, the team assumes a flat, homogeneous and isotropic universe filled with pressureless matter and radiation, then rewrites the modified Friedmann equations and the scalar field equation as an autonomous system of first-order differential equations. This dynamical systems approach, a staple of modern theoretical cosmology, replaces the difficult task of solving the equations outright with the more tractable problem of finding equilibrium states, called critical points, and determining whether they are stable. The dimensionless variables capture the kinetic and potential energy of the scalar field relative to the total energy density, along with rescaled matter and radiation densities. Because of the non-minimal coupling, these rescaled densities are not the physical ones, but they relate to them through the coupling function, a subtlety the authors handle carefully throughout.

A key constraint shapes the whole analysis. For the phase space of solutions to be bounded, and therefore mathematically well-behaved, the coupling parameter omega must be positive; if it were negative, the allowed region of the phase space would become a hyperbola stretching to infinity, and the phantom-like field it describes would require more elaborate compactification techniques. Positivity of omega also guarantees the absence of ghosts and Laplacian instabilities, the pathological modes that render many modified gravity theories unphysical. With these restrictions in place, the authors investigate three widely used potentials in turn: a constant potential, a power-law potential and an exponential potential, visualising the flow of cosmic histories with bounded phase portraits.

For the constant potential, the analysis uncovers a sequence of critical points that mirrors the broad architecture of LambdaCDM. A radiation-dominated saddle point gives way to a mixed matter-and-scalar-field epoch, and the universe finally settles into a scalar-field-dominated attractor that can drive accelerated expansion, provided omega exceeds half of xi squared, a condition easily satisfied when deviations from general relativity are small. Intriguingly, the matter-dominated point in this model always carries some scalar field contribution, a feature absent from simpler quintessence models, though tuning the couplings can suppress it. Numerical evolutions show that choosing omega at least an order of magnitude larger than xi shifts matter-radiation equality to a redshift of order one thousand, matching observations, whereas smaller values of omega push equality to redshifts around one hundred, which is ruled out by data.

The power-law potential, written in a form that keeps the key dynamical variable constant, produces richer behaviour with two additional saddle points at which the scalar field mixes with radiation or matter. These mixed epochs can even exhibit transient acceleration for a narrow range of negative exponents, though the authors caution that this does not guarantee a consistent later cosmology. For the illustrative inverse-square case, viability requires the matter-dominated point to be a saddle and the late-time attractor to exist, which constrains the ratio of omega to xi squared to lie either below one half or above six and a half. Notably, the team finds no tracking behaviour of the kind familiar from minimally coupled quintessence, which they attribute to the non-minimal coupling, since tracking solutions were never designed for such theories. They also verify that their critical points reduce exactly to those of standard Brans-Dicke cosmology under an appropriate limit, a reassuring consistency check.

The exponential potential brings the most surprising result. Because of the non-minimal coupling, the slope variable of the potential becomes dynamical rather than fixed, and this generates an entirely new critical point: a fully scalar-field-dominated, dark and accelerating attractor that cannot exist in standard quintessence with an exponential potential. This point is stable for all relevant values of the coupling and serves as a genuine late-time dark energy attractor. However, the model also exhibits an early epoch of kination, in which the kinetic energy of the scalar field dominates, which does not appear in the standard cosmic history and must be suppressed through careful choice of initial conditions. Fitting the physical density parameters to a LambdaCDM-like evolution suggests that couplings of order one require an essentially flat potential, while smaller couplings allow a gentle exponential slope of order one tenth.

Across all three potentials, a consistent picture emerges: viable cosmologies exist, but they demand that deviations from general relativity be modest, with the coupling xi comfortably below one, and that the parameters satisfy specific stability and acceleration conditions summarised by the authors in a compact table. The dynamical behaviour differs noticeably from other scalar-tensor models despite superficially similar variables, underlining how the combination of non-minimal coupling and non-canonical kinetics reshapes the cosmic phase portrait. None of the models provides an early inflationary starting point from the dynamics alone, meaning these scenarios are best understood as descriptions of the late universe, though the exponential case offers a saddle point that could act as an inflationary transition for suitable parameters.

The study is a background-level analysis, and the authors are explicit about what comes next. Cosmological perturbations in these models must be evolved and confronted with observations of structure growth and the cosmic microwave background before any claim can be made that they ease the Hubble tension or the clustering tension. The team also flags the Palatini formalism, an alternative variational approach in which Brans-Dicke parameters can take much lower values, as a promising direction. Still, the message of the work is striking: a sixty-year-old alternative to Einstein’s gravity, generalised with a non-canonical twist, can reproduce the entire background history of our universe, from radiation domination through matter domination to today’s accelerated expansion, and may yet hold clues to the anomalies that currently trouble precision cosmology.

Subject of Research: Cosmological dynamics of a non-canonical generalised Brans-Dicke scalar-tensor theory as an alternative to LambdaCDM dark energy

Article Title: Cosmological dynamics of a non-canonical generalised Brans-Dicke theory

Article References: Debono, M., Said, J. L., & Farrugia, G. (2026). Cosmological dynamics of a non-canonical generalised Brans-Dicke theory. General Relativity and Gravitation, 58(10), Article 116. https://doi.org/10.1007/s10714-026-03616-0

Image Credits: AI Generated

DOI: 10.1007/s10714-026-03616-0

Keywords: Brans-Dicke theory, dark energy, cosmology, dynamical systems, scalar-tensor gravity, Hubble tension, LambdaCDM, modified gravity, quintessence, critical points, non-canonical scalar field, General Relativity and Gravitation

The work, carried out by Yerlan Myrzakulov of L. N. Gumilyov Eurasian National University in Astana, Saddam Hussain of Zhejiang University of Technology in Hangzhou, and Mohd Shahalam of Integral University in Lucknow, focuses on a class of models known as quintessence. In quintessence scenarios, dark energy is not the unchanging cosmological constant of the standard Lambda-CDM model but a dynamical scalar field that rolls down a potential energy landscape over cosmic time. The specific potential studied here is a double exponential, first proposed in the context of a recently introduced scenario called Quintessence-Driven Slow-Contraction Cold Dark Matter, or Q-SC-CDM. In that scenario, dark energy gradually loses strength, cosmic acceleration eventually ceases, and the universe drifts into a period of slow contraction that terminates in a Big Crunch.

The authors’ key innovation is to add a non-minimal coupling between the scalar field and gravity. In standard quintessence, the field interacts with gravity only through its energy and pressure. In the non-minimally coupled version, the field multiplies the Ricci scalar, the mathematical object that encodes spacetime curvature, through a coupling function of the form F(phi) = 1 minus kappa squared xi phi squared. The dimensionless parameter xi controls the strength of this coupling. When xi equals zero, ordinary minimal coupling is recovered; the value xi equal to one sixth corresponds to the theoretically motivated conformal coupling; and large values of xi can radically reshape cosmic evolution. This kind of coupling has a distinguished pedigree, tracing back to Brans-Dicke theory in 1961, and it has recently gained attention because a strongly coupled Higgs field can drive cosmic inflation.

Technically, the coupling changes everything. The Friedmann equation governing the expansion rate acquires extra terms proportional to xi, including contributions proportional to the Hubble parameter times the field and its derivative, and the field’s equation of motion picks up a curvature-driven term xi R phi. The researchers solved these coupled equations numerically on a spatially flat Friedmann-Lemaitre-Robertson-Walker background, starting from a matter-dominated early universe in which the field is frozen by Hubble damping. Their numerical experiments, sweeping the coupling across the range from 0.01 to 1 and varying the potential parameters alpha and beta, revealed a striking bifurcation. For couplings below roughly 0.2, and for sufficiently large potential parameters, the universe passes through decelerated expansion, then acceleration, and finally a slow-contraction phase in which the Hubble parameter falls to zero and the scale factor collapses. For couplings above roughly 0.2, the contraction phase vanishes entirely: the universe accelerates today and settles into de Sitter expansion, an exponentially expanding state with a constant Hubble rate, in the far future.

But a model that avoids the Big Crunch is worthless if it cannot reproduce the universe we actually observe. The crucial test is whether the model produces a sufficiently long matter-dominated epoch at high redshift, the era in which galaxies, clusters, and the cosmic microwave background patterns formed. The authors found that many parameter combinations producing a future slow-contraction phase fail this test badly: the matter era is cut short, the field’s fractional energy density turns negative or exceeds unity at high redshift, or the dynamics diverge. Such solutions, however dramatic their future, are physically non-viable. This observation alone begins to undermine the Q-SC-CDM picture, since the very trajectories that lead to a Big Crunch cannot sustain a consistent past.

To make the case quantitative, the team confronted three versions of the model with a formidable battery of data. The datasets included 32 cosmic chronometer measurements of the Hubble parameter spanning redshifts from 0.07 to nearly 2; the Pantheon+ compilation of 1701 Type Ia supernovae; the Dark Energy Survey five-year sample of 1829 supernovae; the Union 3 compilation of 2087 supernovae; and seven baryon acoustic oscillation data points from the second data release of the Dark Energy Spectroscopic Instrument, DESI, with the sound horizon at the drag epoch treated as a free parameter. Using Markov chain Monte Carlo sampling with the dynesty nested sampler, they computed posterior distributions and Bayesian evidence for each model relative to Lambda-CDM, interpreting evidence differences on the revised Jeffreys scale.

The results were unambiguous. Model I, with the coupling fixed at xi equal to 0.3, yields a Hubble constant near 70 kilometers per second per megaparsec and a matter density around 0.295, but its Bayesian evidence falls far below that of Lambda-CDM across every dataset combination. Model II, with xi equal to 0.5, fares even worse: it prefers a matter density of only about 0.255, produces a sound horizon near 145 megaparsecs, and shows the largest deviations in the Om(z) diagnostic, a low-redshift test that flags dynamical dark energy by measuring departures from constancy in a normalized combination of Hubble measurements. Both fixed-coupling models predict Hubble constants above 70 kilometers per second per megaparsec, which might superficially look like a route to easing the Hubble tension, but the authors caution that the severe evidence penalty means this should not be read as a genuine resolution of that conflict.

The most interesting case is Model III, in which the coupling is allowed to float between 0.01 and 0.7. Here the data consistently pull the coupling down to a small value, xi of roughly 0.12, close to the conformal value of one sixth. The best-fit evolution closely mimics Lambda-CDM: the matter density stays near unity at high redshift, the effective equation of state remains at zero through the matter era and saturates at minus one in the future, and the potential-related combination y squared minus A squared stays positive throughout, meaning the effective potential never turns negative. Because a negative potential is precisely what triggers recollapse in these models, the observationally favored parameter space simply never enters the Big Crunch regime. Even for this modest coupling, well below the threshold of 0.2 that mechanically eliminates contraction, the data steer the model into the corner of parameter space where the potential parameters are small enough that slow contraction never occurs.

An independent phase-space analysis confirms the picture. Rewriting the field equations as a first-order autonomous system, the authors identified five stationary points and classified them by the eigenvalues of the linearized flow. Points A1, A2, and A4 are saddle points corresponding to transient matter-dominated or non-accelerating phases; A3 is a non-hyperbolic point with de Sitter characteristics; and A5 is a one-parameter family of equilibrium points that acts as a stable attractor. For the observationally preferred parameters, with alpha near ten to the minus four, beta near 0.26, and xi near 0.12, the three transverse eigenvalues at A5 are negative, and phase portraits show trajectories from across the phase space converging onto the equilibrium curve. At that attractor, the field freezes, the Hubble derivative vanishes, the dark energy density parameter approaches one, and the effective equation of state approaches minus one, exactly the de Sitter endpoint implied by the observational fits.

The study leaves the standard cosmological model standing, at least at the level of background evolution: every variant considered was strongly disfavored relative to Lambda-CDM, with evidence differences exceeding five on the logarithmic scale. Yet the work delivers a substantive negative result with real physical content. The coupling strengths that would spare the universe a Big Crunch through brute force are ruled out by data, while the couplings the data prefer steer the universe away from collapse more subtly, by selecting potential shapes that never go negative. The framework also remains phantom-free, with a positive kinetic term that avoids ghost instabilities even as DESI measurements hint at evolving dark energy in other models. The authors note that the real test lies ahead: extending the analysis to cosmological perturbations, confronting the model with full cosmic microwave background data and large-scale structure surveys, and checking its impact on the matter clustering amplitude sigma eight. For now, the verdict is reassuring. The same observations that map the expansion history of the cosmos appear to be quietly vetoing the universe’s most dramatic possible ending, pointing instead toward a future of gentle, eternal expansion.

Subject of Research: Nonminimally coupled quintessence dark energy models and observational constraints on the future fate of the universe

Article Title: Nonminimally coupled quintessence with double exponential potential: observational evidence against big crunch singularity

Article References: Myrzakulov, Y., Hussain, S., & Shahalam, M. (2026). Nonminimally coupled quintessence with double exponential potential: observational evidence against big crunch singularity. The European Physical Journal C, 86(9), Article 1055. https://doi.org/10.1140/epjc/s10052-026-16306-7

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16306-7

Keywords: dark energy, quintessence, nonminimal coupling, Big Crunch singularity, de Sitter expansion, cosmology, scalar field, Bayesian evidence, DESI, supernovae, Lambda-CDM, phase-space analysis

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 technical foundation of the work is a carefully chosen seed metric. The authors begin with a static, spherically symmetric Schwarzschild-type line element whose metric function f(r) gathers every energy contribution into a single expression: the black hole’s mass, its electric charge squared, its magnetic charge squared, a term proportional to the cosmological constant Λ, and a quintessence term of the form α times r to the power minus one minus three omega. Here α measures the strength of the dark energy field and the parameter ω, constrained to lie between −1 and −1/3, encodes how the pressure of that field relates to its energy density. This functional form descends from the classic 2003 solution of V. V. Kiselev, which first showed how the Einstein equations can be solved exactly for a black hole surrounded by quintessence. Crucially, in the new construction the cosmological constant and the quintessence are treated as genuine external matter sources distributed throughout spacetime, not merely as geometric corrections appended to a vacuum solution.

To set this static configuration spinning, de Oliveira and Schmidt employed the Newman–Janis algorithm, a mathematical recipe invented in 1965 that transforms non-rotating black hole metrics into their rotating counterparts. The algorithm is powerful but notorious: different choices in its intermediate steps can lead to ambiguous results. The authors sidestepped this pitfall by adopting a unique and general complexification rule for the radial coordinate, one that systematically converts every power of r in the metric function into its rotating equivalent. Under this rule, powers of r are rewritten in terms of the quantity Σ, which equals the radial coordinate squared plus the square of the spin parameter a times the square of the cosine of the polar angle. The result is a fully rotating line element in Boyer–Lindquist coordinates, whose structure function Δr combines the mass, both charges, the quintessence contribution, and a fourth-power term in r carrying the cosmological constant. Notably, Δr contains no coupling terms between Λ and the spin parameter a—a departure from the standard Kerr–Newman–AdS geometry, where such couplings are introduced precisely to keep the vacuum Einstein equations satisfied.

The payoff of that choice is conceptual clarity about what is holding the spacetime together. Because the cosmological constant and quintessence enter as matter from the very beginning, the Einstein tensor of the new solution is nonzero in ways the vacuum AdS case is not. The authors computed the full stress–energy tensor by evaluating the Ricci tensor and the Einstein tensor for the rotating metric, and they found components that scale with each energy ingredient: pieces proportional to α times ω that encode the quintessence fluid, pieces proportional to the combined charge squared, and pieces proportional to Λ. They also derived the electromagnetic field tensor of the dyonic configuration. Starting from the static solution, where Maxwell’s equations yield a radial electric field of Q_e over r squared and a radial magnetic field of Q_m over r squared, they used the Newman–Penrose formalism with null tetrad vectors to carry the field into the rotating geometry. The resulting Maxwell scalar takes the compact form involving the complex combination Q_e plus i times Q_m divided by two Σ, and the electromagnetic tensor, by symmetry, turns out to be entirely independent of Λ and α.

With the geometry in hand, the researchers mapped its most dramatic features: the horizons and the ergosphere. Setting Δr equal to zero yields the event horizon condition, a fourth-degree equation in the radial coordinate whose roots include an inner horizon, the outer Schwarzschild-like horizon, a cosmological horizon, and a quintessence horizon. Numerical analysis across representative parameter values showed that only three of the four roots are ever real, the fourth remaining complex. The trends are physically intuitive: increasing either the cosmological constant or the quintessence parameter shrinks the real horizon radii, meaning the dark energy content of the universe effectively compresses the region of no return. The ergosphere—the region outside the horizon where nothing can remain static because the spinning black hole drags spacetime itself around with it—was found by setting the time-time metric component to zero. In the equatorial plane this reduces to static circumferences that depend on Λ, α, ω, and the charges but are, curiously, independent of the spin parameter a. Larger quintessence strength again produces smaller static circumferences.

Singularities, the points where the mathematical description of spacetime breaks down entirely, were probed using the Kretschmann scalar, a coordinate-invariant measure of curvature built from the Riemann tensor. For this solution the Kretschmann scalar takes the form of a complicated parameter-dependent function divided by Σ raised to the twelfth power. It diverges only when Σ vanishes, which happens at r equals zero on the equatorial plane—exactly the ring-like singular structure familiar from the Kerr and Kerr–Newman solutions. The dark energy ingredients, however dramatic their effect on the horizon structure, leave this singularity region unchanged. The Ricci scalar tells a subtler story. In the limits where both spin and quintessence vanish, it reduces to minus four Λ, the textbook vacuum value for a spacetime containing only a cosmological constant. But whenever α is nonzero, even in the limit where the cosmological constant is absent, the scalar curvature departs from the vacuum result—a direct fingerprint of treating dark energy as matter rather than geometry.

To see how matter would actually move in this environment, the authors studied a unit-mass test particle on a circular orbit in the equatorial plane, observed from a locally non-rotating frame, a reference system in which the particle carries no angular momentum. From the geodesic equations they extracted expressions for the particle’s angular momentum, mechanical energy, angular velocity, and rotational velocity. Every energy term in the model leaves its mark on these quantities, and in the limit where both α and Λ go to zero the standard Kerr–Newman results are recovered. Intriguingly, even when the black hole itself has no spin, the locally non-rotating frame analysis yields nonzero angular and rotational velocities, a reminder that this reference frame is itself defined by the geometry rather than by the black hole’s rotation alone.

The final piece of the analysis concerned the energy conditions, the classical sanity checks that general relativity imposes on any legitimate matter source. Because the rotating metric contains off-diagonal terms, the usual shortcut expressions cannot be applied; instead, the authors constructed an orthonormal tetrad basis in the locally non-rotating frame and expressed the stress–energy tensor in terms of an energy density ε and three principal pressures. Their findings delineate when this exotic spacetime can be supported by ordinary matter and when it cannot. When the null energy condition is violated, every other energy condition falls with it. Divergences in the energy density and pressures occur at the event horizon and at a second locus defined by the vanishing of ξ—but that second locus always lies inside the horizon, so it does not affect the exterior region of physical interest. For concrete parameter choices, such as ω equal to −2/3 with modest charge, mass, quintessence strength, and cosmological constant, there exists an exterior band of radii where the weak energy condition holds, meaning any observer measures positive energy density and the geometry is supported by non-exotic, classical matter. Outside those bands, the spacetime demands exotic ingredients—precisely the quintessence field of the model, or potentially quantum effects.

The new solution is a rich family rather than a single object. Switching off parameters recovers a parade of familiar geometries: without spin, the dyonic cosmological black hole surrounded by quintessence; without quintessence, the rotating charged cosmological black hole; without a cosmological constant, the dyonic Kerr–Newman black hole in quintessence; with both dark energy terms removed, the Kasuya–Kerr–Newman solution; and with both charges zeroed, the rotating cosmological black hole surrounded by quintessence. The work thus gives theorists a versatile laboratory for exploring how dark energy reshapes black hole shadows, accretion disks, quasinormal ringing, and gravitational lensing—observables that upcoming surveys and the next generation of gravitational wave detectors may one day constrain. For now, the result stands as an exact, matter-inclusive answer to a deceptively simple question: what does a spinning, charged black hole look like in a universe that refuses to sit still?

Subject of Research: An exact general relativistic solution for a rotating dyonic black hole surrounded by quintessence dark energy

Article Title: Dyonic rotating cosmological black hole surrounded by quintessence

Article References: Dyonic rotating cosmological black hole surrounded by quintessence. (n.d.). https://doi.org/10.1140/epjc/s10052-026-16309-4

Image Credits: AI Generated

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

Keywords: black holes, general relativity, quintessence, dark energy, Kerr–Newman solution, Newman–Janis algorithm, event horizon, ergosphere, energy conditions, cosmological constant, dyonic black hole, theoretical physics