Alphabeta Math
Pipeline-generated
How statement and proof provenance work

The first chip identifies the source of the statement or construction; the second identifies the source of its local proof or verification.

  • Literature-sourced: the exact statement appears in a cited source; only wording and notation differ.
  • AI-adapted: a semantically identical restatement of literature-sourced material, modulo indexing, notation, and boundary cases adopted by the library.
  • AI-generated: a genuinely novel statement formulated by AI, with no source for the claim itself.

These labels describe origin, not correctness: citations and verification chips remain separate evidence.

Smooth static spherical geometry and state domain

1 · Prerequisites

2 · Summary

This is a mathematical page about a specified smooth, static, spherically symmetric metric and its state profiles. The exact product manifold, areal-radius form, signature, lapse, diagonal isotropic source tensor, mass function and formal SI dimensions are defined in Smooth static spherical geometry and state domain. The countable-choice convention is only the ACω inherited from the cited smooth Lorentzian and tensor-bundle suppliers; no physical postulate or empirical claim is made.

The reduction lemma Complete spherical Einstein and conservation reduction computes the signed Einstein components from the metric and proves both directions between the full smooth component equations and the mass, lapse and pressure equations on r>0. Its converse uses the contracted Bianchi identity only after the mass, lapse and conservation equations hold. It does not assert an equation for anisotropic, rotating, nonstatic, weak, or distributional sources.

The center result Even radial coefficients give a C² Cartesian center supplies the explicit even-radial-to-Cartesian C2 calculation. The local theorem Regular-center TOV local existence and dependence constructs a positive-lapse, positive-b center by a Volterra integral contraction, gives uniform local dependence on central pressure, and permits continuation only while the positive-b EOS domain is compactly retained. The exact constant-density result Exact regular constant-density solution with verified matching checks the pressure, lapse, mass and vacuum match for the family-specific range 0<u<8/9; it proves neither a general Buchdahl theorem nor realistic causal transport.

Finally, Proper spherical entropy and particle constrained functionals defines proper spherical entropy, particle number and mass-energy on a supplied C2 local state chart. The theorem Constrained entropy criticality, Tolman–Klein equivalence and TOV consequence proves the first-variation equivalence with the Tolman/Klein products and derives the TOV pressure equation under its fixed-radius, positive-b, state-chart and boundary-normalization hypotheses. It proves criticality, not a maximum, existence theorem, or dynamical stability.

3 · Logical flowchart

4 · Definitions, postulates, theorems and proofs

DefinitionDefinition: AI-adaptedProof: Not applicablejudge pass (gpt-6-luna)Open item page →

Smooth static spherical geometry and state domain

Definition

Assume the countable-choice convention of The Axiom of Countable Choice (ACω) for the canonical smooth tangent- and tensor-bundle structures in Smooth time-oriented Lorentzian four-manifolds: the geometric objects of the model. Fix positive formal parameters c>0 and G>0 with SI dimensions [c]=m s−1 and [G]=m3 kg−1 s−2. In this mathematical definition they are fixed parameters, not derived quantities, and no numerical or empirical value is asserted. Let I0⊆R be a nonempty open interval for a length coordinate x0, and define the corresponding time coordinate t:=x0/c∈I0/c, measured in seconds. Let r be the coordinate on (0,∞), and let (S2,γ) be the unit round sphere with the induced smooth round metric of The round metric on the sphere as an induced metric. On the product smooth four-manifold M=I0×(0,∞)×S2, using the canonical product smooth structure of Products of smooth manifolds have a canonical product smooth structure, define, for ϕ,b∈C∞((0,∞),R) with b(r)>0, g=−e2ϕ(r)(dx0)2+b(r)−1dr2+r2γ. Set N=∂x0 as time orientation. The signature convention is (−,+,+,+), x0 and r are measured in metres, the line element has units m2, and ϕ and b are dimensionless. The sphere at fixed (x0,r) has induced metric r2γ. In spherical coordinates its area element is r2sin⁡θ dθ dφ by The round metric on the sphere as an induced metric, so ∫02π∫0πr2sin⁡θ dθ dφ=4πr2; hence r is the areal radius. The flow of N is translation in x0 and preserves g because its coefficients are independent of x0; here a Killing field means a vector field whose local flow preserves the metric. Moreover, g(N,⋅)=−e2ϕdx0, so N is normal to the hypersurfaces x0=constant. In this definition, the metric is static in the precise sense that it has a timelike Killing field normal to such hypersurfaces; N is timelike because g(N,N)=−e2ϕ<0.

A static isotropic diagonal source field is a smooth tensor field of type (1,1) in the sense of Tensor fields of type (r,s) on a Lorentzian four-manifold and their components, defined globally using the product splitting TM≅TI0⊕T(0,∞)⊕TS2 by T(v0,vr,vS)=(−ϵ(r)v0,p(r)vr,p(r)vS). In every standard spherical chart its mixed components are Tab=diag⁡(−ϵ(r),p(r),p(r),p(r)), where ϵ,p∈C∞((0,∞),R) represent energy density and isotropic pressure, respectively, under the formal physical reading, with [ϵ]=[p]=J m−3. The associated covariant tensor has components Tab:=gacTcb; it is smooth and symmetric because both component matrices are diagonal, as defined for covariant 2-tensor fields in Smooth covariant 2-tensor fields on a smooth manifold. Define the mass function M:(0,∞)→R by M(r):=c2r2G(1−b(r)), which has units kg. No Einstein equation, equation of state, central regularity, physical applicability, or empirical claim is part of this definition. The Einstein tensor and component conventions used later are those of The Riemann curvature operator and the derived Ricci, scalar and Einstein tensors and Tensor fields of type (r,s) on a Lorentzian four-manifold and their components.

Remarks

The coordinate domain excludes r=0; a regular center is a separate conclusion and is defined by the center regularity result and the local-existence theorem on this page. Countable Choice is inherited only through the stated canonical smooth-bundle suppliers; the coordinate formulas below add no stronger choice principle.

LemmaStatement: AI-adaptedProof: AI-generatedjudge pass (gpt-6-luna)Open item page →

Even radial coefficients give a C² Cartesian center

Statement

Let δ>0, let f,h∈C2((−δ,δ),R) be even, and let Bδ(0)⊆R3 carry its standard Euclidean coordinates x=(x1,x2,x3) and radius r=∣x∣ from The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn. Then F(x)=f(∣x∣) is C2 on Bδ(0), and each component Hij(x)=δij+xixjh(∣x∣) is C2 there. In particular, if an areal spherical spatial metric has b(r)=1−r2q(r) with even q∈C2((−δ,δ)) and b>0, then its Cartesian coefficients gij(x)=δij+xixjq(∣x∣)b(∣x∣) extend C2 across the center. This result is a coordinate regularity statement; it asserts no field equation or physical claim.

Facts & Assumptions

Given: δ>0, even C2 functions f,h,q on (−δ,δ), and the Euclidean norm on R3.

[F2]

If g∈C2((−δ,δ)) is even, then g′(0)=0, the derivative definition gives g′(r)/r→g′′(0) as r→0, and Peano's Taylor theorem gives g(r)=g(0)+12g′′(0)r2+o(r2) (Peano's form: the normalized Taylor remainder tends to zero).

[F3]

