Dark Energy Stars: Models and Implications

• Dark Energy Stars (DESs) are hypothetical compact objects whose interiors are defined by a dark-energy-like equation of state, replacing ordinary nuclear matter.
• They are modeled through various approaches including anisotropic solutions, Chaplygin-type fluids, and quantum droplet frameworks to counteract gravitational collapse.
• DES models link relativistic stellar structure with dark energy phenomena, providing potential explanations for unusual mass-radius relations and dark matter contributions.

Searching arXiv for the cited dark energy star papers to ground the article in current records. Dark energy stars (DESs) are hypothetical compact objects whose interior stress-energy is governed, wholly or partly, by a dark-energy-like equation of state rather than by ordinary nuclear matter alone. In the arXiv literature, the term spans several related constructions: horizonless ultracompact stars with a core satisfying $p=-\rho$, Chaplygin-type self-bound compact stars in general relativity, anisotropic gravastar-like configurations, and, in a distinct cosmological proposal, quantum droplets associated with vacuum energy and dark matter (Beltracchi et al., 2018). Across these variants, the unifying idea is that sufficiently negative pressure can provide an effective repulsive contribution that counteracts collapse, alters the compactness-redshift relation, and may produce regular interiors without the standard black-hole end state (Banerjee et al., 2019).

1. Conceptual scope and historical placement

The DES concept occupies the intersection of relativistic stellar structure, dark-energy phenomenology, and black-hole alternatives. One line of work treats DESs as finite-size astrophysical objects with an interior equation of state typical of dark energy, especially a central region with $p_r=p_T=-\rho=\text{constant}$, placing them in the same broad family as false vacuum bubbles, vacuum nonsingular black holes, and gravastars (Beltracchi et al., 2018). In that usage, a DES is primarily a compact-object solution of Einstein’s equations with strong negative pressure in the interior and, frequently, anisotropy in the transition region between the core and the outer layers.

A second line of work models DESs as static relativistic stars supported by dark-energy-inspired fluids, especially generalized, extended, or modified Chaplygin equations of state. In these constructions the object is usually treated as a compact star with a finite radius $R$, a Schwarzschild exterior, and a nonzero surface density when $p(R)=0$, so the star is self-bound rather than crust-terminated in the ordinary hadronic sense (Panotopoulos et al., 2020).

A third, conceptually distinct proposal identifies DESs with gravitationally stable quantum droplets whose interiors resemble vacuum with much higher vacuum energy density. In that framework DESs are not merely stellar alternatives but constituents of dark matter and intermediate states in an early-universe scenario linking a high-vacuum-energy cosmic seed to a Friedmann-like universe of radiation plus residual dark matter (Chapline, 2010).

These usages are not equivalent. Some papers study DESs as equilibrium compact stars in standard GR, some as exact anisotropic solutions sourced by a phantom field, and some as effective objects in modified gravity. This suggests that “dark energy star” functions in current literature less as a single sharply delimited object class than as a family of negative-pressure compact configurations.

2. Matter models and equilibrium equations

Most DES models assume static spherical symmetry with line element

$$ds^2=-e^{\nu(r)}dt^2+e^{\lambda(r)}dr^2+r^2(d\theta^2+\sin^2\theta\,d\phi^2),$$

or equivalent Schwarzschild-like parametrizations, and employ either isotropic perfect fluids or anisotropic fluids with $p_r\neq p_t$ (Banerjee et al., 2019). The basic equilibrium system is the Tolman–Oppenheimer–Volkoff structure, generalized when necessary to include anisotropy,

$$\frac{dm}{dr}=4\pi r^2\rho, \qquad \frac{dp_r}{dr}=-(\rho+p_r)\left(\frac{m}{r^2}+4\pi r p_r\right)\left(1-\frac{2m}{r}\right)^{-1}+\frac{2\sigma}{r},$$

with $\sigma\equiv p_t-p_r$ in anisotropic models (Pretel, 2023).

Several matter prescriptions recur in the literature:

Paper	Matter model	Characteristic feature
(Banerjee et al., 2019)	$p_r=\omega\rho$, $\omega<0$	Finch–Skea exact anisotropic interior
(Panotopoulos et al., 2020)	$p_r=p_T=-\rho=\text{constant}$0	isotropic generalized Chaplygin DES
(Panotopoulos et al., 2021)	$p_r=p_T=-\rho=\text{constant}$1	isotropic extended Chaplygin DES
(Pretel, 2023)	$p_r=p_T=-\rho=\text{constant}$2	anisotropic Chaplygin-type DES
(Das et al., 24 Mar 2026)	$p_r=p_T=-\rho=\text{constant}$3	modified Chaplygin DES and universal relations

