Tolman IV: Analytic Stellar Interior Model

• Tolman IV solution is an exact static, spherically symmetric model for relativistic stellar interiors with an isotropic perfect fluid.
• It serves as a benchmark for testing equilibrium, thermodynamic acceptability, and as a seed for anisotropic and modified gravity extensions.
• Analytic expressions for density, pressure, and temperature enable matching with the exterior Schwarzschild metric and application to observed compact stars.

Tolman IV solution is an exact static, spherically symmetric interior solution of Einstein’s equations for an isotropic perfect fluid, and it remains a standard analytic model of relativistic stellar interiors. In current literature it serves several distinct roles: a benchmark for regular finite-radius perfect-fluid configurations, an explicit test case for extending physical acceptability to thermodynamic acceptability, a fully analytic model for compact objects with observed masses and radii, and a seed geometry for anisotropic, charged, braneworld, bimetric, and matter-geometry-coupled deformations (Demissenova et al., 21 May 2026, Panotopoulos, 2024, Andrade et al., 2021).

1. Canonical general-relativistic form

In the standard static spherically symmetric line element,

$$ds^2=-e^{\nu(r)}dt^2 + e^{\lambda(r)}dr^2 + r^2(d\theta^2+\sin^2\theta\, d\phi^2),$$

the Tolman IV interior is specified by

$$e^{\nu}=B^2\left(1+\frac{r^2}{A^2}\right), \qquad e^{-\lambda}= \frac{\left(1-r^2/C^2\right)\left(1+r^2/A^2\right)}{1+2r^2/A^2},$$

with perfect-fluid stress tensor

$$T_{ab}=(\rho+p)u_a u_b + pg_{ab}, \qquad u^a=e^{-\nu/2}\delta^a_0 .$$

For the parametrization used in the thermodynamic analysis, the matter variables are

$$\rho(r)=\frac{3A^2+2r^2}{8\pi(A^2+2r^2)^2}, \qquad p(r)=\frac{C^2-A^2-3r^2}{8\pi(A^2+2r^2)(C^2-r^2)},$$

with central values

$$\rho(0)=\frac{3}{8\pi A^2}, \qquad p(0)=\frac{C^2-A^2}{8\pi A^2C^2},$$

so the central pressure is positive if $C^2>A^2$ (Demissenova et al., 21 May 2026).

The stellar surface is defined by $p(R)=0$, which gives

$C^2=A^2+3R^2.$

Matching to the exterior Schwarzschild solution requires

$$e^{\nu(R)}=1-\frac{2M}{R}, \qquad e^{-\lambda(R)}=1-\frac{2M}{R},$$

hence

$$B^2=\frac{A^2}{A^2+3R^2}, \qquad M=\frac{3R^3}{A^2+3R^2}.$$

In that form the solution has a finite radius, a finite mass, and a smooth match to vacuum (Demissenova et al., 21 May 2026).

A stellar-modeling formulation uses the same metric potentials with constants $$e^{\nu}=B^2\left(1+\frac{r^2}{A^2}\right), \qquad e^{-\lambda}= \frac{\left(1-r^2/C^2\right)\left(1+r^2/A^2\right)}{1+2r^2/A^2},$$0, $$e^{\nu}=B^2\left(1+\frac{r^2}{A^2}\right), \qquad e^{-\lambda}= \frac{\left(1-r^2/C^2\right)\left(1+r^2/A^2\right)}{1+2r^2/A^2},$$1, and $$e^{\nu}=B^2\left(1+\frac{r^2}{A^2}\right), \qquad e^{-\lambda}= \frac{\left(1-r^2/C^2\right)\left(1+r^2/A^2\right)}{1+2r^2/A^2},$$2,

$$e^{\nu}=B^2\left(1+\frac{r^2}{A^2}\right), \qquad e^{-\lambda}= \frac{\left(1-r^2/C^2\right)\left(1+r^2/A^2\right)}{1+2r^2/A^2},$$3

and fixes the three constants by imposing

$$e^{\nu}=B^2\left(1+\frac{r^2}{A^2}\right), \qquad e^{-\lambda}= \frac{\left(1-r^2/C^2\right)\left(1+r^2/A^2\right)}{1+2r^2/A^2},$$4

