---
title: Frozen Star Gravitational Configurations
url: https://www.emergentmind.com/topics/frozen-star
type: topic
---

# Frozen Star Gravitational Configurations

Frozen star denotes a family of gravitational configurations united by an extreme redshift effect but not by a single underlying model. In the older relativistic-collapse literature, the term referred to the Oppenheimer–Snyder result that a collapsing surface asymptotically approaches the Schwarzschild radius in the external Schwarzschild time coordinate, appearing to “freeze” there [1003.1359]. In recent compact-object literature, the same term is used for ultracompact configurations in which the metric redshift factor approaches zero at a finite radius, so that matter and clocks near that radius appear frozen to infinity, while no true event horizon forms. This latter usage includes horizonless solutions in Einstein–Hayward–scalar theory [2312.07224], anisotropic-fluid and string-fluid constructions [2404.15985, 2109.10017], higher-curvature boson and Proca stars [2406.08813, 2512.24584], dark-matter-sourced compact stars [2601.15415], and neutron-star sequences with frozen endpoints [2512.23525, 2509.09338].

## 1. Historical and terminological development

The classical “frozen star” paradox originates in the external description of gravitational collapse. In the Oppenheimer–Snyder dust model, the exterior is Schwarzschild and the surface radius satisfies an integral for the Schwarzschild coordinate time that diverges as the radius approaches \(R_H=2M\); the comoving observer sees collapse complete in finite proper time, but the external observer never sees the surface cross the Schwarzschild radius [1003.1359]. Liu and Zhang generalized this setup to shells of finite thickness falling onto a pre-existing black hole and, crucially, included the effect of changing total gravitating mass. In that treatment the event horizon grows and can swallow shell elements in finite external time, so the enduring “frozen star” picture does not survive once the back-reaction of infalling matter and outer shells is included [1003.1359].

A different historical use appears in Janis–Newman–Winicour spacetimes. There the metric has an infinite-redshift surface at \(r=r_0\), with
\[
z+1=f(r_e)^{-\alpha/2},\qquad f(r)=1-\frac{r_0}{r},
\]
so that clocks near \(r_0\) appear frozen to infinity, even though \(r=r_0\) is a scalar singularity rather than a regular horizon [1607.08219]. Cosmological generalizations retain this infinite-redshift property in expanding or contracting backgrounds [1607.08219].

Modern compact-object work uses the term in a more specific sense: a frozen star is typically a regular or nonsingular ultracompact object that mimics a black hole from the exterior while remaining horizonless. This suggests that the term has bifurcated into two conceptually distinct usages: a coordinate-time description of collapse, and a class of stationary or limiting black-hole mimickers.

## 2. Geometric definition and horizon criterion

A common metric ansatz in contemporary frozen-star constructions is
\[
ds^2=-N(r)\sigma^2(r)\,dt^2+\frac{dr^2}{N(r)}+r^2d\Omega^2,
\]
or an equivalent Schwarzschild-like form with \(f(r)=g(r)\) [2312.07224, 2601.15415, 2512.23525]. The defining geometric feature is that the redshift factor
\[
g_{tt}(r)=-N(r)\sigma^2(r)
\]
or \(-f(r)\) approaches zero arbitrarily closely at some finite radius \(r=r_c\), \(r=r_*\), or \(r=r_{cH}\), without changing sign and without producing a true event horizon [2601.15415, 2512.23525, 2503.16265].

In the Einstein–Hayward–scalar model, an event horizon would appear if \(N(r_h)=0\) for some finite \(r_h>0\), and the frozen Hayward–boson-star analysis restricts attention to star-like solutions with \(N(r)>0\) for all \(r>0\) [2312.07224]. In dark-matter-sourced compact stars, the minimum of \(N(r)\) approaches zero from above, \(\min N(r)\to0^+\), yet the solution remains horizonless unless one exactly tunes to the regular black-hole case [2601.15415]. In four-dimensional non-polynomial gravities, the frozen endpoint appears when the metric function develops a double root extremely close to the stellar surface, with \(f(r_*)\to0\) and \(g(r_*)=f(r_*)\to0\) [2512.23525]. In five-dimensional higher-derivative Proca stars, the frozen state develops a critical horizon at finite radius \(r_c\), where \(-g_{tt}\) and \(1/g_{rr}\) approach zero, but the configuration contains neither curvature singularities nor event horizons [2512.24584].