In the exact Finch–Skea model the interior metric potential is fixed by

$p_r=p_T=-\rho=\text{constant}$4

which yields

$p_r=p_T=-\rho=\text{constant}$5

together with a nontrivial tangential pressure $p_r=p_T=-\rho=\text{constant}$6 and anisotropy $p_r=p_T=-\rho=\text{constant}$7. The center is regular, with $p_r=p_T=-\rho=\text{constant}$8 and $p_r=p_T=-\rho=\text{constant}$9 (Banerjee et al., 2019).

In Chaplygin-type models the surface is defined by $R$0, but the density remains finite. For the generalized or extended Chaplygin forms this gives

$R$1

a standard indicator of self-bound behavior rather than a density profile tapering continuously to zero at the surface (Panotopoulos et al., 2020). In the modified Chaplygin gas model,

$R$2

and causality is enforced through $R$3 (Das et al., 24 Mar 2026).

3. Internal geometry, anisotropy, and formation channels

Anisotropy is central to many DES constructions. In the exact anisotropic Finch–Skea model, the sign of the local gravitational acceleration

$R$4

determines whether gravity is attractive or repulsive. For $R$5, $R$6, and the interior behaves in the repulsive regime usually identified as DES-like (Banerjee et al., 2019). In the same interval the anisotropy becomes positive, $R$7, so the anisotropic force is outward-directed and contributes to support against collapse.

In the anisotropic Chaplygin-type GR model, pressure anisotropy is introduced through the Horvat prescription

$R$8

which vanishes at the center and disappears from hydrostatic equilibrium in the Newtonian limit. Within the Iyer et al. classification, isotropic DESs in this setup are ordinary compact stars, whereas sufficiently large positive $R$9 can drive configurations into the ultra-compact regime (Pretel, 2023).

A distinct anisotropy mechanism arises in the phantom-field exact solution generated from the Schwarzschild interior solution. There the phantom scalar contributes only to the radial sector,

$p(R)=0$0

so the anisotropy factor

$p(R)=0$1

is positive and repulsive. However, the paper emphasizes that full strong-energy-condition violation throughout the star occurs only at the Buchdahl limit $p(R)=0$2, so the phantom component does not dominate the whole interior for $p(R)=0$3 (Sakti et al., 2021).

DES formation has also been modeled dynamically rather than imposed as a static ansatz. In the time-dependent collapse solution, a spherical configuration evolves from an ordinary positive-pressure precursor to a final state containing a dark-energy core with $p(R)=0$4. The transition necessarily passes through an anisotropic inversion zone, because a continuous pressure profile cannot jump directly from $p(R)=0$5 outside to $p(R)=0$6 in the core. In that model the mass remains fixed, no thin shells are introduced, and collapse halts before the surface reaches $p(R)=0$7, so the spacetime develops neither a singularity nor an event horizon (Beltracchi et al., 2018).

Matching to an exterior Schwarzschild geometry is standard in stellar DES models. The exact Finch–Skea construction performs the junction at $p(R)=0$8 using the Israel–Lanczos thin-shell formalism, with surface stress tensor

$p(R)=0$9

where $$ds^2=-e^{\nu(r)}dt^2+e^{\lambda(r)}dr^2+r^2(d\theta^2+\sin^2\theta\,d\phi^2),$$0 is the surface energy density and $$ds^2=-e^{\nu(r)}dt^2+e^{\lambda(r)}dr^2+r^2(d\theta^2+\sin^2\theta\,d\phi^2),$$1 the surface pressure (Banerjee et al., 2019).

4. Stability, energy conditions, and oscillation spectra

DES stability is strongly model-dependent. In the generalized Chaplygin isotropic model, causality, adiabatic stability, and the standard energy conditions are satisfied throughout the star for the studied parameter sets. The sound speed obeys

$$ds^2=-e^{\nu(r)}dt^2+e^{\lambda(r)}dr^2+r^2(d\theta^2+\sin^2\theta\,d\phi^2),$$2