once the stellar mass $$e^{\nu}=B^2\left(1+\frac{r^2}{A^2}\right), \qquad e^{-\lambda}= \frac{\left(1-r^2/C^2\right)\left(1+r^2/A^2\right)}{1+2r^2/A^2},$$5 and radius $$e^{\nu}=B^2\left(1+\frac{r^2}{A^2}\right), \qquad e^{-\lambda}= \frac{\left(1-r^2/C^2\right)\left(1+r^2/A^2\right)}{1+2r^2/A^2},$$6 are known observationally (Panotopoulos, 2024).

2. Regularity, acceptability, and a standard caveat

Tolman IV is repeatedly singled out as a physically acceptable relativistic interior because it has regular metric functions at the center, finite $$e^{\nu}=B^2\left(1+\frac{r^2}{A^2}\right), \qquad e^{-\lambda}= \frac{\left(1-r^2/C^2\right)\left(1+r^2/A^2\right)}{1+2r^2/A^2},$$7 and $$e^{\nu}=B^2\left(1+\frac{r^2}{A^2}\right), \qquad e^{-\lambda}= \frac{\left(1-r^2/C^2\right)\left(1+r^2/A^2\right)}{1+2r^2/A^2},$$8, positive density and pressure for suitable parameters, monotonic decrease of density and pressure outward, and a finite boundary where $$e^{\nu}=B^2\left(1+\frac{r^2}{A^2}\right), \qquad e^{-\lambda}= \frac{\left(1-r^2/C^2\right)\left(1+r^2/A^2\right)}{1+2r^2/A^2},$$9 (Demissenova et al., 21 May 2026).

A more explicit causality check is available in the compact-star modeling literature, where the squared sound speed is

$$T_{ab}=(\rho+p)u_a u_b + pg_{ab}, \qquad u^a=e^{-\nu/2}\delta^a_0 .$$0

which satisfies

$$T_{ab}=(\rho+p)u_a u_b + pg_{ab}, \qquad u^a=e^{-\nu/2}\delta^a_0 .$$1

throughout the interior. The same analysis also lists the standard energy conditions,

$$T_{ab}=(\rho+p)u_a u_b + pg_{ab}, \qquad u^a=e^{-\nu/2}\delta^a_0 .$$2

together with the regularity condition $$T_{ab}=(\rho+p)u_a u_b + pg_{ab}, \qquad u^a=e^{-\nu/2}\delta^a_0 .$$3, as criteria that Tolman IV meets in the isotropic setting (Panotopoulos, 2024).

The principal caveat emphasized in the recent stellar-modeling treatment is not a breakdown of the interior geometry but a surface effect in the one-layer model. The relativistic adiabatic index,

$$T_{ab}=(\rho+p)u_a u_b + pg_{ab}, \qquad u^a=e^{-\nu/2}\delta^a_0 .$$4

diverges at the stellar surface because $$T_{ab}=(\rho+p)u_a u_b + pg_{ab}, \qquad u^a=e^{-\nu/2}\delta^a_0 .$$5 while $$T_{ab}=(\rho+p)u_a u_b + pg_{ab}, \qquad u^a=e^{-\nu/2}\delta^a_0 .$$6. That divergence is explicitly interpreted not as an interior pathology but as a consequence of using a one-layer configuration with nonzero surface density; the proposed resolution is the inclusion of a thin crust or envelope with a polytropic equation of state so that both $$T_{ab}=(\rho+p)u_a u_b + pg_{ab}, \qquad u^a=e^{-\nu/2}\delta^a_0 .$$7 and $$T_{ab}=(\rho+p)u_a u_b + pg_{ab}, \qquad u^a=e^{-\nu/2}\delta^a_0 .$$8 at the physical surface (Panotopoulos, 2024).

3. Thermodynamic acceptability and the Tolman equilibrium law

A major recent development is the extension of the Delgaty–Lake physical acceptability program to thermodynamic acceptability. In addition to regularity, positivity of matter variables, finite radius, and causal sound propagation, the proposed thermodynamic criteria require a positive and regular entropy density, finite total entropy, positive temperature, compatibility with the relativistic Gibbs relation, satisfaction of the Tolman equilibrium law, regularity of thermodynamic quantities, and monotonic outward behavior of thermodynamic variables (Demissenova et al., 21 May 2026).

