Spherical Expanding Shell Model

• The spherical expanding shell model is a framework describing dynamic, evolving spherical layers whose radius, density, and physical state change over time in various physical contexts.
• It employs diverse modeling techniques ranging from the Gross–Pitaevskii equation in Bose–Einstein condensates to relativistic mechanics and radiation-diffusion formulations.
• This model enables practical insights into phenomena such as supernova remnant expansion, gravitational shell fragmentation, and pressure-assisted instabilities validated by simulations.

Searching arXiv for the cited shell-model papers to ground the article in current arXiv records. A spherical expanding shell model is a family of formulations in which the dominant degrees of freedom are attached to a spherical shell or shell-like layer whose radius, thickness, density, phase, or radiative state evolves in time. In the literature considered here, this includes a thin shell-shaped Bose–Einstein condensate released from a radially shifted harmonic trap, a relativistic dividing shell between collapsing and expanding regions, multilayer thin shells separated by vacuum, expansion-free cavities bounded by anisotropic matter, pressure-confined gaseous shells that fragment gravitationally, and directly imaged or reconstructed shells in supernovae and the nearby interstellar medium (Tononi, 22 Apr 2026, Mimoso et al., 2013, Gaspar et al., 2011, Prisco et al., 2011, Wunsch et al., 2010, Witt et al., 2015, Bialy et al., 2021).

1. Shared geometric structure

Across these realizations, the shell is either a thin layer concentrated near a radius $R$, or a finite-thickness region bounded by inner and outer radii. In the shell-shaped Bose–Einstein condensate, the confining potential is

$U(r)=\frac{1}{2}m\omega_r^2(r-R)^2,$

and for $R\gg l_r$, with radial oscillator length

$l_r=\sqrt{\frac{\hbar}{m\omega_r}},$

the condensate is localized in a narrow region around $r\simeq R$. In radiation-diffusion models, the shell occupies $r\ge r_0$ and a moving radiative front is written as $r=f(t)$. In relativistic shell mechanics, the shell is a timelike hypersurface with physical radius $r(\tau)$, while in static Newtonian and GR shell families it is a compact or infinite-thickness mass distribution with a central hole or a pair of separated radial maxima (Tononi, 22 Apr 2026, Smith et al., 11 Feb 2025, Gaspar et al., 2011, Vogt et al., 2010).

Realization	Shell definition	Characteristic relation
Shell-shaped BEC	Thin spherical shell around $r\simeq R$	$$i\hbar \partial_t\Psi = \left[-\frac{\hbar^2}{2m}\nabla^2 + g|\Psi|^2\right]\Psi$$
Relativistic dividing shell	Surface separating $U(r)=\frac{1}{2}m\omega_r^2(r-R)^2,$0 and $U(r)=\frac{1}{2}m\omega_r^2(r-R)^2,$1	$U(r)=\frac{1}{2}m\omega_r^2(r-R)^2,$2 and $U(r)=\frac{1}{2}m\omega_r^2(r-R)^2,$3
Radiation-diffusion shell	Cold shell $U(r)=\frac{1}{2}m\omega_r^2(r-R)^2,$4 with front $U(r)=\frac{1}{2}m\omega_r^2(r-R)^2,$5	$U(r)=\frac{1}{2}m\omega_r^2(r-R)^2,$6

A common structural motif is the separation between shell geometry and internal physics. The geometry supplies the spherical Laplacian, shell area growth, or junction conditions; the internal physics supplies the evolution law, such as the Gross–Pitaevskii equation, a generalized Tolman–Oppenheimer–Volkoff balance, diffusion, ionization balance, or thin-layer momentum conservation. Taken together, these works suggest that the shell formulation is most natural when the dominant dynamics is organized by a distinguished spherical layer rather than by a fully volumetric flow.

2. Quantum and superfluid realization

In the shell-shaped Bose–Einstein condensate realization, the initial wavefunction factorizes as

$U(r)=\frac{1}{2}m\omega_r^2(r-R)^2,$7

with

$U(r)=\frac{1}{2}m\omega_r^2(r-R)^2,$8

and uniform surface density

$U(r)=\frac{1}{2}m\omega_r^2(r-R)^2,$9

The thin-shell approximation requires $R\gg l_r$0 and

$R\gg l_r$1

so that the radial degree of freedom is frozen and the relevant dynamics initially resides on the spherical surface. After trap release, the system evolves under the 3D time-dependent Gross–Pitaevskii equation

$R\gg l_r$2

