Quantum Oppenheimer-Snyder Black Holes

• Quantum Oppenheimer-Snyder black holes are models that replace the classical dust collapse singularity with a quantum bounce or smeared region using various quantization schemes.
• They feature modified exterior metrics, such as a deformed Schwarzschild solution with r⁻⁴ corrections, which impact horizon structure and gravitational-wave signatures.
• These approaches provide insights into singularity resolution, remnant formation, and information recovery, with model-dependent implications for black hole thermodynamics and evaporation.

Quantum Oppenheimer-Snyder black holes are quantum-gravity-motivated reworkings of the classical Oppenheimer-Snyder collapse model in which the dust-ball singularity is replaced by a bounce, a bound quantum state, or a smeared quantum region, and in many constructions the Schwarzschild exterior is deformed by a leading $r^{-4}$ correction. The literature does not contain a single unique model under this name. Instead, it includes loop-quantum-cosmology junction constructions with an effective exterior metric $f(r)=1-\frac{2M}{r}+\frac{\alpha M^2}{r^4}$, polymer quantizations of Lemaître-Tolman-Bondi dust with transient apparent horizons and an outgoing shock, and reduced or affine quantizations in which collapse becomes a Schrödinger or coherent-state problem with a nonzero minimum radius or a discrete spectrum (Lewandowski et al., 2022, Husain et al., 2022, Corda et al., 2021, Piechocki et al., 2020, Góźdź et al., 2023).

1. Classical starting point and principal quantization routes

The classical Oppenheimer-Snyder model describes a homogeneous pressureless dust ball whose interior is FLRW and whose exterior is Schwarzschild. Several quantum versions preserve this starting point but quantize different structures: the interior Friedmann dynamics, the marginally trapped dust geometry, or a reduced radial degree of freedom. As a result, “Quantum Oppenheimer-Snyder black hole” denotes a family of models rather than a single geometry.

Construction	Basic object	Characteristic result
LQC/APS junction	$H^2=\frac{8\pi G}{3}\rho(1-\rho/\rho_c)$, $f(r)=1-\frac{2M}{r}+\frac{\alpha M^2}{r^4}$	bounce, mass gap, deformed Schwarzschild exterior
Polymer LTB	$$\dot\beta+\frac{1}{2x^2}\partial_x[x^3\sin^2\beta]=0$$	transient apparent horizons, outgoing shock, lifetime $\propto M^2$
Rosen/Schrödinger	$V(r)=-GM\mu/r$, $E_n=-\mu(GM)^2/(2\hbar^2n^2)$	gravitational hydrogen atom
Affine or integral quantization	coherent states on the half-plane	bounce and smeared curvature invariants

In the LQC line, the interior scale factor satisfies a modified Friedmann equation with critical density $\rho_c$, and matching to a static exterior yields the quantum-deformed Schwarzschild metric (Lewandowski et al., 2022). In the polymer line, the marginally trapped LTB dynamics is quantized by the “improved-dynamics” prescription $b\to \sin[\mu(x)b]/\mu(x)$, which produces a non-singular bounce and a finite trapped region (Husain et al., 2022). In Rosen’s minisuperspace quantization, the collapse degree of freedom becomes a hydrogen-like radial Schrödinger problem with discrete bound states (Corda et al., 2021). In affine coherent-state and integral quantizations, both the comoving and stationary observer descriptions admit quantum-corrected bounces, though the stationary-observer bounce depends strongly on quantization ambiguities (Piechocki et al., 2020, Góźdź et al., 2023).

2. LQC-derived qOS geometry and horizon structure

A central qOS construction begins with the semiclassical APS interior

$f(r)=1-\frac{2M}{r}+\frac{\alpha M^2}{r^4}$0

with

$f(r)=1-\frac{2M}{r}+\frac{\alpha M^2}{r^4}$1

Matching across the dust boundary by the Darmois-Israel conditions forces the exterior lapse to become

$f(r)=1-\frac{2M}{r}+\frac{\alpha M^2}{r^4}$2

so the exterior line element is

$f(r)=1-\frac{2M}{r}+\frac{\alpha M^2}{r^4}$3

The classical Schwarzschild solution is recovered for $f(r)=1-\frac{2M}{r}+\frac{\alpha M^2}{r^4}$4 (Lewandowski et al., 2022).

The horizon equation

$f(r)=1-\frac{2M}{r}+\frac{\alpha M^2}{r^4}$5

implies a lower bound on black-hole formation. Real positive roots exist only if

$f(r)=1-\frac{2M}{r}+\frac{\alpha M^2}{r^4}$6