In Euclidean polar coordinates on R3∖{0}, the Euclidean metric is dr2+r2γ with γ the unit round sphere metric (The round metric on the sphere as an induced metric).

Proof

1.1F1F2algebra

For r>0, the chain rule and product rule in [F1] give ∂iF=f′(r)xir,∂j∂iF=(f′′(r)−f′(r)r)xixjr2+f′(r)rδij. Evenness implies f′(0)=0. By [F2], f′(r)/r→f′′(0) and f(r)=f(0)+12f′′(0)r2+o(r2) as r→0. Thus F is differentiable at 0 with gradient 0 and Hessian f′′(0)I3. Since the displayed first and second derivatives tend to those values as x→0, they extend continuously to 0, proving F∈C2.

1.2F1F2algebra

For r>0, ∂kHij=δikxjh(r)+δjkxih(r)+xixjh′(r)xkr. Evenness gives h′(0)=0, so [F2] yields h′(r)/r→h′′(0) and h′′(r)−h′(r)/r→0. Differentiating the displayed formula once more, the terms without a factor xi or xj are (δikδjℓ+δjkδiℓ)h(r); every remaining term is a sum of terms containing either one factor xi or xj times h′(r)xℓ/r, or two such factors times h′′(r)−h′(r)/r or h′(r)/r. Since h′ is O(r) and h′′, h′/r are bounded near 0, these remaining terms tend to 0. The Peano expansion in [F2] gives Hij(x)=δij+h(0)xixj+o(r2), hence derivative 0 and second derivative (δikδjℓ+δiℓδjk)h(0) at 0. Thus the first and second derivatives extend continuously and Hij is C2.

2.1F1F3step 1.1step 1.2algebra∎

For r>0, write dr2=(xi dxi)2/r2 and use the polar identity dr2+r2γ=∑idxi2 from [F3]. Therefore b−1dr2+r2γ=∑idxi2+b−1−1r2(xi dxi)2=∑i,j(δij+xixjq(r)b(r))dxi dxj, because b=1−r2q. The quotient q/b is even C2 since b>0, so step 1.2 proves the spatial extension. The time coefficient −e2ϕ(r) extends C2 by step 1.1 whenever ϕ is even C2. This proves the metric clause.

LemmaStatement: AI-adaptedProof: AI-generatedjudge pass (gpt-6-luna)Open item page →

Complete spherical Einstein and conservation reduction

Statement

Assume ZF+ACω with ACω exactly as in The Axiom of Countable Choice (ACω). Let J⊆(0,∞) be an open interval; let the state domain, metric, mixed source tensor, mass function, units, and signature be those of Smooth static spherical geometry and state domain. Write κ=8πG/c4, and let Gab be the mixed Einstein tensor defined from the signed curvature convention in The Riemann curvature operator and the derived Ricci, scalar and Einstein tensors. Suppose ϕ,b,ϵ,p are smooth on J and b>0. As a conditional mathematical system, assume the component equations Gab=κTab. Then on J, M′=4πr2ϵc2,ϕ′=GM/c2+4πGr3p/c4r2b,p′=−G(ϵ+p)(M+4πr3p/c2)c2r2b. Conversely, if these three ordinary differential equations hold with M=c2r(1−b)/(2G), then every mixed Einstein component satisfies Gab=κTab. The identity ∇aGab=0 used for this converse is the geometric theorem of The second Bianchi identity and the vanishing divergence of the Einstein tensor, not a field equation or empirical premise. This lemma is a pure conditional calculation and makes no claim that the model describes nature.

Facts & Assumptions

Given: ZF+ACω as declared in the Statement; an open interval J⊂(0,∞); smooth ϕ,b,ϵ,p with b>0; the static spherical data in Smooth static spherical geometry and state domain; and, for the forward direction, the conditional equations Gab=κTab.

[F1]

The static spherical metric, diagonal mixed source, signature (−,+,+,+), mass definition M=c2r(1−b)/(2G), and formal SI dimensions are those in Smooth static spherical geometry and state domain.

[F2]

With curvature convention R(X,Y)Z=∇X∇YZ−∇Y∇XZ−∇[X,Y]Z, the Christoffel symbols are Γabc=12gad(∂bgdc+∂cgbd−∂dgbc) (Existence and uniqueness of the metric-selected connection of a nondegenerate metric, and its extension to tensor fields), the coordinate curvature is Rabcd=∂cΓadb−∂dΓacb+ΓaceΓedb−ΓadeΓecb (Tensoriality of the curvature, its coordinate formula and its point-normal metric formula), and Gab=Ric⁡ab−12Sgab with Ric⁡bd=Rabad (The Riemann curvature operator and the derived Ricci, scalar and Einstein tensors).

[F3]

The contracted Bianchi identity is ∇aGab=0 for every index b (The second Bianchi identity and the vanishing divergence of the Einstein tensor).

[A1]

The field equation Gab=κTab is a mathematical hypothesis for the forward implication only; it is not adopted or proved as a claim about nature.

[A2]

The smooth tangent- and tensor-bundle construction in the geometric supplier chain assumes ACω as recorded in The Axiom of Countable Choice (ACω); no stronger choice axiom is assumed.

Proof

1.1A2F1F2

In a regular spherical chart away from the angular coordinate poles, the metric in [F1] has the following nonzero Christoffel symbols, with lower-index symmetry understood: Γ00r=ϕ′,Γr00=be2ϕϕ′,Γrrr=−b′2b,Γrθθ=−rb,Γrφφ=−rbsin⁡2θ,Γθrθ=Γφrφ=1r,Γθφφ=−sin⁡θcos⁡θ,Γφθφ=cot⁡θ. They follow by inserting the diagonal metric entries (−e2ϕ,b−1,r2,r2sin⁡2θ) and their derivatives into the Christoffel formula in [F2]. The identities obtained below depend only on r and are tensor-component identities, so their validity extends across the two angular coordinate poles by smoothness.

2.1A2step 1.1F2algebra

Substituting the symbols of step 1.1 in the curvature formula of [F2], contracting as in [F2], and simplifying gives Ric⁡00=−b(ϕ′′+ϕ′2+2ϕ′r+b′ϕ′2b),Ric⁡rr=−b(ϕ′′+ϕ′2+b′ϕ′2b)−b′r,Ric⁡θθ=Ric⁡φφ=1−br2−b′2r−bϕ′r. Hence G00=−(1−b)/r2+b′/r, Grr=−(1−b)/r2+2bϕ′/r, and Gθθ=Gφφ=b(ϕ′′+ϕ′2+ϕ′r)+b′2(ϕ′+1r). These equalities use only the displayed finite component formula, the listed Christoffels, and the contractions defining G; no spherical field equation is cited.

2.2A2F1step 1.1algebra

For the diagonal tensor Tab of [F1], its radial divergence is ∇aTar=p′+ϕ′(ϵ+p). Indeed, in ∂aTar+ΓaacTcr−ΓcarTac, the derivative term is p′, the trace term contributes p(ϕ′−b′/(2b)+2/r), and the last term contributes ϵϕ′+pb′/(2b)−2p/r. Adding the three contributions gives p′+(ϵ+p)ϕ′. The other three components of ∇aTab vanish: for b=0 there is no time dependence and every mixed time-radial Christoffel is zero; for b=θ,φ, angular derivatives vanish and the only polar trace term cot⁡θ p cancels the matching diagonal connection term, while the φ component has no nonzero trace term.