The vortex-free case is obtained when the vortex–antivortex separation $R\gg l_r$3, so that $R\gg l_r$4. The resulting expansion is spherically symmetric and develops a central density peak at the origin together with concentric interference ripples. The central peak comes from inward-moving portions of the shell focusing and interfering at $R\gg l_r$5, while the ripples are radial interference fringes formed during subsequent dispersive evolution. In the simulations highlighted in the paper, the parameters are $R\gg l_r$6, $R\gg l_r$7, $R\gg l_r$8, and $R\gg l_r$9 (Tononi, 22 Apr 2026).

The inclusion of vorticity is implemented through a vortex–antivortex dipole on the sphere. Its phase field is

$l_r=\sqrt{\frac{\hbar}{m\omega_r}},$0

with vortices placed at

$l_r=\sqrt{\frac{\hbar}{m\omega_r}},$1

so that the angular separation is $l_r=\sqrt{\frac{\hbar}{m\omega_r}},$2. The velocity field is

$l_r=\sqrt{\frac{\hbar}{m\omega_r}},$3

Because the sphere is closed, net vorticity vanishes and the minimal topological excitation is a vortex–antivortex dipole.

The main dynamical consequence is curvature-dependent symmetry breaking during free expansion. Vortex cores expand into density holes, spherical symmetry is broken for $l_r=\sqrt{\frac{\hbar}{m\omega_r}},$4, and the cloud aspect ratio

$l_r=\sqrt{\frac{\hbar}{m\omega_r}},$5

shows a non-monotonic dependence on $l_r=\sqrt{\frac{\hbar}{m\omega_r}},$6: it decreases below unity for small $l_r=\sqrt{\frac{\hbar}{m\omega_r}},$7, increases above unity for $l_r=\sqrt{\frac{\hbar}{m\omega_r}},$8, and crosses over near $l_r=\sqrt{\frac{\hbar}{m\omega_r}},$9. A complementary sphericity measure,

$r\simeq R$0

decreases with increasing $r\simeq R$1 and is nearly constant in time, $r\simeq R$2. The paper interprets this as weak nonlinear coupling between angular momentum modes and persistence of the initially imprinted angular structure during expansion (Tononi, 22 Apr 2026).

3. Relativistic and gravitational formulations

In relativistic work, the shell may be a dividing surface in a continuous spacetime, a thin matter layer between vacuum regions, or a finite-thickness anisotropic fluid distribution. In generalized Painlevé–Gullstrand coordinates, the dividing shell in anisotropic spherically symmetric spacetimes is characterized by the areal radius $r\simeq R$3, the Misner–Sharp mass $r\simeq R$4, and the conditions

$r\simeq R$5

together with the generalized equilibrium condition

$r\simeq R$6

A dividing shell separating $r\simeq R$7 from $r\simeq R$8 is defined by $r\simeq R$9 and $r\ge r_0$0; the latter is the generalized TOV balance in the presence of anisotropic stresses (Mimoso et al., 2013).

A more mechanical realization uses multilayer, infinitely thin spherical shells separated by vacuum Schwarzschild regions. For a single shell with inner mass $r\ge r_0$1, gravitational shell mass $r\ge r_0$2, and rest mass $r\ge r_0$3, the shell motion obeys

$r\ge r_0$4

With a linear equation of state $r\ge r_0$5, one has

$r\ge r_0$6

In a many-shell system, crossings are treated as totally transparent, and repeated crossings can lead to mass inflation, with asymptotic scaling

$r\ge r_0$7

This framework realizes an expanding or collapsing thick shell as the continuum limit of many concentric thin layers (Gaspar et al., 2011).

A third relativistic shell picture imposes vanishing expansion scalar, $r\ge r_0$8, on an anisotropic fluid surrounding a central Minkowski cavity and matched to an outer Schwarzschild region. The interior metric is

$r\ge r_0$9

with

$r=f(t)$0

after a radial gauge choice. The shell is then the finite-thickness anisotropic matter region between the inner vacuum cavity and the exterior vacuum, with matching conditions $r=f(t)$1, $r=f(t)$2, $r=f(t)$3, and $r=f(t)$4 for Darmois junctions (Prisco et al., 2011).

More phenomenological expanding-shell mechanics also exists. In a thin-layer approximation with porosity parameter $r=f(t)$5, a classical shell in a uniform medium obeys

$r=f(t)$6

while the relativistic porosity model gives

$r=f(t)$7

