Higher-D Quantum Oppenheimer–Snyder Model

• The paper introduces a higher-dimensional quantum Oppenheimer–Snyder framework where LQC-inspired effective Friedmann dynamics replace the classical singularity with a finite-radius bounce.
• It details how junction conditions yield a quantum-corrected external metric exhibiting modified horizon structure, quasinormal mode spectra, and altered thermodynamic properties.
• The study extends the framework to incorporate a cosmological constant, providing insights into AdS holography and phase-space thermodynamics in quantum gravitational collapse.

The higher-dimensional quantum Oppenheimer–Snyder model is a $D=d+1$ dimensional generalization of Oppenheimer–Snyder collapse in which a homogeneous, isotropic, pressureless dust interior is matched to a static, spherically symmetric exterior, and the interior dynamics is modified by effective higher-dimensional Loop Quantum Cosmology (LQC). In this framework, the junction conditions determine a quantum-corrected exterior metric, the collapse reaches a finite-radius bounce rather than a classical singularity, and the resulting black-hole geometry acquires modified horizon structure, quasinormal spectra, and thermodynamics. A later extension incorporates a cosmological constant and studies the AdS case in extended phase-space thermodynamics. Related AdS/CFT work uses an Oppenheimer–Snyder collapse in asymptotically AdS space as a holographic model of thermalization, but there the term “quantum” refers to the dual strongly coupled quantum field theory rather than to LQC-inspired quantum geometry corrections (Shi et al., 2024, Jiang et al., 6 Jun 2026, Taanila, 2015).

1. Classical higher-dimensional Oppenheimer–Snyder framework

The underlying classical model describes the gravitational collapse of a spherically symmetric, homogeneous, pressureless dust star. In $D=d+1$ spacetime dimensions with $d\ge 3$, the interior is taken to be a spatially flat FRW region and the exterior classical vacuum solution is the Schwarzschild–Tangherlini metric. The interior line element is

$$ds^2_{\rm in} = -dT^2 + a^2(T)\left[dR^2 + R^2 d\Omega^2_{d-1}\right],$$

with stellar surface at comoving radius $R=R_0$ and physical radius

$\tilde{R}(T) = a(T)R_0.$

For pressureless dust, $T^\mu_{\ \nu}=\rho\,u^\mu u_\nu$ with $p=0$, and conservation gives

$\rho(T)=\rho_0\left(\frac{a_0}{a}\right)^d.$

The physical volume of the $d$-ball of radius $D=d+1$0 is

$D=d+1$1

and the paper introduces

$D=d+1$2

so that $D=d+1$3. The classical Friedmann equation is

$D=d+1$4

which at the stellar surface becomes

$D=d+1$5

The corresponding classical exterior is

$D=d+1$6

with horizon radius determined by

$D=d+1$7

This higher-dimensional OS construction preserves the defining ingredients of the four-dimensional model—homogeneity, isotropy, and dust—while replacing the exterior Schwarzschild geometry by its Tangherlini generalization (Shi et al., 2024).

2. LQC-inspired effective dynamics and singularity avoidance

The quantum model replaces the classical interior Friedmann equation by the higher-dimensional effective LQC equation

$D=d+1$8

with quantum parameter

$D=d+1$9

and critical density

$d\ge 3$0

Using $d\ge 3$1 and $d\ge 3$2, this becomes

$d\ge 3$3

The effective dynamics implies a bounce when $d\ge 3$4, equivalently when $d\ge 3$5. The bounce radius is

$d\ge 3$6

and in terms of the scale factor the minimum is fixed by

$d\ge 3$7

Since $d\ge 3$8 never reaches zero, the interior energy density is bounded by $d\ge 3$9 and the classical collapse singularity is avoided.

The existence of the bounce does not depend on $$ds^2_{\rm in} = -dT^2 + a^2(T)\left[dR^2 + R^2 d\Omega^2_{d-1}\right],$$0 beyond the scaling of $$ds^2_{\rm in} = -dT^2 + a^2(T)\left[dR^2 + R^2 d\Omega^2_{d-1}\right],$$1 and $$ds^2_{\rm in} = -dT^2 + a^2(T)\left[dR^2 + R^2 d\Omega^2_{d-1}\right],$$2. In this sense, the higher-dimensional quantum OS model is not merely a higher-dimensional restatement of classical collapse; the nonperturbative quantum geometry input enters directly through the modified Friedmann equation and forces a turning point at finite radius (Shi et al., 2024).

3. Junction conditions and the quantum-corrected exterior geometry

The interior and exterior are matched by Israel–Darmois junction conditions. Continuity of the induced metric on the stellar surface $$ds^2_{\rm in} = -dT^2 + a^2(T)\left[dR^2 + R^2 d\Omega^2_{d-1}\right],$$3 gives