3.1A1A2F1step 2.1algebra

The 00 equation in [A1] is −(1−b)/r2+b′/r=−κϵ. Multiplying by r2 yields 1−b−rb′=κϵr2. Differentiating M=c2r(1−b)/(2G) gives M′=c2(1−b−rb′)/(2G)=4πr2ϵ/c2.

3.2A1A2F1step 2.1algebra

The rr equation in [A1] is −(1−b)/r2+2bϕ′/r=κp. Solving for ϕ′ and substituting 1−b=2GM/(c2r) and κ=8πG/c4 gives ϕ′=(1−b)+κpr22br=GM/c2+4πGr3p/c4r2b.

4.1A1A2F3step 3.2step 2.2algebra

By [F3] and [A1], ∇aTar=κ−1∇aGar=0. Step 2.2 therefore gives p′=−(ϵ+p)ϕ′. Substitution of step 3.2 yields the displayed TOV pressure equation.

5.1A2F1F3step 1.1step 2.1step 2.2step 3.1step 3.2algebra∎

Conversely, the mass and lapse equations give the 00 and rr equations by reversing steps 3.1–2.2. Define Dab=Gab−κTab. Then D00=Drr=0, and its only possible nonzero entries are the equal angular diagonal entries Dθθ=Dφφ, by step 2.1. The pressure equation and step 2.2 give ∇aTab=0 for every b, so [F3] gives ∇aDab=0. In particular, the radial component is ∇aDar=−ΓθθrDθθ−ΓφφrDφφ=−2Dθθ/r. Since r>0, it follows that both angular entries vanish. Thus D=0 and the full mixed field equation holds.

Remarks

  • The claim is conditional mathematics. Its ACω assumption is inherited exactly from the smooth tangent/tensor-bundle and curvature suppliers; no full Axiom of Choice is invoked, and the coordinate algebra performs no selection.
  • The direct component computation verifies all independent Einstein components on r>0. The converse uses Bianchi only after the mass, lapse, and conservation equations hold.
  • The equation is confined to a smooth isotropic static spherical domain with b>0; it says nothing about anisotropic, rotating, nonstatic, weak, or distributional sources.
DefinitionDefinition: AI-adaptedProof: Not applicableOpen item page →

Proper spherical entropy and particle constrained functionals

Definition

Assume ZF+ACω as in The Axiom of Countable Choice (ACω). Fix R>0, an open state chart U⊂(0,∞)2, and a function s∈C2(U,R) with sϵ>0. The coordinates on U are energy density ϵ in J m−3 and particle number density n in m−3; the entropy density has units [s]=J K−1 m−3. Define the local temperature, chemical potential, and pressure functions by Θ(ϵ,n):=1sϵ(ϵ,n),μ(ϵ,n):=−Θ(ϵ,n)sn(ϵ,n),p(ϵ,n):=Θs+μn−ϵ. These are C1 functions on U; their formal SI dimensions are K, J, and J m−3, respectively. By definition ϵ+p=Θs+μn, and differentiation gives dϵ=Θ ds+μ dn and dp=s dΘ+n dμ on the chart.

Let ϵ,n be restrictions to [0,R] of even C2 functions on an open interval containing [−R,R], positive on [0,R], with image in a compact subset of U. This evenness is the center regularity convention. With K:=4πG/c4, define M(r):=4πc2∫0rϵ(t)t2 dt,b(r):=1−2GM(r)c2r(r>0),b(0):=1. Equivalently, b(r)=1−r2q(r) on [−R,R], where q(r):=2K∫01s2ϵ(rs) ds. The compact-parameter differentiation theorem Leibniz's rule on a compact rectangle: an interior parameter derivative with a continuous extension may be passed through a Riemann integral makes q even C2. Therefore Even radial coefficients give a C² Cartesian center gives the corresponding C2 Cartesian spatial metric coefficients. Assume b(r)≥b∗>0 on [0,R]. The displayed integrals are Riemann integrals in the Darboux convention of The lower and upper Darboux integrals of a bounded f on [a,b] as sup⁡PL(f,P) and inf⁡PU(f,P), Darboux integrability as their equality, and the notation ∫abf; their continuous integrands are integrable by A continuous function on [a,b] is Riemann integrable, by Heine-Cantor and Riemann's criterion, and b−1/2 is continuous on this positive range by Continuity and derivatives of positive-base real powers. Use the spatial metric h=b(r)−1dr2+r2γ on (0,R)×S2, where γ is the unit round sphere metric used in Smooth static spherical geometry and state domain. Here h has C2 coefficients and the Cartesian center extension just established; no globally smooth static state is required. Define its volume density in a spherical coordinate chart by det⁡(hij) dr dθ dφ. Since det⁡(hij)=b−1r4sin⁡2θ, integration over the unit sphere gives the radial density 4πr2b(r)−1/2dr. This uses the same spatial metric formula as the smooth static domain, at the present C2 regularity. Define the proper entropy, proper particle number, and spherical mass-energy functionals S[ϵ,n]:=4π∫0Rs(ϵ(r),n(r)) r2b(r)−1/2 dr,Nb[ϵ,n]:=4π∫0Rn(r) r2b(r)−1/2 dr,E[ϵ,n]:=c2M(R)=4π∫0Rϵ(r)r2 dr. Their formal dimensions are [S]=J K−1, [Nb]=1, and [E]=J. These are mathematical functionals on the supplied state chart and positive-b profile domain; they assert no entropy law, equilibrium, maximum, dynamics, or empirical fact.

Remarks

The energy functional is the spherical mass-energy constraint c2M(R), not a generic additive gravitational energy on arbitrary spacetimes. No relation between this state chart and any particular material equation of state is assumed beyond the displayed definitions.

TheoremStatement: AI-adaptedProof: AI-generatedjudge pass (gpt-6-luna)Open item page →

Regular-center TOV local existence and dependence

Statement

Assume ZF+ACω, where ACω is exactly The Axiom of Countable Choice (ACω) and is inherited from the smooth Lorentzian/tensor suppliers. Let ϵ:Ip→(0,∞) be C∞ on an open interval Ip⊆(0,∞), let Pc⊂Ip be a nonempty compact interval, and choose δ>0 such that the closed δ-neighbourhood of Pc lies in Ip. Fix ϕc∈R and use the units and signature of Smooth static spherical geometry and state domain. Put K=4πG/c4, with formal units m J−1; v has units m−2. There is R0>0, common to all pc∈Pc, such that for each pc there is a unique pair of even continuous functions p,v on [−R0,R0] with ∣p(r)−pc∣≤δ and satisfying v(r)=K∫01s2ϵ(p(rs)) ds,p(r)=pc−∫0rt(ϵ(p(t))+p(t))(v(t)+Kp(t))1−2t2v(t) dt. They are C∞, p(0)=pc, v(0)=Kϵ(pc)/3, p′(0)=v′(0)=0, and b(r):=1−2r2v(r)≥1/2 on [0,R0]. Define M(r)=c2Gr3v(r),ϕ(r)=ϕc+∫0rtv(t)+Kp(t)1−2t2v(t) dt(∣r∣≤R0). Then (M,p,ϕ) satisfies the mass, pressure, and lapse equations of Complete spherical Einstein and conservation reduction for 0<r≤R0. The metric extends across r=0 as a C2 Lorentzian metric in Cartesian spatial coordinates, with b(0)=1 and lapse eϕc>0; this is the regular-center conclusion in this item. On the common interval, p,v,ϕ depend continuously in the uniform norm on pc∈Pc. For r>0, the solution continues past any finite endpoint at which its graph remains in a compact subset of {(r,p,v):r>0, p∈Ip, 1−2r2v>0}. No assertion is made that every equation of state reaches a zero-pressure surface, avoids a horizon globally, or gives a dynamically stable solution.