Other shell constructions remain static but are relevant as equilibrium benchmarks: conformastatic Majumdar–Papapetrou thick shells sourced by anisotropic matter, static Newtonian and GR finite-thickness shell families built by superposition and Kelvin inversion, and a cosmological shell model in which the visible universe is treated as an expanding spherical shell with radius $r=f(t)$8 Gpc and horizon arc length $r=f(t)$9 Gpc (Zaninetti, 2011, García-Reyes, 2016, Vogt et al., 2010, Vlahovic, 2012).

4. Radiation, ionization, and line formation

In radiation-diffusion theory, the shell is a cold spherical domain outside an inner radius $r(\tau)$0, driven by a hot boundary at $r(\tau)$1 and bounded by a moving Marshak front $r(\tau)$2. Under grey diffusion, LTE, constant density, and power-law constitutive relations

$r(\tau)$3

the spherical diffusion equation is

$r(\tau)$4

The extension of the Hammer–Rosen model to spherical geometry shows that geometric divergence slows the wavefront relative to planar geometry, and that the effect grows as the inner radius $r(\tau)$5 decreases. For the materials tested, the error in front position at 3 ns is $r(\tau)$6 for $r(\tau)$7 cm, whereas in $r(\tau)$8 at $r(\tau)$9 cm a planar model predicts a front about $r\simeq R$0 faster than the small-radius spherical solution (Smith et al., 11 Feb 2025).

In classical D-type H II regions, the shell is the swept-up neutral layer between a shock radius $r\simeq R$1 and an ionization front radius $r\simeq R$2. The thin-shell momentum equation is

$r\simeq R$3

while the late-time stalled configuration is set by the stagnation radius

$r\simeq R$4

A thick-shell approximation distinguishes shock and ionization-front radii, and a matched asymptotic equation captures both the early thin-shell expansion and the late-time damped oscillations about $r\simeq R$5, with the paper stating that this solution agrees very well with the numerical solution at all times and is superior to previously published solutions (Williams et al., 2018).

In radiative-transfer models for young planetary nebulae, the shell is a hollow expanding H I layer surrounding a point-like He II source. The neutral density follows

$r\simeq R$6

the radial H I column is

$r\simeq R$7

and expansion speeds of $r\simeq R$8 are considered. Raman-scattered He II $r\simeq R$9 line formation is solved with a grid-based Monte Carlo code that tracks H I density variation along each photon path. For a monochromatic He II source, the emergent profiles show an asymmetric double peak, a tertiary red peak, and a significant red tail; the peak separation corresponds to the expansion velocity, tertiary red peaks arise from multiple Rayleigh reflections at the inner surface of the hollow shell, and Raman conversion efficiency increases strongly with expansion speed because the scattering cross section rises sharply near resonance (Choi et al., 2019).

5. Instability, observables, and shell inference

In self-gravitating gas shells, finite thickness changes the instability spectrum. For a shell of surface density

$$i\hbar \partial_t\Psi = \left[-\frac{\hbar^2}{2m}\nabla^2 + g|\Psi|^2\right]\Psi$$0

the uniform-density half-thickness is

$$i\hbar \partial_t\Psi = \left[-\frac{\hbar^2}{2m}\nabla^2 + g|\Psi|^2\right]\Psi$$1

Modeling fragments as uniform oblate spheroids leads to the pressure-assisted gravitational instability dispersion relation, and the paper concludes that if the confining pressure is low, only large fragments are unstable, whereas if the confining pressure is high, fragments smaller than predicted by standard thin-shell analysis become unstable. The resulting dispersion relation is reported to be in good agreement with 3D hydrodynamic simulations (Wunsch et al., 2010).

Direct shell fitting is central in radio supernova work. For SN 2011dh, VLBI visibility fitting with a spherical shell model at $$i\hbar \partial_t\Psi = \left[-\frac{\hbar^2}{2m}\nabla^2 + g|\Psi|^2\right]\Psi$$2 d and 8.4 GHz yields an outer angular radius of $$i\hbar \partial_t\Psi = \left[-\frac{\hbar^2}{2m}\nabla^2 + g|\Psi|^2\right]\Psi$$3. At 7.8 Mpc this corresponds to a linear radius of $$i\hbar \partial_t\Psi = \left[-\frac{\hbar^2}{2m}\nabla^2 + g|\Psi|^2\right]\Psi$$4 cm and an average expansion velocity since explosion of $$i\hbar \partial_t\Psi = \left[-\frac{\hbar^2}{2m}\nabla^2 + g|\Psi|^2\right]\Psi$$5. Combining VLBI radii with synchrotron-self-absorbed spectral radii gives a power-law evolution

