p2.py — src/27-fields

1# framed-rational status: MIXED -- exact framed-rational core; continuum constructs are confined to clearly-labelled [approx] / degenerate-idealisation comparisons (Tier 2/3/5 or measured-data), never to an exact claim.2"""3P2: the order-one coefficient of the bare electromagnetic coupling.45Claim: the coefficient is exactly 1, alpha_bare = 1/(4pi), with NO open order-one6factor, because (i) channel unity makes the EM connection stiffness beta equal to7the one per-cell capacity kappa, (ii) that unit capacity is the action quantum8hbar, and (iii) charge = winding index = one phase quantum (EM-1), so in natural9units the elementary charge is 1 and alpha_bare = e^2/4pi = 1/4pi.1011We back this with the saturation-regime computation (the EM analogue of gravity's12exact static profile): the capacity-bounded Gauss law 4 pi r^2 kappa sin(E) = q13gives the weak-field Coulomb coefficient 1/(4 pi kappa), an EM saturation core14r* = sqrt(q/4 pi kappa), and a FINITE charge self-energy -- the classical15self-energy divergence cut off by the same capacity bound that gives gravity its16slip core.  The 4pi is the 3-D Gauss solid angle (continuum-comparison reading).17"""18import math1920print("=" * 68)21print("P2  the bare electromagnetic coupling coefficient")22print("=" * 68)2324# ---- (1) the coefficient from channel unity -------------------------------25print("\n(1) coefficient from channel unity + charge=action=phase quantum")26print("    gravity:  G_shell      = 1/(4 pi kappa),  kappa = 1   (unit capacity)")27print("    EM:       Coulomb coef = 1/(4 pi beta),   beta  = kappa = 1  (channel unity)")28print("    charge quantum = action quantum = one phase quantum => e = 1 (natural units)")29alpha_bare = 1/(4*math.pi)30print(f"    => alpha_bare = e^2/(4 pi) = 1/(4 pi) = {alpha_bare:.6f}   (coefficient = 1)")3132# ---- (2) saturation-regime profile: weak-field coefficient ----------------33# capacity-bounded Gauss law:  4 pi r^2 kappa sin(E(r)) = q   (kappa=1, q=1)34# weak field:  E(r) ~ q/(4 pi r^2),  potential A0(r) = integral_r^inf E dr'35def E_field(r, q=1.0, kappa=1.0):36    x = q/(4*math.pi*kappa*r**2)37    return math.asin(x) if x <= 1.0 else None          # undefined inside the core38def A0(r, q=1.0, kappa=1.0, R=1e4, n=200000):39    # potential = integral_r^R of E dr'  (R large)40    h = (R-r)/n; s = 0.041    for i in range(n):42        rr = r + (i+0.5)*h; e = E_field(rr, q, kappa)43        s += e*h44    return s45print("\n(2) saturation profile  4 pi r^2 sin(E)=q  -> weak-field Coulomb coefficient")46for r in (5.0, 10.0, 20.0):47    a0 = A0(r); coef = a0*r48    print(f"    r={r:5.1f}:  r*A0(r) = {coef:.5f}   (-> 1/4pi = {alpha_bare:.5f})")4950# ---- (3) the EM saturation core and a finite self-energy ------------------51q = kappa = 1.052rstar = math.sqrt(q/(4*math.pi*kappa))                 # where sin(E)=153print(f"\n(3) EM saturation core  r* = sqrt(q/4 pi kappa) = {rstar:.5f}  (shell units)")54# self-energy with the capacity bound: U = int_{r*}^inf (1/2) E^2 4 pi r^2 dr (finite),55# vs continuum int_0^inf which diverges at r->0.56def self_energy(rmin, R=1e4, n=400000):57    h=(R-rmin)/n; s=0.058    for i in range(n):59        rr=rmin+(i+0.5)*h; e=E_field(rr,q,kappa)60        if e is None: continue61        s += 0.5*e*e*4*math.pi*rr*rr*h62    return s63U_cut = self_energy(rstar)64# continuum (no cutoff, no saturation): (1/2)(q/4pi r^2)^2 4pi r^2 = q^2/(8 pi r^2)65# integral_{eps}^inf = q^2/(8 pi eps)  -> diverges as eps->066print(f"    finite self-energy U(>r*) = {U_cut:.5f}   (continuum int_0 DIVERGES)")67print(f"    continuum q^2/(8 pi r*)   = {q*q/(8*math.pi*rstar):.5f}  (same order, finite)")68print("    => the capacity bound cuts off the classical self-energy divergence,")69print("       the EM analogue of gravity's horizonless slip core.")7071# ---- (4) honest separation from the running ------------------------------72print("\n(4) the residual to the measured value is NOT a bare-coefficient ambiguity")73print(f"    1/alpha_bare = 4 pi = {4*math.pi:.4f}      (coefficient fixed = 1)")74print(f"    1/alpha(0)   = 137.036                  (running P1 + EW mixing)")75print("    P2 closes the bare coefficient; the gap is P1 (running) + the Weinberg")76print("    mixing, both separate problems.  EM's bare coupling is MORE constrained")77print("    than gravity's open order-one constants (a0's 2pi, the entropy c_S'),")78print("    because channel unity locks it to hbar.")79print("=" * 68)80