Facts & Assumptions

Given: The state domain and SI conventions of Smooth static spherical geometry and state domain, ACω of The Axiom of Countable Choice (ACω), the smooth positive function ϵ on Ip, the compact central-pressure interval Pc, the margin δ, and ϕc.

[F1]

For a smooth isotropic spherical profile, the mass and lapse equations and the pressure equation are exactly those in Complete spherical Einstein and conservation reduction. With GM/c2=r3v they become v′=(Kϵ(p)−3v)/r, p′=−r(ϵ(p)+p)(v+Kp)/(1−2r2v), and ϕ′=r(v+Kp)/(1−2r2v) for r>0.

[F2]

If even C2 functions f,h are used in f(∣x∣) and δij+xixjh(∣x∣), these extend C2 to the Cartesian center; for b=1−r2q the spatial metric coefficients are δij+xixjq/b (Even radial coefficients give a C² Cartesian center).

[F3]

A continuous function on a compact interval is Riemann integrable, and the first fundamental theorem differentiates its variable-endpoint integral at points of continuity (A continuous function on [a,b] is Riemann integrable, by Heine-Cantor and Riemann's criterion, The first fundamental theorem: if f is integrable on [a,b] and continuous at c, then F′(c)=f(c); in particular a continuous f has F as a primitive). A continuous parameter integrand with continuous parameter derivative on a compact rectangle can be differentiated through its Riemann integral (Leibniz's rule on a compact rectangle: an interior parameter derivative with a continuous extension may be passed through a Riemann integral).

[F4]

If a locally Lipschitz ODE solution remains in a compact subset of its open time-state domain near a finite endpoint, it extends beyond that endpoint (Compact continuation for locally Lipschitz ODEs: a solution whose graph stays in a compact subset of the domain near a finite endpoint extends past it clause 1).

[F6]

Closed bounded subsets of Rn are compact, and every continuous real-valued function on a nonempty compact subset attains finite extrema (For a nonempty subset of Rn with n≥1, compactness, closedness and boundedness, pseudocompactness, and attainment of extrema by every continuous real-valued function are equivalent).

[F8]

The second FTC evaluates an integral of a continuous derivative; in particular ∫01s2 ds=1/3 and ∫0rt dt=r2/2 (The second fundamental theorem: if G is differentiable on [a,b] with G′=f and f is integrable, then ∫abf=G(b)−G(a)).

[F9]

Integral monotonicity and nonnegativity hold for integrable functions on an ordered compact interval; in particular, if g≥0 and −g≤f≤g, then ∣∫f∣≤∫g by comparing with g and −g (If f≤g on [a,b] and both are integrable then ∫abf≤∫abg; and m(b−a)≤∫abf≤M(b−a)). Linearity permits subtraction and scalar multiplication under the integral (Integrable functions on [a,b] form a set closed under sums and scalar multiples, and ∫ab(λf+μg)=λ∫abf+μ∫abg). Reversing endpoints changes only the sign, and an equal-endpoint integral is zero (The integral with oriented limits: ∫aaf:=0 and ∫baf:=−∫abf), so absolute estimates use the interval with increasingly ordered endpoints.

[F10]

A continuous real function on a compact real interval is uniformly continuous (Heine-Cantor in R: a continuous real function on a compact subset of R is uniformly continuous).

[A1]

The countable-choice use is the canonical smooth tangent/tensor-bundle structure inherited through The Axiom of Countable Choice (ACω); no full Axiom of Choice is used. The integral contraction and its iteration below make no selection from an infinite family.

Proof

technique · Choose one common closed pressure interval and uniform contraction bounds for every $p_c\in P_c$; then bootstrap the fixed point to the stated center regularity
1.1givenF1F5F6F8choosealgebra

Write Pc=[p−,p+] and set P=[p−−δ,p++δ]; the margin hypothesis puts P in Ip, and [F6] makes it compact. Thus ϵmax⁡=max⁡Pϵ and Lϵ=max⁡P∣ϵ′∣ exist and are finite. Put V=Kϵmax⁡/3, and first choose R1≤1 with 2R12V≤1/2. The set Q=[−R1,R1]×P×[0,V] is closed and bounded in R3, hence compact by [F6]. On Q, the denominator of H(r,x,y)=(ϵ(x)+x)(y+Kx)1−2r2y is at least 1/2, so H,Hx,Hy are continuous there by [F5]; [F6] supplies finite bounds C for ∣H∣, Cx for ∣Hx∣, and Cy for ∣Hy∣. Set C1=Cx+KLϵCy/3. Choose R0≤R1 so that CR02/2≤δ and C1R02/2≤1/2.

2.1step 1.1givenF3F5F7F8F9F10algebra

For each pc∈Pc, let Bpc be the closed sup-norm ball of radius δ about the constant function pc in C([−R0,R0]). It is complete: a sup-norm Cauchy sequence has pointwise limits in the closed real interval P, the convergence is uniform by the Cauchy bound, and its uniform limit is continuous. For q∈Bpc define v[q](r)=K∫01s2ϵ(q(rs)) ds,(Tpcq)(r)=pc−∫0rtH(t,q(t),v[q](t)) dt. The function ϵ∘q is uniformly continuous on [−R0,R0] by [F10]. Because ∣rs−r~s∣≤∣r−r~∣ for s∈[0,1], [F9] bounds ∣v[q](r)−v[q](r~)∣ by K/3 times sup⁡s∈[0,1]∣ϵ(q(rs))−ϵ(q(r~s))∣, which tends to zero; hence v[q] is continuous. The integrand defining Tpcq is then continuous, so [F3] makes that map continuous as well. By [F9] and [F8], 0≤v[q]≤Kϵmax⁡/3=V; hence 1−2r2v[q]≥1/2. The mean value theorem [F7], linearity and [F9] applied to the bounds ∣s2(ϵ(q(rs))−ϵ(q~(rs)))∣≤s2Lϵ∥q−q~∥∞ give ∥v[q]−v[q~]∥∞≤KLϵ∥q−q~∥∞/3. Applying [F7] successively in the x and y variables gives ∣H(r,x,y)−H(r,x~,y~)∣≤Cx∣x−x~∣+Cy∣y−y~∣ for state pairs in P×[0,V]. For each endpoint r, [F9] bounds the absolute integral of tH or of the difference t(H−H~) by the corresponding constant times the integral of ∣t∣ over the interval between 0 and r. This last integral is r2/2 by [F8]; reversing endpoints changes only the integral sign, and the r=0 case is zero. Hence the bounds in step 1.1, [F9] and [F8] give ∥Tpcq−pc∥∞≤CR02/2≤δ and ∥Tpcq−Tpcq~∥∞≤(C1R02/2)∥q−q~∥∞≤12∥q−q~∥∞. Thus Tpc maps the ball to itself and is a contraction.