The starting thermodynamic identities are

$$T_{ab}=(\rho+p)u_a u_b + pg_{ab}, \qquad u^a=e^{-\nu/2}\delta^a_0 .$$9

For vanishing chemical potential $$\rho(r)=\frac{3A^2+2r^2}{8\pi(A^2+2r^2)^2}, \qquad p(r)=\frac{C^2-A^2-3r^2}{8\pi(A^2+2r^2)(C^2-r^2)},$$0,

$$\rho(r)=\frac{3A^2+2r^2}{8\pi(A^2+2r^2)^2}, \qquad p(r)=\frac{C^2-A^2-3r^2}{8\pi(A^2+2r^2)(C^2-r^2)},$$1

The total entropy of a static spherical star is then

$$\rho(r)=\frac{3A^2+2r^2}{8\pi(A^2+2r^2)^2}, \qquad p(r)=\frac{C^2-A^2-3r^2}{8\pi(A^2+2r^2)(C^2-r^2)},$$2

Thermal equilibrium in curved spacetime is imposed through the Tolman law,

$$\rho(r)=\frac{3A^2+2r^2}{8\pi(A^2+2r^2)^2}, \qquad p(r)=\frac{C^2-A^2-3r^2}{8\pi(A^2+2r^2)(C^2-r^2)},$$3

For Tolman IV this gives

$$\rho(r)=\frac{3A^2+2r^2}{8\pi(A^2+2r^2)^2}, \qquad p(r)=\frac{C^2-A^2-3r^2}{8\pi(A^2+2r^2)(C^2-r^2)},$$4

which is monotonic decreasing with $$\rho(r)=\frac{3A^2+2r^2}{8\pi(A^2+2r^2)^2}, \qquad p(r)=\frac{C^2-A^2-3r^2}{8\pi(A^2+2r^2)(C^2-r^2)},$$5. The reported physical interpretation is precise: this profile implies no heat flow, vanishing entropy production, and genuine global thermal equilibrium. When the Tolman IV metric and matter variables are substituted into the entropy integral, the resulting entropy functional has an integrand that is finite throughout the interior, positive entropy density when pressure and density are positive, and finite, positive, regular total entropy for physically reasonable parameter choices. In the normalized case $$\rho(r)=\frac{3A^2+2r^2}{8\pi(A^2+2r^2)^2}, \qquad p(r)=\frac{C^2-A^2-3r^2}{8\pi(A^2+2r^2)(C^2-r^2)},$$6, the entropy is finite, positive, monotonically increasing with radius, and asymptotically behaves like $$\rho(r)=\frac{3A^2+2r^2}{8\pi(A^2+2r^2)^2}, \qquad p(r)=\frac{C^2-A^2-3r^2}{8\pi(A^2+2r^2)(C^2-r^2)},$$7 for sufficiently large radii (Demissenova et al., 21 May 2026).

The same study also considers an alternative local ideal-gas prescription,

$$\rho(r)=\frac{3A^2+2r^2}{8\pi(A^2+2r^2)^2}, \qquad p(r)=\frac{C^2-A^2-3r^2}{8\pi(A^2+2r^2)(C^2-r^2)},$$8

which yields a local temperature profile that is locally meaningful but does not satisfy the Tolman equilibrium condition globally, vanishes at $$\rho(r)=\frac{3A^2+2r^2}{8\pi(A^2+2r^2)^2}, \qquad p(r)=\frac{C^2-A^2-3r^2}{8\pi(A^2+2r^2)(C^2-r^2)},$$9 because $$\rho(0)=\frac{3}{8\pi A^2}, \qquad p(0)=\frac{C^2-A^2}{8\pi A^2C^2},$$0, and produces a non-equilibrium entropy functional with a non-integrable singularity at the boundary. The boundary divergence is explicitly interpreted as incompatibility with equilibrium rather than as a pathology of the spacetime geometry. Under the Tolman temperature choice, by contrast, Tolman IV is judged thermodynamically acceptable in the equilibrium sense (Demissenova et al., 21 May 2026).

4. Use as an analytic model for observed compact stars

