Timelike Thin-Shell Evolution in Gravitational Collapse: Geometric and Thermodynamic Perspectives in Classical General Relativity

1. Introduction

Thermodynamic concepts—such as entropy, temperature, and energy flow—have deeply influenced how we understand gravity [1], especially in the context of black holes. Since the pioneering work of Bekenstein, Hawking, and others [2,3,4], geometric quantities like horizon area have been linked to entropy and information content. At the same time, it has long been appreciated that self-gravitating systems exhibit negative specific heat, a feature already visible in Newtonian theory [5,6]. These thermodynamic insights are often associated with equilibrium or stationary horizons—but their role in dynamical, collapsing systems remains less concretely understood.

This paper investigates whether classical general relativity (GR), without quantum input, already supports such thermodynamic behavior during the late stages of gravitational collapse. Specifically, we examine a simple but instructive model: a timelike thin shell separates a regular constant-curvature (de Sitter) interior from a Schwarzschild or Schwarzschild–de Sitter exterior. Thin-shell techniques provide an exact classical framework for this setup, based on Israel’s junction condition and the Birkhoff–Schwarzschild geometry [7,8,9].

After a short formation stage, we focus on a post-transient regime with negligible inflow, where the exterior mass is constant. The junction condition can then be reformulated as an effective potential for the shell radius $R\left(\tau \right)$, allowing a transparent analysis of its dynamics. Within this framework, we derive three core results:

(i)

a closed-form sufficient threshold for outward shell evolution,

(ii)

the boundedness of scalar curvature invariants across the entire spacetime domain, and

(iii)

a simple, falsifiable, mass-scaled frequency bound for near-shell spectral features observable at infinity.

Although purely classical, the setup naturally supports a thermodynamic interpretation: area growth, energy-driven expansion, and behavior consistent with negative heat capacity. This interpretation is discussed separately and used to connect the geometric findings to coarse-graining intuitions—without treating thermodynamics as a fundamental input.

1.1. Scope

The analysis is restricted to classical GR with spherical symmetry, timelike thin shells, and static patches with $f_{\pm }>0$. The lightlike (null) limit [10] is not considered. Quantum-information aspects such as evaporation, Page time, or firewalls lie outside the present scope; see [11,12,13] for background.

1.2. Organization

Section 2 introduces the geometric framework and conservation laws, leading to the effective-potential form of the junction condition. Section 3 derives the sufficient threshold for outward shell evolution. Section 4 proves curvature boundedness. Section 5 introduces the mass-scaled redshift bound for near-shell frequencies. Section 6 discusses physical interpretation, observational falsifiability, and future directions.

2. Geometric Setup and Methods

This section establishes the classical geometric framework used to model late-time gravitational collapse in spherical symmetry. The setup consists of two vacuum regions—an interior and an exterior—joined across a timelike thin shell. The entire analysis remains within general relativity (GR), using exact solutions and matching conditions.

2.1. Exterior (Vacuum, Spherical Symmetry)

By Birkhoff’s theorem, the vacuum exterior of a spherically symmetric source must be Schwarzschild or Schwarzschild–de Sitter (SdS):

$$d{s}_{+}^2=-f_{+}\left(r\right)d{t}_{+}^2+f_{+}^{-1}\left(r\right)d{r}^2+r^{2}d{\Omega }^2,f_{+}\left(r\right)=1-{\displaystyle \frac{2M}{r}}-{\displaystyle \frac{{\mathsf{\Lambda }}_{+}{r}^2}{3}},$$

(1)

where M is the exterior mass parameter (the ADM mass if ${\mathsf{\Lambda }}_{+}=0$). The analysis is confined to the static region where $f_{+}\left(r\right)>0$, consistent with timelike shell evolution.

2.2. Interior (Regular Constant Curvature)

The interior is modeled as a regular de Sitter region with constant positive curvature:

$$d{s}_{-}^2=-(1-H^2{r}^2)d{t}_{-}^2+{(1-H^2{r}^2)}^{-1}d{r}^2+r^{2}d{\Omega }^2,H^2={\mathsf{\Lambda }}_{-}/3>0,$$