3.1step 2.1constructalgebra

Iterate canonically from q0(r)=pc, qj+1=Tpcqj. The differences satisfy ∥qj+1−qj∥∞≤2−j∥q1−q0∥∞, so the geometric tail estimate makes (qj) Cauchy in the complete ball. Its uniform limit p is continuous and remains in the ball. The contraction estimate makes Tpc continuous, hence Tpcp=p; any two fixed points have distance at most half their distance and therefore coincide. If q is even, then v[q] is even, H(r,q(r),v[q](r)) is even, and its product with r is odd, so Tpcq is even. Every iterate and hence p is even; the fixed point equations give the two displayed integral identities.

4.1step 3.1F1F3F5algebra

The first FTC in [F3] makes p continuously differentiable and yields p′(r)=−rH(r,p(r),v[p](r)), including p′(0)=0. The parameter-integral theorem [F3], together with the Ck composition and algebra rules [F5], applied repeatedly on [−R0,R0]×[0,1], shows that whenever p∈Ck, v[p] is Ck because (r,s)↦s2ϵ(p(rs)) has continuous r-derivatives through order k. The denominator is bounded below by 1/2, so H(r,p(r),v[p](r)) is Ck; the integral equation then gives p∈Ck+1. Induction proves p,v are smooth and even. Rewriting the v integral as r3v(r)=K∫0rt2ϵ(p(t)) dt, then applying the product rule and the first FTC gives v′=(Kϵ(p)−3v)/r for r>0; the pressure equation is the defining derivative of p.

4.2A1step 1.1step 3.1F1F5F6F7F8F9algebra

The same R0 and contraction factor work for every central value in Pc. Subtracting the fixed-point equations for pa,pb gives ∥pa−pb∥∞≤∣a−b∣+12∥pa−pb∥∞, hence ∥pa−pb∥∞≤2∣a−b∣. The integral formula for v, [F9] and [F8] give ∥va−vb∥∞≤KLϵ∥pa−pb∥∞/3. The lapse integrand J(r,x,y)=r(y+Kx)/(1−2r2y) and its x,y derivatives are continuous and bounded on the same compact set Q by [F6]. Applying [F7] successively in x,y bounds its state increment by a constant times ∣pa−pb∣+∣va−vb∣; integrating over 0≤r≤R0 with [F9] yields ∥ϕa−ϕb∥∞≤CJR0(∥pa−pb∥∞+∥va−vb∥∞) for a finite CJ. Thus the solution map is continuous in the sup norm. If another regular-center solution stays in the same pressure ball, integrating its mass and pressure equations yields the same two fixed-point equations; uniqueness in step 3.1 proves it is this solution.

5.1A1step 4.1F1F2F5algebra

Set b=1−2r2v and define ϕ by the displayed integral. Then b≥1/2, ϕ is smooth and even, and ϕ′(0)=0. The product rule and the equation for v give M′=4πr2ϵ(p)/c2 from M=c2r3v/G; also GM/c2=r3v, so ϕ′=r(v+Kp)/b is the lapse equation in [F1]. Apply [F2] with q=2v and h=q/b to obtain the C2 Cartesian metric extension. Its center coefficients are g00(0)=−e2ϕc and gij(0)=δij, so it is Lorentzian at the center.

6.1A1step 4.1F1F4F5F6∎

For r>0, let U={(r,p,v):r>0, p∈Ip, 1−2r2v>0} and define the vector field F(r,p,v)=(−r(ϵ(p)+p)(v+Kp)1−2r2v, Kϵ(p)−3vr). By [F5], F is C1 on the open set U. Around each point of U choose a closed time-state cylinder [r−,r+]×B‾((p∗,v∗),η)⊂U. It is a closed bounded subset of R3, hence compact by [F6]. The state derivative of F is continuous there, and its operator norm is continuous, so [F6] gives a finite uniform bound L. For each fixed r∈[r−,r+], the state ball is convex and open; the mean-value inequality [thm-mean-value-inequality-for-total-derivatives] gives ∥F(r,z)−F(r,z~)∥2≤L∥z−z~∥2 on it. Thus F is locally Lipschitz in the state variable, locally uniformly in time, as required by [F4]. If the solution graph remains in a compact subset of U near a finite endpoint, [F4] extends (p,v) past that endpoint. The integral defining ϕ then also extends. This is only the stated compact-interior continuation criterion; it gives no global surface or stability claim.

Remarks

The integral contraction resolves the coordinate singularity at r=0 directly. The ordinary ODE continuation theorem is used only on r>0, where the vector field is smooth. The proof is conditional mathematics under the stated EOS and ACω convention; it does not assert that an observed star follows this model.

TheoremStatement: AI-adaptedProof: AI-generatedOpen item page →

Exact regular constant-density solution with verified matching

Statement

Assume ZF+ACω exactly as in The Axiom of Countable Choice (ACω) and use the geometry and units of Smooth static spherical geometry and state domain. Let ρ0>0 be a constant mass density in kg m−3, R>0 in metres, ϵ0=ρ0c2 in J m−3, MR=4πρ0R3/3 in kilograms, and u=2GMRc2R∈(0,8/9),a=1−u,B(r)=1−ur2/R2. For 0≤r≤R, define M(r)=MR(r/R)3,b(r)=B(r)2,p(r)=ϵ0B(r)−a3a−B(r),eϕ(r)=3a−B(r)2. Set ϵ(r)=ϵ0 for 0≤r<R and ϵ(r)=0 for r≥R. For r≥R, also set M(r)=MR, b(r)=1−2GMR/(c2r), p(r)=0, and eϕ(r)=b(r). Then p(R)=0, b>0, the interior lapse is positive, and the central pressure is finite and equals ϵ0(1−a)/(3a−1). The interior and exterior fields solve the spherical Einstein and conservation system of Complete spherical Einstein and conservation reduction on their open pieces; the interior equations extend to the center by continuity. The interior metric has a C2 Cartesian center; the interior and exterior join at r=R as a C1 metric that is piecewise C2, with matching first derivatives at the interface. For this interface assertion, a radial distribution on an open radial-coordinate interval U is a linear functional T on Cc∞(U), the space of smooth real functions with compact support in U, such that for every compact K⊂U there exist CK≥0 and an integer NK≥0 satisfying ∣T[ψ]∣≤CKmax⁡0≤j≤NKsup⁡K∣ψ(j)∣ whenever supp⁡ψ⊆K, the Dirac distribution is δR[ψ]=ψ(R), and a delta-supported surface stress means a source term SabδR. The prescribed piecewise source has no such term. The proof also shows that 0<u<8/9 is necessary and sufficient for these algebraic formulas to give real interior data with b>0, finite nonnegative pressure, and positive lapse throughout [0,R]. This range is specific to the constant-density family and is not a general Buchdahl theorem.

Facts & Assumptions

Given: ρ0>0, R>0, u∈(0,8/9), ϵ0=ρ0c2, the metric convention of Smooth static spherical geometry and state domain, ACω of The Axiom of Countable Choice (ACω), and the reduction and radial regularity results below.

