The research team, led by M. Castillo Alarcón of the Center for Research and Advanced Studies of the National Polytechnic Institute in Mexico City, together with Leonel Bixano, I. A. Sarmiento-Alvarado, Claudio Salas-Pérez and the veteran dark matter theorist Tonatiuh Matos, set out to correct a persistent weakness in the literature. Most existing models of black holes wrapped in dark matter halos begin by assuming a relationship between the two fundamental functions that define a static, spherically symmetric spacetime: the temporal metric function A(r), which governs gravitational redshift and time dilation, and the radial metric function B(r), which governs spatial curvature along the radial direction. In the pure Schwarzschild vacuum solution these functions are exact reciprocals of one another, satisfying A(r)B(r) = 1. Many authors simply carry that reciprocal relation over to halo systems as a convenient ansatz, without checking whether the Einstein field equations actually permit it once matter is present.

Castillo Alarcón and colleagues refused to make that assumption. They started from the most general static, spherically symmetric line element, with A(r) and B(r) treated as fully independent functions, and modeled the dark matter halo as an anisotropic fluid, one whose radial pressure differs from its tangential pressure. This flexibility matters because dark matter need not behave like a perfect fluid; different microscopic candidates, from cold collisionless particles to ultralight scalar fields, can generate different stress-energy configurations. By feeding three widely used halo density profiles directly into the Einstein equations, the Burkert profile favored by observations of dwarf and low-surface-brightness galaxies, the Einasto profile that fits simulated halos so well, and the multistate scalar field dark matter model developed over decades by Matos and collaborators, the team solved for the metric functions, the radial pressure and the tangential pressure simultaneously, without imposing any of them by hand.

The construction works as follows. The radial function B(r) is determined by the total gravitational mass enclosed within a given radius, which combines the black hole mass with the accumulated halo mass. Crucially, the authors assume that dark matter vanishes inside a critical radius r_in, equal to three times the Schwarzschild radius, because any dark matter orbiting closer than the innermost stable circular orbit of a massive particle would be captured by the black hole. Inside that boundary the spacetime is exactly Schwarzschild, in accordance with Birkhoff’s theorem, though the surrounding halo imprints a uniform gravitational redshift on the inner vacuum region that cannot be removed globally. Outside the boundary, the temporal function A(r) is built from a deformation potential that measures precisely how the halo warps the clock rate relative to the isolated black hole case.

The radial pressure of the halo is reconstructed from the observed circular velocities of stars, using the exact relativistic relation between circular geodesics and the gravitational potential. The tangential pressure then follows from the conservation of the stress-energy tensor, the anisotropic generalization of the Tolman-Oppenheimer-Volkoff equation that governs hydrostatic equilibrium in relativistic stars. With the full system in hand, the authors verified that all three density profiles satisfy the strong energy condition, meaning the halo source is physically reasonable rather than exotic in a pathological sense. They also demonstrated that the transition between the inner vacuum region and the outer halo at r_in satisfies the Darmois-Israel junction conditions: the induced metric and the extrinsic curvature match perfectly on both sides, so no thin shell of exotic surface matter, no surface pressure and no surface tension appear at the boundary. The halo begins smoothly, not as a membrane.

The headline result emerged from both numerical computation and a semi-analytical proof. Even though no reciprocal relation was ever assumed, the product A(r)B(r) equals one plus corrections of order ten to the minus six throughout the halo. The reason is a small dimensionless parameter that the authors call the halo strength, defined as eight pi times the gravitational constant times the central halo density times the square of the scale radius, divided by the square of the speed of light. For Milky Way parameters this factor is roughly a millionth. The Einstein equations show that the variation of the product A(r)B(r) is controlled entirely by this parameter, with the radial pressure contributing only a second-order correction and the tangential pressure not entering the relation at all. In other words, the reciprocal structure so often assumed in the literature is not an arbitrary guess; it is a derived consequence of the sheer weakness of the halo’s relativistic pull.

This finding has immediate practical value. Shadow calculations, gravitational lensing studies, quasinormal mode analyses and accretion disk models that adopt A(r) = 1/B(r) for halo systems can now do so with a theoretical license rather than a leap of faith. The authors note that the result is robust even if the sharp inner cutoff of the halo is replaced by a smooth taper, and that the three density profiles, despite their very different origins as phenomenological fits, simulation-motivated shapes and microscopic field theory constructions, produce nearly identical metric functions. From the perspective of light rays skimming the black hole, the team’s null geodesic calculations show that one essentially cannot tell whether the black hole is embedded in a halo or floating in perfect vacuum; the differences only become visible at distances far from the hole.