A distinct but closely related implementation appears in the anisotropic-fluid frozen-star model of Brustein and collaborators, where the interior is regularized by a small parameter \(\varepsilon\) and
\[
ds^2_{\rm int}=-\varepsilon^2 dt^2+\frac{dr^2}{\varepsilon^2}+r^2d\Omega^2,
\qquad \varepsilon^2\ll1,
\]
so the entire interior has an arbitrarily large but finite redshift and matches to a Schwarzschild exterior at \(r=R_s\) [2404.15985, 2301.09712].

## 3. Source sectors and model classes

No single matter model defines the frozen-star concept. Recent realizations span anisotropic fluids, bosonic and fermionic fields, nonlinear electromagnetic sectors, string fluids, dark-matter halos, and higher-curvature gravities.

In the anisotropic-fluid constructions, the central equation of state is maximally negative radial pressure,
\[
p_r=-\rho,
\]
often together with \(p_\perp=0\) or perturbatively small deviations from it [2109.10017, 2301.09712, 2404.15985]. In these models Einstein’s equations admit a nearly null interior redshift factor, and the condition \(\rho+p_r=0\) is repeatedly identified as the mechanism behind both the near-horizon behavior and the freeze-out of linear fluctuations [2109.10017].

The string-fluid realization identifies the matter source with the string fluid resulting from the decay of an unstable \(D\)-brane or brane–antibrane system at the end of open-string tachyon condensation [2404.15985]. The effective Lagrangian can be recast into a Born–Infeld form, the interior electric displacement obeys
\[
D(r)=\frac{\alpha'}{4Gr^2},
\]
and the stress tensor has
\[
\rho(r)=\frac{1-\varepsilon^2}{8\pi G r^2},\qquad p_r(r)=-\rho(r),\qquad p_\perp(r)=0
\]
inside the star [2404.15985].

A second large class couples matter fields to regular-black-hole sectors of Bardeen or Hayward type. The Einstein–Hayward–scalar theory minimally couples a complex, massive, free scalar field to Hayward nonlinear electrodynamics and yields ground and excited frozen Hayward–boson stars [2312.07224]. Closely related Bardeen constructions include frozen Bardeen–boson stars, frozen Bardeen–Dirac stars, and frozen Bardeen–Proca stars, all obtained by adding scalar, Dirac, or Proca matter to the nonlinear electromagnetic background so that an event horizon is disrupted and replaced by a critical surface or critical horizon [2503.14984, 2312.07400, 2503.16265].

Further generalizations place the phenomenon in modified gravity and ordinary-matter settings. Four-dimensional non-polynomial gravities produce frozen neutron stars at the end of the neutron-star sequence for the BSk19, SLy4, and AP4 equations of state [2512.23525]. Einstein–nonlinear-electrodynamics models with Bardeen or Hayward magnetic monopoles yield frozen neutron stars at a critical magnetic charge \(q_c\) [2509.09338]. Dark-matter halo models produce horizonless compact stars that approach a frozen state when the condition \(P_r=-\rho\) is relaxed to broader anisotropic families [2601.15415]. Five-dimensional boson and Proca stars under infinite towers of higher-curvature corrections also admit frozen states when the coupling exceeds a threshold [2406.08813, 2512.24584].

## 4. Limiting behavior and representative solution families

The frozen limit usually arises by tuning a control parameter to a critical regime. In boson and Proca stars, the relevant parameter is often the matter frequency \(\omega\to0\). In magnetically charged compact stars, it is frequently a critical magnetic charge. In higher-curvature or non-polynomial gravities, it is a coupling strength together with the central density.

In the Einstein–Hayward–scalar system, each radial excitation \(n=0,1,2,\dots\) has a branch of solutions labeled by \(\omega\). For \(q<q_{\rm crit}(n)\), the mass–frequency curve is multivalued and spirals as in ordinary boson stars. As \(q\to q_{\rm crit}(n)\), the spiral opens, \(\omega_{\min}\to0\), and a single-valued branch extends to \(\omega=0\) [2312.07224]. In that limit the scalar profile develops a very steep wall at a critical radius \(r_c\), \(N(r_c)\sim0\), and \(g_{tt}(r<r_c)\approx0\), yet no true horizon forms [2312.07224]. Numerically, the Hayward case gives
\[
q_{\rm crit}^{\rm ground}\approx0.494,\qquad
q_{\rm crit}^{\rm 1st\,excited}\approx0.514,
\]
and, remarkably, the critical radius \(r_c\) and critical mass \(M_{\rm crit}\) at \(\omega=0\) are independent of the radial node number \(n\) for fixed \(q_{\rm crit}\) [2312.07224].

Bardeen-based models show analogous limiting behavior with model-specific details. Frozen Bardeen–Dirac stars exist for \(q_c<q<q_b\), and for \(s=0.3\) the Dirac construction gives \(q_c\approx0.746\); the metric outside the critical horizon tracks that of an extremal Bardeen black hole, while the interior remains nonsingular [2312.07400]. In the Bardeen–Proca case, sufficiently large magnetic charge allows \(\omega\to0\), \(N(r_{cH})\to0\), \(\sigma(r_{cH})\to0\), and \(g_{tt}\to0\) inside \(r_{cH}\), while the ADM mass remains finite and the Noether charge tends to zero [2503.16265]. In the Bardeen–boson-star null-geodesic analysis, the frozen configurations are the \(q>q_c\) solutions whose metric becomes very similar to that of black holes near the critical horizon, while remaining fully horizonless and possessing complete null geodesics [2503.14984].

In non-polynomial gravity, increasing the modification parameter \(\alpha\) enlarges both neutron-star radius and mass, and a frozen state emerges when the metric functions approach zero extremely close to the stellar surface [2512.23525]. The allowed observational window is roughly \(\alpha\in[10^7,10^8]\,{\rm m}^2\), and within this range frozen-star solutions appear for BSk19, SLy4, and AP4 [2512.23525]. In the model \(h=\psi/(1-\alpha^2\psi^2)\), the quoted thresholds are
\[
\alpha_{\min}\sim3.4\times10^7\,{\rm m}^2\ {\rm for\ BSk19},\quad
\alpha_{\min}\sim4.4\times10^7\,{\rm m}^2\ {\rm for\ SLy4},\quad
\alpha_{\min}\sim6.1\times10^7\,{\rm m}^2\ {\rm for\ AP4},
\]
above which frozen neutron stars appear within the causal range of central densities [2512.23525]. In dark-matter-sourced compact stars, the frozen regime is reached when the minimum of \(N(r)\) approaches zero from above; for the Einasto-profile Case I, the critical densities are \(\rho_0^{\rm crit}\approx0.68391\) for \(n=1/2\), \(\rho_0^{\rm crit}\approx0.56899\) for \(n=3/4\), and \(\rho_0^{\rm crit}\approx0.40977\) for \(n=1\) [2601.15415].

## 5. Perturbations, thermodynamics, and formation

One prominent claim in the anisotropic-fluid frozen-star program is that the exact frozen configuration is dynamically bald. In the matched interior–crust–exterior construction with \(p_r=-\rho\), linearized analysis gives \(H_2=0\), \(\delta\rho=0\), and, after imposing the inner boundary conditions, \(H_0=0\), so that every metric perturbation and every matter perturbation vanishes identically at linear order [2109.10017]. This is the strongest available formulation of “freezing out” in the perturbative sense.

When the equation of state is perturbed away from maximal negativity, the star “defrosts.” In the Cowling treatment of internal fluid modes, a small parameter \(\varepsilon\) yields a leading-order dispersion relation
\[
\omega(\varepsilon,\ell)=\sqrt{\frac{2\varepsilon}{3R^2}\,\bigl[\ell(\ell+1)-4\bigr]},
\]
so the real frequencies vanish in the strict frozen limit \(\varepsilon\to0\) [2304.04984]. With even-parity metric perturbations included, the defrosting parameter \(\gamma\) produces non-radial oscillations whose sound speed is non-relativistic and proportional to \(\gamma\), while the lifetime is parametrically long and proportional to \(1/\gamma^2\) [2410.00493]. These results were presented as a starting point for calculating the spectrum of emitted gravitational waves from an excited frozen star [2410.00493].

Thermodynamically, frozen stars have been argued to reproduce black-hole behavior to leading order. Using a Euclidean-action calculation, the entropy is
\[
S=\frac{A}{4}+{\cal O}(\varepsilon^2,\lambda),
\]
and the temperature is
\[
T_{FS}=\frac{1}{8\pi M}+{\cal O}(\varepsilon^2,\lambda),
\]
so both are perturbatively close to the Bekenstein–Hawking values [2310.11572]. The same work estimated a transition probability \(e^{-A/4}\) for a collapsing shell to tunnel into a frozen star, and then argued that the phase-space factor \(e^{+A/4}\) makes the total formation probability of order unity [2310.11572]. A later tunneling analysis sharpened that conclusion, stating that the transition probability is unity, up to negligible corrections, because the Euclidean-action difference is approximately zero and the determinant prefactor is unity [2508.02100].

## 6. Exterior mimicry, observables, and unresolved issues

Across the modern literature, frozen stars are presented as black-hole mimickers. In the string-fluid model, the exterior is exactly Schwarzschild and for \(\varepsilon^2\lesssim10^{-22}\) it is stated to be indistinguishable from Schwarzschild under Event Horizon Telescope bounds [2404.15985]. In frozen Hayward, Bardeen, boson, Dirac, Proca, and higher-curvature models, the exterior geometry is repeatedly said to coincide with, or become nearly indistinguishable from, the corresponding extremal black hole outside the critical radius [2312.07224, 2312.07400, 2512.24584].

The proposed discriminants are correspondingly subtle. Several models predict light-ring structures resembling black holes but not identical to them. Frozen Bardeen–Proca stars have two light rings, one stable inside and one unstable outside the critical horizon, and these move further apart as \(\omega\) decreases [2503.16265]. Frozen Bardeen–Dirac stars also have an outer unstable light ring, but inside the critical horizon the photon velocity becomes very close to zero, leading to a “light ball” [2312.07400]. In frozen Bardeen–boson stars, photons undergo sharp deflections near the critical horizon, travel nearly straight inside it, and can spend a very long coordinate time in the interior as seen from infinity [2503.14984].

Other proposed observables include gravitational-wave echoes, late-time ringdown modifications, nonzero tidal deformability, and accretion signatures. The dark-matter compact-star model states that frozen stars mimic black-hole shadows, ringdown spectra, and strong lensing, yet may exhibit reflectivity, echoes, and slight shifts in quasinormal frequencies; it also states that linearized axial perturbations are stable for viable parameter ranges [2601.15415]. The non-polynomial-gravity neutron-star analysis identifies late-time ringdown, echoes, tidal deformability \(\Lambda\), and accretion-disk and iron-\(K\alpha\) diagnostics as possible tests, while emphasizing that the frozen configurations remain allowed by present NICER and GW170817 constraints [2512.23525].

Several central issues remain open in the cited literature. The five-dimensional higher-derivative Proca-star work states explicitly that linear stability remains to be fully studied [2512.24584]. The frozen-neutron-star analysis likewise notes that a rigorous stability analysis under radial perturbations is essential and that electromagnetic or gravitational-wave imprints remain to be thoroughly investigated [2509.09338]. This suggests that the current status of frozen stars is strongest at the level of existence proofs, limiting geometries, and kinematical phenomenology, and less settled at the level of full dynamical selection, nonlinear stability, and observational identifiability.

Source: https://www.emergentmind.com/topics/frozen-star