Analytical Spherical Symmetry: Models & Applications

• General Analytical Spherically Symmetric Model is a methodological framework that exploits perfect spherical symmetry to reduce multi-dimensional phenomena to tractable radial differential equations.
• It delivers closed-form or semi-closed-form solutions across diverse fields such as photoacoustics, black-hole accretion, light propagation, and stellar interior modeling.
• While the model offers precise analytical benchmarks, its reliance on strict symmetry and idealized assumptions limits its realism in more complex, heterogeneous systems.

The expression general analytical spherically symmetric model has been used for several distinct but structurally related constructions: a unified photoacoustic forward model for arbitrary radial initial data, parametrized black-hole accretion frameworks, exact Hubble-flow distortions around a mass concentration, fast ray-mapping formulas in static spherically symmetric spacetimes, exact stellar-interior solutions, and gauge-invariant perturbation formalisms (Li et al., 12 Jan 2026, Yang et al., 2020, Baushev, 2020, Claros et al., 5 Nov 2025). This suggests that the phrase functions less as the name of a single canonical theory than as a methodological descriptor for models in which spherical symmetry reduces the governing problem to radial ODEs, shell integrals, or low-dimensional dynamical systems, while preserving enough structure to obtain closed-form or semi-closed-form results.

1. Shared meaning and analytical architecture

Across the cited literatures, spherical symmetry is imposed at the level of the unknown fields rather than only at the level of observables. In wave problems this means $p_0(\mathbf r)=p_0(r)$; in gravitation it means metrics of the form $ds^2=-A(r)dt^2+B(r)dr^2+r^2d\Omega^2$, $$ds^2=-N^2(r)dt^2+\frac{B^2(r)}{N^2(r)}dr^2+r^2d\Omega^2$$, or related generalizations; in Finsler geometry it means $F=r\,\phi(u,s,v,t)$; and in perturbation theory it means a static spherical background with all perturbations decomposed into spherical harmonics (Li et al., 12 Jan 2026, Claros et al., 5 Nov 2025, Yang et al., 2020, Khani-Moghaddam et al., 2024, Liu et al., 2022). The analytical gain is always the same: angular dependence is either integrated out exactly or encoded in harmonics, leaving a radial problem with explicit coefficients.

A second common feature is that these models are usually exact within a sharply delimited regime. The photoacoustic model assumes linear acoustics, an infinite homogeneous lossless medium, instantaneous heating, and spherical symmetry; the accretion models assume steady, radial, test-fluid flow; the light-propagation formulas assume static spherically symmetric metrics and direct photon trajectories; the stellar models assume staticity and spherical symmetry; and the perturbation framework assumes a static spherical background (Li et al., 12 Jan 2026, Yang et al., 2020, Claros et al., 5 Nov 2025, Pant, 2010, Liu et al., 2022). Analytical tractability is therefore obtained by symmetry reduction plus restrictive constitutive assumptions.

Domain	Representative ansatz	Principal analytical output
Photoacoustics	$p(r,t)$, $p_0(r)$	Closed-form pressure field
Black-hole accretion	$N(r), B(r)$ or $A(r),B(r),C(r)$	First integrals, sonic conditions
Light propagation	$A(r),B(r)$	Emission-angle/image-plane mapping
Stellar interiors	$e^{\nu(r)}, e^{\lambda(r)}$ or spheroidal geometry	Exact mass-radius and stability relations
Perturbations	$ds^2=-A(r)dt^2+B(r)dr^2+r^2d\Omega^2$0 background	Gauge-invariant master equations
Finsler geometry	$ds^2=-A(r)dt^2+B(r)dr^2+r^2d\Omega^2$1	Curvature classification and PDE constraints

A common misconception is that spherical symmetry forces the Schwarzschild gauge $ds^2=-A(r)dt^2+B(r)dr^2+r^2d\Omega^2$2. Several of these works explicitly avoid that restriction: general accretion uses either $ds^2=-A(r)dt^2+B(r)dr^2+r^2d\Omega^2$3 or $ds^2=-A(r)dt^2+B(r)dr^2+r^2d\Omega^2$4 independently, light propagation treats arbitrary $ds^2=-A(r)dt^2+B(r)dr^2+r^2d\Omega^2$5 and $ds^2=-A(r)dt^2+B(r)dr^2+r^2d\Omega^2$6, and perturbation theory is developed for the fully general static form $ds^2=-A(r)dt^2+B(r)dr^2+r^2d\Omega^2$7 (Bahamonde et al., 2015, Claros et al., 5 Nov 2025, Liu et al., 2022).