(2)

so that all interior curvature invariants are finite, e.g., the Kretschmann scalar $K_{-}=8{\mathsf{\Lambda }}_{-}^2/3$ [4].

2.3. Timelike Thin Shell and Junction Condition

The two spacetime regions are joined across a timelike shell at areal radius $R\left(\tau \right)$, where $\tau$ is the proper time of observers comoving with the shell. The shell is endowed with surface stress–energy tensor ${S^a}_b=\mathrm{diag}(-\sigma ,p,p)$, and its motion obeys the Israel junction condition [7,9]:

$$\sqrt{f_{+}\left(R\right)+{\dot{R}}^2}-\sqrt{f_{-}\left(R\right)+{\dot{R}}^2}=\kappa \left(R\right),\kappa \left(R\right)=4\pi \sigma \left(R\right)R.$$

(3)

Additionally, the shell satisfies surface-energy conservation,

$${\displaystyle \frac{d\left(\sigma A\right)}{d\tau }}+p{\displaystyle \frac{dA}{d\tau }}=\mathsf{\Phi }A,A=4\pi R^2,$$

(4)

where $\mathsf{\Phi }$ is the net normal energy flux across the shell.

2.4. Post-Transient Zero-Inflow Regime

After an initial formation stage, the analysis focuses on a late-time regime in which energy inflow is negligible. This is enforced by the condition

$${T^r}_t{|}_{\mathrm{shell}}=0⟹dM\simeq 0,R_S=2M=\mathrm{const},$$

(5)

ensuring that the exterior mass parameter remains fixed and attention shifts to the shell dynamics governed by $R\left(\tau \right)$ and the interior scale ${\mathsf{\Lambda }}_{-}$.

2.5. Domain and Conventions

All analysis is restricted to static patches where $f_{\pm }\left(R\right)>0$, ensuring that the shell remains timelike. Natural units are used throughout: $G=c=\hslash =k_B=1$, and the metric signature is $(-+++)$.

The Israel condition can be recast in an effective-potential form by solving for $\dot{R}$ in (3):

$${\dot{R}}^2+V\left(R\right)=0,V\left(R\right)=f_{-}\left(R\right)-{\displaystyle \frac{{\left(f_{+}\left(R\right)-f_{-}\left(R\right)-\kappa {\left(R\right)}^2\right)}^2}{4\kappa {\left(R\right)}^2}}.$$

(6)

Turning points occur where $V\left(R_{*}\right)=0$, and the sign of $V^{\prime }\left(R_{*}\right)$ determines the local evolution tendency.

To control sign behavior, a linear surface equation of state $p=w\sigma$ with $w>-1/2$ and vanishing flux $\mathsf{\Phi }=0$ is assumed. Then energy conservation (4) reduces to

$${\displaystyle \frac{d\sigma }{dR}}=-{\displaystyle \frac{2(\sigma +p)}{R}}⟹{\kappa }^{\prime }\left(R\right)={\displaystyle \frac{d}{dR}}\left(4\pi \sigma R\right)=4\pi \sigma (-1-2w)<0,$$

(7)

ensuring monotonic behavior of the effective potential in physically relevant cases.

3. A Sufficient Outward-Evolution Threshold

Balancing interior de Sitter acceleration ($\sim H^{2}R$) against exterior attraction ($\sim M/R^2$) singles out

$$R_{\mathrm{thr}}:={\left({\displaystyle \frac{3M}{{\mathsf{\Lambda }}_{-}}}\right)}^{1/3}.$$

(8)

Proposition 1

(Sufficient outward evolution at zero inflow). Assume $\mathsf{\Phi }=0$ (so that M is time-independent), finite surface stresses with $p=w\sigma$ and $w>-1/2$, and $f_{\pm }\left(R\right)>0$ (static patches). If for some ${\tau }_0$ the shell satisfies $R\left({\tau }_0\right)\ge R_{\mathrm{thr}}$, then $\dot{R}\left({\tau }_0\right)\ge 0$ and the shell evolves outward on a finite interval beyond ${\tau }_0$.

