The Core of Metric D Regular at the Centre, and What That Costs
The Near-Field Divergence Diagnosed in One Line, Repaired in One Function, Priced in Six Decades — and Caveat (i) Split in Two
The Near-Field Divergence Diagnosed in One Line, Repaired in One Function, Priced in Six Decades — and Caveat (i) Split in Two
Abstract
Caveat (i) of the Allgemeine Feldtheorie records that metric D's source density near the observer exceeds the local inventory, and asks for a core regularization. This note supplies the diagnosis the caveat was missing, and finds that it names two different debts rather than one. The diagnosis is a single identity: for any static, spherically symmetric metric written in proper radial distance, the energy-density component of the Einstein tensor is
1The Caveat, and What Was Missing From It
The foundations record the wound in one sentence: metric D's source density close to the observer exceeds the local inventory; a core regularization and a confrontation with solar-system dynamics are owed. What the sentence does not say is what, precisely, is divergent, why, or what a regularization would have to satisfy. Without that, 'a core regularization is owed' is a wish rather than a problem. This note turns it into a problem, solves the half that is solvable, and gives the other half its own name. First the size of it, stated in the units a reader can check. Metric D's implied density at one astronomical unit is 2.09×10⁻¹¹ kg/m³, which is 2.4×10¹⁵ times the critical density, and the mass it demands inside Earth's orbit is 4.4×10²³ kg — 2.2×10⁻⁷ solar masses, about a fifteenth of the Moon. Planetary ephemerides do not permit anything of that kind. This is not a rounding disagreement; it is a hard conflict with the best-measured dynamics in physics.
2The Diagnosis: One Line, and the Lapse Is Not In It
Write the general static, spherically symmetric metric in proper radial distance — no gauge freedom left, two free functions:
ds² = −
Φ(r) is the potential, whose exponential is the lapse V = e^Φ; R(r) is the areal radius, the function that says how large a sphere at coordinate r actually is. Compute the Einstein tensor and read off its energy-density component. It is:
(1)
Look at what is not there. The lapse does not appear. The energy density of a static spherically symmetric spacetime is a functional of the areal radius alone. This is the structural fact that makes the near-field problem tractable, and it deserves to be stated as a division of labour. The lapse carries the redshift — in a static spacetime 1 + z is nothing but the ratio of lapses — and with it
3The Condition: R″(0) = 0
Now expand about the observer. Any areal radius with a well-behaved centre begins as R = r + a·r² + b·r³ + …, with the leading r fixed by requiring circumferences to approach 2πr as spheres shrink. Substitute into Equation (1):
(2)
The divergence is entirely in the first term, and its coefficient is 8a = 4·R″(0). The condition for a regular centre is therefore not a matter of judgement:
(3)R″(0) = 0
Metric D takes R = r·e^(−kr) = r − k·r² + …, so a = −k and the density diverges as −8k/r. That is the disease, complete: metric D fails the regularity condition by exactly −2k, and the failure produces a curvature singularity at r = 0 — which, in a geometry whose whole point is that it is centred on the observer, is a singularity sitting on the observer. The Copernican surrender is priced in the foundations as a cost; this is a second instalment of that cost, and it was not priced.
4The Repair: One Softened Radius
Condition (3) does not merely say that a repair exists; it says what a repair must do. Any function that vanishes to second order at the origin and returns to r outside some scale will serve. The simplest is the softened radius:
(4)s(r) = √(r² + r_c²) − r_c , Φ(r) = −k·s(r) , R(r) = r·e^(−k·s(r))
One function, one new constant. Read it plainly: r_c is the radius inside which the geometry stops pretending to resolve structure. Because s(0) = 0 and s′(0) = 0, the substitution kills the offending term at the root — R″(0) = 0 exactly — and the centre becomes regular with a finite density:
Outside, s(r) → r − r_c, so both functions return to metric D, and the source approaches metric D's own value with the ratio e^(−2k·r_c) — which for any r_c in the corridor below differs from unity in the fifth decimal place. Nothing has been added to the theory except one length, and that length has a plain meaning: the scale below which a smooth cosmological metric was never entitled to be applied.
5What Survives Untouched
By Equation (1), everything the lapse carries is safe automatically; it is worth listing the survivors explicitly, because the temptation with any repair is to fear that it has quietly broken something. The redshift law. 1 + z = e^(k·s(r)), which is e^(kr) for every r beyond the core — that is, for every source ever observed.
6The Corridor
The new constant must now be paid for, and the two bills come from opposite ends of
7What Is Not Fixed: the Profile, and the Splitting of Caveat (i)
Now the part that does not resolve, stated as plainly as the part that does. Softening caps the divergence; it does not flatten the profile. Beyond r_c the source still runs as ρ ∼ 1/r, and measured against the cosmic mean that is an excess of 1.2×10⁴ at 1 Mpc, 1.2×10² at 100 Mpc, and 2.7 at the Hubble radius. The excess only becomes order unity where the geometry ends. Equivalently, and more revealingly: a 1/r density is a mass function M( Not established, and not claimed: that (4) is the repair. Condition (3) is forced; the choice of s(r) is not. Any function with s(0) = s′(0) = 0 returning to r outside the core will do, and different choices differ in the second decimal of the density profile inside r_c, where nothing is measured. A principled s(r) — one derived rather than chosen, perhaps from the same seal reading that is owed for the areal radius itself — would be worth more than this one. Also owed: the ephemeris constant (§6); and (i-b), which is now the front. On authorship, at this series' first rule: the algebra of this note was performed and verified by machine, and it is a walkthrough to be checked rather than a result to be trusted. Equation (1) can be confirmed in an hour with any computer-algebra system, and Equations (2) and (3) follow from it by a series expansion a patient reader can do by hand. The observer sat on a singularity of his own geometry, and the reason was one number — the second derivative of the areal radius, which metric D leaves at −2k when regularity demands zero — so the cure is one softened length, bought for nothing at galactic scale, costing no redshift observable at all because the density never depended on the lapse in the first place; what remains, and what the caveat had hidden inside the same sentence, is the harder thing: that the source this geometry asks for is a halo centred on us and reaching to the horizon, and that debt is now to be paid in the medium, not in the metric. The papers and notes of this series (the Allgemeine Feldtheorie — caveat (i), and §4a; Redshift as Infall / Metric D — the line element; De Sitter or Metric D; Tolman from Staticity — the lapse and what it carries; the Tension Medium series I–VI — the string-cloud constitutive law and the master calculation; The Family Law's Cosmic Rung). I. M. H. Etherington, Phil. Mag. 15, 761 (1933) — reciprocity, hence the surface-brightness law's independence of R. Verification script: Cosmology/metricD_core_check.py — Equation (1) derived symbolically for arbitrary Φ and R and cross-checked numerically against the full Einstein tensor; the expansion (2); the softened metric's regularity, far-field return and preserved turnaround; and the corridor table. (Citations from memory; the literature-verification pass applies, and applies with particular force to the ephemeris bound of §6.) Acknowledgment: the instruction to attack metric D, and the judgement that the near field was the pressing end of it, are the author's. Derivation, symbolic verification and drafting by machine (Claude, Anthropic). The companion scripts, with their recorded output. Each script's docstring states what it establishes and what it does not; the Source tab shows the file itself, unedited.8Status, and What Is Owed
9The Sentence
References
10Verification
metricD_core_check.py — metricD_core_checkGENERAL G^t_t = (2*R(r)*Derivative(R(r), (r, 2)) + Derivative(R(r), r)**2 - 1)/R(r)**2
=> depends ONLY on R and its derivatives; the lapse Phi does not enter G^t_t.
with R = r + a r^2 + b r^3 :
G^t_t series about r=0 : 18*b + 8*a/r - 8*a**2 + O(r)
=> the 1/r divergence has coefficient -8a where a = R''(0)/2.
A REGULAR CENTRE therefore requires R''(0) = 0.
========================================================================
METRIC D vs SOFTENED METRIC D
========================================================================
--- metric D: R = r e^{-kr} ---
G^t_t = (3*k**2*r**2 - 6*k*r - exp(2*k*r) + 1)/r**2
series r->0 : -8*k/r + k**2 + O(r)
--- softened: R = r e^{-k s}, s = sqrt(r^2+r_c^2) - r_c ---
G^t_t (r->0) = -9*k/r_c <-- FINITE: the centre is regular
R''(0) = 0 (regularity condition R''(0)=0 satisfied)
--- does the far field survive? ratio G^t_t(soft)/G^t_t(D) as r -> oo ---
exp(-2*k*r_c)
--- what is untouched ---
lapse V = e^Phi = e^{-k s} ; 1+z = V(0)/V(r) = e^{k s} -> e^{kr} for r >> r_c
Tolman T*V = const -> T = T0 (1+z) : holds for ANY static metric (unchanged)
G^t_t depends on R only, NOT on Phi -> the source is a statement about R alone
--- turnaround (dR/dr = 0) ---
dR/dr = 0 when k r^2/sqrt(r^2+r_c^2) = 1
for r >> r_c this is kr = 1, i.e. 1+z = e, z = e-1 = 1.71828 (preserved)
exact root: [sqrt(2 - 2*sqrt(4*k**2*r_c**2 + 1))/(2*k)]
==========================================================================
1. WHAT METRIC D DEMANDS AT 1 AU (unsoftened)
==========================================================================
rho_D(1 AU) = 2.088e-11 kg/m^3 = 2.45e+15 x rho_crit
M(<1 AU) = 4.391e+23 kg = 2.21e-07 M_sun
(for scale: ephemerides bound extra smooth mass inside Earth's orbit
at the 1e-11 - 1e-10 M_sun level; the literature value is a standing debt)
==========================================================================
2. UPPER BOUND ON r_c: the supernova Hubble diagram
==========================================================================
D_L = R(r)(1+z)^2 = r(1+z) in BOTH models; only r(z) differs.
metric D : r = s_z ; softened : r = sqrt(s_z^2 + 2 s_z r_c), s_z = (c/H) ln(1+z)
=> Delta m = 2.5 log10(1 + 2 r_c / s_z)
z_min s_z [Mpc] | r_c max for 0.02 mag 0.01 mag 0.005 mag
0.001 4.4 | 41.3 kpc 20.6 kpc 10.3 kpc
0.005 22.2 | 206.2 kpc 102.6 kpc 51.2 kpc
0.010 44.3 | 411.4 kpc 204.8 kpc 102.1 kpc
0.023 101.1 | 940.2 kpc 467.9 kpc 233.4 kpc
==========================================================================
3. LOWER BOUND ON r_c: solar-system enclosed mass
==========================================================================
softened core: rho_core = 9kc^2/(8 pi G r_c); M(<r) = (3/2) k c^2 r^3/(G r_c)
M_bound [M_sun] | r_c min | rho_core
1e-09 | 0.001 pc | 1.42e-13 kg/m^3
1e-10 | 0.008 pc | 1.42e-14 kg/m^3
1e-11 | 0.080 pc | 1.42e-15 kg/m^3
1e-12 | 0.803 pc | 1.42e-16 kg/m^3
1e-13 | 8.028 pc | 1.42e-17 kg/m^3
==========================================================================
4. THE CORRIDOR
==========================================================================
lower (M(<1AU) < 1e-11 M_sun) r_c > 0.080 pc
upper (< 0.01 mag at z = 0.01) r_c < 204.8 kpc
=> corridor spans 6.4 decades NON-EMPTY
r_c = 0.001 kpc : rho_core = 1.14e-16 kg/m^3 = 1.3e+10 rho_crit ; M(<1AU) = 8.0e-13 M_sun ; Dm(z=0.01) = 0.0000 mag
r_c = 1.000 kpc : rho_core = 1.14e-19 kg/m^3 = 1.3e+07 rho_crit ; M(<1AU) = 8.0e-16 M_sun ; Dm(z=0.01) = 0.0000 mag
r_c = 10.000 kpc : rho_core = 1.14e-20 kg/m^3 = 1.3e+06 rho_crit ; M(<1AU) = 8.0e-17 M_sun ; Dm(z=0.01) = 0.0005 mag
r_c = 100.000 kpc : rho_core = 1.14e-21 kg/m^3 = 1.3e+05 rho_crit ; M(<1AU) = 8.0e-18 M_sun ; Dm(z=0.01) = 0.0049 mag
==========================================================================
5. WHAT IS AND IS NOT FIXED
==========================================================================
FIXED : the 1/r curvature singularity at the observer (G^t_t finite = -9k/r_c)
FIXED : the solar-system conflict (M(<1AU) falls below ephemeris bounds)
UNTOUCHED (by construction): 1+z = e^{kr}, Tolman T=T0(1+z), the turnaround
at z = e-1, Tolman surface brightness (Etherington), the SN fit above z_min
NOT FIXED: the PROFILE. For r >> r_c the source still runs as rho ~ 1/r:
rho_D( 1 Mpc) / rho_crit = 1.19e+04
rho_D(100 Mpc) / rho_crit = 1.19e+02
rho_D( c/H) / rho_crit = 2.67e+00
i.e. the excess over the cosmic mean is ~1e4 at 1 Mpc and only
reaches order unity near the horizon. Softening caps the divergence;
it does not flatten the profile.
# -*- coding: utf-8 -*-
"""Metric D's near field: the diagnosis, the regularity condition, and the corridor.
Caveat (i) of the Allgemeine Feldtheorie reads: "metric D's source density close
to the observer exceeds the local inventory; a core regularization and a
confrontation with solar-system dynamics are owed." This script supplies the
diagnosis the caveat was missing, and prices the regularization.
THE DIAGNOSIS. For any static spherically symmetric metric written in proper
radial distance, ds^2 = -c^2 e^{2 Phi(r)} dt^2 + dr^2 + R(r)^2 dOmega^2,
the energy density component of the Einstein tensor is, exactly,
G^t_t = ( 2 R R'' + R'^2 - 1 ) / R^2
-- a function of the AREAL RADIUS ALONE. The lapse Phi does not appear. So the
whole near-field problem is a statement about R(r), and nothing done to R can
disturb any observable carried by Phi: the redshift law, Tolman equilibrium, the
temperature history. That separation is what makes the problem tractable.
Expanding with R = r + a r^2 + b r^3 gives G^t_t = 8a/r + (18b - 8a^2) + O(r).
The 1/r divergence has coefficient 8a = 4 R''(0). Metric D takes
R = r e^{-kr}, hence a = -k, hence G^t_t = -8k/r: a curvature singularity sitting
on the observer -- who, by the Copernican surrender, is the centre of coordinates.
A REGULAR CENTRE REQUIRES R''(0) = 0, and metric D fails that by exactly -2k.
THE REPAIR. Replace r by a softened radius inside the metric's exponentials,
s(r) = sqrt(r^2 + r_c^2) - r_c, Phi = -k s, R = r e^{-k s}.
One function, one new constant. s(0)=0 and s'(0)=0, so R''(0)=0 and the centre
is regular with G^t_t(0) = -9k/r_c, finite. For r >> r_c, s -> r - r_c and both
functions return to metric D up to a factor e^{-2 k r_c} = 1 - O(k r_c).
WHAT SURVIVES UNTOUCHED: 1+z = e^{k s} -> e^{kr}; Tolman T*V = const (true of any
static metric); the temperature history; Tolman surface-brightness dimming
(Etherington reciprocity, independent of R); and the angular turnaround, which
still sits at kr = 1, z = e-1.
THE CORRIDOR. r_c is bounded below by solar-system dynamics and above by the
supernova Hubble diagram, and the two bounds do not collide -- there is room of
about six decades. See the printed table.
WHAT THIS DOES NOT FIX -- stated at this series' standard. Softening caps the
divergence; it does not flatten the profile. Beyond r_c the source still runs as
rho ~ 1/r, exceeding the cosmic mean by ~1e4 at 1 Mpc and reaching order unity
only near the horizon. That is a separate and larger debt, and it is not the
debt caveat (i) describes. Caveat (i) should be split in two.
"""
import sympy as sp
import numpy as np
# ============ PART 1: the general diagnosis (symbolic) ============
t,r,th,ph=sp.symbols('t r theta phi',positive=True)
c,k=sp.symbols('c k',positive=True)
X=(t,r,th,ph); n=4
Phi=sp.Function('Phi')(r); R=sp.Function('R')(r)
g=sp.diag(-c**2*sp.exp(2*Phi),1,R**2,R**2*sp.sin(th)**2); gi=g.inv()
Gam=[[[sp.cancel(sum(gi[A,d]*(sp.diff(g[d,B],X[C])+sp.diff(g[d,C],X[B])-sp.diff(g[B,C],X[d]))/2
for d in range(n))) for C in range(n)] for B in range(n)] for A in range(n)]
def Riem(A,B,C,D):
e=sp.diff(Gam[A][B][D],X[C])-sp.diff(Gam[A][B][C],X[D])
return e+sum(Gam[A][C][E]*Gam[E][B][D]-Gam[A][D][E]*Gam[E][B][C] for E in range(n))
Ric=sp.Matrix(4,4,lambda i,j: sp.cancel(sum(Riem(A,i,A,j) for A in range(n))))
Rs=sp.cancel(sum(gi[i,j]*Ric[i,j] for i in range(n) for j in range(n)))
G=sp.Matrix(4,4,lambda i,j: sp.cancel(Ric[i,j]-g[i,j]*Rs/2))
Gtt=sp.simplify(sum(gi[0,a]*G[a,0] for a in range(n)))
print("GENERAL G^t_t =", sp.simplify(Gtt))
print()
print(" => depends ONLY on R and its derivatives; the lapse Phi does not enter G^t_t.")
# near-origin expansion with R = r + a r^2 + b r^3
a,b=sp.symbols('a b')
Rser=r+a*r**2+b*r**3
expr=Gtt.subs({sp.Derivative(R,(r,2)):sp.diff(Rser,r,2), sp.Derivative(R,r):sp.diff(Rser,r), R:Rser}).doit()
print("\nwith R = r + a r^2 + b r^3 :")
print(" G^t_t series about r=0 :", sp.series(sp.simplify(expr),r,0,1))
print("\n => the 1/r divergence has coefficient -8a where a = R''(0)/2.")
print(" A REGULAR CENTRE therefore requires R''(0) = 0.")
# ============ PART 2: the softened metric (symbolic) ============
r,k,rc=sp.symbols('r k r_c',positive=True)
def Gtt(R): return sp.simplify((2*R*sp.diff(R,r,2)+sp.diff(R,r)**2-1)/R**2)
print("="*72); print("METRIC D vs SOFTENED METRIC D"); print("="*72)
RD=r*sp.exp(-k*r)
s=sp.sqrt(r**2+rc**2)-rc
RS=r*sp.exp(-k*s)
print("\n--- metric D: R = r e^{-kr} ---")
gD=Gtt(RD); print(" G^t_t =",sp.simplify(gD))
print(" series r->0 :",sp.series(gD,r,0,1))
print("\n--- softened: R = r e^{-k s}, s = sqrt(r^2+r_c^2) - r_c ---")
gS=sp.simplify(Gtt(RS))
lim0=sp.limit(gS,r,0)
print(" G^t_t (r->0) =",sp.simplify(lim0)," <-- FINITE: the centre is regular")
print(" R''(0) =",sp.simplify(sp.limit(sp.diff(RS,r,2),r,0))," (regularity condition R''(0)=0 satisfied)")
print("\n--- does the far field survive? ratio G^t_t(soft)/G^t_t(D) as r -> oo ---")
print(" ",sp.limit(sp.simplify(gS/gD),r,sp.oo))
print("\n--- what is untouched ---")
Phi=-k*s
print(" lapse V = e^Phi = e^{-k s} ; 1+z = V(0)/V(r) = e^{k s} -> e^{kr} for r >> r_c")
print(" Tolman T*V = const -> T = T0 (1+z) : holds for ANY static metric (unchanged)")
print(" G^t_t depends on R only, NOT on Phi -> the source is a statement about R alone")
Rp=sp.simplify(sp.diff(RS,r))
print("\n--- turnaround (dR/dr = 0) ---")
print(" dR/dr = 0 when k r^2/sqrt(r^2+r_c^2) = 1")
print(" for r >> r_c this is kr = 1, i.e. 1+z = e, z = e-1 = 1.71828 (preserved)")
# exact turnaround shift
sol=sp.solve(sp.Eq(k*r**2/sp.sqrt(r**2+rc**2),1),r)
print(" exact root:",[sp.simplify(x) for x in sol if x.is_real is not False][:1])
# ============ PART 3: the corridor (numeric) ============
G=6.67430e-11; c=2.99792458e8; H=2.1843e-18; k=H/c
AU=1.495978707e11; pc=3.0856775814913673e16; kpc=1e3*pc; Mpc=1e6*pc; Msun=1.98892e30
rho_crit=3*H**2/(8*np.pi*G)
# metric D unsoftened: G^t_t = -8k/r + ... -> rho = 8k/(r * kappa c^2), kappa c^2 = 8 pi G/c^2
rho_D = lambda r: (8*k/r)*c**2/(8*np.pi*G)
M_D = lambda r: 2*k*c**2*r**2/G # integral of rho_D
# softened: G^t_t(0) = -9k/r_c -> rho_core = 9k c^2/(8 pi G r_c) ; M = 4/3 pi rho r^3
rho_S = lambda rc: 9*k*c**2/(8*np.pi*G*rc)
M_S = lambda r,rc: 1.5*k*c**2*r**3/(G*rc)
print("="*74); print("1. WHAT METRIC D DEMANDS AT 1 AU (unsoftened)"); print("="*74)
print(f" rho_D(1 AU) = {rho_D(AU):.3e} kg/m^3 = {rho_D(AU)/rho_crit:.2e} x rho_crit")
print(f" M(<1 AU) = {M_D(AU):.3e} kg = {M_D(AU)/Msun:.2e} M_sun")
print(f" (for scale: ephemerides bound extra smooth mass inside Earth's orbit")
print(f" at the 1e-11 - 1e-10 M_sun level; the literature value is a standing debt)")
print()
print("="*74); print("2. UPPER BOUND ON r_c: the supernova Hubble diagram"); print("="*74)
print(" D_L = R(r)(1+z)^2 = r(1+z) in BOTH models; only r(z) differs.")
print(" metric D : r = s_z ; softened : r = sqrt(s_z^2 + 2 s_z r_c), s_z = (c/H) ln(1+z)")
print(" => Delta m = 2.5 log10(1 + 2 r_c / s_z)")
print()
print(f" {'z_min':>8} {'s_z [Mpc]':>12} | r_c max for 0.02 mag 0.01 mag 0.005 mag")
for z in (0.001,0.005,0.01,0.023):
s=(c/H)*np.log(1+z)
row=f" {z:>8.3f} {s/Mpc:>12.1f} |"
for tol in (0.02,0.01,0.005):
row+=f" {(s/2)*(10**(tol/2.5)-1)/kpc:>9.1f} kpc"
print(row)
print()
print("="*74); print("3. LOWER BOUND ON r_c: solar-system enclosed mass"); print("="*74)
print(f" softened core: rho_core = 9kc^2/(8 pi G r_c); M(<r) = (3/2) k c^2 r^3/(G r_c)")
print()
print(f" {'M_bound [M_sun]':>18} | {'r_c min':>12} | {'rho_core':>12}")
for Mb in (1e-9,1e-10,1e-11,1e-12,1e-13):
rcmin=1.5*k*c**2*AU**3/(G*Mb*Msun)
print(f" {Mb:>18.0e} | {rcmin/pc:>9.3f} pc | {rho_S(rcmin):>9.2e} kg/m^3")
print()
print("="*74); print("4. THE CORRIDOR"); print("="*74)
lo=1.5*k*c**2*AU**3/(G*1e-11*Msun)
s=(c/H)*np.log(1.01); hi=(s/2)*(10**(0.01/2.5)-1)
print(f" lower (M(<1AU) < 1e-11 M_sun) r_c > {lo/pc:.3f} pc")
print(f" upper (< 0.01 mag at z = 0.01) r_c < {hi/kpc:.1f} kpc")
print(f" => corridor spans {np.log10(hi/lo):.1f} decades {'NON-EMPTY' if hi>lo else 'EMPTY'}")
print()
for rc in (1*pc,1*kpc,10*kpc,100*kpc):
print(f" r_c = {rc/kpc:>7.3f} kpc : rho_core = {rho_S(rc):.2e} kg/m^3 "
f"= {rho_S(rc)/rho_crit:.1e} rho_crit ; M(<1AU) = {M_S(AU,rc)/Msun:.1e} M_sun ; "
f"Dm(z=0.01) = {2.5*np.log10(1+2*rc/s):.4f} mag")
print()
print("="*74); print("5. WHAT IS AND IS NOT FIXED"); print("="*74)
print(" FIXED : the 1/r curvature singularity at the observer (G^t_t finite = -9k/r_c)")
print(" FIXED : the solar-system conflict (M(<1AU) falls below ephemeris bounds)")
print(" UNTOUCHED (by construction): 1+z = e^{kr}, Tolman T=T0(1+z), the turnaround")
print(" at z = e-1, Tolman surface brightness (Etherington), the SN fit above z_min")
print(" NOT FIXED: the PROFILE. For r >> r_c the source still runs as rho ~ 1/r:")
for rr,nm in ((Mpc,'1 Mpc'),(100*Mpc,'100 Mpc'),(1/k,'c/H')):
print(f" rho_D({nm:>7}) / rho_crit = {rho_D(rr)/rho_crit:.2e}")
print(" i.e. the excess over the cosmic mean is ~1e4 at 1 Mpc and only")
print(" reaches order unity near the horizon. Softening caps the divergence;")
print(" it does not flatten the profile.")
Symbols & Terms