2. Radial wave models in photoacoustics

In photoacoustics, the term denotes a unified analytical solution for the acoustic pressure generated by a spherically symmetric initial pressure distribution under linear acoustics, constant sound speed $ds^2=-A(r)dt^2+B(r)dr^2+r^2d\Omega^2$8, no attenuation, no dispersion, no boundaries, an infinite medium, and instantaneous heating (Li et al., 12 Jan 2026). Starting from the standard wave equation and a Green’s-function representation, the angular integral can be performed exactly, yielding a shell-integral form and then the compact result

$ds^2=-A(r)dt^2+B(r)dr^2+r^2d\Omega^2$9

This formula is presented as the general spherically symmetric analytical model because it applies to any radial initial profile $$ds^2=-N^2(r)dt^2+\frac{B^2(r)}{N^2(r)}dr^2+r^2d\Omega^2$$0 for all $$ds^2=-N^2(r)dt^2+\frac{B^2(r)}{N^2(r)}dr^2+r^2d\Omega^2$$1 and $$ds^2=-N^2(r)dt^2+\frac{B^2(r)}{N^2(r)}dr^2+r^2d\Omega^2$$2 (Li et al., 12 Jan 2026).

Its interpretation is geometrically transparent. The term $$ds^2=-N^2(r)dt^2+\frac{B^2(r)}{N^2(r)}dr^2+r^2d\Omega^2$$3 is identified with a converging wave, while $$ds^2=-N^2(r)dt^2+\frac{B^2(r)}{N^2(r)}dr^2+r^2d\Omega^2$$4 is identified with a diverging wave. The absolute value handles transmission through the center. The prefactor $$ds^2=-N^2(r)dt^2+\frac{B^2(r)}{N^2(r)}dr^2+r^2d\Omega^2$$5 encodes spherical spreading. The same derivation also yields the shell-integral representation

$$ds^2=-N^2(r)dt^2+\frac{B^2(r)}{N^2(r)}dr^2+r^2d\Omega^2$$6

which is useful near $$ds^2=-N^2(r)dt^2+\frac{B^2(r)}{N^2(r)}dr^2+r^2d\Omega^2$$7 where the closed form is only apparently singular (Li et al., 12 Jan 2026).

The model is then specialized to several common source profiles: a uniform sphere, Gaussian, exponential, and power-law initial pressures. In the far field, the converging term becomes negligible near the arrival time $$ds^2=-N^2(r)dt^2+\frac{B^2(r)}{N^2(r)}dr^2+r^2d\Omega^2$$8, and the signal reduces to

$$ds^2=-N^2(r)dt^2+\frac{B^2(r)}{N^2(r)}dr^2+r^2d\Omega^2$$9

so the measured waveform directly mirrors the radial source profile up to the factor $F=r\,\phi(u,s,v,t)$0 (Li et al., 12 Jan 2026). The stated applications are photoacoustic imaging system design, signal interpretation, and benchmarking of numerical solvers; the implementation is provided in the repository SlingBAG_Ultra (Li et al., 12 Jan 2026).

3. Gravitational transport, accretion, and radial observables

In relativistic accretion theory, a general analytical spherically symmetric model typically means a background metric plus fluid equations reduced to a pair of first integrals and a critical-point condition. One such formulation uses the Rezzolla–Zhidenko parametrization of a generic static spherically symmetric black-hole spacetime,

$F=r\,\phi(u,s,v,t)$1

with deviations from Schwarzschild encoded by $F=r\,\phi(u,s,v,t)$2 (Yang et al., 2020). For stationary radial flow of a perfect fluid, baryon conservation and energy–momentum conservation give

$F=r\,\phi(u,s,v,t)$3

and the sonic point satisfies $F=r\,\phi(u,s,v,t)$4 together with a metric-dependent regularity condition (Yang et al., 2020). A particularly sharp result is that for ideal photon gas the sonic radius obeys

$F=r\,\phi(u,s,v,t)$5

which coincides exactly with the photon-sphere condition, so the sonic radius equals the photon-sphere radius for any metric of this class (Yang et al., 2020).

A more metric-agnostic accretion framework starts from

$F=r\,\phi(u,s,v,t)$6