Proof.(sketch). Differentiate (6) and evaluate at a point with $H^2{R}^3=M$ (for ${\mathsf{\Lambda }}_{+}=0$). One has

$$\begin{array}{cccc}\hfill f_{+}^{\prime }\left(R\right)& ={\displaystyle \frac{2M}{R^2}},\hfill & \hfill f_{-}^{\prime }\left(R\right)& =-2H^{2}R,\hfill \\ \hfill \Delta f& :=f_{+}-f_{-}=-{\displaystyle \frac{2M}{R}}+H^2{R}^2,\hfill & \hfill \Delta f^{\prime }& ={\displaystyle \frac{2M}{R^2}}+2H^{2}R=4H^{2}R>0,\hfill \end{array}$$

using $H^2{R}^3=M$. Writing $V^{\prime }=f_{-}^{\prime }-{\displaystyle \frac{\partial }{\partial R}}\left[{\displaystyle \frac{{(\Delta f-{\kappa }^2)}^2}{4{\kappa }^2}}\right]$ and using ${\kappa }^{\prime }\left(R\right)<0$ from (7), one finds $V^{\prime }\left(R_{\mathrm{thr}}\right)<0$ within the static domain $f_{\pm }>0$. Hence $-V$ increases for $R>R_{\mathrm{thr}}$ in a neighborhood of $R_{\mathrm{thr}}$, implying outward evolution on a finite interval when $\dot{R}\left({\tau }_0\right)\ge 0$. □

Remark 1

(Dimensionless form and static-patch consistency). Let $\chi ={\mathsf{\Lambda }}_{-}{M}^2$. Then

$${\displaystyle \frac{R_{\mathrm{thr}}}{R_S}}={\displaystyle \frac{1}{2}}{\left({\displaystyle \frac{3}{\chi }}\right)}^{1/3}.$$

Timelike consistency requires $R_{\mathrm{thr}}>R_S$ and $H{R}_{\mathrm{thr}}<1$ (and, if ${\mathsf{\Lambda }}_{+}>0$, $R_{\mathrm{thr}}$ below the cosmological horizon), i.e. $f_{\pm }\left(R\right)>0$ throughout [8].

Beyond providing control over the shell’s motion, the static patch framework also ensures that curvature remains bounded. This regularity is addressed next.

4. Boundedness of Curvature Scalars in the Covered Domain

A central question in collapse scenarios is whether curvature invariants diverge somewhere in the constructed spacetime, signaling a singularity. In this section, we show that the entire region covered by the thin-shell construction remains regular. This follows directly from the known properties of the chosen interior and exterior geometries, together with the well-controlled shell junction.

Proposition 2

(Bounded curvature). Let $R\left(\tau \right)\ge R_{min}>0$, with $f_{\pm }\left(R\right)>0$ (i.e., the shell remains within the static patches), and assume that the shell surface stresses $\sigma \left(\tau \right),p\left(\tau \right)$ remain finite. Then all scalar curvature invariants are finite throughout the spacetime region covered by the construction.

Proof.(sketch). The interior region is a regular constant-curvature spacetime, specifically de Sitter, in which all scalar invariants are finite by construction. For example, the Kretschmann scalar is

$$K_{-}=R_{\mu \nu \rho \sigma }{R}^{\mu \nu \rho \sigma }={\displaystyle \frac{8}{3}}{\mathsf{\Lambda }}_{-}^2.$$

In the exterior, which is Schwarzschild–de Sitter (SdS), the curvature is also fully determined and remains finite as long as the radius stays above the Schwarzschild horizon. The Kretschmann scalar is

$$K_{+}\left(r\right)={\displaystyle \frac{48M^2}{r^6}}+{\displaystyle \frac{8}{3}}{\mathsf{\Lambda }}_{+}^2,$$

which is manifestly bounded for $r\ge R_{min}>0$. □

At the shell, the metric is continuous, and any curvature singularities would arise only from distributional sources due to the jump in extrinsic curvature. However, as long as the shell energy density $\sigma$ and pressure p remain finite, the distributional curvature remains controlled. This is a standard result from the Israel junction formalism [7,9].

