A new preprint presents an analytical model of a star with a mixed core: ordinary matter and a second incompressible fluid interpreted as dark matter are present together, while the outer envelope contains ordinary matter alone. The two fluids are conserved separately and interact through the common spacetime geometry. The result is a controlled mathematical setting for asking how a core-confined second component relates to the internal structure and static limits of a star.
The study does not analyze observations or a sample of stars. Its objects are parameterized mathematical configurations: static, spherical arrangements of two fluids with constant energy densities in the core, joined to an ordinary-matter envelope. That distinction matters because the paper is testing the behavior of a model, not reporting evidence that dark matter has been found inside an observed star.
A boundary-defined central pressure
The model’s central pressure is not chosen as an independent starting value. It is fixed by the conditions imposed at the core-envelope interface and at the stellar surface. At the core boundary, the dark-matter pressure is set to zero; ordinary-matter pressure is matched across the interface; and ordinary-matter pressure is set to zero at the surface.
Those matching conditions also connect the geometry on the two sides of the interface. The temporal part of the metric remains continuous there, and the envelope geometry carries the dark component’s gravitational contribution through the matching calculation. In practical terms, the outer layer is ordinary-matter-only in composition, but it remains linked to what the core contributes to the shared spacetime geometry.
The authors use the point where the central pressure diverges to mark a Buchdahl-like critical boundary. The numerical analysis identifies a region with finite positive central pressure and real metric functions, while the critical curve marks central-pressure divergence. It is therefore a boundary of the static solutions examined by the model.
The one-fluid check and the critical curve
The construction passes a basic internal check: in the one-fluid limit, it recovers the Schwarzschild constant-density star and the critical compactness 2M/R = 8/9. Within this framework, that recovery ties the two-fluid calculation to the corresponding constant-density one-fluid case.
For the mixed-core cases, the critical condition is evaluated numerically after requiring the metric functions and the relevant integral to remain real and well defined. The reported critical branch is approximately linear; that description is an approximation rather than an exact divergence condition. The approximation improves as α increases across the values examined: α = 2, α = 2.5, and α = 3.33, with the relative dark-matter density parameter f examined over 0 < f < 1.
The numerical sequences show another consistent trend. At fixed α, the critical inner compactness Ci is lower at larger relative dark-matter density f. The result describes how the model’s pressure-divergence boundary varies across the internal-composition parameters used in the calculation.
The same outside, different inside
The mass-radius relation gives the model another way to distinguish global properties from internal structure. For fixed ordinary-matter density ρo, fixed f, and fixed α, the relation scales cubically with the stellar radius. It also contains the second-fluid contribution through the factor f α−3.
That dependence leads to a central interpretive point: different combinations of Ci, f, and α can produce the same global compactness while corresponding to different internal matter distributions. A single overall compactness in this model therefore does not uniquely identify how the two components are arranged inside the star.
The paper presents this degeneracy, together with the pressure and matching results, as a reason to use the construction as an analytical benchmark. It isolates the gravitational effects of a centrally concentrated second component in a form that can be compared with more detailed models later.
What the calculation cannot establish
The critical curve should not be read as a stability limit. In the paper, it marks divergence of the central pressure; it is not a demonstrated dynamical-stability or marginal-stability boundary. Determining an actual stability boundary would require a coupled radial-perturbation analysis, which is outside the calculation reported here.
The model is also deliberately idealized. It uses constant energy densities, isotropic perfect fluids, static spherical symmetry, and an abrupt termination of the dark component at the core boundary. The resulting mass-radius curves are not presented as realistic neutron-star sequences, and the framework does not supply microscopic equations of state.
For the same reason, the preprint offers no observational predictions, empirical validation, or result for tidal deformability. The open question is whether its qualitative trends survive when realistic microscopic equations of state, coupled stability calculations, and extensions that examine observables are added.
A benchmark, not a detection
The document is an arXiv version 1 preprint dated 20 Aug 2026. Its contribution is a mathematical reference case: a two-fluid star whose core composition, interface conditions, pressure behavior, and mass-radius scaling can be followed analytically and numerically. Its conclusions describe that model’s static solutions, not the properties of observed stars or evidence for dark matter inside them.
Paper data and sources
Original title: An Analytical Two Incompressible Fluid Star with a Mixed Ordinary Dark Matter Core and an Ordinary Matter Envelope
Authors: Milko Estrada, Santiago Esteban Perez Bergliaffa
Journal/Repository: arXiv
Status: Preprint, not yet peer-reviewed
First online: 2026-08-20
DOI: Not available
Original paper · Full text