[F1]

The metric signature, source components, curvature sign, and SI units are those of Smooth static spherical geometry and state domain; the spherical Einstein/conservation system is equivalent on r>0 to the mass, lapse, and pressure equations in Complete spherical Einstein and conservation reduction.

[F2]

For b=1−r2q with even C2 q and b>0, the spatial metric has the Cartesian form δij+xixjq/b, and even radial lapse coefficients extend C2 (Even radial coefficients give a C² Cartesian center).

[F3]

The real-power derivative formula gives (x)′=1/(2x) on (0,∞); repeated differentiation gives smoothness of positive-base powers (Continuity and derivatives of positive-base real powers).

[F6]

A smooth real function with nonzero derivative at a point has a smooth local inverse there (Choice-free smooth inverse function theorem in Euclidean space).

[A1]

For the branch-range claim, for any real u>0 consider the same algebraic formulas wherever their square roots are real; a regular positive branch means that the interior has b>0, finite p≥0, and positive lapse at every radius in [0,R]. No empirical stellar model or causal transport law is assumed.

[A2]

The geometric suppliers use ACω as specified in The Axiom of Countable Choice (ACω); no full Axiom of Choice is used.

Proof

1.1givenA1A2F3algebra

Since 0<u<8/9, a=1−u>1/3. On 0≤r≤R, a≤B(r)≤1, so b=B2>0 and 3a−B(r)≥3a−1>0. Thus the lapse and pressure are finite and the lapse is positive; p(R)=0, p(r)≥0, and p(0)=ϵ0(1−a)/(3a−1).

2.1givenA1step 1.1algebra

To check sharpness, first take 8/9≤u<1. At u=8/9, a=1/3 and B(0)=1=3a, so the pressure denominator vanishes at the center while its numerator is positive, and eϕ(0)=0. If 8/9<u<1, then a<1/3 and a<3a<1=B(0), whereas B(R)=a; continuity gives r∗∈(0,R) with B(r∗)=3a, where the pressure denominator vanishes and the lapse is zero. If u=1, then b(R)=0; if u>1, B(R) is not real. Thus no u≥8/9 gives a regular positive branch, while step 1.1 proves that every 0<u<8/9 does.

2.2A2step 1.1F2F3algebra

The mass formula gives M′=3MRr2/R3=4πρ0r2=4πr2ϵ0/c2. The functions B,ϕ,p are even and smooth near r=0 because 1−ur2/R2>0 and 3a−B>0. Writing b=1−r2q with q=u/R2, [F2] gives a C2 Cartesian center for the interior metric. The Einstein equations hold at the center by continuity of the curvature and source components from r>0.

2.3A2step 1.1F1F3F4algebra

For 0<r<R, the chain, product, and quotient rules [F4], together with [F3], give B′=−urR2B,ϕ′=urR2B(3a−B),(pϵ0)′=−2aurR2B(3a−B)2. Also 1+p/ϵ0=2a/(3a−B), so p′=−(ϵ0+p)ϕ′. Using u=8πGϵ0R2/(3c4) and p/ϵ0=(B−a)/(3a−B), direct simplification yields ϕ′=4πGr(ϵ0/3+p)c4b=GM/c2+4πGr3p/c4r2b, where M=MRr3/R3. Substituting the displayed lapse into p′=−(ϵ0+p)ϕ′ is the TOV pressure equation in [F1], so the interior mass, lapse, and pressure equations all hold.

2.4A2F1F4step 1.1algebra

For r≥R, M′=0, ϵ=p=0, and b=1−uR/r>0 since u<1. The chain and quotient rules [F4] give ϕ′=12b′/b=uR/(2r2b)=GMR/(c2r2b), so the vacuum mass and lapse equations hold and the pressure equation is 0=0.

3.1A2F1F3F4F5F6step 1.1step 2.2step 2.3step 2.4algebra∎

Let b−(r)=1−ur2/R2 on the interior side and b+(r)=1−uR/r on the exterior side. Both formulas are smooth and positive on separate open neighborhoods of R. By [F3]–[F4], their negative one-half powers are smooth there. For a common value ℓR, define ℓ±(r)=ℓR+∫Rrb±(s)−1/2 ds. By [F5] these functions have derivatives dℓ±/dr=b±−1/2>0; repeated differentiation of these identities makes ℓ± smooth. By [F6], after shrinking the neighborhoods, each has a smooth local inverse r±(ℓ) satisfying dr±/dℓ=b±(r±(ℓ)). Use ℓ<ℓR for the interior and ℓ>ℓR for the exterior as the matching chart, so r±(ℓR)=R and dr±/dℓ(ℓR)=a on both sides. At r=R, [F1] gives b−(R)=b+(R)=a2 and both lapse values equal a. The lapse derivatives with respect to r agree by steps 2.3–2.4: (eϕ)−′(R)=u/(2Ra)=(b+)′(R). Hence in the matching chart the coefficients g00=−e2ϕ and r±(ℓ)2γ have matching values and first derivatives, while gℓℓ=1 on both sides. Thus the metric is C1 and piecewise C2 in that chart. For a piecewise C2 coefficient f(ℓ) with one-sided values and derivatives, define Df[ψ]:=−∫f(ℓ)ψ′(ℓ) dℓ for ψ∈Cc∞, and define δℓR[ψ]=ψ(ℓR) and δℓR′[ψ]=−ψ′(ℓR). These functionals satisfy the stated test bounds: for nonempty compact K⊂U, choose a compact interval [a,b]⊂U containing K and let M bound ∣f∣ there on both sides of the interface. For supp⁡ψ⊆K, ∣∫fψ∣≤(b−a)Msup⁡K∣ψ∣, ∣Df[ψ]∣≤(b−a)Msup⁡K∣ψ′∣, and ∣D2f[ψ]∣=∣∫fψ′′∣≤(b−a)Msup⁡K∣ψ′′∣; hence the orders 0,1,2 suffice. Also ∣δℓR[ψ]∣≤sup⁡K∣ψ∣ and ∣δℓR′[ψ]∣≤sup⁡K∣ψ′∣, with both evaluations zero if ℓR∉K; empty K admits only the zero test function. Linearity follows from these formulas. Integration by parts on the two sides gives D2f=fpiecewise′′+[f′]δℓR+[f]δℓR′, where [f]=f(ℓR+)−f(ℓR−) and similarly for [f′]. The matched metric coefficients and first derivatives make both jumps zero. The Christoffel coefficients are formed from the metric, its inverse, and its first derivatives; since the metric is nondegenerate, they are continuous and piecewise C1. With the curvature convention of [F1], its coordinate formula on the smooth sides is Rabcd=∂cΓadb−∂dΓacb+ΓaceΓedb−ΓadeΓecb; in this interface calculation take those derivatives in the distributional sense just defined. Since the Christoffel coefficients are continuous, the curvature formula produces no surface-supported term. The source fields are piecewise ordinary functions and likewise contain no delta-supported surface stress.

Remarks

The u<8/9 pressure restriction and endpoint divergence belong to this explicit uniform-density branch. The areal radius is used as a smooth coordinate on each side; the matching chart across r=R is the proper-radial coordinate ℓ constructed in the proof. The argument neither proves the general Buchdahl theorem nor supplies a finite-sound-speed matter dynamics for the incompressible equation of state.