the adiabatic index satisfies $$ds^2=-e^{\nu(r)}dt^2+e^{\lambda(r)}dr^2+r^2(d\theta^2+\sin^2\theta\,d\phi^2),$$3, and the ten lowest radial oscillation modes have positive frequencies, implying dynamical stability against radial collapse (Panotopoulos et al., 2020).

The anisotropic Chaplygin-type GR analysis sharpens this picture by explicitly relating the turning point of the equilibrium sequence to radial instability. The conventional criterion

$$ds^2=-e^{\nu(r)}dt^2+e^{\lambda(r)}dr^2+r^2(d\theta^2+\sin^2\theta\,d\phi^2),$$4

selects the stable branch, and at the maximum-mass configuration the squared frequency of the fundamental radial mode satisfies

$$ds^2=-e^{\nu(r)}dt^2+e^{\lambda(r)}dr^2+r^2(d\theta^2+\sin^2\theta\,d\phi^2),$$5

This reproduces the standard turning-point interpretation familiar from neutron-star theory, but now for anisotropic DESs (Pretel, 2023).

The exact Finch–Skea model evaluates viability through NEC, WEC, SEC, and DEC, as well as compactness and surface redshift. For the representative configuration Vela X-1 with $$ds^2=-e^{\nu(r)}dt^2+e^{\lambda(r)}dr^2+r^2(d\theta^2+\sin^2\theta\,d\phi^2),$$6, the plotted combinations satisfy these energy conditions throughout the interior, the compactness obeys the Buchdahl bound $$ds^2=-e^{\nu(r)}dt^2+e^{\lambda(r)}dr^2+r^2(d\theta^2+\sin^2\theta\,d\phi^2),$$7, and the surface redshift remains finite, typically below unity (Banerjee et al., 2019).

By contrast, the phantom-field DES yields a more cautionary result. Although it exhibits dark-energy-like features and positive anisotropy, it violates the causality conditions

$$ds^2=-e^{\nu(r)}dt^2+e^{\lambda(r)}dr^2+r^2(d\theta^2+\sin^2\theta\,d\phi^2),$$8

and is not stable against gravitational cracking, assessed through

$$ds^2=-e^{\nu(r)}dt^2+e^{\lambda(r)}dr^2+r^2(d\theta^2+\sin^2\theta\,d\phi^2),$$9

The paper therefore treats this DES as an ultra-compact toy model rather than as a physically stable compact star (Sakti et al., 2021).

The same instability analysis intersects directly with echo phenomenology. For ultra-compact configurations with $p_r\neq p_t$0, the phantom-field DES admits a photon sphere and an echo time

$p_r\neq p_t$1

The phantom contribution increases the echo time and lowers the echo frequency relative to the constant-density Schwarzschild interior solution, which the authors attribute to a deeper effective potential well (Sakti et al., 2021).

5. Rotation, tidal response, and universal relations

Slow rotation in DESs is usually treated within Hartle or Hartle–Thorne perturbation theory. In the extended Chaplygin isotropic model, the metric acquires a frame-dragging term

$p_r\neq p_t$2

and, after harmonic decomposition, only the $p_r\neq p_t$3 sector survives asymptotically. The resulting moment of inertia is

$p_r\neq p_t$4

with $p_r\neq p_t$5. The reported trends are that $p_r\neq p_t$6 increases with stellar mass, grows faster in the non-rotating reference sequence, and is smaller for the rotating star than for the non-rotating one at fixed mass (Panotopoulos et al., 2021).

The fully anisotropic Hartle–Thorne study with modified Chaplygin fluid extends this by keeping both monopole and quadrupole deformations. There the moment of inertia becomes

$p_r\neq p_t$7

and the quadrupole moment is written as

$p_r\neq p_t$8