$$i\hbar \partial_t\Psi = \left[-\frac{\hbar^2}{2m}\nabla^2 + g|\Psi|^2\right]\Psi$$6

which implies almost free expansion over $$i\hbar \partial_t\Psi = \left[-\frac{\hbar^2}{2m}\nabla^2 + g|\Psi|^2\right]\Psi$$7 d to 453 d. The shell is almost circular in outline, although the paper allows the possibility of some asymmetry in brightness around the ridge (Witt et al., 2015).

Three-dimensional dust reconstruction provides an observational counterpart at interstellar scales. The Per–Tau Shell is identified as an extended near-spherical shell $$i\hbar \partial_t\Psi = \left[-\frac{\hbar^2}{2m}\nabla^2 + g|\Psi|^2\right]\Psi$$8 pc in diameter, centered at $$i\hbar \partial_t\Psi = \left[-\frac{\hbar^2}{2m}\nabla^2 + g|\Psi|^2\right]\Psi$$9 pc, with shell radius $U(r)=\frac{1}{2}m\omega_r^2(r-R)^2,$00 pc. Its age is estimated to be $U(r)=\frac{1}{2}m\omega_r^2(r-R)^2,$01 Myrs, based on standard spherical shell scalings and a characteristic turbulent speed $U(r)=\frac{1}{2}m\omega_r^2(r-R)^2,$02. The paper interprets it as a large expanding shell formed by previous stellar and supernova feedback, within which the swept-up ISM condensed to form the shell itself and the Perseus and Taurus molecular clouds (Bialy et al., 2021).

A different observational use of shell modeling appears in SN 2024ggi. There, a single spherical shell predicts only concave photospheric-radius evolution $U(r)=\frac{1}{2}m\omega_r^2(r-R)^2,$03, whereas the observed evolution changes from concave to convex. The paper therefore introduces a two-shell toy model composed of an outer low-mass spherical shell and an inner elongated massive shell: during the first week the spherical shell forms the photosphere, and later the elongated shell dominates it. This establishes a specific limitation of one-component spherical shells in supernova photospheric modeling (Shiran et al., 18 Dec 2025).

6. Limitations and extensions

The shell formalism is powerful precisely because it compresses complex dynamics into a geometrically privileged layer, but each implementation rests on restrictive assumptions. In the Bose–Einstein condensate problem, the analysis assumes the thin-shell approximation, weak effective 2D interactions, and zero-temperature mean-field evolution under the 3D GPE; stronger interactions, thicker shells, many-vortex configurations, and finite-temperature effects are left open. In radiation diffusion, the heat wave is assumed supersonic, so hydrodynamics is neglected. In multilayer relativistic shells, collisions are totally transparent and equations of state are simple. Several shell families in GR are explicitly static, so they provide equilibrium configurations rather than expanding dynamics (Tononi, 22 Apr 2026, Smith et al., 11 Feb 2025, Gaspar et al., 2011, García-Reyes, 2016, Vogt et al., 2010).

Spherical symmetry is also often only approximate. The Per–Tau Shell is described as near-spherical but also as an “ovaloid” shell with a broad density peak of width $U(r)=\frac{1}{2}m\omega_r^2(r-R)^2,$04 pc; SN 2011dh is almost circular in outline rather than exactly so; and SN 2024ggi requires an inner elongated shell in addition to an outer spherical one. In the quantum-superfluid context, curvature itself generates effects absent in flat geometries, such as the non-monotonic dependence of the aspect ratio on vortex dipole separation. This suggests that the shell model is most informative when used as a controlled reduction rather than as an assertion of exact symmetry (Tononi, 22 Apr 2026, Witt et al., 2015, Bialy et al., 2021, Shiran et al., 18 Dec 2025).

The surveyed literature also indicates clear extension paths. The shell-shaped BEC framework explicitly points to cylinders, cones, tori and surfaces of revolution, and ellipsoids as natural generalizations of the spherical case. The diffusion model is extended in parallel to cylindrical geometry. Static shell families in Newtonian gravity and conformastatic GR provide analytical shell profiles that could serve as initial data or equilibrium benchmarks for time-dependent problems. Taken together, these developments suggest that the spherical expanding shell model functions both as a concrete model class and as a prototype for a broader theory of curved, expanding shell-like media (Tononi, 22 Apr 2026, Smith et al., 11 Feb 2025, Vogt et al., 2010, García-Reyes, 2016).