$$ds^2_{\rm in} = -dT^2 + a^2(T)\left[dR^2 + R^2 d\Omega^2_{d-1}\right],$$4

while continuity of the extrinsic curvature yields $$ds^2_{\rm in} = -dT^2 + a^2(T)\left[dR^2 + R^2 d\Omega^2_{d-1}\right],$$5 and, for the angular components,

$$ds^2_{\rm in} = -dT^2 + a^2(T)\left[dR^2 + R^2 d\Omega^2_{d-1}\right],$$6

with $$ds^2_{\rm in} = -dT^2 + a^2(T)\left[dR^2 + R^2 d\Omega^2_{d-1}\right],$$7 a conserved energy along $$ds^2_{\rm in} = -dT^2 + a^2(T)\left[dR^2 + R^2 d\Omega^2_{d-1}\right],$$8. Choosing the time coordinate so that $$ds^2_{\rm in} = -dT^2 + a^2(T)\left[dR^2 + R^2 d\Omega^2_{d-1}\right],$$9 implies $R=R_0$0, and combining the matching with $R=R_0$1 gives

$R=R_0$2

Substituting the classical interior Friedmann equation reproduces the Schwarzschild–Tangherlini metric, while substituting the LQC-corrected equation yields

$R=R_0$3

and hence

$R=R_0$4

Horizons are roots of $R=R_0$5. The function $R=R_0$6 has a unique minimum at

$R=R_0$7

and the critical mass $R=R_0$8 is defined by $R=R_0$9. For $\tilde{R}(T) = a(T)R_0.$0 there are two horizons $\tilde{R}(T) = a(T)R_0.$1; for $\tilde{R}(T) = a(T)R_0.$2 there is a double horizon; for $\tilde{R}(T) = a(T)R_0.$3 no horizon forms. The bounce radius and horizon structure together imply that, when $\tilde{R}(T) = a(T)R_0.$4, the collapsing star crosses $\tilde{R}(T) = a(T)R_0.$5, reaches $\tilde{R}(T) = a(T)R_0.$6, bounces, re-expands, and exits through a white-hole asymptotic region in the extended spacetime (Shi et al., 2024).

With a cosmological constant, the same matching logic is retained. The classical exterior becomes Schwarzschild-(A)dS,

$\tilde{R}(T) = a(T)R_0.$7

with

$\tilde{R}(T) = a(T)R_0.$8

The interior quantum dynamics is modified to

$\tilde{R}(T) = a(T)R_0.$9

and the effective exterior is written as

$T^\mu_{\ \nu}=\rho\,u^\mu u_\nu$0

The junction conditions are identical to the classical case; the only change is $T^\mu_{\ \nu}=\rho\,u^\mu u_\nu$1. In the $T^\mu_{\ \nu}=\rho\,u^\mu u_\nu$2 limit, the Schwarzschild-(A)dS result is recovered (Jiang et al., 6 Jun 2026).

4. Scalar perturbations and quasinormal spectra

The quantum-corrected exterior supports a standard perturbation analysis for a massless scalar field. Using the separation

$T^\mu_{\ \nu}=\rho\,u^\mu u_\nu$3

the radial equation in the metric $T^\mu_{\ \nu}=\rho\,u^\mu u_\nu$4 is

$T^\mu_{\ \nu}=\rho\,u^\mu u_\nu$5

Introducing the tortoise coordinate $T^\mu_{\ \nu}=\rho\,u^\mu u_\nu$6 through $T^\mu_{\ \nu}=\rho\,u^\mu u_\nu$7 gives the Schrödinger-like form

$T^\mu_{\ \nu}=\rho\,u^\mu u_\nu$8

with effective potential

$T^\mu_{\ \nu}=\rho\,u^\mu u_\nu$9

Quasinormal-mode boundary conditions are purely ingoing at the horizon, $p=0$0 as $p=0$1, and purely outgoing at infinity, $p=0$2 as $p=0$3.

The spectra were computed with a 13th-order WKB method with Padé improvement and checked against time-domain finite element evolutions. The reported trends are systematic: increasing $p=0$4 at fixed $p=0$5 reduces both $p=0$6 and $p=0$7; increasing $p=0$8 at fixed $p=0$9 increases both $\rho(T)=\rho_0\left(\frac{a_0}{a}\right)^d.$0 and $\rho(T)=\rho_0\left(\frac{a_0}{a}\right)^d.$1; relative to Schwarzschild–Tangherlini, quantum corrections always reduce $\rho(T)=\rho_0\left(\frac{a_0}{a}\right)^d.$2, while $\rho(T)=\rho_0\left(\frac{a_0}{a}\right)^d.$3 increases in lower dimensions and decreases in higher dimensions. For the illustrative fundamental mode with $\rho(T)=\rho_0\left(\frac{a_0}{a}\right)^d.$4, $\rho(T)=\rho_0\left(\frac{a_0}{a}\right)^d.$5, and $\rho(T)=\rho_0\left(\frac{a_0}{a}\right)^d.$6 in the paper’s units, the values are