Tolman IV has been used as a fully analytic isotropic model for compact objects with measured masses and radii. A recent application considers PSR J0740+6620 and HESS J1731-347, with the Tolman IV constants fixed entirely by the matching conditions (Panotopoulos, 2024).

For PSR J0740+6620, the observational inputs are

$$\rho(0)=\frac{3}{8\pi A^2}, \qquad p(0)=\frac{C^2-A^2}{8\pi A^2C^2},$$1

which yield

$$\rho(0)=\frac{3}{8\pi A^2}, \qquad p(0)=\frac{C^2-A^2}{8\pi A^2C^2},$$2

The reported stability values are

$$\rho(0)=\frac{3}{8\pi A^2}, \qquad p(0)=\frac{C^2-A^2}{8\pi A^2C^2},$$3

For HESS J1731-347, the corresponding data are

$$\rho(0)=\frac{3}{8\pi A^2}, \qquad p(0)=\frac{C^2-A^2}{8\pi A^2C^2},$$4

giving

$$\rho(0)=\frac{3}{8\pi A^2}, \qquad p(0)=\frac{C^2-A^2}{8\pi A^2C^2},$$5

with

$$\rho(0)=\frac{3}{8\pi A^2}, \qquad p(0)=\frac{C^2-A^2}{8\pi A^2C^2},$$6

In both cases the model is reported to have a regular center, finite central pressure and density, monotonic decrease of $$\rho(0)=\frac{3}{8\pi A^2}, \qquad p(0)=\frac{C^2-A^2}{8\pi A^2C^2},$$7 and $$\rho(0)=\frac{3}{8\pi A^2}, \qquad p(0)=\frac{C^2-A^2}{8\pi A^2C^2},$$8, causal sound propagation, satisfaction of the energy conditions, and $$\rho(0)=\frac{3}{8\pi A^2}, \qquad p(0)=\frac{C^2-A^2}{8\pi A^2C^2},$$9. In the same work, Tolman IV is contrasted with the Kohler–Chao solution, whose pressure and density never vanish at a finite radius and therefore do not define a proper stellar surface. Within that comparison, Tolman IV is the analytically tractable interior that supports hydrostatic equilibrium at a finite boundary (Panotopoulos, 2024).

5. Seed geometry for anisotropy and complexity-driven constructions

Tolman IV also functions as a seed solution in anisotropic model building. In the gravitational-decoupling framework with minimal geometric deformation, the isotropic Tolman IV configuration is used to extract a prototype Herrera complexity factor,

$C^2>A^2$0

which for the Tolman IV seed becomes

$C^2>A^2$1

That expression is then generalized to the “like-Tolman IV” condition

$C^2>A^2$2

and imposed as the supplementary condition that closes the decoupling system (Andrade et al., 2021).

In this construction the Tolman IV seed does not simply reappear unchanged. Rather, the complexity factor becomes the engine that generates anisotropy while preserving the Tolman IV seed geometry in deformed form. The resulting Tolman IV-based anisotropic model is reported to be physically acceptable: $C^2>A^2$3, $C^2>A^2$4, and $C^2>A^2$5 are finite at the center and decrease monotonically, $C^2>A^2$6, $C^2>A^2$7 for $C^2>A^2$8, the dominant energy condition holds, the radial and tangential sound speeds remain subluminal, and the redshift decreases outward. At the same time, the Schwarzschild matching conditions imply the compactness restriction

$C^2>A^2$9

which is stricter than Buchdahl’s isotropic bound $p(R)=0$0, so the interval

$p(R)=0$1

is forbidden in that Tolman IV-based anisotropic extension. When tested against the compactness values of SMC X-1 and Cen X-3, the model remains physically well behaved but is not the closest density-ratio fit among the four seed choices studied (Andrade et al., 2021).

A plausible implication is that Tolman IV also sits naturally inside the covariant anisotropic Tolman–Oppenheimer–Volkoff generating and reconstruction frameworks for regular isotropic stars, even though it is not written out explicitly there. The anisotropic covariant formalism is designed to deform known isotropic interiors with finite-radius surfaces satisfying $p(R)=0$2, and Tolman IV fits that class of regular isotropic seeds (Carloni et al., 2017).

6. Extensions in modified gravity, charged matter, and local interaction physics