For $f(r)=1-\frac{2M}{r}+\frac{\alpha M^2}{r^4}$7 no horizon forms; for $f(r)=1-\frac{2M}{r}+\frac{\alpha M^2}{r^4}$8 the configuration is extremal; for $f(r)=1-\frac{2M}{r}+\frac{\alpha M^2}{r^4}$9 there are two distinct horizons $H^2=\frac{8\pi G}{3}\rho(1-\rho/\rho_c)$0 (Lewandowski et al., 2022).

This exterior admits an effective stress tensor $H^2=\frac{8\pi G}{3}\rho(1-\rho/\rho_c)$1 with local energy density

$H^2=\frac{8\pi G}{3}\rho(1-\rho/\rho_c)$2

a purely quantum contribution that vanishes as $H^2=\frac{8\pi G}{3}\rho(1-\rho/\rho_c)$3. The model explicitly states that this contribution may play the role of “dark matter” around astrophysical black holes (Lewandowski et al., 2022).

3. Bounce, causal structure, and the fate of horizons

In the static qOS extension, the classical singularity is replaced by a bounce at

$H^2=\frac{8\pi G}{3}\rho(1-\rho/\rho_c)$4

The maximal analytic extension parallels Reissner-Nordström: an infinite chain of static regions $H^2=\frac{8\pi G}{3}\rho(1-\rho/\rho_c)$5, trapped or anti-trapped regions $H^2=\frac{8\pi G}{3}\rho(1-\rho/\rho_c)$6, and inner wormhole regions $H^2=\frac{8\pi G}{3}\rho(1-\rho/\rho_c)$7, with alternating black-hole and white-hole sectors. In the collapse interpretation, the dust world-tube starts in a past asymptotically flat $H^2=\frac{8\pi G}{3}\rho(1-\rho/\rho_c)$8 region, crosses the outer horizon, bounces at $H^2=\frac{8\pi G}{3}\rho(1-\rho/\rho_c)$9, and emerges into a future $f(r)=1-\frac{2M}{r}+\frac{\alpha M^2}{r^4}$0 region via a white-hole horizon (Lewandowski et al., 2022).

The polymer LTB treatment produces a different causal narrative. For Oppenheimer-Snyder initial data, the dust boundary evolves as

$f(r)=1-\frac{2M}{r}+\frac{\alpha M^2}{r^4}$1

with minimum radius $f(r)=1-\frac{2M}{r}+\frac{\alpha M^2}{r^4}$2 at the bounce. The model tracks an outer horizon near $f(r)=1-\frac{2M}{r}+\frac{\alpha M^2}{r^4}$3 and an inner dynamical horizon that forms with it, then evolves until the trapped region disappears. After the bounce, a discontinuity in the effective field produces an outgoing shock wave, and the black-hole lifetime is found to be proportional to the square its mass (Husain et al., 2022).

The information-loss implications are correspondingly model dependent. In the polymer construction, the singularity is replaced by a space-time bounce and the event horizon by transient apparent horizons; the model therefore states that there is in principle a causal channel for information to escape and that Hawking emission would be cut off well before the Page time (Husain et al., 2022). In the static LQC junction model, by contrast, the maximal extension still contains outer and inner Killing horizons and a white-hole sector (Lewandowski et al., 2022). A plausible implication is that singularity resolution is common across these constructions, while the global causal structure is not.

The same LQC framework also admits the quantum Swiss Cheese model, obtained by reversing the junction: the dust ball is removed from the APS universe and replaced by the static qOS exterior. Its Penrose diagram contains an infinite sequence of black-hole and white-hole tunnels, and the “current” universe lies in an asymptotically flat region bordering a white-hole horizon, realizing a version of black-hole cosmology or fecund universes (Lewandowski et al., 2022).

4. Schrödinger, shell, and coherent-state formulations

Rosen’s quantization of Oppenheimer-Snyder collapse leads to a hydrogen-like radial Hamiltonian

$f(r)=1-\frac{2M}{r}+\frac{\alpha M^2}{r^4}$4

with spectrum

$f(r)=1-\frac{2M}{r}+\frac{\alpha M^2}{r^4}$5

and, in the two-particle picture $f(r)=1-\frac{2M}{r}+\frac{\alpha M^2}{r^4}$6,

$f(r)=1-\frac{2M}{r}+\frac{\alpha M^2}{r^4}$7

The associated Bohr radius

$f(r)=1-\frac{2M}{r}+\frac{\alpha M^2}{r^4}$8

places the ground state at the Planck scale. In this description, the classical $f(r)=1-\frac{2M}{r}+\frac{\alpha M^2}{r^4}$9 singularity is replaced by a non-singular two-particle system whose components are labeled “nucleus” and “electron,” with the “electron” interpreted as a horizon oscillation mode by de Broglie’s hypothesis (Corda et al., 2021).