For the angular frequencies considered, anisotropy affects $p_r\neq p_t$9, radius, $$\frac{dm}{dr}=4\pi r^2\rho, \qquad \frac{dp_r}{dr}=-(\rho+p_r)\left(\frac{m}{r^2}+4\pi r p_r\right)\left(1-\frac{2m}{r}\right)^{-1}+\frac{2\sigma}{r},$$0, $$\frac{dm}{dr}=4\pi r^2\rho, \qquad \frac{dp_r}{dr}=-(\rho+p_r)\left(\frac{m}{r^2}+4\pi r p_r\right)\left(1-\frac{2m}{r}\right)^{-1}+\frac{2\sigma}{r},$$1, $$\frac{dm}{dr}=4\pi r^2\rho, \qquad \frac{dp_r}{dr}=-(\rho+p_r)\left(\frac{m}{r^2}+4\pi r p_r\right)\left(1-\frac{2m}{r}\right)^{-1}+\frac{2\sigma}{r},$$2, $$\frac{dm}{dr}=4\pi r^2\rho, \qquad \frac{dp_r}{dr}=-(\rho+p_r)\left(\frac{m}{r^2}+4\pi r p_r\right)\left(1-\frac{2m}{r}\right)^{-1}+\frac{2\sigma}{r},$$3, and $$\frac{dm}{dr}=4\pi r^2\rho, \qquad \frac{dp_r}{dr}=-(\rho+p_r)\left(\frac{m}{r^2}+4\pi r p_r\right)\left(1-\frac{2m}{r}\right)^{-1}+\frac{2\sigma}{r},$$4 more strongly than rotation itself; moreover, larger anisotropic strength moves $$\frac{dm}{dr}=4\pi r^2\rho, \qquad \frac{dp_r}{dr}=-(\rho+p_r)\left(\frac{m}{r^2}+4\pi r p_r\right)\left(1-\frac{2m}{r}\right)^{-1}+\frac{2\sigma}{r},$$5 closer to the Kerr value $$\frac{dm}{dr}=4\pi r^2\rho, \qquad \frac{dp_r}{dr}=-(\rho+p_r)\left(\frac{m}{r^2}+4\pi r p_r\right)\left(1-\frac{2m}{r}\right)^{-1}+\frac{2\sigma}{r},$$6 (Jyothilakshmi et al., 2024).

Tidal response is another major discriminator. In the isotropic generalized Chaplygin model,

$$\frac{dm}{dr}=4\pi r^2\rho, \qquad \frac{dp_r}{dr}=-(\rho+p_r)\left(\frac{m}{r^2}+4\pi r p_r\right)\left(1-\frac{2m}{r}\right)^{-1}+\frac{2\sigma}{r},$$7

and the dimensionless deformability decreases as compactness increases. The paper concludes that DESs occupy a distinct region in the deformability-versus-compactness plane and may therefore be distinguishable from ordinary neutron stars or quark stars if binary measurements become sufficiently precise (Panotopoulos et al., 2020).

The anisotropic Chaplygin-type GR model shows that the Love number $$\frac{dm}{dr}=4\pi r^2\rho, \qquad \frac{dp_r}{dr}=-(\rho+p_r)\left(\frac{m}{r^2}+4\pi r p_r\right)\left(1-\frac{2m}{r}\right)^{-1}+\frac{2\sigma}{r},$$8 rises with compactness up to a maximum and then decreases, while positive anisotropy generally reduces the maximum Love number and modifies the surface redshift and moment of inertia most strongly in the high-mass branch (Pretel, 2023).

A 2026 modified Chaplygin gas analysis systematizes these macroscopic observables into the $$\frac{dm}{dr}=4\pi r^2\rho, \qquad \frac{dp_r}{dr}=-(\rho+p_r)\left(\frac{m}{r^2}+4\pi r p_r\right)\left(1-\frac{2m}{r}\right)^{-1}+\frac{2\sigma}{r},$$9-$\sigma\equiv p_t-p_r$0-$\sigma\equiv p_t-p_r$1-$\sigma\equiv p_t-p_r$2 universal relations. Using

$\sigma\equiv p_t-p_r$3

the paper finds that DESs obey very tight relations among compactness, moment of inertia, tidal deformability, and $\sigma\equiv p_t-p_r$4-mode frequency, but that these relations are very similar to those of quark stars. The key new claim is that gravitational binding energy,

$\sigma\equiv p_t-p_r$5

breaks this degeneracy through the $\sigma\equiv p_t-p_r$6, $\sigma\equiv p_t-p_r$7, and $\sigma\equiv p_t-p_r$8 relations (Das et al., 24 Mar 2026).

6. Astrophysical applications and observational status

DES models have been compared with several categories of data: heavy pulsar masses, gravitational-wave tidal constraints, NICER-like radius measurements, and the existence of compact objects in the mass gap. The exact Finch–Skea model was motivated partly by objects whose masses and radii were stated to be difficult to reconcile with standard neutron-star equations of state, including PSR J1416-2230, Vela X-1, 4U 1608-52, Her X-1, and PSR J1903+327, with the overall sequence reported to be consistent with masses near $\sigma\equiv p_t-p_r$9 and with the Buchdahl limit (Banerjee et al., 2019).