and derives general first integrals for a perfect fluid with barotropic equation of state, together with critical-point formulas such as

$F=r\,\phi(u,s,v,t)$7

This makes spherical accretion comparable across black holes, naked singularities, and other compact objects without rederiving the conservation equations case by case (Bahamonde et al., 2015).

The same analytical logic appears in cosmological radial flows. For the Hubble stream perturbed by a spherically symmetric mass concentration, shell dynamics reduce to a Friedmann-like equation with effective density parameters $F=r\,\phi(u,s,v,t)$8, $F=r\,\phi(u,s,v,t)$9, and $p(r,t)$0, leading to exact implicit relations for the present-day velocity $p(r,t)$1 in terms of the present radius $p(r,t)$2 (Baushev, 2020). The dimensionless parameters

$p(r,t)$3

organize both the expanding and infalling branches. This is a different physical problem from accretion, but it exemplifies the same reduction of spherical dynamics to a one-dimensional analytical model (Baushev, 2020).

For null transport, the recent light-propagation formalism for generic static spherically symmetric metrics

$p(r,t)$4

recasts the bending integral into a small set of metric-dependent radial integrals $p(r,t)$5 and then derives analytical relations between $p(r,t)$6 and $p(r,t)$7, including the linear generalized Beloborodov formula

$p(r,t)$8

The full approximation is then used for fast computation of accretion-disk images, polarization, and pulsar luminosity curves in Johannsen–Psaltis, Rezzolla–Zhidenko, and EMDA spacetimes (Claros et al., 5 Nov 2025).

4. Static compact objects, interiors, and shells

In stellar structure, the phrase often refers to exact interior solutions of Einstein’s equations under static spherical symmetry. One example is a parametric class of perfect-fluid balls in canonical coordinates,

$p(r,t)$9

constructed by introducing $p_0(r)$0 and $p_0(r)$1, reducing Tolman’s equation to a linear ODE, and then imposing an integrability ansatz that yields explicit forms for $p_0(r)$2, $p_0(r)$3, pressure, and density (Pant, 2010). The resulting models have infinite central pressure and density, but the central ratio $p_0(r)$4 remains finite and less than one over the physically admissible parameter range (Pant, 2010).

A more geometry-driven construction uses a static, spherically symmetric interior whose constant-time slices form a three-spheroid. In that model the metric coefficient $p_0(r)$5 is fixed by spheroidal geometry, $p_0(r)$6 follows from a second-order ODE, and the density ratio

$p_0(r)$7

controls the mass-radius relation (R et al., 6 Jan 2026). The model satisfies the standard energy conditions and hydrostatic equilibrium, and the radial stability analysis selects $p_0(r)$8 as the value required for dynamical stability (R et al., 6 Jan 2026).

Thin-shell traversable wormholes provide a dynamic shell-based version of the same analytical program. Two arbitrary static spherically symmetric bulks are joined across a timelike shell with radius $p_0(r)$9, and the Lanczos junction conditions yield the shell stress tensor $N(r), B(r)$0 (Garcia et al., 2011). The throat dynamics reduce to a mechanical equation

$N(r), B(r)$1

and linear stability of a static throat $N(r), B(r)$2 is equivalent to $N(r), B(r)$3 (Garcia et al., 2011). The broad conclusion is that stability is equivalent to suitable properties of the exotic material on the throat, so the analytical model converts a geometric matching problem into an inequality for the shell equation of state (Garcia et al., 2011).

5. Generalized geometries and modified-gravity realizations

In Finsler geometry, a general spherically symmetric analytical model is not a spacetime metric but a Finsler structure

$N(r), B(r)$4

with a fixed constant $N(r), B(r)$5-form $N(r), B(r)$6 and a smooth function $N(r), B(r)$7 (Khani-Moghaddam et al., 2024). Within this class, every such metric is proved to be semi C-reducible, and the condition for vanishing mean stretch curvature is expressed as a PDE involving the antisymmetric combinations $N(r), B(r)$8, $N(r), B(r)$9, and $A(r),B(r),C(r)$0 (Khani-Moghaddam et al., 2024). Here spherical symmetry is broadened into an invariant ansatz for non-Riemannian curvature tensors.