TheoremStatement: AI-adaptedProof: AI-generatedOpen item page →

Constrained entropy criticality, Tolman–Klein equivalence and TOV consequence

Statement

Assume ZF+ACω as in The Axiom of Countable Choice (ACω). Fix R>0, a state chart and positive even C2 profiles (ϵ,n) with compact range and b≥b∗>0 as in Proper spherical entropy and particle constrained functionals. Let Θ,μ,p,S,Nb,E be the functions and functionals defined there. Define the lapse by ϕ(r)=12log⁡b(R)−∫rRGM(t)/c2+4πGt3p(t)/c4t2b(t) dt, whose integrand extends continuously to 0 because M(t)=O(t3), so it solves the lapse equation of Complete spherical Einstein and conservation reduction with boundary normalization ϕ(R)=12log⁡b(R). An admissible profile curve is a map λ↦(ϵλ,nλ) on an open interval about 0, C1 in the sup-norm space C0([0,R])2, whose values are even C2 profiles in the stated compactly-contained state chart and positive-b domain, and whose derivative at 0 is a smooth pair compactly supported in (0,R) (extended evenly at the center). Define constrained stationarity to mean that S has zero derivative at the base profile along every such curve on which both E and Nb are constant. Then constrained stationarity is equivalent to the existence of constants α∈R and β∈R such that μ(r)Θ(r)=α,Θ(r)eϕ(r)=ΘK,μ(r)eϕ(r)=μK are constant on (0,R). In that case the profile satisfies p′(r)=−(ϵ(r)+p(r))ϕ′(r)=−G(ϵ(r)+p(r))(M(r)+4πr3p(r)/c2)c2r2b(r), the Tolman and Klein products are ΘK=Θ(R)b(R) and μK=μ(R)b(R), and the two multipliers are α=μ/Θ and β=1/ΘK. The conclusion is a first-variation criticality statement only; it proves neither a local or global entropy maximum nor existence or dynamical stability of a physical fluid.

Facts & Assumptions

Given: The state chart, positive profiles, fixed radius, uniform positive lower bound on b, and functionals of Proper spherical entropy and particle constrained functionals, together with the conditional mass/lapse equations of Complete spherical Einstein and conservation reduction.

[F1]

The static geometry, signature, and formal SI convention are those of Smooth static spherical geometry and state domain. The functional derivatives are computed on the fixed-radius profile domain with b=1−2Kr∫0rϵ(t)t2dt, K=4πG/c4, S=4π∫0Rsr2b−1/2dr, Nb=4π∫0Rnr2b−1/2dr, and E=4π∫0Rϵr2dr (Proper spherical entropy and particle constrained functionals).

[F2]

The local state chart obeys dϵ=Θ ds+μ dn and ϵ+p=Θs+μn by the definitions in Proper spherical entropy and particle constrained functionals.

[F3]

From the defining integrals and b=1−2GM/(c2r) in Proper spherical entropy and particle constrained functionals, continuity and boundedness of ϵ and p on [0,R] follow from A continuous function on [a,b] is Riemann integrable, by Heine-Cantor and Riemann's criterion, and comparison of the mass integrand with a constant times r2 using If f≤g on [a,b] and both are integrable then ∫abf≤∫abg; and m(b−a)≤∫abf≤M(b−a) gives M(r)=O(r3). The first FTC gives M′(r)=4πr2ϵ(r)/c2 for 0<r<R and hence b′(r)=2GM(r)/(c2r2)−2Krϵ(r). The lapse integrand in the Statement extends continuously to r=0 by M(r)=O(r3), bounded p, and b≥b∗>0; applying the first FTC to that integral gives ϕ′(r)=(GM/c2+4πGr3p/c4)/(r2b) for 0<r<R. These are the same mass and lapse expressions appearing in Complete spherical Einstein and conservation reduction, but the identities here are derived directly at the stated profile regularity, without applying that lemma's stronger smooth-field hypotheses.

[F4]

A C1 map F:U⊆Rm+n→Rn with F(a,b)=0 and invertible derivative in its n correction variables has a local C1 implicit correction map (The Euclidean implicit function theorem with derivative formula).

[F5]

For any nonempty interval there is a smooth nonnegative bump supported in a smaller interval and positive on its interior (A smooth bump between concentric Euclidean balls).

[F6]