The Tolman IV geometry has been transplanted into several non-GR settings, usually by preserving the analytic structure of the GR interior and attributing deviations to an extra source, nonminimal coupling, or local curvature correction.

Context	Tolman IV role	Reported effect
Randall–Sundrum braneworld (Ovalle et al., 2013)	GR seed for minimal geometric deformation	Positive radial deformation $p(R)=0$3 lowers effective mass and compactness
Bigravity (Singh et al., 2020)	Physical interior metric on a de Sitter background	Lower $p(R)=0$4 gives a stiffer EoS; $p(R)=0$5 recovers the GR-like limit
$p(R)=0$6 (Bhar et al., 2021)	Exact Tolman IV perfect-fluid interior for LMC X-4	Larger $p(R)=0$7 lowers mass, compactness, and redshift, but raises sound speed and adiabatic index
Charged $p(R)=0$8 (Sharif et al., 2023)	Exact interior ansatz with MIT bag EOS and Reissner–Nordström matching	Only Her X-I and 4U 1820-30 remain stable for both $p(R)=0$9
Anisotropic non-minimally coupled theory (Sharif et al., 2024)	Tolman IV spacetime plus MIT bag EOS for six compact stars	Model I is physically acceptable for both $C^2=A^2+3R^2.$0; model II is stable only for $C^2=A^2+3R^2.$1
Curvature-corrected Yukawa interaction (Zamperlini et al., 26 Mar 2026)	Regular interior metric for local inertial-frame quantum corrections	Radial symmetry is preserved locally; shifts are of order $C^2=A^2+3R^2.$2 for solution IV

In the braneworld construction, Tolman IV is retained as the perfect-fluid seed while only the radial metric component is deformed, yielding an exact analytic interior on the brane. The deformation is positive, the effective mass

$C^2=A^2+3R^2.$3

is reduced, and the broader conclusion is that bulk gravitons soften stellar configurations by reducing compactness (Ovalle et al., 2013).

In $C^2=A^2+3R^2.$4 gravity, the Tolman IV interior remains regular and singularity-free for the LMC X-4 application. The paper reports that increasing $C^2=A^2+3R^2.$5 decreases $C^2=A^2+3R^2.$6, increases $C^2=A^2+3R^2.$7 and $C^2=A^2+3R^2.$8, lowers $C^2=A^2+3R^2.$9, $$e^{\nu(R)}=1-\frac{2M}{R}, \qquad e^{-\lambda(R)}=1-\frac{2M}{R},$$0, $$e^{\nu(R)}=1-\frac{2M}{R}, \qquad e^{-\lambda(R)}=1-\frac{2M}{R},$$1, compactness, and surface redshift, and increases both the sound speed and the adiabatic index. For the tabulated LMC X-4 model, $$e^{\nu(R)}=1-\frac{2M}{R}, \qquad e^{-\lambda(R)}=1-\frac{2M}{R},$$2 rises from $$e^{\nu(R)}=1-\frac{2M}{R}, \qquad e^{-\lambda(R)}=1-\frac{2M}{R},$$3 at $$e^{\nu(R)}=1-\frac{2M}{R}, \qquad e^{-\lambda(R)}=1-\frac{2M}{R},$$4 to $$e^{\nu(R)}=1-\frac{2M}{R}, \qquad e^{-\lambda(R)}=1-\frac{2M}{R},$$5 at $$e^{\nu(R)}=1-\frac{2M}{R}, \qquad e^{-\lambda(R)}=1-\frac{2M}{R},$$6, while the compactness decreases from $$e^{\nu(R)}=1-\frac{2M}{R}, \qquad e^{-\lambda(R)}=1-\frac{2M}{R},$$7 to $$e^{\nu(R)}=1-\frac{2M}{R}, \qquad e^{-\lambda(R)}=1-\frac{2M}{R},$$8 (Bhar et al., 2021). In bigravity, the same Tolman IV interior is coupled to a constant-curvature de Sitter background; lower $$e^{\nu(R)}=1-\frac{2M}{R}, \qquad e^{-\lambda(R)}=1-\frac{2M}{R},$$9 yields a stiffer EoS, whereas increasing $$B^2=\frac{A^2}{A^2+3R^2}, \qquad M=\frac{3R^3}{A^2+3R^2}.$$0 weakens the coupling and moves the model toward the Minkowski-background or GR-like limit (Singh et al., 2020).