To test whether the model survives contact with reality, the researchers confronted it with two of the most precise measurements in modern astrophysics. The first is the pericenter precession of the star S2, which whips around Sagittarius A* every sixteen years on a highly eccentric orbit, passing within about fourteen hundred Schwarzschild radii of the black hole. The GRAVITY Collaboration has measured the relativistic precession of this orbit to remarkable accuracy. Integrating the geodesic equations in the halo spacetime for all three density profiles, the team obtained a precession of 12.1687 arcminutes per orbit, with a fitted scaling factor of 1.11 plus or minus 0.21, comfortably inside the observational range reported by GRAVITY.

The second test is the shadow of Sagittarius A* itself. Using the general formula for the angular radius of a shadow in a static, spherically symmetric spacetime, and adopting the measured distance to the galactic center of 8178 parsecs, the authors computed an angular diameter of 50.1746 microarcseconds for all three halo profiles. The Event Horizon Telescope reports 51.8 plus or minus 2.3 microarcseconds, so the prediction sits squarely within the error bars. Two independent observables, one probing the orbital dynamics of matter far from the hole and the other probing the bending of light near the photon sphere, are simultaneously satisfied. The dark matter halo, at least with the densities inferred for the Milky Way, does not push the spacetime outside what astronomers have already measured.

The study does leave honest caveats. The construction delivers a self-consistent static equilibrium but does not prove dynamical stability under radial perturbations; because the radial pressure depends nonlocally on the enclosed halo mass, a proper stability analysis would require solving the linearized Einstein equations with specified constitutive relations, a task the authors explicitly defer to future work. The treatment of the region inside the innermost stable orbit, where captured dark matter would accumulate, is also acknowledged as a simplification. Nevertheless, the paper’s central contribution stands: the near-reciprocity of the metric functions in black hole-halo systems is now a proven property rooted in the feeble relativistic strength of galactic dark matter, quantified by a single small parameter. As horizon-scale interferometry sharpens its resolution and stellar orbit monitoring extends to fainter stars closer to the event horizon, the framework offers theorists a rigorous foundation for asking the next question: how small a deviation from Schwarzschild would it take for the Event Horizon Telescope, or its successors, to finally see the dark matter that surrounds the shadow?

Subject of Research: Exact general relativistic modeling of a Schwarzschild black hole surrounded by a dark matter halo

Article Title: Determination of the metric of a Schwarzschild black hole surrounded by a halo of dark matter

Article References: Determination of the metric of a Schwarzschild black hole surrounded by a halo of dark matter. (n.d.). https://doi.org/10.1140/epjc/s10052-026-16375-8

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16375-8

Keywords: black holes, dark matter, general relativity, Schwarzschild metric, dark matter halo, Sagittarius A*, Event Horizon Telescope, GRAVITY Collaboration, S2 star, anisotropic fluid, Einstein field equations, shadow radius

The central subject of the new perspective is the shock development problem: the question of what happens to a compressible fluid after a shock wave has begun to form. Shocks are among the most familiar singularities in nature, responsible for the crack of a sonic boom and the abrupt pressure fronts that propagate through air after an explosion. Mathematically, however, they pose a profound challenge. The equations of inviscid compressible fluid dynamics, the Euler equations, describe smooth flow up to a moment when the gradient of certain quantities blows up. Before that moment, the solution is classical and differentiable; after it, discontinuities appear and the equations in their classical form cease to make sense. Understanding how to continue the solution past this breakdown point, in a way that is both physically meaningful and mathematically rigorous, had remained open for more than a century and a half.

Christodoulou’s most recent monograph, The Shock Development Problem, published by the EMS Publishing House in 2019, finally resolved this question. The new article reviews the mathematical methods of that monograph, and in doing so it reveals an unexpected intellectual lineage. The machinery used to track the formation and propagation of shocks was not invented in isolation; it is, in essential respects, an adaptation of the geometric framework that Christodoulou had earlier built together with Sergiu Klainerman in their landmark 1993 work, The Global Nonlinear Stability of the Minkowski Space, published in the Princeton Mathematical Series. That result demonstrated, using purely geometric methods, that Minkowski spacetime, the flat spacetime of special relativity, is stable against small nonlinear perturbations, thereby establishing the nonlinear stability of the simplest solution of Einstein’s field equations.

The bridge between these two domains lies in what mathematicians call Lorentzian geometry, the geometry of spacetimes in which the metric has a signature that distinguishes timelike directions from spacelike ones. In general relativity, the fundamental objects of study are Lorentzian manifolds, and the natural way to analyze them is through foliations, families of hypersurfaces slicing the manifold, and through the behavior of null surfaces, the analogues of light fronts. Christodoulou realized that even when one studies a compressible fluid in flat spacetime, a problem that superficially has nothing to do with gravity, the same Lorentzian geometric structures are the correct tools. The characteristic surfaces along which information propagates in a fluid, and along which shocks in particular form, are naturally viewed as hypersurfaces in a Lorentzian metric constructed from the fluid equations themselves. This reframing allows the powerful estimates and construction techniques of geometric analysis to be applied to fluid mechanics.