Finite-dimensional componentwise algebra and compositions of C1 maps are C1 (Ck Euclidean maps are closed under componentwise algebra and composition). Continuous parameter derivatives pass through compact profile integrals (Leibniz's rule on a compact rectangle: an interior parameter derivative with a continuous extension may be passed through a Riemann integral). The resulting parameter-derivative integrals are continuous because their integrands are uniformly continuous on compact parameter rectangles (For a nonempty subset of Rn with n≥1, compactness, closedness and boundedness, pseudocompactness, and attainment of extrema by every continuous real-valued function are equivalent, Heine-Cantor: a continuous map from a compact metric space to any metric space is uniformly continuous) and the integral error is bounded by interval length times the uniform error (Uniformly close integrable functions have integrals differing by at most the interval length times their uniform error); hence the finite-parameter constraint maps below are C1. Integration by parts holds for differentiable factors with integrable derivatives (If u,v are differentiable on [a,b] with u′,v′ integrable, then ∫abuv′=u(b)v(b)−u(a)v(a)−∫abu′v).

[F8]

Each continuous integrand on a compact interval is Riemann integrable in the Darboux convention (A continuous function on [a,b] is Riemann integrable, by Heine-Cantor and Riemann's criterion).

[F9]

For an integrable real function f on a compact interval, the derivative of x↦∫axf(t) dt equals f(x) at every point where f is continuous (The first fundamental theorem: if f is integrable on [a,b] and continuous at c, then F′(c)=f(c); in particular a continuous f has F as a primitive). Combined with [F8], this applies to the continuous integrands used below.

[F10]

If a continuous nonnegative function on a compact interval is not identically zero, its integral is strictly positive: monotonicity makes the integral nonnegative, while A continuous f≥0 on [a,b] with ∫abf=0 is identically 0 rules out equality to zero (If f≤g on [a,b] and both are integrable then ∫abf≤∫abg; and m(b−a)≤∫abf≤M(b−a)).

[A1]

The countable-choice assumption is inherited from the geometric supplier chain through The Axiom of Countable Choice (ACω); the finite-dimensional correction theorem used in the variational argument is choice-free.

Proof

1.1givenF1F6F7F8algebra

For a smooth compactly supported variation (ϵ˙,n˙) in (0,R), differentiating the mass definition and b=1−2GM/(c2r) gives E˙=4π∫0Rϵ˙r2dr,b˙(r)=−2Kr∫0rt2ϵ˙(t)dt,ddλ∣0bλ(r)−1/2=Krb(r)−3/2∫0rt2ϵ˙(t)dt. The compact-range and positive-b hypotheses give uniform differentiability under the integrals by [F6]. Small multiples of these variations remain admissible because the profile range is compactly contained in the open chart and b≥b∗>0.

1.2givenF1F6F7F8algebra

These first variations also apply to every admissible sup-norm C1 profile curve, with the same tangent at zero. Indeed the linear operator Lf(r):=−2Kr−1∫0rt2f(t)dt for r>0, with Lf(0)=0, maps C0([0,R]) to itself and satisfies ∥Lf∥∞≤(2KR2/3)∥f∥∞. Applying this bound to difference quotients shows that bλ=1+Lϵλ is C1 in the sup norm with derivative Lϵλ′. Evaluation of a sup-norm C1 curve and of its derivative gives jointly continuous functions of (λ,r): continuity follows by splitting the change into a sup-norm parameter change and a spatial change at a fixed parameter. On a sufficiently small closed parameter interval the profiles stay in a common compact subset of U and bλ≥b∗/2, by sup-norm continuity. Consequently the integrands defining S,Nb,E and their parameter derivatives are jointly continuous, using the C1 state functions and the positive-base power rule. The compact-rectangle differentiation theorem in [F6] therefore applies directly to these curves, so the first variations depend only on the tangent at zero, with no assumption of affine dependence on λ.

2.1F1F2F5F10step 1.1algebra

Using sϵ=1/Θ and sn=−μ/Θ, the first variations are S˙=4π∫0Rr2[b−1/2(ϵ˙Θ−μΘn˙)+s ddλ∣0bλ−1/2]dr,N˙b=4π∫0Rr2[b−1/2n˙+n ddλ∣0bλ−1/2]dr. For a nonnegative bump ηn≢0, the continuous integrand ηnr2b−1/2 is nonnegative and nonzero on (0,R), so [F10] gives N˙b=4π∫ηnr2b−1/2dr>0 for the variation (0,ηn), whose E˙=0. For another nonnegative bump ηϵ≢0, [F10] similarly gives E˙=4π∫ηϵr2dr>0. Hence the derivative map (ϵ˙,n˙)↦(E˙,N˙b) is onto R2.

3.1givenF4F6step 2.1algebra

Fix any tangent variation (ϵ˙,n˙) with E˙=N˙b=0, and choose the two bump directions from step 2.1 whose constraint derivative columns are linearly independent. Define G(λ,aϵ,an) to be the pair of constraint changes (E−E[ϵ,n],Nb−Nb[ϵ,n]) along the profiles (ϵ+λϵ˙+aϵηϵ, n+λn˙+anηn). The compact-in-chart and positive-b hypotheses give an open parameter neighbourhood where these profiles remain admissible; [F6] makes G a C1 map there. Its derivative in (aϵ,an) at (0,0,0) has the two independent columns from step 2.1 and is invertible. Apply the C1 implicit function theorem [F4] with path parameter λ and correction variables (aϵ,an). It supplies C1 corrections aϵ(λ),an(λ) with aϵ(0)=an(0)=0 that keep both constraints fixed. Differentiating the constraint equations at zero and using E˙=N˙b=0 gives the invertible correction matrix times (aϵ′(0),an′(0))=0, so both correction derivatives vanish. The correction path remains admissible for small λ, and its derivative at zero is the given smooth compactly supported tangent. Therefore constrained stationarity implies S˙=0 for every tangent variation in the kernel of (E˙,N˙b).

3.2F1F2F6F8F9step 1.1step 2.1algebra

Put q=b−1/2 and L=S+αNb−βE. Substituting the variations from steps 1.1 and 2.1 and integrating by parts in the single radial variable gives L˙=4π∫0Rr2[q(α−μΘ)n˙+{qΘ+K∫rRt(s+αn)b(t)3/2dt−β}ϵ˙]dr. For the nonlocal term, set Q(r)=∫0rt2ϵ˙(t)dt and H(r)=K∫rRt(s+αn)b−3/2dt. By [F8]–[F9], Q,H are differentiable with Q′(r)=r2ϵ˙(r) and H′(r)=−Kr(s+αn)b−3/2. Also H(R)=Q(0)=0, so integration by parts gives K∫0Rr(s+αn)b−3/2Q(r)dr=∫0Rr2H(r)ϵ˙(r)dr. Thus the displayed formula has no omitted boundary term.

4.1step 2.1step 3.1algebra

Because the constraint derivative is onto and S˙ vanishes on its kernel, finite-dimensional linear algebra gives constants α,β such that S˙+αN˙b−βE˙=0 for every compactly supported variation. Conversely, if such multipliers exist, then S˙=0 for every variation with E˙=N˙b=0. The multiplier argument uses only the two bump directions and finite-dimensional algebra.

4.2F2F5F10step 3.2

The two bracket coefficients in step 3.2 are continuous. Each can be tested against every smooth compactly supported variation; if a coefficient were positive or negative at an interior point, continuity and [F5] would give a nonnegative bump supported where it keeps that sign. The resulting continuous integrand, or its negative in the negative case, is nonnegative and not identically zero, so [F10] makes its integral strictly positive; this contradicts the vanishing first variation. Since q>0, this gives μ/Θ=α and qΘ+K∫rRt(s+αn)b(t)3/2dt=β. By the Euler identity in [F2], s+αn=(ϵ+p)/Θ.

5.1A1F1F2F3F7F9step 4.2algebra

Differentiate the last identity on (0,R); [F9] differentiates its variable-endpoint integral, and [F7] differentiates q=b−1/2, giving Θ′Θ=−b′2b−Kr(ϵ+p)b. The mass and b′ identities in [F3] yield b′=2GM/(c2r2)−2Krϵ, so Θ′Θ=−GM/(c2r2)+Krpb=−ϕ′. Thus Θeϕ is constant. Since μ/Θ=α, μeϕ=αΘeϕ is constant as well. At r=R, the boundary normalization makes eϕ(R)=b(R), and the energy bracket has integral zero there; hence β=1/(Θ(R)b(R))=1/ΘK.

6.1A1F2F3F7step 5.1algebra

Differentiating ϵ+p=Θs+μn using the calculus rules [F7] and substituting dϵ=Θds+μdn gives dp=s dΘ+n dμ. Since μ/Θ=α, μ′=αΘ′, so p′=sΘ′+nμ′=(s+αn)Θ′=ϵ+pΘΘ′=−(ϵ+p)ϕ′. Substituting the lapse equation from [F3] gives the TOV pressure equation.

7.1A1F2F3step 1.2step 3.2step 4.2step 5.1∎

Conversely, suppose the stated Tolman and Klein products are constant with the boundary normalization. Then μ/Θ is constant. The Euler and differential identities in [F2] imply p′=−(ϵ+p)ϕ′, so the TOV pressure equation follows from [F3]. Reversing the algebra in step 5.1 shows that the energy bracket in step 4.2 has derivative zero; its boundary value at R is β=1/ΘK, so it equals β throughout. The number bracket vanishes because μ/Θ=α. Therefore step 3.2 gives L˙=0 for every variation, and any constraint-preserving variation has S˙=0. This proves constrained stationarity and both directions.

Remarks

The variations establish a first-order critical-point equivalence only. Neither the source derivation nor this proof establishes a maximum, a second-variation sign, existence of a critical profile, a general equation of state, or dynamical stability. In the formal SI bookkeeping, α=μ/Θ has units J K−1, ΘK has units K, μK has units J, and β=1/ΘK has units K−1. The boundary normalization is a clock convention for this fixed-radius calculation, not a vacuum-matching assumption.

5 · Examples, experiments and counterexamples

None yet.

Sources