Charged matter-geometry-coupled models likewise employ the Tolman IV ansatz as the exact interior geometry. In the $$B^2=\frac{A^2}{A^2+3R^2}, \qquad M=\frac{3R^3}{A^2+3R^2}.$$1 analysis, Tolman IV is combined with the MIT bag EOS $$B^2=\frac{A^2}{A^2+3R^2}, \qquad M=\frac{3R^3}{A^2+3R^2}.$$2 and Reissner–Nordström matching; charge reduces the effective density and anisotropy and modifies compactness and redshift, while only Her X-I and 4U 1820-30 remain stable for both signs of $$B^2=\frac{A^2}{A^2+3R^2}, \qquad M=\frac{3R^3}{A^2+3R^2}.$$3 (Sharif et al., 2023). In the related non-minimally coupled study based on $$B^2=\frac{A^2}{A^2+3R^2}, \qquad M=\frac{3R^3}{A^2+3R^2}.$$4 and $$B^2=\frac{A^2}{A^2+3R^2}, \qquad M=\frac{3R^3}{A^2+3R^2}.$$5, Tolman IV again provides a regular interior with analytically determined constants $$B^2=\frac{A^2}{A^2+3R^2}, \qquad M=\frac{3R^3}{A^2+3R^2}.$$6; model I yields physically acceptable structures for all six stars and both $$B^2=\frac{A^2}{A^2+3R^2}, \qquad M=\frac{3R^3}{A^2+3R^2}.$$7 values, whereas model II is stable only for $$B^2=\frac{A^2}{A^2+3R^2}, \qquad M=\frac{3R^3}{A^2+3R^2}.$$8 (Sharif et al., 2024).

A distinct line of work uses a Tolman IV-type gravitational potential rather than the standard isotropic perfect-fluid solution. In the Einstein–Maxwell quark-star model with quadratic EOS $$B^2=\frac{A^2}{A^2+3R^2}, \qquad M=\frac{3R^3}{A^2+3R^2}.$$9, the ansatz

$$e^{\nu}=B^2\left(1+\frac{r^2}{A^2}\right), \qquad e^{-\lambda}= \frac{\left(1-r^2/C^2\right)\left(1+r^2/A^2\right)}{1+2r^2/A^2},$$00

produces regular central metric potentials and causal radial sound speed, but in the charged case the charge density is singular at the center (Malaver, 2015).

Tolman IV has also been used outside fluid-structure acceptability in order to probe local interaction physics. For curvature-corrected Yukawa exchange in compact stars, the local orthonormal Ricci components satisfy $$e^{\nu}=B^2\left(1+\frac{r^2}{A^2}\right), \qquad e^{-\lambda}= \frac{\left(1-r^2/C^2\right)\left(1+r^2/A^2\right)}{1+2r^2/A^2},$$01, so the correction remains radially symmetric in the local inertial frame; the reported energy shifts for realistic stars are extremely small, of order $$e^{\nu}=B^2\left(1+\frac{r^2}{A^2}\right), \qquad e^{-\lambda}= \frac{\left(1-r^2/C^2\right)\left(1+r^2/A^2\right)}{1+2r^2/A^2},$$02 for solution IV (Zamperlini et al., 26 Mar 2026).

The name should be distinguished from the broader family of Tolman metrics and from Tolman–Oppenheimer–Volkoff equations. An $$e^{\nu}=B^2\left(1+\frac{r^2}{A^2}\right), \qquad e^{-\lambda}= \frac{\left(1-r^2/C^2\right)\left(1+r^2/A^2\right)}{1+2r^2/A^2},$$03-gravity kinetic-theory paper derives a TOV-like relation for stationary Vlasov matter,

$$e^{\nu}=B^2\left(1+\frac{r^2}{A^2}\right), \qquad e^{-\lambda}= \frac{\left(1-r^2/C^2\right)\left(1+r^2/A^2\right)}{1+2r^2/A^2},$$04

but it does not discuss Tolman IV specifically (Jain et al., 2015). Within the contemporary literature, “Tolman IV solution” denotes the particular regular interior geometry that continues to function as an analytic benchmark for equilibrium, acceptability, deformation, and compact-star phenomenology.