The isotropic generalized Chaplygin study obtains neutron-star-like mass-radius curves, with radii around $p_r=\omega\rho$0–$p_r=\omega\rho$1 km for masses near $p_r=\omega\rho$2–$p_r=\omega\rho$3, and compactness around $p_r=\omega\rho$4, safely below $p_r=\omega\rho$5. These configurations are presented as regular, physically admissible, and potentially distinguishable through tidal measurements and radial seismology rather than through bulk mass-radius data alone (Panotopoulos et al., 2020).

The anisotropic Chaplygin-type GR model was explicitly applied to GW190814. In the isotropic case, some models with $p_r=\omega\rho$6 and $p_r=\omega\rho$7 were found consistent with the secondary component. With anisotropy the parameter space broadens, and the paper states that the GW190814 secondary can be consistently described as a stable anisotropic DES, especially for $p_r=\omega\rho$8 and $p_r=\omega\rho$9 or $\omega<0$0 in the chosen $\omega<0$1-range (Pretel, 2023).

The slow-rotation anisotropic study with Bowers–Liang prescription further reports consistency with GW170817, GW190814, and massive pulsars such as PSR J2215+5135 and PSR J0740+6620. In that analysis all three equation-of-state sets satisfy the GW170817 tidal bound $\omega<0$2, while $\omega<0$3 configurations satisfy GW190814 across all sets (Jyothilakshmi et al., 2024).

The modified Chaplygin universal-relation study uses the GW170817 deformability constraint to infer canonical DES properties. Adopting the model-independent bound $\omega<0$4, it obtains

$\omega<0$5

$\omega<0$6

These inferences are specific to the DES model family under consideration, not model-independent statements about compact stars in general (Das et al., 24 Mar 2026).

7. Cosmological and modified-gravity extensions, and persistent controversies

Not all DES literature is confined to compact-star phenomenology in GR. In modified Rastall teleparallel gravity, DESs are modeled as compact objects containing ordinary baryonic matter plus a dark-energy sector with

$\omega<0$7

and the teleparallel-Rastall structure is encoded through

$\omega<0$8

Within this framework the dark-energy sector has $\omega<0$9, $p_r=p_T=-\rho=\text{constant}$00, and $p_r=p_T=-\rho=\text{constant}$01, while the effective total configuration is reported to satisfy equilibrium, possess finite redshift and admissible compactness, and remain stable except for a noted instability at $p_r=p_T=-\rho=\text{constant}$02 (Ditta et al., 2024).

The cosmological DES proposal is more radical. There, a finite “cosmic seed” with vacuum energy density $p_r=p_T=-\rho=\text{constant}$03 becomes unstable because it exceeds the de Sitter horizon scale and fractures through near-horizon quantum critical fluctuations into primordial DES-like droplets of characteristic mass

$p_r=p_T=-\rho=\text{constant}$04

These droplets later coalesce, transfer most of their mass-energy into radiation, and leave residual dark matter clumps with characteristic mass scale

$p_r=p_T=-\rho=\text{constant}$05

The same scenario is proposed to account jointly for the present dark matter density, a CMB origin near $p_r=p_T=-\rho=\text{constant}$06, and fluctuation amplitudes connected to $p_r=p_T=-\rho=\text{constant}$07 after renormalization of initially large fluctuations (Chapline, 2010).

The principal controversies surrounding DESs arise from this heterogeneity of definitions and outcomes. Stability claims are not universal: some GR Chaplygin-type models satisfy causality, energy conditions, and radial-mode stability, whereas the phantom-field exact solution violates causality and cracking criteria and is used only as an ultra-compact toy model (Panotopoulos et al., 2020). Likewise, the intended ontology varies sharply: DESs can be black-hole alternatives, self-bound compact stars, gravastar-like end states of collapse, or quantum droplets constituting dark matter (Sakti et al., 2021).

For that reason, the DES literature is best read as a collection of mathematically and physically distinct negative-pressure compact-object programs. Their common theme is the replacement of ordinary high-density matter, or of the black-hole interior, by a dark-energy-like medium. Their differences concern the microphysics, the role of anisotropy, the interpretation of stability, and the degree to which present observations can discriminate them from neutron stars, quark stars, gravastars, or black holes.