The same program computes the qOS entropy as a function of the black-hole principal quantum number in terms of Bekenstein-Hawking entropy and three sub-leading corrections, and states that the coefficient of the Bekenstein-Hawking entropy is reduced to a quarter of the traditional value. It further argues that a rescaled spectrum matches the semi-classical Bohr-like approach for large principal quantum number and that the time evolution of the system solves the black-hole information paradox (Corda et al., 2021).

A related Schrödinger and Klein-Gordon treatment reduces the collapse to a self-gravitating thin shell or a two-body problem with effective Hamiltonian

$$\dot\beta+\frac{1}{2x^2}\partial_x[x^3\sin^2\beta]=0$$0

and interprets the black hole as a self-interacting, highly excited, spherically symmetric, massive quantum shell generated by matter condensing on the apparent horizon. That work states that black holes have neither horizons nor singularities and that there is neither information loss in black-hole evaporation, nor black-hole complementarity, nor firewall paradox (Corda, 2023). This is one of the strongest examples of model dependence in the qOS literature.

Affine coherent-state quantization yields a comoving lower symbol

$$\dot\beta+\frac{1}{2x^2}\partial_x[x^3\sin^2\beta]=0$$1

with non-singular solution

$$\dot\beta+\frac{1}{2x^2}\partial_x[x^3\sin^2\beta]=0$$2

For the stationary observer, a bounce appears only if the fiducial parameter $$\dot\beta+\frac{1}{2x^2}\partial_x[x^3\sin^2\beta]=0$$3 exceeds a critical $$\dot\beta+\frac{1}{2x^2}\partial_x[x^3\sin^2\beta]=0$$4, and then $$\dot\beta+\frac{1}{2x^2}\partial_x[x^3\sin^2\beta]=0$$5, so the bounce occurs outside the photon sphere (Piechocki et al., 2020). In integral quantization on the affine half-plane, expectation values of curvature monomials remain finite for packets with a hole excised around the classical singularity, and even when the $$\dot\beta+\frac{1}{2x^2}\partial_x[x^3\sin^2\beta]=0$$6 limit diverges, the variance remains strictly positive, so the singularity is never sharply localized but is irreducibly smeared by quantum fluctuations (Góźdź et al., 2023).

5. Perturbations, shadows, and gravitational-wave probes

For the static qOS metric, null geodesics and perturbations preserve much of the Schwarzschild machinery but with a modified lapse. The photon sphere obeys

$$\dot\beta+\frac{1}{2x^2}\partial_x[x^3\sin^2\beta]=0$$7

and to first order in the quantum parameter $$\dot\beta+\frac{1}{2x^2}\partial_x[x^3\sin^2\beta]=0$$8,

$$\dot\beta+\frac{1}{2x^2}\partial_x[x^3\sin^2\beta]=0$$9

The leading correction therefore shrinks both the photon-sphere radius and the shadow radius. In the same geometry, scalar and vector perturbations are governed by a Regge-Wheeler-Zerilli equation with effective potential

$\propto M^2$0

and the quasinormal frequencies satisfy $\propto M^2$1, indicating linear stability (Yang et al., 2022).

The decrease in damping rate persists in broader perturbative studies. For scalar, electromagnetic, and Dirac test fields, the qOS quasinormal spectrum remains close to Schwarzschild in $\propto M^2$2, while $\propto M^2$3 is significantly reduced; in the eikonal limit the null-geodesic–QNM correspondence continues to hold (Skvortsova, 2024). In the qOS–de Sitter extension,

$\propto M^2$4

increasing $\propto M^2$5 shrinks $\propto M^2$6, reduces the critical impact parameter $\propto M^2$7, and reduces the shadow size $\propto M^2$8 (Luo et al., 2024).

Extreme-mass-ratio inspirals provide a separate observational channel. In the static eccentric analysis, the dimensionless parameter is $\propto M^2$9; the quantum correction slows the evolution of the semi-latus rectum and eccentricity, generates a secular phase advance of several radians after $V(r)=-GM\mu/r$0 yr of inspiral even for $V(r)=-GM\mu/r$1, and yields mismatches larger than the distinguishability threshold in a LISA-style analysis (Yang et al., 29 Sep 2025). In the rotating qOS spacetime,

$V(r)=-GM\mu/r$2

the quantum correction also produces detectable EMRI dephasing and loss of waveform faithfulness, but larger spin $V(r)=-GM\mu/r$3 suppresses the effect. For the example $V(r)=-GM\mu/r$4, $V(r)=-GM\mu/r$5, $V(r)=-GM\mu/r$6, $V(r)=-GM\mu/r$7, the cumulative dephasing after one year is $V(r)=-GM\mu/r$8 rad for $V(r)=-GM\mu/r$9, depending on spin (Zhang et al., 25 Apr 2026).