$\rho(T)=\rho_0\left(\frac{a_0}{a}\right)^d.$7

and

$\rho(T)=\rho_0\left(\frac{a_0}{a}\right)^d.$8

These results identify the perturbative imprint of the quantum-corrected OS background on ringdown observables (Shi et al., 2024).

5. Thermodynamics, remnants, and extended phase space

For the asymptotically flat quantum-corrected black holes, the Hawking temperature is determined by

$\rho(T)=\rho_0\left(\frac{a_0}{a}\right)^d.$9

In $d$0 ($d$1), the paper gives

$d$2

and

$d$3

Assuming the first law $d$4, the entropy is obtained by integration. The leading large-$d$5 results are

$d$6

$d$7

$d$8

The Hawking temperature rises to a peak and then decreases to zero at the extremal radius, and the heat capacity exhibits an extra divergence that signals an additional phase transition introduced by quantum corrections. The free energy $d$9 is larger than the classical value for small $D=d+1$00 and smaller otherwise (Shi et al., 2024).

In the AdS extension, the cosmological constant is treated as pressure,

$D=d+1$01

with thermodynamic volume

$D=d+1$02

Entropy is defined by

$D=d+1$03

where $D=d+1$04 is set by $D=d+1$05 so that $D=d+1$06. At small $D=d+1$07, the paper gives, for example,

$D=d+1$08

$D=d+1$09

$D=d+1$10

A notable result is that the entropy is independent of $D=d+1$11 for any $D=d+1$12; quantum corrections enter only through $D=d+1$13.

The small-AdS-black-hole regime is thermally regularized. Classically,

$D=d+1$14

so the temperature diverges. In the quantum-corrected model, the square-root structure $D=d+1$15 with $D=d+1$16 implies instead

$D=d+1$17

In five dimensions, the heat capacity

$D=d+1$18

has a new divergence at small $D=d+1$19, where $D=d+1$20 changes sign and a thermally stable near-remnant branch appears. The Gibbs free energy $D=d+1$21 shows a swallow-tail for $D=d+1$22, and the critical exponents are the mean-field values $D=d+1$23, $D=d+1$24, $D=d+1$25, and $D=d+1$26, independent of $D=d+1$27 (Jiang et al., 6 Jun 2026).

6. AdS/holographic line, interpretive distinctions, and limitations

A separate but related line of work studies Oppenheimer–Snyder collapse in asymptotically AdS space as a holographic model of thermalization. There the bulk spacetime is $D=d+1$28-dimensional, the boundary CFT lives in $D=d+1$29 dimensions, the exterior is AdS-Schwarzschild,

$D=d+1$30

and the interior is a hyperbolic FRW dust region,

$D=d+1$31

The junction conditions imply

$D=d+1$32

with $D=d+1$33. Thermalization is probed by the geodesic approximation for heavy operators,

$D=d+1$34

The paper finds that the OS collapse trajectory is always faster than the pressureless thin-shell model with the same asymptotic mass and initial radius, and that the renormalized geodesic length approaches equilibrium faster in OS collapse. For large angular separations and intermediate times, the geodesic length can be multivalued, and choosing the continuous branch yields kinks with discontinuous derivative (Taanila, 2015).

These papers therefore use two distinct notions of “quantum.” In the holographic AdS setting, the bulk collapse is classical and the quantum content lies in the dual nonequilibrium boundary field theory and in the geodesic approximation to boundary correlators. In the higher-dimensional quantum OS model and its cosmological-constant extension, the quantum input is instead the LQC-inspired effective Friedmann dynamics of the interior, encoded through $D=d+1$35, $D=d+1$36, $D=d+1$37, and $D=d+1$38, and then transferred to the exterior through the junction conditions.

The scope and limitations are correspondingly specific. The LQC-based constructions assume a homogeneous, isotropic, pressureless dust interior, exact spherical symmetry, a static and spherically symmetric exterior, and quantum corrections encoded solely through the interior effective Friedmann equation; anisotropies, inhomogeneities, and quantum backreaction in the exterior are neglected. The perturbative analysis is restricted to scalar perturbations, and the global fate of the inner Cauchy horizon remains open. In the AdS extension, thermodynamic expressions are computed for small $D=d+1$39 with $D=d+1$40 treated as fixed, and promoting $D=d+1$41 to a thermodynamic variable would require a generalized first law and Smarr relation. Open questions identified in the literature include gravitational perturbations beyond the scalar sector, detailed Hawking radiation spectra and greybody factors, ensemble analyses with fluctuating $D=d+1$42, and the dynamical endpoint of the bounce in AdS, including white-hole transition versus remnant scenarios (Shi et al., 2024, Jiang et al., 6 Jun 2026).