This is not the first time Christodoulou has shown such a connection. In 2007 he published The Formation of Shocks in 3-Dimensional Fluids, also in the EMS Monographs in Mathematics, in which he proved that shocks are an inevitable feature of smooth, arbitrarily small initial data for compressible Euler flow in three dimensions. The result was striking because it showed that the nonlinearity of the fluid equations is so persistent that no smooth initial condition, however small and however carefully chosen, can avoid developing infinite gradients in finite time. Two years later, in 2009, he applied analogous methods to the relativistic setting in The Formation of Black Holes in General Relativity, proving that trapped surfaces, the hallmark boundaries of black holes, can form in the evolution of spacetimes from smooth data that is nowhere trapped initially. The structural parallels between these two problems, shock formation in fluids and trapped surface formation in gravitation, are no accident; both concern the focusing of characteristic cones and the breakdown of smooth evolution.

The new perspective article addresses the step beyond formation: what happens past the point where the shock first appears. This is the shock development problem proper. In the 2007 monograph, Christodoulou characterized the geometry of the boundary of the region of smooth flow at the exact moment of shock formation, showing that this boundary possesses a singular structure, often described in terms of a crease set, where the smooth solution meets the incipient discontinuity. The question then becomes how to define and construct the solution for later times, when the shock surface separates regions of fluid in different states. The answer required constructing the shock front itself as a free boundary, tracking its evolution with estimates of a precision rarely attempted in fluid mechanics, and verifying that the resulting weak solution satisfies the physical conservation laws in the correct distributional sense.

Technically, the framework demands an intricate interplay between hyperbolic estimates and differential geometry. One must control the derivatives of the fluid variables up to very high order on the smooth side of the shock, establish that the shock surface remains regular in an appropriate sense, and match the two across the discontinuity. The Lorentzian viewpoint is indispensable here because the shock surface and the characteristic hypersurfaces emanating from it form a geometric system whose curvature properties dictate where and how the solution can be continued. Quantities analogous to the curvature components studied in general relativity play the role of error terms that must be estimated; the deformation tensors of the foliations play the role of connection coefficients. In this way, a problem of classical fluid dynamics is recast as a problem in the geometric analysis of a Lorentzian manifold, and the hard-won lessons of decades spent on the Einstein equations become directly transferable.

The significance of this unification extends beyond fluid mechanics itself. The methods developed for the shock development problem have already stimulated further research in the mathematical community. Among the works building on these foundations is a 2022 preprint by Leon Abbrescia and Jared Speck investigating the emergence of the singular boundary from the crease in three-dimensional compressible Euler flow, which takes up precisely the geometric structures that Christodoulou’s monograph isolated. Christodoulou himself also provided a shorter account of the shock development program in a 2022 article in the Journal of Mathematical Physics, offering a compact summary of the ideas now expanded in the present review. The appearance of this new perspective in General Relativity and Gravitation, a journal at the heart of the relativity community, underscores the message that the boundaries between gravitational physics and fluid dynamics are far more porous than traditional disciplinary divisions suggest.

There is also a broader conceptual payoff. The fact that the same geometric techniques resolve questions in both general relativity and classical fluid mechanics points toward a unified theory of singularities in hyperbolic systems, the class of equations, including the wave equation, the Einstein equations, and the Euler equations, whose solutions propagate at finite speed and whose nonlinearities tend to concentrate energy into ever thinner regions. Whether the singularity is a black hole horizon forming in a collapsing spacetime or a pressure discontinuity racing through air, the mathematical pathology is structurally related: the characteristic cones focus, gradients blow up, and the smooth description of the system terminates. Geometric analysis in the Lorentzian setting provides a language in which this commonality can be made precise and, more importantly, exploited.

For Christodoulou, whose career has been marked by results of extraordinary technical depth, from his early work on black hole dynamics and the memory effect to his proofs of global stability and singularity formation, the new article serves as both a retrospective and a roadmap. It presents the shock development problem not as an isolated curiosity of fluid mechanics but as a chapter in a larger story about the geometry of propagation, one that began with the stability of spacetime itself and continues wherever waves and fronts evolve in a Lorentzian world. Readers of the journal, whether they come to it from relativity or from the theory of partial differential equations, will find in it a demonstration of how the deepest tools of one field can illuminate the oldest open problems of another.