Modified gravity supplies several further meanings. In the $A(r),B(r),C(r)$1-$A(r),B(r),C(r)$2 model, the most general spherically symmetric vacuum solution is obtained by a full Hamiltonian reduction of the ADM action with preferred foliation; the constraint algebra generates a CMC or maximal slicing condition as a tertiary constraint, and the resulting solutions are generally non-static, incompatible with asymptotic flatness, and parametrized by a conserved mass together with foliation data such as the mean curvature $A(r),B(r),C(r)$3 and a transverse-traceless parameter $A(r),B(r),C(r)$4 (Loll et al., 2017). The four-dimensional Ricci scalar becomes

$A(r),B(r),C(r)$5

so for $A(r),B(r),C(r)$6 and $A(r),B(r),C(r)$7 the spacetime is generically not Ricci-flat (Loll et al., 2017).

In Cotton gravity, the general static spherically symmetric vacuum solution is written in terms of two metric functions $A(r),B(r),C(r)$8 and $A(r),B(r),C(r)$9, with

$A(r),B(r)$0

and $A(r),B(r)$1 containing two additional integration constants (Gogberashvili et al., 2023). The resulting geometry has a striking difference from Schwarzschild: curvature invariants acquire a singularity at the photon sphere $A(r),B(r)$2, which the authors argue probably obstructs the formation of stellar Schwarzschild-radius black holes (Gogberashvili et al., 2023). In the weak-field limit, the model also generates a repulsive velocity-squared term in the radial force law, proposed as a dark-energy-like effect (Gogberashvili et al., 2023).

Lorentz-breaking massive gravity yields yet another analytical classification. With the hedgehog Stückelberg ansatz

$A(r),B(r)$3

the requirement $A(r),B(r)$4 reduces to $A(r),B(r)$5, and the admissible potentials form a commutative ring $A(r),B(r)$6 (Li et al., 2015). A subring $A(r),B(r)$7 with $A(r),B(r)$8 yields Schwarzschild, AdS, or dS metrics, while $A(r),B(r)$9 produces new solutions such as furry black holes and logarithmically corrected metrics (Li et al., 2015). This is a rare case where the analytical spherically symmetric model is organized algebraically, at the level of the action rather than only at the level of solutions.

6. Perturbative completions, applications, and recurring limitations

A general analytical spherically symmetric model is often most useful when it supports a perturbation theory. For a static background

$e^{\nu(r)}, e^{\lambda(r)}$0

metric perturbations can be decomposed in the A–K basis, gauge-transformation laws can be written explicitly, and gauge-invariant master variables can be constructed in both the Detweiler easy gauge and the Regge–Wheeler gauge (Liu et al., 2022). For $e^{\nu(r)}, e^{\lambda(r)}$1, one obtains odd- and even-parity master equations whose coefficients are explicit functionals of $e^{\nu(r)}, e^{\lambda(r)}$2 and $e^{\nu(r)}, e^{\lambda(r)}$3, reducing to the Regge–Wheeler and Zerilli equations in the Schwarzschild limit (Liu et al., 2022). This provides a reusable analytical perturbation toolkit for non-vacuum black holes, stars, and Effective-One-Body backgrounds.

The practical uses emphasized across these works are rapid forward modeling, parameter inference, benchmarking, and analytical control of observables. In photoacoustics, the closed form is used for ultrafast forward simulation and signal interpretation (Li et al., 12 Jan 2026). In light propagation, the analytic ray-mapping formulas enable fast disk imaging and polarization calculations in non-Schwarzschild metrics (Claros et al., 5 Nov 2025). In stellar and shell models, exact solutions clarify which energy conditions, matching conditions, or stability inequalities are structurally required (R et al., 6 Jan 2026, Garcia et al., 2011).

Several limitations recur. Exactness is usually contingent on strong restrictions: perfect spherical symmetry; staticity or stationarity; homogeneous, lossless, or nondispersive media; test-fluid or test-field approximations; or highly specific constitutive assumptions (Li et al., 12 Jan 2026, Yang et al., 2020, Claros et al., 5 Nov 2025, Pant, 2010). Analytical does not necessarily mean globally regular: some stellar models have infinite central density, Cotton gravity develops a photosphere singularity, and thin-shell wormholes require exotic surface matter (Pant, 2010, Gogberashvili et al., 2023, Garcia et al., 2011). This suggests that the enduring value of the general analytical spherically symmetric model lies less in universal realism than in its role as a controllable benchmark: it isolates the consequences of symmetry, identifies which quantities are fixed by geometry and which by matter or closure relations, and supplies closed or nearly closed expressions against which more general numerical models can be tested.