Moreover, because we restrict to static patches with $f_{\pm }\left(R\right)>0$, the shell does not cross into regions containing horizons or trapped surfaces. This ensures that the construction remains within a regular coordinate patch, and the entire covered spacetime is free of curvature singularities [8].

Beyond geometric regularity, classical timelike shell models can also constrain persistent late-time features. In particular, a mass-scaled bound emerges from combining shell locality with redshift effects.

5. A Mass-Scaled Frequency Bound for Near-Shell Modes

Late-time observations often probe the presence of quasi-stationary or slowly decaying features near compact objects. In the present context, such features may be modeled as near-shell modes, and their observability at infinity is governed by redshift and confinement effects.

A minimal localization criterion for such modes is $k_{\mathrm{loc}}R\gtrsim \xi$, with $\xi =\mathcal{O}\left(1\right)$, ensuring that the mode is stored near the shell with meaningful structure. The Tolman redshift relation then yields the observed frequency at infinity:

$${\omega }_{\infty }=\sqrt{f_{+}\left(R\right)}{\omega }_{\mathrm{loc}},$$

where ${\omega }_{\mathrm{loc}}$ is the local proper frequency near the shell [4]. A full derivation and interpretation of this relation in the thin-shell context is given in Appendix A.

Assuming ${\omega }_{\mathrm{loc}}\simeq k_{\mathrm{loc}}$ and defining the redshifted cyclic frequency $f_c={\omega }_{\infty }/\left(2\pi \right)$, we obtain the general mass-scaled expression:

$$f_c{R}_S={\displaystyle \frac{\xi }{2\pi }}{\displaystyle \frac{\sqrt{1-{\displaystyle \frac{2M}{R}}-{\displaystyle \frac{{\mathsf{\Lambda }}_{+}{R}^2}{3}}}}{R/R_S}}.$$

(9)

For a Schwarzschild exterior (${\mathsf{\Lambda }}_{+}=0$), the right-hand side is maximized at $R=\frac{3}{2}{R}_S$, yielding the compact upper bound:

$$\begin{array}{|c|}\hline f_c{R}_S\le {\displaystyle \frac{\xi }{3\sqrt{3}\pi }}\\ \hline\end{array}.$$

(10)

This inequality provides a falsifiable, dimensionless constraint on persistent near-shell spectral features in the absence of exterior inflow and under classical GR assumptions. A robust observational violation of Equation (10) would falsify the near-shell storage scenario as described.

6. Discussion and Outlook

6.1. Classical Content

Within classical GR and spherical symmetry, with a timelike thin shell evolving in static patches $f_{\pm }>0$, a constant-curvature (de Sitter) interior, and a post-transient zero-inflow regime ($dM\simeq 0$), the analysis establishes three conservative statements. First, a closed-form sufficient threshold for outward evolution, Equation (8), emerges from balancing interior de Sitter acceleration against exterior attraction. Second, scalar curvature invariants remain bounded throughout the spacetime region covered by the construction (Proposition 2), ensuring regularity without intersecting exterior trapping surfaces. Third, a mass-scaled redshift bound limits persistent near-shell frequencies as observed at infinity, Equation (10). These results derive directly from the Israel junction geometry and standard shell conservation laws, requiring no inputs beyond classical GR.

6.2. Energy Conditions

The de Sitter interior saturates the null and weak energy conditions but violates the strong energy condition, consistent with vacuum energy properties. The thin shell, governed by a linear surface equation of state $p=w\sigma$ with $w>-1/2$, satisfies the null, weak, and dominant energy conditions for physically admissible configurations ($\sigma >0$, $\left|w\right|\le 1$). Under these assumptions, the standard singularity theorems [14] are not activated within the domain, as the model avoids the focusing conditions that typically drive geodesic incompleteness [8]. This highlights how controlled violations of energy conditions can yield regular classical outcomes in collapse scenarios [15].

6.3. Thermodynamic Interpretation