6. Thermodynamics, evaporation, remnants, and primordial-black-hole applications

Several qOS papers treat Hawking emission semiclassically on the deformed exterior

$E_n=-\mu(GM)^2/(2\hbar^2n^2)$0

For minimally coupled massless scalar radiation, one study derives

$E_n=-\mu(GM)^2/(2\hbar^2n^2)$1

and concludes that the positive $E_n=-\mu(GM)^2/(2\hbar^2n^2)$2 term slows the late-time evaporation and eventually halts it at finite mass, leaving a stable remnant. A higher-order WKB quasinormal-mode analysis finds $E_n=-\mu(GM)^2/(2\hbar^2n^2)$3 for $E_n=-\mu(GM)^2/(2\hbar^2n^2)$4, supporting linear stability of the remnant. For nonminimal coupling, the fate depends on $E_n=-\mu(GM)^2/(2\hbar^2n^2)$5: the paper states that $E_n=-\mu(GM)^2/(2\hbar^2n^2)$6 leads to accelerated late-time evaporation with no remnant, while $E_n=-\mu(GM)^2/(2\hbar^2n^2)$7 again produces a stable remnant (Tan et al., 30 Dec 2025).

The tunneling approach gives a complementary entropy formula. In Painlevé-Gullstrand coordinates, massive-particle tunneling yields an emission rate $E_n=-\mu(GM)^2/(2\hbar^2n^2)$8 and an entropy

$E_n=-\mu(GM)^2/(2\hbar^2n^2)$9

so the qOS entropy contains the standard area term plus a logarithmic correction with positive coefficient (Tan, 9 Feb 2025).

Primordial-black-hole applications exploit the same thermodynamic deformation. In the qOS PBH study, the Hawking temperature

$\rho_c$0

is always lower than the Schwarzschild temperature, and the net primary photon spectrum is suppressed by up to $\rho_c$1–$\rho_c$2 orders of magnitude as $\rho_c$3. Folding this through the extragalactic $\rho_c$4-ray background from HEAO-1, COMPTEL, and EGRET relaxes the usual constraints and broadens the allowed mass window in which qOS PBHs can make all the dark matter. For the benchmark $\rho_c$5, the paper quotes

$\rho_c$6

to be compared with the narrower Schwarzschild window (Wang et al., 1 Apr 2026).

7. Extensions, variants, and model-dependent controversies

The qOS framework has been generalized in several directions. The charged qOS model embeds the quartic correction in nonlinear electrodynamics and matches the APS interior to a charged thin shell. The shell has equilibrium radius

$\rho_c$7

and small perturbations around $\rho_c$8 are oscillatory, so the shell is stable in the sense of the paper (Mazharimousavi, 12 Feb 2025). Higher-dimensional qOS collapse replaces the four-dimensional exterior by

$\rho_c$9

with bounce radius

$b\to \sin[\mu(x)b]/\mu(x)$0

and yields dimension-dependent quasinormal shifts together with an extra phase transition in the heat capacity (Shi et al., 2024).

Additional matter sectors further enrich the phenomenology. The Lorentz-term LQC version modifies the near-horizon geometry, slows quasinormal damping relative both to Schwarzschild and to the earlier Euclidean-term-only qOS model, and yields an entropy of the form

$b\to \sin[\mu(x)b]/\mu(x)$1

together with an extra small-radius phase transition in the heat capacity (Ou et al., 2 Aug 2025). Quintessence, string clouds, and perfect-fluid dark matter generate more elaborate lapse functions, modify photon spheres, shadow sizes, QNMs, greybody factors, and thermodynamic stability, and can either enhance or suppress observational signatures depending on the parameter set (Ahmed et al., 5 Aug 2025, Ahmed et al., 26 Feb 2026, Sajjad et al., 27 Apr 2026).

The main controversy is not whether the classical singularity is altered, but how the replacement should be interpreted. In the LQC junction model, $b\to \sin[\mu(x)b]/\mu(x)$2 gives a genuine two-horizon geometry with a Reissner-Nordström-like maximal extension (Lewandowski et al., 2022). In the polymer LTB treatment, the event horizon is replaced by transient apparent horizons and an outgoing shock (Husain et al., 2022). In the affine coherent-state model, the distant-observer bounce lies outside the photon sphere unless the minimal model is modified (Piechocki et al., 2020). In the shell-based Schrödinger-Klein-Gordon program, the paper explicitly states that black holes have neither horizons nor singularities (Corda, 2023). This suggests that singularity resolution is the most recurrent feature across the qOS literature, whereas the status of the horizon, remnant formation, and information recovery remains strongly dependent on the chosen quantization scheme and effective description.