Thermodynamic language serves as an interpretive framework, distinct from the geometric derivations. In the static-patch, zero-inflow setup with fixed exterior mass M, the shell radius $R\left(\tau \right)$ naturally links to a geometric area functional $S\left[R\right]:=\alpha R^2$ (where $\alpha$ is an arbitrary scale). Outward evolution implies $dS>0$ without external inflow, mirroring the phenomenology of self-gravitating systems with negative specific heat capacity ($C<0$). In this scenario, larger shell configurations correspond to higher “entropy” states, consistent with the idea that energy loss promotes expansion rather than collapse. Crucially, $C<0$ allows for local entropy production via $dS\sim dR$ during the evolution, potentially reaching an entropic maximum before geometric horizon formation. No statistical entropy, temperature, or microstates are invoked; the analogy complements the horizon-free junction geometry, offering a classical lens on late-time stability [16].

6.4. Scope and Limitations

The results are confined to classical GR, spherical symmetry, timelike shells in static patches, finite surface stresses (with $p=w\sigma$, $w>-1/2$), and exclusion of the null limit [10]. The threshold in Equation (8) is sufficient, not necessary. These restrictions ensure analytic control but limit applicability to highly idealized situations. Generalizations involving strong time dependence, rotation, or anisotropies are left for future work.

6.5. Observational Handle and Falsifiability

The mass–scaled inequality in Equation (10) furnishes a null test: given an independent estimate of M, any long–lived near–shell spectral feature measured at infinity should respect the bound. A robust, systematic violation would falsify the near–shell storage picture under the stated assumptions (timelike shell in static patches, negligible inflow, time–independent exterior mass). The bound also provides a practical template for targeted searches formulated directly in the mass–scaled variable $f_c{R}_S$.

6.6. Outlook

Natural extensions could test the robustness of the classical findings while linking to observational strategies. Key directions include: (i) incorporating the null limit [10] to model early, ultrarelativistic collapse phases, and exploring how lightlike shell dynamics might affect the outward-evolution threshold (Equation (8)) and entropy production ($dS\sim dR$); (ii) studying perturbative stability against nonspherical fluctuations, particularly regarding their impact on the boundedness of curvature invariants; and (iii) allowing for time-dependent exterior mass $M\left(t\right)$ due to inflow or radiation, to explore how mass evolution modifies the frequency bound (Equation (10)) and potentially introduces observational signatures. These avenues could connect the junction framework to broader spacetime symmetries, astrophysical realism, and alternative gravity theories.

Appendix D.1. Schwarzschild Exterior (Lambda+=0)

Let $x:=R/M>2$ so that $R/R_S=x/2$. Define

$$g\left(x\right):={\displaystyle \frac{\sqrt{1-{\displaystyle \frac{2}{x}}}}{x/2}}={\displaystyle \frac{2}{x}}\sqrt{1-{\displaystyle \frac{2}{x}}},x>2.$$

(A4)

Then (A3) reads $f_c{R}_S\gtrsim {\displaystyle \frac{\xi }{2\pi }}g\left(x\right)$. A straightforward maximization gives $g^{\prime }\left(x\right)=0\iff x=3$, and

$$g\left(3\right)={\displaystyle \frac{2}{3}}\sqrt{1-{\displaystyle \frac{2}{3}}}={\displaystyle \frac{2}{3\sqrt{3}}}.$$

Therefore,

$$f_c{R}_S\lesssim {\displaystyle \frac{\xi }{2\pi }}g\left(3\right)={\displaystyle \frac{\xi }{3\sqrt{3}\pi }},$$

(A5)

which is the bound (10) in the main text.

Appendix D.2. Including a Cosmological Term

For ${\mathsf{\Lambda }}_{+}>0$, the same recipe applies on the static patch $f_{+}\left(R\right)>0$. Since the extra $-{\mathsf{\Lambda }}_{+}{R}^2/3$ lowers $\sqrt{f_{+}\left(R\right)}$ at large R, the Schwarzschild value is an upper envelope; any ${\mathsf{\Lambda }}_{+}>0$ tightens the bound.

Appendix D.3. Interpretation

Given an independent mass estimate M, any persistent late–time feature localized near the shell and measured at infinity must obey the mass–scaled inequality above. A robust violation would falsify the near–shell storage picture within the assumptions of the main text.